如何解决我应该将时不变变量放在 glmm.PQL 函数中的什么位置以执行时间相关的负二项式回归?
我是一名医生,最近为了应用负二项式回归,我处理和清理了一个数据。
首先,我尝试在 R 中使用包 MASS 的函数 glm.nb
,但在确保模型将实现数据对于一个独特的参与者(一组观察中可能存在的相关性)。
然后,我意识到我可以使用包 MASS 的 glmmPQL
或包 lme4 的 glmer
并在其中使用家庭负二项式家庭链接。
问题是我想知道我可以在模型的哪个部分嵌入偏移量(治疗天数的对数)以及我应该如何插入一个 id(例如性别和df 中的基线年龄)?
我最近的尝试是:
(glmmPQL (event ~ treatment + offset (log(person.time)),random= list (id=~1,gender=~1,baseline.age=~1),family= negative.binomial (theta=1.75),data=df ))
遇到与内存相关的错误(可能是因为代码错误)。 数据示例:
df<-data.frame(id=rep(1:3,each=4),treatment=sample(c(0,1),12,replace = T),event=sample(c(0,person.time=sample(c(15,31,30),age=rep(c(65,58,74),gender=rep(c("m","f","m"),each=4))
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。