The LikelihoodRatioStatistic function computes the likelihood ratio statistic of a random variable R evaluated on a dataset V.

The first parameter R can be a distribution (see Statistics[Distribution]), a random variable, or an algebraic expression involving random variables (see Statistics[RandomVariable]).

The second parameter V can be an Array of data samples or a symbol representing an Array of data samples (in which case the option samplesize must be specified).

The options argument can contain one or more of the options shown below.

samplesize = deduce or posint -- If set to 'deduce' (default) the function attempts to automatically determine the number of data samples provided in V.  This parameter must be specified if the number of samples in V is not immediately detectable.

bounds = none or range -- If this option specifies a range, then the function will attempt to calculate the maximum likelihood estimate within the given bounds when constructing the likelihood ratio statistic.  Otherwise the function will attempt to calculate the maximum likelihood estimate over the entire real line (default).

params = list(name=realcons) -- This option specifies parameters that should be evaluated prior to calculating the likelihood ratio statistic.  This function will only work if exactly one indeterminate is present in the final expression - hence all remaining indeterminates must be evaluated prior to generating the final expression.

ignore = truefalse -- This option is used to specify how to handle non-numeric data. If ignore is set to true all non-numeric items in V will be ignored.

weights = rtable -- Vector of weights (one-dimensional rtable). If weights are given, the LikelihoodRatioStatistic function will scale each data point to have given weight. Note that the weights provided must have type realcons and the results are floating-point, even if the problem is specified with exact values. Both the data array and the weights array must have the same number of elements.

Note that for numeric solutions this function uses NLPSolve which is not sensitive to assumptions made with the assume command.  For symbolic solutions this function uses solve instead and ignores the bounds option.

$\mathrm{with}\left(\mathrm{Statistics}\right)\:$

$\mathrm{LikelihoodRatioStatistic}\left(\mathrm{Binomial}\left(10\,\mathrm{\theta }\right)\,A\,\mathrm{samplesize}=1\right)$

$2\ln \left(\left(\genfrac{}{}{0ex}{}{10}{A_1}\right){10}^{-A_1}{A}_1^{A_1}{\left(1-\frac{A_1}{10}\right)}^{10-A_1}\right)-2\ln \left(\left(\genfrac{}{}{0ex}{}{10}{A_1}\right){\mathrm{\theta }}^{A_1}{\left(1-\mathrm{\theta }\right)}^{10-A_1}\right)$

$S≔\mathrm{Sample}\left(\mathrm{Binomial}\left(10\,0.4\right)\,1000\right)\:$

$\mathrm{LikelihoodRatioStatistic}\left(\mathrm{Binomial}\left(10\,\mathrm{\theta }\right)\,S\,\mathrm{bounds}=0..1\right)$

$-13401.4225072572-12140.\ln \left(1.-\mathrm{\theta }\right)-7860.\ln \left(\mathrm{\theta }\right)$

$\mathrm{assume}\left(0<\mathrm{\sigma }\right)$

$\mathrm{LikelihoodRatioStatistic}\left(\mathrm{Normal}\left(0\&comma;\mathrm{\sigma }\right)\&comma;A\&comma;\mathrm{samplesize}=2\right)$

$\frac{2\ln \left(2\right){\mathrm{\&sigma;~}}^2-2\ln \left(A_1^2+A_2^2\right){\mathrm{\&sigma;~}}^2+4\ln \left(\mathrm{\&sigma;~}\right){\mathrm{\&sigma;~}}^2-2{\mathrm{\&sigma;~}}^2+A_1^2+A_2^2}{{\mathrm{\&sigma;~}}^2}$