QIC and quasi-Likelihood for GEE

Calculate quasi-likelihood under the independence model criterion () for Generalized Estimating Equations.

QIC(object, ..., typeR = FALSE)
QICu(object, ..., typeR = FALSE)
quasiLik(object, ...)

Arguments

object

a fitted model object of class "gee", "geepack", "geem", "wgee", or "yags".

...

for and , optionally more fitted model objects.

typeR

logical, whether to calculate (R). (R) is based on quasi-likelihood of a working correlation \(R\) model. Defaults to FALSE, and (I) based on independence model is returned.

Value

If just one object is provided, returns a numeric value with the corresponding ; if more than one object are provided, returns a data.frame with rows corresponding to the objects and one column representing or .

Note

This implementation is based partly on (revised) code from packages yags (R-Forge) and ape.

References

Pan, W. 2001 Akaike's Information Criterion in Generalized Estimating Equations. Biometrics 57, 120–125

Hardin J. W., Hilbe, J. M. 2003 Generalized Estimating Equations. Chapman & Hall/CRC

Author

Kamil Bartoń

See also

Methods exist for geegee (package gee), geeglmgeepack (geepack), geemgeeM (geeM), wgee (wgeesel, the package's QIC.gee function is used), and yags (yags on R-Forge). There is also a QIC function in packages MESS and geepack, returning some extra information (such as CIC and QICc). yags and compar.gee from package ape both provide QIC values.

Examples