Calculate the modified GEE variance estimator proposed by Morel, Bokossa and Neerchal (2003).

1 | ```
GEE.var.mbn(formula,id,family=gaussian,data,corstr="independence",d=2,r=1)
``` |

`formula` |
an object of class "formula" (or one that can be coerced to that class): a formula expression as for other regression models to be fitted, of the form response ~ predictors. The details of formula specification can be seen in glm() and gee(). |

`id` |
a vector which identifies the clusters. The length of id should be the same as the total number of observations. Data is assumed to be sorted so that observations on a cluster are contiguous rows for all entities in the formula. |

`family` |
a family object: a list of functions and expressions for defining link and variance functions to be used in the model. This can be a character string naming a family function including "gaussian", "binomial", and "poisson". See family for details of family functions in glm() and gee(). |

`data` |
an optional data frame, list or environment (or object coercible by as.data.frame to a data frame) containing the variables in the model. If not found in data, the variables are taken from environment(formula), typically the environment from which glm is called. |

`corstr` |
a character string specifying working correlation structure: "independence", "AR-M","exchangeable", "unstructured" are possible. |

`d` |
adjustment factor in the variance estimator. The default value is 2. Details see Morel, Bokossa and Neerchal (2003). |

`r` |
adjustment factor in the variance estimator. The default value is 1. Details see Morel, Bokossa and Neerchal (2003). |

Fit the model by GEE and calculate the variance estimator proposed by Morel, Bokossa and Neerchal (2003).

`cov.beta ` |
estimate of the variance for beta |

`cov.var ` |
estimate of the variance-covariance matrix for the variance estimator |

Ming Wang <mwang@phs.psu.edu>

Maintainer: Zheng Li <zheng.li@outlook.com>

Morel JG, Bokossa MC, and Neerchal NK. Small sample correction for the variance of GEE estimators. Biometrical Journal 2003;45(4): 395-409.

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 | ```
data(dental)
data_alt <- reshape(dental, direction="long", timevar="Time",
varying=names(dental)[3:6], v.names="response", times=c(8,10,12,14))
data_alt <- data_alt[order(data_alt$subject),]
data_alt$gender <- as.numeric(data_alt$gender)
data_alt$Time <- sqrt(data_alt$Time)
formula <- response~Time+gender
mbn.ind <- GEE.var.mbn(formula,id="subject",family=gaussian,
data_alt,corstr="independence",d=2,r=1) ##Independence correlation structure;
mbn.exch <- GEE.var.mbn(formula,id="subject",family=gaussian,
data_alt,corstr="exchangeable",d=2,r=1) ##Exchangeable correlation structure;
mbn.ar1 <- GEE.var.mbn(formula,id="subject",family=gaussian,
data_alt,corstr="AR-M",d=2,r=1) ##AR-1 correlation structure;
mbn.unstr <- GEE.var.mbn(formula,id="subject",family=gaussian,
data_alt,corstr="unstructured",d=2,r=1) ##Unstructured correlation structur;
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.