threeML.classicMLE.likelihood_ratio_test module¶
-
class
threeML.classicMLE.likelihood_ratio_test.
LikelihoodRatioTest
(joint_likelihood_instance0: threeML.classicMLE.joint_likelihood.JointLikelihood, joint_likelihood_instance1: threeML.classicMLE.joint_likelihood.JointLikelihood)¶ Bases:
object
-
property
TS_distribution
¶
-
by_mc
(n_iterations=1000, continue_on_failure=False, save_pha=False)¶ Compute the Likelihood Ratio Test by generating Monte Carlo datasets and fitting the current models on them. The fraction of synthetic datasets which have a value for the TS larger or equal to the observed one gives the null-hypothesis probability (i.e., the probability that the observed TS is obtained by chance from the null hypothesis)
- Parameters
n_iterations – number of MC iterations to perform (default: 1000)
continue_of_failure – whether to continue in the case a fit fails (False by default)
save_pha – Saves pha files for reading into XSPEC as a cross check. Currently only supports OGIP data. This can become slow! (False by default)
- Returns
tuple (null. hyp. probability, TSs, frame with all results, frame with all likelihood values)
-
get_models
(id)¶
-
get_simulated_data
(id: int)¶
-
property
null_hypothesis_probability
¶
-
plot_TS_distribution
(show_chi2=True, scale=1.0, **hist_kwargs)¶ - Parameters
show_chi2 –
scale –
hist_kwargs –
- Returns
-
property
reference_TS
¶
-
property