如何解决JM 的jointModel 抛出错误,左截断数据
我在依赖时间的 cox 模型中遇到以下错误。不确定是什么问题。它既不适用于教科书中的示例,也不适用于我的数据。
jointModel 中的错误(lmeFit.pbc,tdCox.pbc,timeVar = "year",method = "spline-PH-aGH") : 在 coxph() 中使用参数 'model = TRUE' 和 cluster()。
library("JM")
library("lcmm")
pbc <- pbc2[c("id","serBilir","drug","year","years","status2","spiders")]
pbc$start <- pbc$year
splitID <- split(pbc[c("start","years")],pbc$id)
pbc$stop <- unlist(lapply(splitID,function (d) c(d$start[-1],d$years[1]) ))
pbc$event <- with(pbc,ave(status2,id,FUN = function (x) c(rep(0,length(x)-1),x[1]) ))
lmeFit.pbc <- lme(log(serBilir) ~ drug * (year + I(year^2)),random = ~ year + I(year^2) | id,data = pbc2)
tdCox.pbc <- coxph(Surv(start,stop,event) ~ drug * spiders + cluster(id),data = pbc,x = TRUE,model = TRUE)
jointFit8.pbc <- jointModel(lmeFit.pbc,method = "spline-PH-aGH")
summary(jointFit8.pbc)
版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。