Доброе утро! Коллеги, кто-нибудь имел дело с функцией Surv , которая

Question

Доброе утро! Коллеги, кто-нибудь имел дело с функцией Surv , которая

используется для построения кривых Каплана-Майера (выживаемости).

Вопрос такой: Имею данные формата
hour alive lab
0 0 1_16 - CSE
12 5 1_16 - CSE
24 6 1_16 - CSE
36 9 1_16 - CSE
48 12 1_16 - CSE
60 25 1_16 - CSE
72 25 1_16 - CSE
84 30 1_16 - CSE

где первая колонка - это время, в котором проводилось измерение, вторая - число погибших экспериментальных животных, а третья - генотип.

Вопрос - как мне эти данные загнать в
Surv
- обьект? Если я делаю вот так:
Surv(time = tmp$hour, event = tmp$alive)

то получаю предупреждение, и обьект, забитый NA
В Surv(time = tmp$hour, event = tmp$alive) :
Invalid status value, converted to NA

Гугл не помог(
Сможете подсказать? Большое спасибо!

#programming #r #russian

0

20.07.2020

1 ответов

9 просмотров

Ivan Struzhkov · Accepted Answer

Ivan Struzhkov

в документации четко сказанно что грузить в event The status indicator, normally 0=alive, 1=dead. Other choices are T/F (TRUE = death) or 1/2 (2=death). For interval censored data, the status indicator is 0=right censored, 1= event at time, 2=left censored, 3=interval censored. https://www.rdocumentation.org/packages/survival/versions/2.11-4/topics/Surv

0

20.07.2020

169 похожих чатов

Доброе утро! Коллеги, кто-нибудь имел дело с функцией Surv , которая

1 ответов

Похожие вопросы