# RJC: Rotnitzky–Jewell's Criterion for Generalized Estimating... In glmtoolbox: Set of Tools to Data Analysis using Generalized Linear Models

## 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.

 ``` 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) ```