Hi,
I'm running GEE using geepack.
I set corstr = "ar1" as below:
> m.ar <- geeglm(L ~ O + A,
+ data = firstgrouptxt, id = id,
+ family = binomial, corstr = "ar1")
> summary(m.ar)
Call:
geeglm(formula = L ~ O + A, family = binomial,
data = firstgrouptxt, id = id, corstr = "ar1")
Coefficients:
Estimate Std.err Wald Pr(>|W|)
(Intercept) -2.62516 0.21154 154.001 <2e-16 ***
ontask 0.00498 0.12143 0.002 0.9673
attachmentB 0.73216 0.35381 4.282 0.0385 *
attachmentC 0.25960 0.33579 0.598 0.4395
---
Signif. codes: 0 *** 0.001 ** 0.01 * 0.05 . 0.1 1
Estimated Scale Parameters:
Estimate Std.err
(Intercept) 1.277 0.3538
Correlation: Structure = ar1 Link = identity
Estimated Correlation Parameters:
Estimate Std.err
alpha 0.978 0.005725
Number of clusters: 49 Maximum cluster size: 533
Then, it shows that :
Correlation: Link = identity
Why is it not Link = logit?
Thank you,
Sachi
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
