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没有列“时间”，“增量”

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

我正在尝试运行以下代码，但收到一条错误消息，提示未找到对象“ delta”。你知道是什么问题吗？

1 个答案: