PAC: Pardo-Alonso's Criterion for Generalized Estimating Equations

In glmtoolbox: Set of Tools to Data Analysis using Generalized Linear Models

Pardo-Alonso's Criterion for Generalized Estimating Equations

Description

Computes the Pardo-Alonso's criterion (PAC) for one or more objects of the class glmgee.

Usage

PAC(..., verbose = TRUE, digits = max(3, getOption("digits") - 2))

Arguments

...	one or several objects of the class glmgee.
verbose	an (optional) logical switch indicating if should the report of results be printed. As default, verbose is set to TRUE.
digits	an (optional) integer indicating the number of digits to print. As default, digits is set to max(3, getOption("digits") - 2).

Value

A data.frame with the values of the PAC for each glmgee object in the input.

References

Pardo M.C., Alonso R. (2019) Working correlation structure selection in GEE analysis. Statistical Papers 60:1447–1467.

Vanegas L.H., Rondon L.M., Paula G.A. (2023) Generalized Estimating Equations using the new R package glmtoolbox. The R Journal 15:105-133.

See Also

QIC, CIC, RJC, AGPC, SGPC, GHYC

Examples

###### Example 1: Effect of ozone-enriched atmosphere on growth of sitka spruces
data(spruces)
mod1 <- size ~ poly(days,4) + treat
fit1 <- glmgee(mod1, id=tree, family=Gamma(log), data=spruces)
fit2 <- update(fit1, corstr="AR-M-dependent")
fit3 <- update(fit1, corstr="Stationary-M-dependent(2)")
fit4 <- update(fit1, corstr="Exchangeable")
PAC(fit1, fit2, fit3, fit4)

###### Example 2: Treatment for severe postnatal depression
data(depression)
mod2 <- depressd ~ visit + group
fit1 <- glmgee(mod2, id=subj, family=binomial(logit), data=depression)
fit2 <- update(fit1, corstr="AR-M-dependent")
fit3 <- update(fit1, corstr="Stationary-M-dependent(2)")
fit4 <- update(fit1, corstr="Exchangeable")
PAC(fit1, fit2, fit3, fit4)

###### Example 3: Treatment for severe postnatal depression (2)
mod3 <- dep ~ visit*group
fit1 <- glmgee(mod3, id=subj, family=gaussian(identity), data=depression)
fit2 <- update(fit1, corstr="AR-M-dependent")
fit3 <- update(fit1, corstr="Exchangeable")
PAC(fit1, fit2, fit3)

glmtoolbox documentation built on Sept. 11, 2024, 7:32 p.m.