No such module: plotutils about skyarea

[farr]

Hi Leo. Plotutils is another package that I have on GitHub: https://github.com/farr/plotutils

[lpsinger]

What about using https://versions.ligo.org/cgit/lalsuite/tree/lalinference/python/lalinference/plot/pp.py here?

[farr]

Yeah, go for it!

Will

On Tue, Oct 18, 2016 at 8:47 PM Leo Singer [email protected] wrote:

What about using
https://versions.ligo.org/cgit/lalsuite/tree/lalinference/python/lalinference/plot/pp.py
here?

—
You are receiving this because you commented.

Reply to this email directly, view it on GitHub
#47 (comment), or mute
the thread
https://github.com/notifications/unsubscribe-auth/AAFfwM_GrXx359eJlliXj4qVi_WMzJaWks5q1SI9gaJpZM4KaDZW
.

[lpsinger]

Oh, OK! Do you have some example inputs for this script?

[lpsinger]

Wait a second. The PPPlot axes class will work for the first call to plot_cumulative_distribution, but not the last four.

Do you mind if I exclude process_areas from the Debian and SL packages, to avoid the import error messages? Is this code used in the pipeline or any of the review tests?

[farr]

Hi Leo,

It may be required for the LI validation tests (diagonal p-p plots in sky
location). Here is an example input set (and corresponding output
areas.dat.bz2).

You could try replacing plot_cumulative_distribution from plotutils with
its code:

def plot_cumulative_distribution(pts, *args, **kwargs):
    """Plots the 1D normalized empirical CDF for the given points.
    Additional arguments are passed to matplotlib's ``plot``.

    """
    pts = np.atleast_1d(pts)
    pts = np.sort(np.concatenate(([0], pts)))

    pp.plot(pts, np.linspace(0, 1, pts.shape[0]), *args, **kwargs)

That shouldn’t cause too much trouble.

Will

On Tue, Oct 18, 2016 at 9:40 PM Leo Singer [email protected]
http://mailto:[email protected] wrote:

Wait a second. The PPPlot axes class will work for the first call to

plot_cumulative_distribution, but not the last four.

Do you mind if I exclude process_areas from the Debian and SL packages,
to avoid the import error messages? Is this code used in the pipeline or
any of the review tests?

—
You are receiving this because you commented.

Reply to this email directly, view it on GitHub
#47 (comment), or mute
the thread
https://github.com/notifications/unsubscribe-auth/AAFfwHmw2yLwDJdC3Wl91BsYP6gWuxnxks5q1S6xgaJpZM4KaDZW
.

​

[lpsinger]

Should there be an attachment somewhere?

[farr]

Oops---that's what happens when you reply to the email rather than editing the form. Here's the attachment.
pa-example.tar.gz

[lpsinger]

I get this error message:

$ process_areas
Traceback (most recent call last):
  File "/Users/leo/local/bin/process_areas", line 87, in <module>
    plt.xscale('log')
  File "/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/matplotlib/pyplot.py", line 1631, in xscale
    gca().set_xscale(*args, **kwargs)
  File "/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/matplotlib/axes/_base.py", line 2833, in set_xscale
    self.autoscale_view(scaley=False)
  File "/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/matplotlib/axes/_base.py", line 2168, in autoscale_view
    x0, x1 = xlocator.view_limits(x0, x1)
  File "/opt/local/Library/Frameworks/Python.framework/Versions/2.7/lib/python2.7/site-packages/matplotlib/ticker.py", line 1617, in view_limits
    "Data has no positive values, and therefore can not be "
ValueError: Data has no positive values, and therefore can not be log-scaled.```

[farr]

Oops. I sent you a data set that doesn't have any injection data (i.e. there are no entries in the searched_area column). You need to run process_areas with the --noinj flag:

farr$ tar xzf pa-example.tar.gz
farr$ process_areas --noinj
/Users/farr/Library/Enthought/Canopy_64bit/User/lib/python2.7/site-packages/matplotlib/font_manager.py:273: UserWarning: Matplotlib is building the font cache using fc-list. This may take a moment.
  warnings.warn('Matplotlib is building the font cache using fc-list. This may take a moment.')
farr$ bzcat areas.dat.bz2 | less