This is a brief introduction to the theory of generalized linear models.
In generalized linear models, the response is assumed to possess a probability distribution of the exponential form. That is, the probability density of the response Y for continuous response variables, or the probability function for discrete responses, can be expressed as
for some functions a, b, and c that determine the specific distribution. For fixed , this is a oneparameter exponential family of distributions. The functions a and c are such that and , where w is a known weight for each observation. A variable representing w in the input data set can be specified in the WEIGHT statement. If no WEIGHT statement is specified, for all observations.
Standard theory for this type of distribution gives expressions for the mean and variance of Y:
where the primes denote derivatives with respect to . If represents the mean of Y, then the variance expressed as a function of the mean is
where V is the variance function.
Probability distributions of the response Y in generalized linear models are usually parameterized in terms of the mean and dispersion parameter instead of the natural parameter . The probability distributions that are available in the GENMOD procedure are shown in the following list. The zeroinflated Poisson and zeroinflated negative binomial distributions are not generalized linear models. However, the zeroinflated distributions are included in PROC GENMOD since they are useful extensions of generalized linear models. See Long (1997) for a discussion of the zeroinflated Poisson and zeroinflated negative binomial distributions. The PROC GENMOD scale parameter and the variance of Y are also shown.
The negative binomial and the zeroinflated negative binomial distributions contain a parameter k, called the negative binomial dispersion parameter. This is not the same as the generalized linear model dispersion , but it is an additional distribution parameter that must be estimated or set to a fixed value.
For the binomial distribution, the response is the binomial proportion . The variance function is , and the binomial trials parameter n is regarded as a weight w.
The density function for the Tweedie distribution when is expressed in terms of the parameters of the compound Poisson distribution. For more information about this representation, see the section Tweedie Distribution For Generalized Linear Models. For , the Tweedie random variable has positive support and its density function can be expressed in terms of stable distributions as defined in Hougaard (1986).
If a weight variable is present, is replaced with , where w is the weight variable.
PROC GENMOD works with a scale parameter that is related to the exponential family dispersion parameter instead of working with itself. The scale parameters are related to the dispersion parameter as shown previously with the probability distribution definitions. Thus, the scale parameter output in the “Analysis of Parameter Estimates” table is related to the exponential family dispersion parameter. If you specify a constant scale parameter with the SCALE= option in the MODEL statement, it is also related to the exponential family dispersion parameter in the same way.
For distributions other than the zeroinflated Poisson or zeroinflated negative binomial, the mean of the response in the ith observation is related to a linear predictor through a monotonic differentiable link function g.
Here, is a fixed known vector of explanatory variables, and is a vector of unknown parameters.
There are two link functions and linear predictors associated with zeroinflated distributions: one for the zero inflation probability , and another for the mean parameter . See the section ZeroInflated Models for more details about zeroinflated distributions.
Loglikelihood functions for the distributions that are available in the procedure are parameterized in terms of the means and the dispersion parameter . Zeroinflated log likelihoods are parameterized in terms two parameters, and . The parameter is the zeroinflation probability, and is a function of the distribution mean. The relationship between the mean of the zeroinflated Poisson and zeroinflated negative binomial distributions and the parameter is defined in the section Response Probability Distributions. The term represents the response for the ith observation, and represents the known dispersion weight. The loglikelihood functions are of the form
where the sum is over the observations. The forms of the individual contributions
are shown in the following list; the parameterizations are expressed in terms of the mean and dispersion parameters.
For the discrete distributions (binomial, multinomial, negative binomial, and Poisson), the functions computed as the sum of the terms are not proper loglikelihood functions, since terms involving binomial coefficients or factorials of the observed counts are dropped from the computation of the log likelihood, and a dispersion parameter is included in the computation. Deletion of factorial terms and inclusion of a dispersion parameter do not affect parameter estimates or their estimated covariances for these distributions, and this is the function used in maximum likelihood estimation. The value of used in computing the reported loglikelihood function is either the final estimated value, or the fixed value, if the dispersion parameter is fixed. Even though it is not a proper loglikelihood function in all cases, the function computed as the sum of the terms is reported in the output as the log likelihood. The proper loglikelihood function is also computed as the sum of the terms in the following list, and it is reported as the full log likelihood in the output.
Normal:
Inverse Gaussian:
Gamma:
Negative binomial:
Poisson:
Binomial:
Multinomial (k categories):
Zeroinflated Poisson:
Zeroinflated negative binomial:
Tweedie:
The GENMOD procedure uses a ridgestabilized NewtonRaphson algorithm to maximize the loglikelihood function with respect to the regression parameters. By default, the procedure also produces maximum likelihood estimates of the scale parameter as defined in the section Response Probability Distributions for the normal, inverse Gaussian, negative binomial, and gamma distributions.
On the rth iteration, the algorithm updates the parameter vector with
where is the Hessian (second derivative) matrix, and is the gradient (first derivative) vector of the loglikelihood function, both evaluated at the current value of the parameter vector. That is,
and
In some cases, the scale parameter is estimated by maximum likelihood. In these cases, elements corresponding to the scale parameter are computed and included in and .
If is the linear predictor for observation i and g is the link function, then , so that is an estimate of the mean of the ith observation, obtained from an estimate of the parameter vector .
The gradient vector and Hessian matrix for the regression parameters are given by
where is the design matrix, is the transpose of the ith row of X, and V is the variance function. The matrix is diagonal with its ith diagonal element
where
The primes denote derivatives of g and V with respect to . The negative of is called the observed information matrix. The expected value of is a diagonal matrix with diagonal values . If you replace with , then the negative of is called the expected information matrix. is the weight matrix for the Fisher scoring method of fitting. Either or can be used in the update equation. The GENMOD procedure uses Fisher scoring for iterations up to the number specified by the SCORING option in the MODEL statement, and it uses the observed information matrix on additional iterations.
The estimated covariance matrix of the parameter estimator is given by
where is the Hessian matrix evaluated using the parameter estimates on the last iteration. Note that the dispersion parameter, whether estimated or specified, is incorporated into . Rows and columns corresponding to aliased parameters are not included in .
The correlation matrix is the normalized covariance matrix. That is, if is an element of , then the corresponding element of the correlation matrix is , where .
Two statistics that are helpful in assessing the goodness of fit of a given generalized linear model are the scaled deviance and Pearson’s chisquare statistic. For a fixed value of the dispersion parameter , the scaled deviance is defined to be twice the difference between the maximum achievable log likelihood and the log likelihood at the maximum likelihood estimates of the regression parameters.
Note that these statistics are not valid for GEE models.
If is the loglikelihood function expressed as a function of the predicted mean values and the vector of response values, then the scaled deviance is defined by
For specific distributions, this can be expressed as
where D is the deviance. The following table displays the deviance for each of the probability distributions available in PROC GENMOD. The deviance cannot be directly calculated for zeroinflated models. Twice the negative of the log likelihood is reported instead of the proper deviance for the zeroinflated Poisson and zeroinflated negative binomial.
Distribution 
Deviance 

