我使用 R 中的 JM 包进行联合模型分析。因此我需要首先拟合线性混合模型和 Cox 生存模型:
lmefit <- lme(DEST ~ Week + FL + FT +
DAC + Age + Gender + MED, random = ~ 1|ID,
data = Long_Set, na.action = na.omit)
coxFit <- coxph(Surv(TT_Event, Event) ~ FT + MED + DAC,
data = Event_Set, x = TRUE)
我的联合模型如下所示:
jointFit <- jointModel(lmefit, coxFit, timeVar = "Week")
现在,当我运行联合模型时,出现以下错误:“fitter 错误(X,Y,strats,offset,init,control,weights = weights,:无法拟合 0 次失败的 Cox 模型”
我不明白这一点,因为我的数据框在 Event 列中既有 1 又有 0,而且混合/cox 模型输出表明事件已检测到 111 个事件。对于 lme(obs 数 = 1205 和 groups = 194)和 cox(obs 数 = 194 和 events = 111)。
我的纵向/lme 部分的数据框看起来像这样(此处并非所有变量都存在):
dfL <- data.frame(ID = c(1, 1, 1, 2, 2, 3), DEST = c(0.4, 1.8, 1.2, 3.2, 3.6, 2.8), Week = c(0, 4, 16, 4, 20, 8), Event = c(1, 1, 1, 0, 0, 1), TT_Event = c(16, 20, 8))
对于生存/考克斯部分看起来非常相似:
dfC <- data.frame(ID = c(1, 2, 3, 4, 5, 6), DEST = c(1.2, 3.6, 2.8, 2.4, 1.9, 3.4), Week = c(16, 20, 8, 36, 24, 32), Event = c(1, 0, 1, 1, 1, 0), TT_Event = c(16, 20, 8, 36, 24, 32))
此外,我没有任何缺失值或其他东西,所以我不知道这里出了什么问题
谢谢 :)