Я пытаюсь запустить приведенный ниже код, но я получаю сообщение об ошибке, в котором говорится, что объект «delta» не найден. Вы знаете, в чем проблема?

Question

Я пытаюсь запустить приведенный ниже код, но я получаю сообщение об ошибке, в котором говорится, что объект 'delta' не найден. Знаете ли вы, в чем проблема?

Вот код:

library(KMsurv)
data()
data(burn)
burn
options(max.print=2772)


library(survival)
KM.fit<-survfit(Surv(time,delta)~type,conf.type="none",data=burn)
summary(KM.fit)

Вот сообщение об ошибке:

 KM.fit<-survfit(Surv(time,delta)~type,conf.type="none",data=burn)

Ошибка в Surv ( время, дельта): объект 'дельта' не найден

Answer 1

У burn нет столбцов 'time', 'delta'

names(burn)
#[1] "Obs" "Z1"  "Z2"  "Z3"  "Z4"  "Z5"  "Z6"  "Z7"  "Z8"  "Z9"  "Z10" "Z11" "T1"  "D1"  "T2"  "D2"  "T3"  "D3" 

Если мы предоставим правильные имена, это должно сработать, т. Е. Рассмотреть возможность проверки lung набора данных

names(lung)
#[1] "inst"      "time"      "status"    "age"       "sex"       
#[6] "ph.ecog"   "ph.karno"  "pat.karno" "meal.cal"  "wt.loss"  

и теперь примените функцию

survfit(Surv(time, status) ~ ph.ecog, data=lung)
#Call: survfit(formula = Surv(time, status) ~ ph.ecog, data = lung)

#   1 observation deleted due to missingness 
#            n events median 0.95LCL 0.95UCL
#ph.ecog=0  63     37    394     348     574
#ph.ecog=1 113     82    306     268     429
#ph.ecog=2  50     44    199     156     288
#ph.ecog=3   1      1    118      NA      NA