README.md

PDF Corner Plots

Python script for plotting probability density functions, stored as importance-weighted lists of samples.

Phil Marshall and Daniel Foreman-Mackey since 12 December, 2011

Installation

setenv PAPPY_DIR ${HOME}/python/pappy
setenv PATH ${PATH}:${PAPPY_DIR}
setenv PYTHONPATH ${PAPPY_DIR}:${PYTHONPATH}

Getting started

CornerPlotter.py --test

This suggests the following examples:

CornerPlotter.py -o test1.png examples/thetas.cpt,blue,shaded
CornerPlotter.py -o test2.png -w 1 -n 2,3,4 examples/localgroup.cpt,red,shaded
CornerPlotter.py -o test3.png --just-2D -w 1 -n 3,4 examples/localgroup.cpt,purple,shaded
CornerPlotter.py -o test4.png --conditional --plot-points -n 2,3,4 examples/localgroup.cpt,green,shaded

Note the data format required: plain text with an informative header, that defines the axis labels and ranges. For example:

# importance, $\log_{10} M_{\rm MW} / M_{\odot}$,  $\log_{10} M_{\rm M31} / M_{\odot}$,  $\log_{10} M_{\rm M31} / M_{\rm MW}$,  $\log_{10} M\prime_{\rm LG} / M_{\odot}$,  
# 0,1,        10.0,14.0,    10.0,14.0,    -3.0,3.0,     10.0,14.0,
1.228662286689141346e-10 1.148898129665888135e+01 1.206671702374781496e+01 5.777357270889337260e-01 1.216860210104861117e+01
9.115455789900011041e-01 1.220497250067228023e+01 1.144655000180502746e+01 -7.584224988672524370e-01 1.227479314997625259e+01
...

Here's the test2.png corner plot: