我应该将时不变变量放在 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 (将#修改为@）