RJC: Rotnitzky–Jewell's Criterion for Generalized Estimating...

Description Usage Arguments Value References See Also Examples

View source: R/geeglm.R

Description

Computes the Rotnitzky–Jewell's criterion (RJC) for one or more objects of the class glmgee.

Usage

1
RJC(..., verbose = TRUE)

Arguments

...

one or several objects of the class glmgee which are obtained from the fit of generalized estimating equations.

verbose

an (optional) logical switch indicating if should the report of results be printed. By default, verbose is set to be TRUE.

Value

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

References

Hin L-Y, Carey V.J., Wang Y-G (2007) Criteria for Working–Correlation–Structure Selection in GEE: Assessment via Simulation. The American Statistician 61, 360–364.

See Also

QIC, CIC, GHYC, AGPC, SGPC

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
## Example 1
mod <- size ~ poly(days,4) + treat
fit1 <- glmgee(mod, id=tree, family=Gamma("log"), data=spruce, corstr="Exchangeable")
fit2 <- update(fit1, corstr="AR-1")
fit3 <- update(fit1, corstr="Stationary-M-dependent(2)")
fit4 <- update(fit1, corstr="Independence")
RJC(fit1, fit2, fit3, fit4)

## Example 2
mod <- dep ~ visit + group
fit1 <- glmgee(mod, id=subj, family=gaussian, corstr="Exchangeable", data=depression)
fit2 <- update(fit1, corstr="AR-1")
fit3 <- update(fit1, corstr="Non-Stationary-M-dependent(2)")
fit4 <- update(fit1, corstr="Independence")
RJC(fit1, fit2, fit3, fit4)

## Example 3
mod <- depressd ~ visit + group
fit1 <- glmgee(mod, id=subj, family=binomial, corstr="Exchangeable", data=depression)
fit2 <- update(fit1, corstr="AR-1")
fit3 <- update(fit1, corstr="Stationary-M-dependent(2)")
fit4 <- update(fit1, corstr="Independence")
RJC(fit1, fit2, fit3, fit4)

glmtoolbox documentation built on June 9, 2021, 9:07 a.m.

Related to RJC in glmtoolbox...