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]]

