Question

我使用来自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")

现在，当我运行联合模型时，我收到以下错误：＆＃34;钳工中的错误（X，Y，strats，偏移，初始化，控制，权重=权重，：不适合Cox模型，0失败＆＃34;

我不明白这一点，因为我的数据帧在Event列中都有1和0，而且mixed / cox模型输出表明事件已经检测到111个事件。对于lme（obs = 1205和group = 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))

此外，我没有任何遗漏的价值观或其他什么，所以我不知道这里出了什么问题

Thnx：）

如何使用R

0 个答案: