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 $\mb {X}$ 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(n,$\pi $) variable is the sum of n independent Bernoulli trials with event probability $\pi $. Each Bernoulli trial results in either an event or a nonevent (with probability $1-\pi $). 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 n of Bernoulli trials. The value of the first variable, events, is the number of events out of n. The values of both events and (trialsevents) 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 43.10 summarizes the options available in the MODEL statement. These are subsequently discussed in detail in alphabetical order by option category.

Table 43.10: 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

OBSWEIGHT=

Specifies the weight variable for the observation level unit

Statistical Computations

ALPHA=$\alpha $

Determines the confidence level ($1-\alpha $)

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

LWEIGHT

Determines how weights are used

NOCENTER

Suppresses centering and scaling of $\mb {X}$ columns during the estimation phase

REFLINP

Specifies a value for the linear predictor

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 $\mb {L}$ 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