Quasilikelihood estimation uses only the first and second moment of the response. In the case of independent data, this requires only a specification of the mean and variance of your data. The GLIMMIX procedure estimates parameters by quasilikelihood, if the following conditions are met:
The response distribution is unknown, because of a userspecified variance function.
There are no Gside random effects.
There are no Rside covariance structures or at most an overdispersion parameter.
Under some mild regularity conditions, the function

known as the log quasilikelihood of the ith observation, has some properties of a loglikelihood function (McCullagh and Nelder, 1989, p. 325). For example, the expected value of its derivative is zero, and the variance of its derivative equals the negative of the expected value of the second derivative. Consequently,

can serve as the score function for estimation. Quasilikelihood estimation takes as the gradient and “Hessian” matrix—with respect to the fixedeffects parameters —the quantities




In this expression, is a matrix of derivatives of with respect to the elements in , and is a diagonal matrix containing variance functions, . Notice that is not the second derivative matrix of . Rather, it is the negative of the expected value of . thus has the form of a “scoring Hessian.”
The GLIMMIX procedure fixes the scale parameter at 1.0 by default. To estimate the parameter, add the statement
random _residual_;
The resulting estimator (McCullagh and Nelder, 1989, p. 328) is

where if the NOREML option is in effect, m = f otherwise, and f is the sum of the frequencies.
See Example 41.4 for an application of quasilikelihood estimation with PROC GLIMMIX.