r - 多変量 survreg によって生成された生存曲線をプロットする方法

Question

この質問は 、多変量モデルをプロットしたいという点で、ここの質問に続きます。

パッケージのリリースノートsurvregの例から生存関数をプロットしたいとしましょう

モデルの場所

survreg(Surv(time, status) ~ ph.ecog + age + strata(sex), lung)

私が理解していることから、次の問題は問題に対処しているように見えますが、データセットを提供せず、1 つの変数のみを使用します。この例では、2 つの連続変数と 1 つの因子があります。

Answer 1

肺年齢の平均と肺 $ph.ecog のモードを使用して、仮説上の個人の予測生存曲線をプロットします (ヘルプ ページの例に従います)。

?predict.survreg
pct <- 1:98/100
ptime <- predict(res, newdata=data.frame(ph.ecog=1, age=62.44737, sex=1),  type='quantile',
                  p=pct, se=TRUE)
 matplot(cbind(ptime$fit, ptime$fit + 2*ptime$se.fit,
                          ptime$fit - 2*ptime$se.fit)/30.5, 1-pct,
         xlab="Months", ylab="Survival", type='l', lty=c(1,2,2), col=1)