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