[程式] SAS Cox model

看板Statistics作者 (國父)時間7年前 (2018/01/28 22:15), 7年前編輯推噓3(304)
留言7則, 4人參與, 7年前最新討論串1/1
[軟體程式類別]: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

01/30 07:38, 7年前 , 1F
Cox是Semi-parametric,無法給出特定時間的存活機率
01/30 07:38, 1F

01/30 13:31, 7年前 , 2F
樓上不要亂回答。
01/30 13:31, 2F

01/31 13:23, 7年前 , 3F
錯了,抱歉
01/31 13:23, 3F

01/31 22:22, 7年前 , 4F
如果是用R有辦法算出來嗎??
01/31 22:22, 4F

02/02 19:16, 7年前 , 5F
baseline out=資料檔名 survival=_all_
02/02 19:16, 5F

02/04 20:32, 7年前 , 6F
加baseline指令仍沒辦法讓codefile提供每人的時間點
02/04 20:32, 6F

02/04 20:33, 7年前 , 7F
估計存活機率值,感謝J大回覆
02/04 20:33, 7F
文章代碼(AID): #1QRThnzh (Statistics)