geeFit is basically a wrapper around glm.fit but also
computes the residuals and their derivatives with respect to the estimated coefficients.geeFit(y, x, link=c("identity","log","logit"), ...)"identity",
"log" or "logit".glm.fity and fitted.values.geeFit estimates coefficients $\theta$ using an estimating equation
$W{R-k^{-1}(\theta^T W)}=0$ where $k$ is the identity,
log or logit link function. glm.fit is used to estimate
the coefficients. The family object used in the glm.fit call depends
on the argument link. It will be gaussian(), quasipoisson() or
binomial() when link is "identity", "log" and "logit"
respectively.geeFit is intended to be used as help function in drFit, ebeFit
and obeFit.
glm.fit