Gosho's Criterion for Generalized Estimating Equations
Computes the Gosho's criterion (DEW) for one or more objects of the class glmgee.
DEW(...)
... |
one or several objects of the class glmgee which are obtained from the fit of generalized estimating equations. |
A matrix with the values of the DEW for each glmgee object provided in the input.
Gosho M., Hamada C., Yoshimura I. (2011) Criterion for the Selection of a Working Correlation Structure in the Generalized Estimating Equation Approach for Longitudinal Balanced Data. Communications in Statistics – Theory and Methods 40, 3839–3856.
Gosho M. (2014) Criteria to Select a Working Correlation Structure in SAS. Journal of Statistical Software, Code Snippets 57.
## Example 1
mod <- size ~ poly(days,4) + treat
fit1 <- glmgee(mod, id=tree, family=Gamma("log"), data=spruce, corstr="AR-1")
fit2 <- update(fit1, corstr="Exchangeable")
fit3 <- update(fit1, corstr="Stationary-M-dependent(2)")
fit4 <- update(fit1, corstr="Independence")
DEW(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="Exchangeable")
fit3 <- update(fit1, corstr="Non-Stationary-M-dependent(2)")
fit4 <- update(fit1, corstr="Independence")
DEW(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="Exchangeable")
fit3 <- update(fit1, corstr="Stationary-M-dependent(2)")
fit4 <- update(fit1, corstr="Independence")
DEW(fit1, fit2, fit3, fit4)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.