The GLIMMIX Procedure
-
Overview
-
Getting Started
-
SyntaxPROC GLIMMIX StatementBY StatementCLASS StatementCODE StatementCONTRAST StatementCOVTEST StatementEFFECT StatementESTIMATE StatementFREQ StatementID StatementLSMEANS StatementLSMESTIMATE StatementMODEL StatementNLOPTIONS StatementOUTPUT StatementPARMS StatementRANDOM StatementSLICE StatementSTORE StatementWEIGHT StatementProgramming StatementsUser-Defined Link or Variance Function
-
DetailsGeneralized Linear Models TheoryGeneralized Linear Mixed Models TheoryGLM Mode or GLMM ModeStatistical Inference for Covariance ParametersDegrees of Freedom MethodsEmpirical Covariance ("Sandwich") EstimatorsExploring and Comparing Covariance MatricesProcessing by SubjectsRadial Smoothing Based on Mixed ModelsOdds and Odds Ratio EstimationParameterization of Generalized Linear Mixed ModelsResponse-Level Ordering and ReferencingComparing the GLIMMIX and MIXED ProceduresSingly or Doubly Iterative FittingDefault Estimation TechniquesDefault OutputNotes on Output StatisticsODS Table NamesODS Graphics
-
ExamplesBinomial Counts in Randomized BlocksMating Experiment with Crossed Random EffectsSmoothing Disease Rates; Standardized Mortality RatiosQuasi-likelihood Estimation for Proportions with Unknown DistributionJoint Modeling of Binary and Count DataRadial Smoothing of Repeated Measures DataIsotonic Contrasts for Ordered AlternativesAdjusted Covariance Matrices of Fixed EffectsTesting Equality of Covariance and Correlation MatricesMultiple Trends Correspond to Multiple Extrema in Profile LikelihoodsMaximum Likelihood in Proportional Odds Model with Random EffectsFitting a Marginal (GEE-Type) ModelResponse Surface Comparisons with Multiplicity AdjustmentsGeneralized Poisson Mixed Model for Overdispersed Count DataComparing Multiple B-SplinesDiallel Experiment with Multimember Random EffectsLinear Inference Based on Summary DataWeighted Multilevel Model for Survey Data
- References
The "Fit Statistics" table provides statistics about the estimated model. The first entry of the table corresponds to the negative of twice the (possibly restricted) log likelihood, log pseudo-likelihood, or log quasi-likelihood. If the estimation method permits the true log likelihood or residual log likelihood, the description of the first entry reads accordingly. Otherwise, the fit statistics are preceded by the words Pseudo- or Quasi-, for Pseudo- and Quasi-Likelihood estimation, respectively.
Note that the (residual) log pseudo-likelihood in a GLMM is the (residual) log likelihood of a linearized model. You should not compare these values across different statistical models, even if the models are nested with respect to fixed and/or G-side random effects. It is possible that between two nested models the larger model has a smaller pseudo-likelihood. For this reason, IC=NONE is the default for GLMMs fit by pseudo-likelihood methods.
See the IC= option of the PROC GLIMMIX statement and Table 44.2 for the definition and computation of the information criteria reported in the "Fit Statistics" table.
For generalized linear models, the GLIMMIX procedure reports Pearson’s chi-square statistic
![\[ X^2 = \sum _{i} \frac{w_ i (y_ i - \widehat{\mu }_ i)^2}{a(\widehat{\mu }_ i)} \]](images/statug_glimmix0914.png){kind=small}
where 
is the variance function evaluated at the estimated mean.
For GLMMs, the procedure typically reports a generalized chi-square statistic,
![\[ X_ g^2 = \widehat{\mb{r}}’\mb{V}(\widehat{\btheta }^*)^{-1} \widehat{\mb{r}} \]](images/statug_glimmix0916.png){kind=small}
so that the ratio of 
or 
and the degrees of freedom produces the usual residual dispersion estimate.
If the R-side scale parameter 
is not extracted from 
, the GLIMMIX procedure computes
![\[ X_ g^2 = \widehat{\mb{r}}’\mb{V}(\widehat{\btheta })^{-1} \widehat{\mb{r}} \]](images/statug_glimmix0919.png){kind=small}
as the generalized chi-square statistic. This is the case, for example, if R-side covariance structures are varied by a GROUP= effect or if the scale parameter is not profiled for an R-side TYPE=CS , TYPE=SP , TYPE=AR , TYPE=TOEP , or TYPE=ARMA covariance structure.
For METHOD= LAPLACE , the generalized chi-square statistic is not reported. Instead, the Pearson statistic for the conditional distribution appears in the "Conditional Fit Statistics" table.
If your model contains smooth components (such as TYPE=RSMOOTH ), then the "Fit Statistics" table also displays the residual degrees of freedom of the smoother. These degrees of freedom are computed as
![\[ \mi{df}_{\mi{smooth},\mi{res}} = f - \mr{trace}(\mb{S}) \]](images/statug_glimmix0920.png){kind=small}
where 
is the "smoother" matrix—that is, the matrix that produces the predicted values on the linked scale.
The ODS name of the "Fit Statistics" table is FitStatistics.