Hi, List,
I am using the package 'gee' downloaded from CRAN several days ago. I
tried to
do a 2 degree freedom contrast but couldn't get any result. It must be so
obvious but I just couldn't get it. Any hint would be greatly appreciated.
Here
is a sample code:
set.seed(123)
mydf <-
data.frame(trm=rep(letters[1:3],each=9),animal=rep(1:9,each=3),res=sample(c(rep(0,12),rep(1,15))))
mygee <-
gee(res~trm,id=animal,data=mydf,family=binomial,corstr='unstructured',
contrasts=list(trm=matrix(c(1,0,-1,1,-1,0),nrow=3)))
I couldn't find the result of contrasts except this:
> mygee$contrasts
$trm
[,1] [,2]
a 1 1
b 0 -1
c -1 0
where can I extract the p-value of this contrast testing though?
Thank you in advance,
Auston
[[alternative HTML version deleted]]
