У меня есть модель Кокса:
res <- coxph(Surv(time,status)~drugName*strata(preventative),data=df)
время - целое число от 0 до N
статус - целое число 0 или 1
drugName - категориальная переменная, например, drugA, drugB и т. д.
профилактический - логическое значение ИСТИНА / ЛОЖЬ
От:
summary(res)
Я вижу результаты для:
drugA
drugB
drugC
drugA:strata(preventative)=TRUE
drugB:strata(preventative)=TRUE
drugC:strata(preventative)=TRUE
Однако я не могу получить результаты для страт, где профилактический = ЛОЖЬ.
Клинический эпидемиолог заверил меня, что это возможно, однако я использую R, а она использует SAS.
Помимо повторного расчета и изменения предупредительного термина ИСТИНА / ЛОЖЬ, как я могу сделать это в R?
РЕДАКТИРОВАТЬ # 1 Согласно ответу Юлия Вайнора.
Я сейчас попробовал:
res <- coxph(Surv(time,status)~drugName*strata(preventative),data=df)
и
res <- coxph(Surv(time,status)~drugName:strata(preventative),data=df)
Соответствующие результаты:
и (новое предложение)
Если новое предложение поможет просто очистить значения категориальной переменной профилактического , есть ли идеи, почему drugNamecitaloprame в первой модели имеет значение, но NA во второй модели?