Normal 

Poisson 

Binomial 

Gamma 

Inverse Gaussian 

Multinomial 

Negative binomial 

Zeroinflated Poisson 

Zeroinflated negative binomial 

In the binomial case, , where is a binomial count and is the binomial number of trials parameter.
In the multinomial case, refers to the observed number of occurrences of the jth category for the ith subpopulation defined by the AGGREGATE= variable, is the total number in the ith subpopulation, and is the category probability.
Pearson’s chisquare statistic is defined as
and the scaled Pearson’s chisquare is .
The scaled version of both of these statistics, under certain regularity conditions, has a limiting chisquare distribution, with degrees of freedom equal to the number of observations minus the number of parameters estimated. The scaled version can be used as an approximate guide to the goodness of fit of a given model. Use caution before applying these statistics to ensure that all the conditions for the asymptotic distributions hold. McCullagh and Nelder (1989) advise that differences in deviances for nested models can be better approximated by chisquare distributions than the deviances can themselves.
In cases where the dispersion parameter is not known, an estimate can be used to obtain an approximation to the scaled deviance and Pearson’s chisquare statistic. One strategy is to fit a model that contains a sufficient number of parameters so that all systematic variation is removed, estimate from this model, and then use this estimate in computing the scaled deviance of submodels. The deviance or Pearson’s chisquare divided by its degrees of freedom is sometimes used as an estimate of the dispersion parameter . For example, since the limiting chisquare distribution of the scaled deviance has degrees of freedom, where n is the number of observations and p is the number of parameters, equating to its mean and solving for yields . Similarly, an estimate of based on Pearson’s chisquare is . Alternatively, a maximum likelihood estimate of can be computed by the procedure, if desired. See the discussion in the section Type 1 Analysis for more about the estimation of the dispersion parameter.
The Akaike information criterion (AIC) is a measure of goodness of model fit that balances model fit against model simplicity. AIC has the form
where p is the number of parameters estimated in the model, and LL is the log likelihood evaluated at the value of the estimated parameters. An alternative form is the corrected AIC given by
where n is the total number of observations used.
The Bayesian information criterion (BIC) is a similar measure. BIC is defined by
See Akaike (1981, 1979) for details of AIC and BIC. See Simonoff (2003) for a discussion of using AIC, AICC, and BIC with generalized linear models. These criteria are useful in selecting among regression models, with smaller values representing better model fit. PROC GENMOD uses the full log likelihoods defined in the section LogLikelihood Functions, with all terms included, for computing all of the criteria.
There are several options available in PROC GENMOD for handling the exponential distribution dispersion parameter. The NOSCALE and SCALE options in the MODEL statement affect the way in which the dispersion parameter is treated. If you specify the SCALE=DEVIANCE option, the dispersion parameter is estimated by the deviance divided by its degrees of freedom. If you specify the SCALE=PEARSON option, the dispersion parameter is estimated by Pearson’s chisquare statistic divided by its degrees of freedom.
Otherwise, values of the SCALE and NOSCALE options and the resultant actions are displayed in the following table.
NOSCALE 
SCALE=value 
Action 

