作者mickey1231 (國父)
看板Statistics
標題[程式] SAS Cox model
時間Sun Jan 28 22:15:10 2018
[軟體程式類別]: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://webptt.com/m.aspx?n=bbs/Statistics/M.1517148913.A.F6B.html
※ 編輯: mickey1231 (49.159.162.155), 01/28/2018 22:21:25
1F：推 Rhomboid: Cox是Semi-parametric，無法給出特定時間的存活機率 01/30 07:38
2F：→ bmka: 樓上不要亂回答。 01/30 13:31
3F：推 Rhomboid: 錯了，抱歉 01/31 13:23
4F：→ mickey1231: 如果是用R有辦法算出來嗎？？ 01/31 22:22
5F：推 jasonfun44: baseline out=資料檔名 survival=_all_ 02/02 19:16
6F：→ mickey1231: 加baseline指令仍沒辦法讓codefile提供每人的時間點 02/04 20:32
7F：→ mickey1231: 估計存活機率值，感謝Ｊ大回覆 02/04 20:33