r 생존 요약이 표시되지 않음

Question

survival 패키지로 생성 된 survfit 개체의 제한된 평균 생존 시간과 표준 오차를 인쇄하려고합니다. 내 전체 데이터에서 다음 예제 발췌를 실행하면 rmean을 요약에 표시 할 수 없습니다 (내 전체 데이터에서 모든 1000 obs를 실행할 때 표시되지 않음).

require(survival) 

dat <- structure(list(treat = c("Treat", "Control", "Treat", "Treat", 
"Control", "Treat", "Treat", "Control", "Control", "Treat", "Treat", 
"Control", "Control", "Control", "Treat", "Control", "Treat", 
"Control", "Control", "Control", "Control", "Control", "Control", 
"Control", "Control", "Control", "Control", "Control", "Control", 
"Control", "Control", "Control", "Control", "Control", "Treat", 
"Control", "Control", "Control", "Treat", "Control", "Control", 
"Control", "Control", "Treat", "Control", "Treat", "Control", 
"Treat", "Treat", "Control"), days = c(40, 90, 90, 90, 
90, 90, 90, 90, 90, 59, 51, 68, 90, 90, 51, 2, 29, 64, 90, 90, 
90, 81, 35, 19, 55, 90, 90, 90, 90, 90, 90, 90, 40, 3, 38, 90, 
32, 90, 90, 67, 42, 90, 90, 90, 66, 90, 90, 90, 90, 90), status = c(1, 
0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 0, 1, 1, 1, 1, 0, 0, 0, 1, 
1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 1, 1, 1, 0, 1, 0, 0, 1, 1, 0, 0, 
0, 1, 0, 0, 0, 0, 0)), .Names = c("treat", "days", "status" 
), row.names = c(NA, 50L), class = "data.frame") 


t.Surv <- Surv(dat$days, dat$status==1) 
print(survfit(t.Surv~treat, data=dat), show.rmean=TRUE) 

# Call: survfit(formula = t.Surv ~ treat, data = dat) 

#    records n.max n.start events median 0.95LCL 0.95UCL 
# treat=Control  35 35  35  13  NA  81  NA 
# treat=Treat  15 15  15  6  NA  51  NA 

Answer 1

인수는 print.rmean하지 show.rmean입니다.

> print(survfit(t.Surv~treat, data=dat), print.rmean=TRUE) 
Call: survfit(formula = t.Surv ~ treat, data = dat) 

       records n.max n.start events *rmean *se(rmean) median 0.95LCL 0.95UCL 
treat=Control  35 35  35  13 73.0  4.51  NA  81  NA 
treat=Treat  15 15  15  6 71.9  5.96  NA  51  NA 
    * restricted mean with upper limit = 90