[程式] SAS Cox model

[mickey1231 (國父)]

[軟體程式類別]:SAS

[程式問題]:Cox proportional hazard model

[軟體熟悉度]:熟悉

[問題敘述]:我用proc phreg建好模型，想把模型套到測試資料與驗證資料驗證模型
           用code file指令將模型係數寫出，套到驗證資料後只有給我一個XBEAT係數
           請問有辦法讓程式幫我跑出各時間點的存活機率值嗎？
           謝謝～


[程式範例]:

data train;
input time y x1 x2 ...;
cards;
3 0 1 2 ...
0 0 2 3 ...
1 1 4 5 ...
5 2 3 1 ...
2 1 3 4 ...
;
run;

proc phreg data=train;
model time*y(0)=x1--x100 /selection=stepwise;
code file="\\cox.sas";
run;

data test2;
set test;
%include "\\cox.sas";
run;







--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 49.159.162.155
※ 文章網址: https://www.ptt.cc/bbs/Statistics/M.1517148913.A.F6B.html
※ 編輯: mickey1231 (49.159.162.155), 01/28/2018 22:21:25
推 Rhomboid: Cox是Semi-parametric，無法給出特定時間的存活機率
→ bmka: 樓上不要亂回答。
推 Rhomboid: 錯了，抱歉
→ mickey1231: 如果是用R有辦法算出來嗎？？
推 jasonfun44: baseline out=資料檔名 survival=_all_
→ mickey1231: 加baseline指令仍沒辦法讓codefile提供每人的時間點
→ mickey1231: 估計存活機率值，感謝Ｊ大回覆