Я делаю модель Кокса. Данные содержат пропущенные значения, поэтому я использовал мышей для вменения пропущенных.
Сначала я сделал вменение:
df <- read.table("File.txt")
p <- quickpred(df, minup = 0.5)
imp <- mice(df, m=5, pred=p)
Затем я сделал свою модель Кокса:
Cox <- with(imp, coxph(Surv(time = Month, event = die) ~ var1 + var2 + var3))
И я объединил это:
u <- pool(Cox)
Coef <- u$pooled$estimate
Но чтобы продолжить мой анализ, мне нужна базовая оценка опасности.
Когда у нас есть данные без вменения, это просто, мы можем использовать:
baz <- basehaz(coxph object, centered = FALSE)
Но вот, если я попытаюсь:
baz <- basehaz(Cox, centered = FALSE)
Это не будет работать, как Кокс объект мира.
На самом деле:
baz <- basehaz(u, centered = FALSE)
не сработало. Что касается объектов списка пула u $, я не думаю, что есть базовая оценка опасности, и если это так, я действительно извиняюсь.
Я пытался найти, как это получить, но не нашел.
Мои вопросы просты:
1) Можем ли мы получить базовую оценку опасности с помощью функции мышей, и если да, то как?
2) Если мы не можем, есть ли решение для этого? Вы думаете, что используете:
Com <- complete(imp, "long")
а затем сделать модель Кокса на этом простом информационном кадре осуществимой и иметь смысл?
Спасибо!