[程式] SAS 在COX model中測試PH assumption
[軟體程式類別]:
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://www.ptt.cc/bbs/Statistics/M.1620668820.A.5F5.html
※ 編輯: pdono12345 (36.230.24.207 臺灣), 05/11/2021 01:48:28
→
05/11 08:05,
3年前
, 1F
05/11 08:05, 1F
→
05/11 08:05,
3年前
, 2F
05/11 08:05, 2F
→
05/11 08:40,
3年前
, 3F
05/11 08:40, 3F
→
05/11 08:40,
3年前
, 4F
05/11 08:40, 4F
→
05/11 21:05,
3年前
, 5F
05/11 21:05, 5F