PROC GLIMMIX: MODEL Statement :: SAS/STAT(R) 9.22 User's Guide

The GLIMMIX Procedure

MODEL Statement

MODEL response <(response-options)> = <fixed-effects> </ model-options> ;

MODEL events/trials = <fixed-effects> < / model-options> ;

The MODEL statement is required and names the dependent variable and the fixed effects. The fixed-effects determine the $\text{[math]}$ matrix of the model (see the section Notation for the Generalized Linear Mixed Model for details). The specification of effects is the same as in the GLM or MIXED procedure. In contrast to PROC GLM, you do not specify random effects in the MODEL statement. However, in contrast to PROC GLM and PROC MIXED, continuous variables on the left and right side of the MODEL statement can be computed through PROC GLIMMIX programming statements.

An intercept is included in the fixed-effects model by default. It can be removed with the NOINT option.

The dependent variable can be specified by using either the response syntax or the events/trials syntax. The events/trials syntax is specific to models for binomial data. A binomial( $\text{[math]}$ , $\text{[math]}$ ) variable is the sum of $\text{[math]}$ independent Bernoulli trials with event probability $\text{[math]}$ . Each Bernoulli trial results in either an event or a nonevent (with probability $\text{[math]}$ ). You use the events/trials syntax to indicate to the GLIMMIX procedure that the Bernoulli outcomes are grouped. The value of the second variable, trials, gives the number $\text{[math]}$ of Bernoulli trials. The value of the first variable, events, is the number of events out of $\text{[math]}$ . The values of both events and (trials $\text{[math]}$ events) must be nonnegative and the value of trials must be positive. Observations for which these conditions are not met are excluded from the analysis. If the events/trials syntax is used, the GLIMMIX procedure defaults to the binomial distribution. The response is then the events variable. The trials variable is accounted in model fitting as an additional weight. If you use the response syntax, the procedure defaults to the normal distribution.

There are two sets of options in the MODEL statement. The response-options determine how the GLIMMIX procedure models probabilities for binary and multinomial data. The model-options control other aspects of model formation and inference. Table 38.6 summarizes important response-options and model-options. These are subsequently discussed in detail in alphabetical order by option category.

Table 38.6 Summary of Important MODEL Statement Options
Option	Description
Response Variable Options
DESCENDING	reverses the order of response categories
EVENT=	specifies the event category in binary models
ORDER=	specifies the sort order for the response variable
REFERENCE=	specifies the reference category in generalized logit models
Model Building
DIST=	specifies the response distribution
LINK=	specifies the link function
NOINT	excludes fixed-effect intercept from model
OFFSET=	specifies the offset variable for linear predictor
Statistical Computations
ALPHA= $\text{[math]}$	determines the confidence level ( $\text{[math]}$ )
CHISQ	requests chi-square tests
DDF=	specifies the denominator degrees of freedom (list)
DDFM=	specifies the method for computing denominator degrees of freedom
HTYPE=	selects the type of hypothesis test
NOCENTER	suppresses centering and scaling of $\text{[math]}$ columns during the estimation phase
ZETA=	tunes sensitivity in computing Type III functions
Statistical Output
CL	displays confidence limits for fixed-effects parameter estimates
CORRB	displays the correlation matrix of fixed-effects parameter estimates
COVB	displays the covariance matrix of fixed-effects parameter estimates
COVBI	displays the inverse covariance matrix of fixed-effects parameter estimates
E, E1, E2, E3	displays the $\text{[math]}$ matrix coefficients
INTERCEPT	adds a row for the intercept to test tables
ODDSRATIO	displays odds ratios and confidence limits
SOLUTION	displays fixed-effects parameter estimates (and scale parameter in GLM models)
STDCOEF	displays standardized coefficients

Top of Page