作者pdono12345 (小孩念機械一定打斷他腿)
看板Statistics
標題[程式] SAS 在COX model中測試PH assumption
時間Tue May 11 01:46:48 2021
[軟體程式類別]:
SAS
[程式問題]:
在 COX model 中測試 PH assumption
[軟體熟悉度]:
新手
[問題敘述]:
想請問大神們:
我用 proc phreg 在跑 PH assumption 時，result 出來很多 variable 的 hazard ratio
都會出現 E 夾雜在數字中，例:HR = 8.61E18，而且有些 HR 都爆大(千,萬)，然後95% CI 其
中一邊都是 0，一邊都是"."。
我的 event 是死亡，老師要我檢查有沒有變項的 event = 0，但都沒有，我也有把某些 event
數太少的變相合併(但不知道每個 cell 至少要大於多少才是正確的)，合併之後在跑一次？
結果還是跟上述問題一樣
在這邊也想請問大大，上述問題會是跟我的sample size太小(組一=72,組二=51)有關嗎？
還是跟event數 數量分佈有關？
舉其中一個變項的event數為例：
https://imgur.com/8cGf7CI
另外，reault 也都會出現以下的 warning : Convergence was not attained in 25
iterations
[程式範例]:
proc phreg data=PHassumption nosummary;
class VAR1(ref=first) VAR2(ref=first) VAR3(ref=first) VAR4(ref=last)
VAR5(ref=first) VAR6(ref=first) VAR7(ref=first) VAR8(ref=first)
VAR9(ref=first) VAR10(ref=first) VAR11(ref=first)/param=ref;
model follow*death(0)= VAR1 VAR2 VAR3 VAR4 VAR5
VAR6 VAR7 VAR8 VAR9 VAR10
t_VAR1_2 t_VAR1_3 t_VAR2_1 t_VAR3_2
t_VAR4_2 t_VAR4_3 t_smokeST1 t_SOP41
/ rl;
t_VAR1_2=(VAR1=2)*logtime;
t_VAR1_3=(VAR1=3)*logtime;
t_VAR2_1=(VAR2=1)*logtime;
t_VAR3_2=(VAR3=2)*logtime;
t_VAR4_2=(VAR4=2016)*logtime;
t_VAR4_3=(VAR4=2015)*logtime;
t_VAR5_1=(VAR5=1)*logtime;
run;
先謝謝各位!
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 36.230.24.207 (臺灣)
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Statistics/M.1620668820.A.5F5.html
※ 編輯: pdono12345 (36.230.24.207 臺灣), 05/11/2021 01:48:28
1F：→ andrew43: 試試逐步加入變數看看。不過樣本數小和這麼多自變數不 05/11 08:05
2F：→ andrew43: OK 05/11 08:05
3F：→ pdono12345: A大的意思是用selection :stepwise的方式嗎，那樣本 05/11 08:40
4F：→ pdono12345: 數小，自變數有限放幾個比較適合？ 05/11 08:40
5F：→ andrew43: 樣本數難說，但你目前樣本數光分成2類我都覺得太少 05/11 21:05