PROC GLIMMIX: Implied Variance Functions :: SAS/STAT(R) 9.22 User's Guide

The GLIMMIX Procedure

Implied Variance Functions

While link functions are not unique for each distribution (see Table 38.9 for the default link functions), the distribution does determine the variance function $\text{[math]}$ . This function expresses the variance of an observation as a function of the mean, apart from weights, frequencies, and additional scale parameters. The implied variance functions $\text{[math]}$ of the GLIMMIX procedure are shown in Table 38.15 for the supported distributions. For the binomial distribution, $\text{[math]}$ denotes the number of trials in the events/trials syntax. For the negative binomial distribution, $\text{[math]}$ denotes the scale parameter. The multiplicative scale parameter $\text{[math]}$ is not included for the other distributions. The last column of the table indicates whether $\text{[math]}$ has a value equal to $\text{[math]}$ for the particular distribution.

Table 38.15 Variance Functions in PROC GLIMMIX
		Variance function
DIST=	Distribution	$\text{[math]}$	$\text{[math]}$
BETA	beta	$\text{[math]}$	No
BINARY	binary	$\text{[math]}$	Yes
BINOMIAL $\text{[math]}$ \| $\text{[math]}$ BIN $\text{[math]}$ \| $\text{[math]}$ B	binomial	$\text{[math]}$	Yes
EXPONENTIAL $\text{[math]}$ \| $\text{[math]}$ EXPO	exponential	$\text{[math]}$	Yes
GAMMA $\text{[math]}$ \| $\text{[math]}$ GAM	gamma	$\text{[math]}$	No
GAUSSIAN $\text{[math]}$ \| $\text{[math]}$ G $\text{[math]}$ \| $\text{[math]}$ NORMAL $\text{[math]}$ \| $\text{[math]}$ N	normal	$\text{[math]}$	No
GEOMETRIC $\text{[math]}$ \| $\text{[math]}$ GEOM	geometric	$\text{[math]}$	Yes
INVGAUSS $\text{[math]}$ \| $\text{[math]}$ IGAUSSIAN $\text{[math]}$ \| $\text{[math]}$ IG	inverse gaussian	$\text{[math]}$	No
LOGNORMAL $\text{[math]}$ \| $\text{[math]}$ LOGN	lognormal	$\text{[math]}$	No
NEGBINOMIAL $\text{[math]}$ \| $\text{[math]}$ NEGBIN $\text{[math]}$ \| $\text{[math]}$ NB	negative binomial	$\text{[math]}$	Yes
POISSON $\text{[math]}$ \| $\text{[math]}$ POI $\text{[math]}$ \| $\text{[math]}$ P	Poisson	$\text{[math]}$	Yes
TCENTRAL $\text{[math]}$ \| $\text{[math]}$ TDIST $\text{[math]}$ \| $\text{[math]}$ T	$\text{[math]}$	$\text{[math]}$	No

To change the variance function, you can use SAS programming statements and the predefined automatic variables, as outlined in the following section. Your definition of a variance function will override the DIST= option and its implied variance function. This has the following implication for parameter estimation with the GLIMMIX procedure. When a user-defined link is available, the distribution of the data is determined from the DIST= option, or the respective default for the type of response. In a GLM, for example, this enables maximum likelihood estimation. If a user-defined variance function is provided, the DIST= option is not honored and the distribution of the data is assumed unknown. In a GLM framework, only quasi-likelihood estimation is then available to estimate the model parameters.

Top of Page