作者tokyo291 (工口工口)
看板R_Language
标题[问题] R 使用不同gee指令会有不收敛的情况
时间Thu Feb 12 16:19:51 2015
[问题类型]:
使用不同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
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 140.116.52.130
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/R_Language/M.1423729194.A.D44.html
※ 编辑: tokyo291 (140.116.52.130), 02/12/2015 17:44:17
1F：→ andrew43: 我没看到指出不收敛。你可以贴出相关讯息吗？ 02/12 18:08
2F：→ andrew43: 我也看到不相同结果。只能猜是有参数预设不同…… 02/12 18:08