get_lrt_results

Synopsis

Access the simulation results of the likelihood ratio test.

Syntax

get_lrt_results()

Description

The get_lrt_results() command returns the likelihood ratio test
results computed by the plot_lrt() command, comparing fits of the
null model to fits of the alternative model using faked data with
Poisson noise. The likelihood ratio based on the observed data is
returned, along with the p-value, used to reject or accept the null model.

The following attributes are returned by get_lrt_results():

samples - parameter samples from all iterations, in the order of
the parameter list (as many as samples as specified in the 'num'
parameter of plot_lrt() are returned; default is 500)

Call the get_lrt_results() command with no arguments to print
to the screen the results of the likelihood ratio test performed by
the plot_lrt() command. View the array of ratios comparing fits
to simulated data done with a simple (null) power-law model
versus the more-complex (alternative) power-law plus
Gaussian model, along with the likelihood ratio and p-value
based on the observed data.

Example 2

lrt = get_lrt_results()

Store the data returned by get_lrt_results() to variable "lrt", in order to access and manipulate individual attributes within a script.