我正在尝试使用R运行多元序数GEE模型。因为我正在使用序数数据(依次为1、2、3),所以我考虑使用两个包:1)geepack中的'ordgee'和2)来自multgee的“ ordLORgee”。我的结果也与之相关,所以我使用了logit链接。
我的数据集包括40,000多个纵向观测值。如果我使用较小的数据集(即500、1000),则来自multgee和geepack的代码可以工作,但是当我尝试使用完整的数据集运行这些功能时,需要3个多小时才能收到一条错误消息,指出我的R连接是迷路了。
我想试着从multgee中找到
ordLORgee(output ~ day + pre2 + pre3 + pre4, id = person, link = 'logit', data = sampledf)
我想通过geepack适应:
ordgee(ordered(output) ~ day + pre2 + pre3 + pre4, id = person, data = final)
谁能告诉我这个问题的来源?理想情况下,我想使用multgee软件包,因为该软件包/ GEE方法相对较新,我无法从在线获得太多帮助。谢谢。