作者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/cn.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