cbcBayesThermoInt.py File Reference

Prototypes
def	cbcBayesThermoInt.extract_temp (filename)
	Extracts the PTMCMC temperature from the header lines of the given file. More...
def	cbcBayesThermoInt.get_mean_logl (filename)
	Returns the mean value of log(L) from the given filename, excluding the first 50% of samples as burn-in. More...

Go to the source code of this file.

Namespaces
namespace	cbcBayesThermoInt

Variables
string	cbcBayesThermoInt.usage
	cbcBayesThermoInt.parser = OptionParser(usage=usage)
	cbcBayesThermoInt.metavar
	cbcBayesThermoInt.default
	cbcBayesThermoInt.help
	cbcBayesThermoInt.options
	cbcBayesThermoInt.args
	cbcBayesThermoInt.betas = np.array([1.0/extract_temp(f) for f in args])
	cbcBayesThermoInt.logls = np.array([get_mean_logl(f) for f in args])
	cbcBayesThermoInt.inds = np.argsort(betas)[::-1]
	cbcBayesThermoInt.ebetas = np.concatenate((betas, [0.0]))
	cbcBayesThermoInt.elogls = np.concatenate((logls, [logls[-1]]))
	cbcBayesThermoInt.ebetas2 = np.concatenate((betas[::2], [0.0]))
	cbcBayesThermoInt.elogls2 = np.concatenate((logls[::2], [logls[::2][-1]]))
	cbcBayesThermoInt.evidence = -trapezoid(elogls, ebetas)
	cbcBayesThermoInt.evidence2 = -trapezoid(elogls2, ebetas2)
	cbcBayesThermoInt.devidence = np.abs(evidence - evidence2)