Take the 2-minute tour ×
Cross Validated is a question and answer site for people interested in statistics, machine learning, data analysis, data mining, and data visualization. It's 100% free, no registration required.

I am learning the generalized estimating equations (GEE) and the geepack R package. There are some questions that I am a little confused.

In a GEE-constructed model, we have $Var(Y_{it})=\phi_{it}\cdot V(\mu_{it})$, where $\phi$ is the scale parameter. We further decompose $Var(Y)$ into $V^{1/2}R(\alpha)V^{1/2}$ where $\alpha$ is the correlation parameter. Three link-models are specified in geepack, for $\mu,\phi,\alpha$, respectively. See this PDF file for details.

(1) In GEE1, can I say we only need to make sure that the mean structure is correctly specified, i.e. the link model $g(\mu)=X\beta$ is correct, while it doesn't matter whether the link models for $\phi$ and $\alpha$ are correctly specified?

(2) By default, the geese.fit function makes the scale value scale.value = 1.0 -- does it say $\phi=1.0$? It is understandable that the default alpha=NULL as people can specify different correlation structures and the program will assign appropriate alpha values accordingly. My question is: how often people try to explicitly model the scale parameter $\phi$?

(3) This question is closely related to (2) about the scale parameter $\phi$. Recall the variance function is $Var(Y_{it})=\phi_{it}\cdot V(\mu_{it})$. In the Gaussian case, we have $\phi=\sigma^2$ and $V(\mu_{it})=1$; in the binomial case, we have $\phi=1$ and $V(\mu_{it})=\mu_{it}(1-\mu_{it})$; in the Poisson case, we have $\phi=1$ and $V(\mu_{it})=\mu_{it}$. Can I say that, in the negative binomial case, $\phi=1$ and $V(\mu_{it})=\mu_{it}+\varphi\mu_{it}^2$? Here $\varphi$ is the NB2 dispersion parameter, NOT the scale parameter $\phi$ in GEE.

Thank you very much!

share|improve this question
add comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Browse other questions tagged or ask your own question.