Как вывести ИСТИННЫЙ / ЛОЖНЫЙ результат страты из регрессии Кокса? - PullRequest
1 голос
/ 28 мая 2019

У меня есть модель Кокса:

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)

Соответствующие результаты:

enter image description here

и (новое предложение)

enter image description here

Если новое предложение поможет просто очистить значения категориальной переменной профилактического , есть ли идеи, почему drugNamecitaloprame в первой модели имеет значение, но NA во второй модели?

1 Ответ

1 голос
/ 28 мая 2019

Ваши текущие результаты уже включают то, что вы хотите: drugA соответствует FALSE, а drugA + drugA:strata(preventative)=TRUE соответствует true.Чтобы получить те же результаты более чистым способом, вместо этого вы можете запустить

coxph(Surv(time, status) ~ drugName:strata(preventative), data = df)
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...