The marginal model is commonly used in analyzing longitudinal data when the populationaveraged effect is of interest. To estimate the regression parameters in the marginal model, Liang and Zeger (1986) proposed the generalized estimating equations method, which is widely used.
Suppose , represent the jth response of the ith subject, which has a vector of covariates . There are measurements on subject i, and the maximum number of measurements per subject is T.
Suppose the responses of the ith subject be with corresponding means . For generalized linear models, the marginal mean of the response is related to a linear predictor through a link function , and the variance of depends on the mean through a variance function .
An estimate of the parameter in the marginal model can be obtained by solving the generalized estimating equations,
where is the working covariance matrix of .
Only the mean and the covariance of are required in the GEE method; a full specification of the joint distribution of the correlated responses is not needed. This is particularly convenient because the joint distribution for noncontinuous responses involves highorder associations and is complicated to specify. Moreover, the regression parameter estimates are consistent even when the working covariance is incorrectly specified. Because of these properties, the GEE method is popular in situations where the marginal effect is of interest and the responses are not continuous. However, the GEE approach can lead to biased estimates when missing responses depend on previous responses. The weighted GEE method, which is described in the section Weighted Generalized Estimating Equations under the MAR Assumption, can provide unbiased estimates.
Suppose is an "working" correlation matrix that is fully specified by the vector of parameters . The covariance matrix of is modeled as
where is an diagonal matrix whose jth diagonal element is and is an diagonal matrix whose jth diagonal is , where is a weight variable that is specified in the WEIGHT statement. If there is no WEIGHT statement, for all i and j. If is the true correlation matrix of , then is the true covariance matrix of .
In practice, the working correlation matrix is usually unknown and must be estimated. It is estimated in the iterative fitting process by using the current value of the parameter vector to compute appropriate functions of the Pearson residual:
If you specify the working correlation matrix as , which is the identity matrix, the GEE reduces to the independence estimating equation.
Table 43.11 shows the working correlation structures that are supported by the GEE procedure and the estimators that are used to estimate the working correlations.
Table 43.11: Working Correlation Structures and Estimators
Working Correlation Structure 
Estimator 

Fixed 

where is the jkth element of a constant, userspecified correlation matrix 
The working correlation is not estimated in this case. 
Independent 


The working correlation is not estimated in this case. 
mdependent 



Exchangeable 



Unstructured 



Autoregressive AR(1) 



The dispersion parameter is estimated by
where is the total number of measurements and p is the number of regression parameters.
The square root of is reported by PROC GEE as the scale parameter in the "Parameter Estimates for Response Model with ModelBased Standard Error" output table. If a fixed scale parameter is specified by using the NOSCALE option in the MODEL statement, then the fixed value is used in estimating the modelbased covariance matrix and standard errors.
The quasilikelihood information criterion (QIC) was developed by Pan (2001) as a modification of Akaike’s information criterion (AIC) to apply to models fit by the GEE approach.
Define the quasilikelihood under the independent working correlation assumption, evaluated with the parameter estimates under the working correlation of interest as
where the quasilikelihood contribution of the jth observation in the ith cluster is defined in the section Quasilikelihood Functions and are the parameter estimates that are obtained by using the GEE approach with the working correlation of interest R.
QIC is defined as
where is the robust covariance estimate and is the inverse of the modelbased covariance estimate under the independent working correlation assumption, evaluated at , which are the parameter estimates that are obtained by using the GEE approach with the working correlation of interest R.
PROC GEE also computes an approximation to , which is defined by Pan (2001) as
where p is the number of regression parameters.
Pan (2001) notes that QIC is appropriate for selecting regression models and working correlations, whereas is appropriate only for selecting regression models.
See McCullagh and Nelder (1989) and Hardin and Hilbe (2003) for discussions of quasilikelihood functions. The contribution of observation j in cluster i to the quasilikelihood function that is evaluated at the regression parameters is expressed by , where is defined in the following list. These definitions are used in the computation of the quasilikelihood information criteria (QIC) for goodness of fit of models that are fit by the GEE approach. The are prior weights, if any, that are specified in the WEIGHT or FREQ statement. Note that the definition of the quasilikelihood for the negative binomial differs from that given in McCullagh and Nelder (1989). The definition used here allows the negative binomial quasilikelihood to approach the Poisson as .
Normal:
Inverse Gaussian:
Gamma:
Negative binomial:
Poisson:
Binomial:
Multinomial (s categories):