The GLIMMIX Procedure

The Likelihood Ratio Test

The likelihood ratio test (LRT) compares the likelihoods of two models where parameter estimates are obtained in two parameter spaces, the space normal upper Omega and the restricted subspace normal upper Omega 0. In the GLIMMIX procedure, the full model defines normal upper Omega and the test-specification in the COVTEST statement determines the null parameter space normal upper Omega 0. The likelihood ratio procedure consists of the following steps (see, for example, Bickel and Doksum 1977, p. 210):

  1. Find the estimate ModifyingAbove bold-italic theta With caret of bold-italic theta element-of normal upper Omega. Compute the likelihood upper L left-parenthesis ModifyingAbove bold-italic theta With caret right-parenthesis.

  2. Find the estimate ModifyingAbove bold-italic theta With caret Subscript 0 of bold-italic theta element-of normal upper Omega 0. Compute the likelihood upper L left-parenthesis ModifyingAbove bold-italic theta With caret Subscript 0 Baseline right-parenthesis.

  3. Form the likelihood ratio

    lamda overbar equals StartFraction upper L left-parenthesis ModifyingAbove bold-italic theta With caret right-parenthesis Over upper L left-parenthesis ModifyingAbove bold-italic theta With caret Subscript 0 Baseline right-parenthesis EndFraction
  4. Find a function f left-parenthesis lamda overbar right-parenthesis that has a known distribution. f left-parenthesis dot right-parenthesis serves as the test statistic for the likelihood ratio test.

Please note the following regarding the implementation of these steps in the COVTEST statement of the GLIMMIX procedure.

  • The function f left-parenthesis dot right-parenthesis in step 4 is always taken to be

    lamda equals 2 log left-brace lamda overbar right-brace

    which is twice the difference between the log likelihoods for the full model and the model under the COVTEST restriction.

  • For METHOD=RSPL and METHOD=RMPL, the test statistic is based on the restricted likelihood.

  • For GLMMs involving pseudo-data, the test statistics are based on the pseudo-likelihood or the restricted pseudo-likelihood and are based on the final pseudo-data.

  • The parameter space normal upper Omega for the full model is typically not an unrestricted space. The GLIMMIX procedure imposes boundary constraints for variance components and scale parameters, for example. The specification of the subspace normal upper Omega 0 must be consistent with these full-model constraints; otherwise the test statistic lamda does not have the needed distribution. You can remove the boundary restrictions with the NOBOUND option in the PROC GLIMMIX statement or the NOBOUND option in the PARMS statement.

Last updated: March 08, 2022