Блог

Я использовал пакет опроса, чтобы нарисовать взвешенный график выживания Каплана Мейера, как это:

 library(survey) data(pbc, package="survival") pbc$randomized lt;- with(pbc, !is.na(trt) amp; trtgt;0) biasmodel lt;- glm(randomized~age*edema,data=pbc) pbc$randprob lt;- fitted(biasmodel)  dpbclt;-svydesign(id=~1, prob=~randprob, strata=~edema, data=subset(pbc,randomized)) s2 lt;-svykm(Surv(time,statusgt;0) ~ sex, design = dpbc)  svyjskm(s2, pval = T, table = T, design = dpbc)  

Теперь я хотел бы рассчитать показатели выживаемости в двух группах, аналогичные тем, которые я получил бы при использовании summary(x, times = c(1:5)) на объекте survfit. Кто-нибудь знает, как я могу извлечь эти значения?

Заранее большое спасибо!