At 04:26 8/04/2013, you wrote:
Hello,
I have the following r-codes for solving a quasilikelihood estimating
equation:
>library(geepack)
>fit<-geese(y~x1+x2+x3,jack=TRUE,id=id,scale.fix=TRUE,data=dat,mean.link =
"logit", corstr="independence")
Now my question is how can I calculate the confidence interval of the
parameters of the above model "fit"?
[[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.
It pays to look at the help page. usually with a regression there is
a summary method
From ?geese and the index page of geepack there is summary.geese
names(summary(m2))
[1] "mean" "correlation"
"scale" "call" "model" "control" "error"
"clusz"
> summary(m2)[[1]]
estimate san.se wald p
(Intercept) 1.3476092 0.1573571 73.3423807 0.00000000
x 0.1118360 0.1159304 0.9306116 0.33470405
trt -0.1068224 0.1936977 0.3041417 0.58129752
x:trt -0.3023841 0.1710601 3.1247883 0.07710989
geepack uses sandwich standard errors (san.se)
summary(m2t)[[2]]
estimate san.se wald p
alpha 0.5978407 0.08107156 54.37934 1.653122e-13
HTH
Duncan
Duncan Mackay
Department of Agronomy and Soil Science
University of New England
Armidale NSW 2351
Email: home: mac...@northnet.com.au
______________________________________________
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.
