Как построить кривую выживаемости, созданную с помощью Survreg (выживаемость в упаковке R)?

Я пытаюсь подогнать и построить модель Вейбулла для данные о выживании. У данных есть только одна ковариата, когорта, которая работает с 2006 по 2010 год. Итак, есть ли идеи о том, что добавить к двум следующим строкам кода, чтобы построить кривую выживаемости когорты 2010 года?

library(survival)
s <- Surv(subSetCdm$dur,subSetCdm$event)
sWei <- survreg(s ~ cohort,dist='weibull',data=subSetCdm)

Достижение того же с Модель Cox PH довольно проста, со следующими строками. Проблема в том, что Survfit () не принимает объекты типа Survreg.

sCox <- coxph(s ~ cohort,data=subSetCdm)
cohort <- factor(c(2010),levels=2006:2010)
sfCox <- survfit(sCox,newdata=data.frame(cohort))
plot(sfCox,col='green')

Используя легкое данных (из пакета выживания), вот что я пытаюсь выполнить.

#create a Surv object
s <- with(lung,Surv(time,status))

#plot kaplan-meier estimate, per sex
fKM <- survfit(s ~ sex,data=lung)
plot(fKM)

#plot Cox PH survival curves, per sex
sCox <- coxph(s ~ as.factor(sex),data=lung)
lines(survfit(sCox,newdata=data.frame(sex=1)),col='green')
lines(survfit(sCox,newdata=data.frame(sex=2)),col='green')

#plot weibull survival curves, per sex, DOES NOT RUN
sWei <- survreg(s ~ as.factor(sex),dist='weibull',data=lung)
lines(survfit(sWei,newdata=data.frame(sex=1)),col='red')
lines(survfit(sWei,newdata=data.frame(sex=2)),col='red')