[问题类型]:
使用不同GEE的指令，得到的结果不同
[软件熟悉度]:
使用者(已经有用R 做过不少作品)
[问题叙述]:
我想使用GEE来分析一笔资料
一开始使用的是geese()的指令(也有用geeglm()结果一样)
分别执行1.无交互作用项 2.有交互作用项
这两个model
但是第2个模型的结果没问题，第1个的结果则为估计不收敛
于是我用SAS和R的另一个GEE分析的指令gee()去比较
结果都是正常，想请问造成此差异可能的原因
以及该如何调整
[程式范例]:
geese()的部分
geefit1<-geese(resp~group+time+I(time^2), family=binomial(link="logit"),
data=dat4, id=id, corstr = "ar1", scale.fix=F)
summary(geefit1)
geefit2<-geese(resp~group+time+I(time^2)+group*time+group*I(time^2),
family=binomial(link="logit"),
data=dat4, id=id, corstr = "ar1", scale.fix=F)
summary(geefit2)
SAS的部分
/*model 1*/
PROC GENMOD data=dat4 DESCEND;
class id group(ref="2") num(ref="1") / param=ref;
MODEL resp=group week week*week /Dist=bin link=logit ;
REPEATED SUBJECT=id /type=AR(1) ;
run;
/*model 2*/
PROC GENMOD data=dat4 DESCEND;
class id group(ref="2") num(ref="1") / param=ref;
MODEL resp=group week week*week group*week group*week*week /Dist=bin
link=logit;
REPEATED SUBJECT=id/type=AR(1);
run;
gee()的部分
a1<-gee(resp~group+time+I(time^2), id=id, data=dat4,
family=binomial(link="logit"), corstr="AR-M")
summary(a1)
a2<-gee(resp~group+time+I(time^2)+group*time+group*I(time^2), id=id,
data=dat4,
family=binomial(link="logit"), corstr="AR-M")
summary(a2)
原始资料(week和time是相同的变量)
http://tinyurl.com/p2yd6jk
[关键字]:
gee geese geeglm