怎样应用R中GEE包用广义估计方程处理纵向重复观测数据，新人求助

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

手头有一个数据集，每个ID有若干条观测（数量各自不同），每条观测我都按照时间排了顺序（times）。我查了这种数据由于不独立，不能用logistic做回归，得用广义估计方程GEE分析，但是根据R语言中geepack包的说明，始终看不太明白，求高手指点程序该怎么写。数据自变量有8个，因变量1个，都已经根据自身数据特点定义为二分类的因子。
下面是的程序：
library(geepack)
mydata <- mydata[order(mydata$id, my$times),]

mf <- as.formula(outcome~age+gender+nation+smoke+F1+F2+F3+F4)

gee1 <- geeglm(mf,data=mydata,id=id,family=binomial, corstr="unstructured")



Error during wrapup: NA/NaN/Inf in 'y'  
我的数据集里8个自变量和1个因变量都不存在缺失值。求指点怎么才能正确的跑出GEE结果，我的qq是170213693，在线等。。

ID	age	nation	gender	smoke	F1	F2	F3	F4	times	outcome
1	1	1	0	0	0	0	0	1	1	1
1	1	1	0	0	0	0	0	1	2	0
1	1	1	0	0	0	1	0	0	3	1
2	0	1	0	1	1	0	0	0	1	1
2	0	1	0	1	1	0	0	0	2	1
3	0	1	1	1	0	1	0	0	1	0
3	0	1	1	1	0	1	0	0	2	0
3	0	1	1	1	0	1	0	0	3	0
3	0	1	1	1	0	0	1	0	4	1
4	1	0	1	0	0	0	0	1	1	0
5	1	0	1	0	0	0	0	1	1	1
5	1	0	1	0	0	0	0	1	2	1
5	1	0	1	0	1	0	0	0	3	0
5	1	0	1	0	1	0	0	0	4	0
6	0	1	0	1	0	1	0	0	1	0
6	0	1	0	1	0	1	0	0	2	0
6	0	1	0	1	0	1	0	0	3	1
6	0	1	0	1	0	1	0	0	4	0
6	0	1	0	1	0	0	1	0	5	1
6	0	1	0	1	0	0	1	0	6	0

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：广义估计方程 重复观测 GEE logistic ogistic R语言 广义估计方程 纵向数据

更正mydata <- mydata[order(mydata$id,mydata$times),]

不要沉啊。。。。

麻烦问一下 您懂了吗？

一场蜕变 发表于 2017-9-13 15:48
麻烦问一下 您懂了吗？

不好意思。。没有

为什么我的R说无法找到geepack这个包呢？

Setehk 发表于 2020-5-28 16:39
为什么我的R说无法找到geepack这个包呢？

很可能网络不佳。

install.packages("geepack")
library(geepack)

??geepack