1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
|
# gee support @(#) tests.q 1.2 96/09/27
library(gee)
source("testgee.dump")
myr <- gee(gsind~x2+x3+x4,id=id,data=testgee)
myr
summary(myr)
summary(gee(gs~x2+x4,id=id,data=testgee,corstr="exchangeable"))
summary(gee(gsind~sin(x2)+log(x3)+x4,id=id,data=testgee))
#summary(gee(gsind~gam(x2~x1)$fitted.values+x4,id=id,data=testgee,subset=id<20))
summary(gee(lgind~x2+x3+x4,id=id,data=testgee,family=binomial))
glm(lgind~x2+x3+x4,data=testgee, start=rep(0,4), # needed in R
family=quasi(link=logit,variance="mu(1-mu)"))
summary(gee(lgind~x2+x3+x4,id=id,data=testgee,family=binomial,scale.fix=TRUE))
#glm(lgind~x2+x3+x4,data=testgee,family=binomial)
summary(gee(lg~x2+x3+x4,id=id,data=testgee,family=binomial,corstr="exchangeable"))
summary(gee(lg~x2+x3+x4,id=id,data=testgee,family=binomial,corstr="stat_M_dep",M=1))
summary(gee(lg~x2+x3+x4,id=id,data=testgee,family=binomial(link=probit)))
# next line does not work in R's glm: needs starting values for quasi
#summary(gee(lg~x2+x3+x4,id=id,data=testgee,family=quasi(link=probit,variance="mu(1-mu)")))
summary(gee(lg~x2+x3+x4, id=id, data=testgee, b=rep(0,4), family=quasi(link=probit,variance="mu(1-mu)")))
logmat <- cbind(testgee[,"lg1"],testgee[,"nn"])
xmat <- cbind(testgee[,"x2"], testgee[,"x3"])
summary(gee(logmat~xmat,id = testgee[,"id"],family=binomial(link=probit)))
sub<-rep(c(FALSE,FALSE,FALSE,FALSE,TRUE),100)
summary(gee(lg~x2+x4,id=id,data=testgee,family=binomial,subset=sub))
R <- c(1.0,0.2,0.3,0.1,
0.2,1.0,0.2,0.3,
0.3,0.2,1.0,0.2,
0.1,0.3,0.2,1.0)
R <- matrix(R,nrow=4)
summary(gee(lg1~x2+x4,id=id,data=testgee,family=poisson(link=log),
subset=sub,corstr="fixed",R=R))
summary(gee(lg1~x2+x4,id=id,data=testgee,family=Gamma(link=log),
subset=sub,corstr="fixed",R=R))
|
