Let L {\displaystyle L} be the likelihood function which depends on a univariate parameter θ {\displaystyle \theta } and let x {\displaystyle x} be the data. The score U ( θ ) {\displaystyle U(\theta )} is defined as
The Fisher information is6
where ƒ is the probability density.
The statistic to test H 0 : θ = θ 0 {\displaystyle {\mathcal {H}}_{0}:\theta =\theta _{0}} is S ( θ 0 ) = U ( θ 0 ) 2 I ( θ 0 ) {\displaystyle S(\theta _{0})={\frac {U(\theta _{0})^{2}}{I(\theta _{0})}}}
which has an asymptotic distribution of χ 1 2 {\displaystyle \chi _{1}^{2}} , when H 0 {\displaystyle {\mathcal {H}}_{0}} is true. While asymptotically identical, calculating the LM statistic using the outer-gradient-product estimator of the Fisher information matrix can lead to bias in small samples.7
Note that some texts use an alternative notation, in which the statistic S ∗ ( θ ) = S ( θ ) {\displaystyle S^{*}(\theta )={\sqrt {S(\theta )}}} is tested against a normal distribution. This approach is equivalent and gives identical results.
where L {\displaystyle L} is the likelihood function, θ 0 {\displaystyle \theta _{0}} is the value of the parameter of interest under the null hypothesis, and C {\displaystyle C} is a constant set depending on the size of the test desired (i.e. the probability of rejecting H 0 {\displaystyle H_{0}} if H 0 {\displaystyle H_{0}} is true; see Type I error).
The score test is the most powerful test for small deviations from H 0 {\displaystyle H_{0}} . To see this, consider testing θ = θ 0 {\displaystyle \theta =\theta _{0}} versus θ = θ 0 + h {\displaystyle \theta =\theta _{0}+h} . By the Neyman–Pearson lemma, the most powerful test has the form
Taking the log of both sides yields
The score test follows making the substitution (by Taylor series expansion)
and identifying the C {\displaystyle C} above with log ( K ) {\displaystyle \log(K)} .
If the null hypothesis is true, the likelihood ratio test, the Wald test, and the Score test are asymptotically equivalent tests of hypotheses.89 When testing nested models, the statistics for each test then converge to a Chi-squared distribution with degrees of freedom equal to the difference in degrees of freedom in the two models. If the null hypothesis is not true, however, the statistics converge to a noncentral chi-squared distribution with possibly different noncentrality parameters.
A more general score test can be derived when there is more than one parameter. Suppose that θ ^ 0 {\displaystyle {\widehat {\theta }}_{0}} is the maximum likelihood estimate of θ {\displaystyle \theta } under the null hypothesis H 0 {\displaystyle H_{0}} while U {\displaystyle U} and I {\displaystyle I} are respectively, the score vector and the Fisher information matrix. Then
asymptotically under H 0 {\displaystyle H_{0}} , where k {\displaystyle k} is the number of constraints imposed by the null hypothesis and
and
This can be used to test H 0 {\displaystyle H_{0}} .
The actual formula for the test statistic depends on which estimator of the Fisher information matrix is being used.10
In many situations, the score statistic reduces to another commonly used statistic.11
In linear regression, the Lagrange multiplier test can be expressed as a function of the F-test.12
When the data follows a normal distribution, the score statistic is the same as the t statistic.
When the data consists of binary observations, the score statistic is the same as the chi-squared statistic in the Pearson's chi-squared test.
Rao, C. Radhakrishna (1948). "Large sample tests of statistical hypotheses concerning several parameters with applications to problems of estimation". Mathematical Proceedings of the Cambridge Philosophical Society. 44 (1): 50–57. Bibcode:1948PCPS...44...50R. doi:10.1017/S0305004100023987. /wiki/Mathematical_Proceedings_of_the_Cambridge_Philosophical_Society ↩
Silvey, S. D. (1959). "The Lagrangian Multiplier Test". Annals of Mathematical Statistics. 30 (2): 389–407. doi:10.1214/aoms/1177706259. JSTOR 2237089. https://doi.org/10.1214%2Faoms%2F1177706259 ↩
Breusch, T. S.; Pagan, A. R. (1980). "The Lagrange Multiplier Test and its Applications to Model Specification in Econometrics". Review of Economic Studies. 47 (1): 239–253. doi:10.2307/2297111. JSTOR 2297111. /wiki/Trevor_S._Breusch ↩
Fahrmeir, Ludwig; Kneib, Thomas; Lang, Stefan; Marx, Brian (2013). Regression : Models, Methods and Applications. Berlin: Springer. pp. 663–664. ISBN 978-3-642-34332-2. 978-3-642-34332-2 ↩
Kennedy, Peter (1998). A Guide to Econometrics (Fourth ed.). Cambridge: MIT Press. p. 68. ISBN 0-262-11235-3. 0-262-11235-3 ↩
Lehmann and Casella, eq. (2.5.16). ↩
Davidson, Russel; MacKinnon, James G. (1983). "Small sample properties of alternative forms of the Lagrange Multiplier test". Economics Letters. 12 (3–4): 269–275. doi:10.1016/0165-1765(83)90048-4. /wiki/Economics_Letters ↩
Engle, Robert F. (1983). "Wald, Likelihood Ratio, and Lagrange Multiplier Tests in Econometrics". In Intriligator, M. D.; Griliches, Z. (eds.). Handbook of Econometrics. Vol. II. Elsevier. pp. 796–801. ISBN 978-0-444-86185-6. 978-0-444-86185-6 ↩
Burzykowski, Andrzej Gałecki, Tomasz (2013). Linear mixed-effects models using R : a step-by-step approach. New York, NY: Springer. ISBN 978-1-4614-3899-1.{{cite book}}: CS1 maint: multiple names: authors list (link) 978-1-4614-3899-1 ↩
Taboga, Marco. "Lectures on Probability Theory and Mathematical Statistics". statlect.com. Retrieved 31 May 2022. https://www.statlect.com/fundamentals-of-statistics/score-test ↩
Cook, T. D.; DeMets, D. L., eds. (2007). Introduction to Statistical Methods for Clinical Trials. Chapman and Hall. pp. 296–297. ISBN 978-1-58488-027-1. 978-1-58488-027-1 ↩
Vandaele, Walter (1981). "Wald, likelihood ratio, and Lagrange multiplier tests as an F test". Economics Letters. 8 (4): 361–365. doi:10.1016/0165-1765(81)90026-4. /wiki/Economics_Letters ↩