2210 ! Chapter 38: The GLIMMIX Procedure Details: GLIMMIX Procedure Generalized Linear Models Theory A generalized linear model consists of the following: a linear predictor Á D x 0 a monotonic mapping between the mean of the data and the linear predictor a response distribution in the exponential family of distributions A density or mass function in this family can be written as f .y/ D exp y b. / C c.y; f . // for some functions b. / and c. /. The parameter  is called the natural (canonical) parameter. The parameter is a scale parameter, and it is not present in all exponential family distributions. See Table 38.15 for a list of distributions for which Á 1. In the case where observations are weighted, the scale parameter is replaced with =w in the preceding density (or mass function), where w is the weight associated with the observation y. The mean and variance of the data are related to the components of the density, EOEY D D b 0 . /, VarOEY D b 00 . /, where primes denote first and second derivatives. If you express  as a function of , the relationship is known as the natural link or the canonical link function. In other words, modeling data with a canonical link assumes that  D x 0 ; the effect contributions are additive on the canonical scale. The second derivative of b. /, expressed as a function of , is the variance func- tion of the generalized linear model, a. / D b 00 .Â. //. Note that because of this relationship, the distribution determines the variance function and the canonical link function. You cannot, however, proceed in the opposite direction. If you provide a user-specified variance function, the GLIMMIX procedure assumes that only the first two moments of the response distribution are known. The full distribution of the data is then unknown and maximum likelihood estimation is not possible. Instead, the GLIMMIX procedure then estimates parameters by quasi-likelihood. Maximum Likelihood The GLIMMIX procedure forms the log likelihoods of generalized linear models as L. ; I y/ D n X i D1 f i l. i ; I y i ; w i / where l. i ; I y i ; w i / is the log likelihood contribution of the i th observation with weight w i and f i is the value of the frequency variable. For the determination of w i and f i , see the WEIGHT and FREQ statements. The individual log likelihood contributions for the various distributions are as follows.