cbcBayesBurstPPAnalysis.py File Reference

Prototypes
def	cbcBayesBurstPPAnalysis.fractional_rank (x, xs)
	Returns the fraction of samples, xs, that fall below the value x. More...
def	cbcBayesBurstPPAnalysis.pp_plot (ps, title=None, outfile=None)
	Generates a p-p plot for the given ps. More...
def	cbcBayesBurstPPAnalysis.pp_kstest_pvalue (ps)
	Returns the K-S p-value for the test of the given ps against a uniform distribution on [0,1]. More...
def	cbcBayesBurstPPAnalysis.read_posterior_samples (f)
	Returns a named numpy array of the posterior samples in the file f. More...
def	cbcBayesBurstPPAnalysis.output_html (outdir, ks_pvalues, injnum, skypp=False)
	Outputs the HTML page summarizing the results. More...

Go to the source code of this file.

Data Structures
class	cbcBayesBurstPPAnalysis.LIGOLWContentHandlerExtractSimBurstTable

Namespaces
namespace	cbcBayesBurstPPAnalysis

Variables
dictionary	cbcBayesBurstPPAnalysis.posterior_name_to_sim_burst_extractor
dictionary	cbcBayesBurstPPAnalysis.posterior_name_to_latex_name
string	cbcBayesBurstPPAnalysis.USAGE
	cbcBayesBurstPPAnalysis.parser = optparse.OptionParser(USAGE)
	cbcBayesBurstPPAnalysis.action
	cbcBayesBurstPPAnalysis.type
	cbcBayesBurstPPAnalysis.dest
	cbcBayesBurstPPAnalysis.help
	cbcBayesBurstPPAnalysis.default
	cbcBayesBurstPPAnalysis.None
	cbcBayesBurstPPAnalysis.options
	cbcBayesBurstPPAnalysis.args
	cbcBayesBurstPPAnalysis.injs = lsctables.SimBurstTable.get_table(utils.load_filename(options.injxml,contenthandler=LIGOLWContentHandlerExtractSimBurstTable))
list	cbcBayesBurstPPAnalysis.parameters = ['frequency', 'quality', 'hrss', 'ra', 'dec', 'psi', 'time', 'alpha','polar_eccentricity']
dictionary	cbcBayesBurstPPAnalysis.pvalues = { }
	cbcBayesBurstPPAnalysis.posfiles = args
int	cbcBayesBurstPPAnalysis.Ninj = 0
def	cbcBayesBurstPPAnalysis.psamples = read_posterior_samples(posfile)
	cbcBayesBurstPPAnalysis.true_params = injs[index]
def	cbcBayesBurstPPAnalysis.samples = psamples[par]
dictionary	cbcBayesBurstPPAnalysis.true_value = posterior_name_to_sim_burst_extractor[par](true_params)
def	cbcBayesBurstPPAnalysis.p = fractional_rank(true_value, samples)
dictionary	cbcBayesBurstPPAnalysis.ks_pvalues = {}
	cbcBayesBurstPPAnalysis.ps
	cbcBayesBurstPPAnalysis.title
	cbcBayesBurstPPAnalysis.outfile
bool	cbcBayesBurstPPAnalysis.skypp = False
int	cbcBayesBurstPPAnalysis.found = 0
	cbcBayesBurstPPAnalysis.inf = os.path.join(os.path.realpath(options.skypp),'p-p.%s'%i)
	cbcBayesBurstPPAnalysis.outf = os.path.join(options.outdir,'sky_p-p.%s'%i)
	cbcBayesBurstPPAnalysis.outdir