各位前輩您好:
想請教如何外推存活率且為利用relative survival的方法。透過S(t)=S*(t)R(t)的關係式
S(t)為all-cause survival probability among cancerpatients
S*(t)為expected survival(可透過國家的生命表估計)
我有參照 Royston & Lambert的書 (Flexible parametric survival analysis usingSTATA: beyond the Cox Model)的Chapter 6.8 外推方式 and Chapter 8的relative survival
為了練習，我使用這本書的data "ew_breast_ch8.dta"，且不放covariates
步驟如下:
(1) fitting the PH model ,using "bhazard" option，語法如下:

use ew_breast_ch8, clear
stset survtime, failure(dead==1)  id(ident)
gen age = int(min(agediag + _t,99))
gen year = int(year(datediag) + _t)
merge m:1 sex region dep year age using popmort_uk, keepusing(rate) keep(match)
stpm2 , scale(hazard) df(1) bhazard(rate) nolog

// Predict at times 0(1)50，外推50年//
range timev 0 50 51
predict rs0_m1, survival zeros timevar(timev)

但我參考 The STATA journal(2009),9:2,p265-290"Further development of flexible parametric models for survivalanalysis"
提到，如果我在stpm2下bhazard(rate)的話，我在預測時(predict指令)，預測的是相對存活(relative survival)，不是存活率。
但，我要外推的是all-cause survival probability, 不是relative survival.
我跑完上面的語法後會出現一個欄位為timev從0到50(以一年為單位切割)，以及0到50年的relative survival的數據，所以我要透過這樣的關係式去外推存活率S(t)=S*(t)R(t)，但還少了0-50年每一年的S*(t)。
我有用"strs"的指令去估計expected survival,S*(t)，但是沒辦法外推，只能估計出追蹤時間內的。

所以想求助各位前輩 如何在STATA估計出0-50年的expected survival probability，以期能夠透過這公式S(t)=S*(t)R(t)外推存活率??
或其實根本不用特別算出expected survival，STATA中stpm2有其他指令可以轉換成survival probability??

謝謝

我是STATA初學者，為了跑relative survival 才開始學STATA，請多指教，感謝

主要是被要求要畫出類似下圖，Y軸是存活率:(出自Therese M-L Andersson,2013)