The GEE Procedure
Generalized Estimating Equations
The marginal model is commonly used in analyzing longitudinal data when the population-averaged 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 high-order 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.
Working Correlation Matrix
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 13 shows the working correlation structures that are supported by the GEE procedure and the estimators that are used to estimate the working correlations.
Table 13: Working Correlation Structures and Estimators
| Working Correlation Structure | Estimator |
|---|---|
| Fixed | |
 where  is the jkth element of a constant, user-specified correlation matrix  |
The working correlation is not estimated in this case. |
| Independent | |
 |
The working correlation is not estimated in this case. |
| m-dependent | |
 |
  |
| Exchangeable | |
 |
  |
| Unstructured | |
 |
 |
| Autoregressive AR(1) | |
  |
  |
Dispersion Parameter
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 Model-Based 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 model-based covariance matrix and standard errors.
Quasi-likelihood Information Criterion
The quasi-likelihood 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 quasi-likelihood under the independent working correlation assumption, evaluated with the parameter estimates under the working correlation of interest as
where the quasi-likelihood contribution of the jth observation in the ith cluster is defined in the section Quasi-likelihood 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 model-based 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.
Quasi-likelihood Functions
See McCullagh and Nelder (1989) and Hardin and Hilbe (2003) for discussions of quasi-likelihood functions. The contribution of observation j in cluster i to the quasi-likelihood 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 quasi-likelihood 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 quasi-likelihood for the negative binomial differs from that given in McCullagh and Nelder (1989). The definition used here allows the negative binomial quasi-likelihood to approach the Poisson as 
.
-
Normal:
-
Inverse Gaussian:
-
Gamma:
-
Negative binomial:
-
Poisson:
-
Binomial:
-
Multinomial (s categories):