Present 
Present 
Scale fixed at value 
Present 
Not present 
Scale fixed at 1 
Not present 
Not present 
Scale estimated by ML 
Not present 
Present 
Scale estimated by ML, 
starting point at value 

Present (negative binomial) 
Not present 
k fixed at 0 
The meaning of the scale parameter displayed in the “Analysis Of Parameter Estimates” table is different for the gamma distribution than for the other distributions. The relation of the scale parameter as used by PROC GENMOD to the exponential family dispersion parameter is displayed in the following table. For the binomial and Poisson distributions, is the overdispersion parameter, as defined in the “Overdispersion” section, which follows.
Distribution 
Scale 

Normal 

Inverse Gaussian 

Gamma 

Binomial 

Poisson 

In the case of the negative binomial distribution, PROC GENMOD reports the “dispersion” parameter estimated by maximum likelihood. This is the negative binomial parameter k defined in the section Response Probability Distributions.
Overdispersion is a phenomenon that sometimes occurs in data that are modeled with the binomial or Poisson distributions. If the estimate of dispersion after fitting, as measured by the deviance or Pearson’s chisquare, divided by the degrees of freedom, is not near 1, then the data might be overdispersed if the dispersion estimate is greater than 1 or underdispersed if the dispersion estimate is less than 1. A simple way to model this situation is to allow the variance functions of these distributions to have a multiplicative overdispersion factor :
Binomial:
Poisson:
An alternative method to allow for overdispersion in the Poisson distribution is to fit a negative binomial distribution, where , instead of the Poisson. The parameter k can be estimated by maximum likelihood, thus allowing for overdispersion of a specific form. This is different from the multiplicative overdispersion factor , which can accommodate many forms of overdispersion.
The models are fit in the usual way, and the parameter estimates are not affected by the value of . The covariance matrix, however, is multiplied by , and the scaled deviance and log likelihoods used in likelihood ratio tests are divided by . The profile likelihood function used in computing confidence intervals is also divided by . If you specify a WEIGHT statement, is divided by the value of the WEIGHT variable for each observation. This has the effect of multiplying the contributions of the loglikelihood function, the gradient, and the Hessian by the value of the WEIGHT variable for each observation.
The SCALE= option in the MODEL statement enables you to specify a value of for the binomial and Poisson distributions. If you specify the SCALE=DEVIANCE option in the MODEL statement, the procedure uses the deviance divided by degrees of freedom as an estimate of , and all statistics are adjusted appropriately. You can use Pearson’s chisquare instead of the deviance by specifying the SCALE=PEARSON option.
The function obtained by dividing a loglikelihood function for the binomial or Poisson distribution by a dispersion parameter is not a legitimate loglikelihood function. It is an example of a quasilikelihood function. Most of the asymptotic theory for log likelihoods also applies to quasilikelihoods, which justifies computing standard errors and likelihood ratio statistics by using quasilikelihoods instead of proper log likelihoods. For details on quasilikelihood functions, see McCullagh and Nelder (1989, Chapter 9), McCullagh (1983); Hardin and Hilbe (2003).
Although the estimate of the dispersion parameter is often used to indicate overdispersion or underdispersion, this estimate might also indicate other problems such as an incorrectly specified model or outliers in the data. You should carefully assess whether this type of model is appropriate for your data.