Re: [問題]Hazard ratio與發生密度
各位，以下是我分析所採用的SAS語法:
PROC SURVIVAL data=vo.dead93wt filetype=sas design=WOR;
NEST STRATA CITYN;
TOTCNT CITYCNT _MINUS1_;
WEIGHT ewt;
SUBGROUP plant_diet1_M1;
LEVEL 2 ;
reflevel plant_diet1_M1=2 ;
EVENT death;
MODEL follow_t=plant_diet1_M1;
TEST ADJWALDF;
PRINT BETA="BETA" SEBETA="S.E." DEFT="DESIGN EFFECT"
T_BETA="T:BETA=0" P_BETA="P-VALUE" HR LOWHR UPHR
DF="DF" ADJWALDF="ADJUSTED WALD F" ADJWALDP="P-VALUE"
/SEBETAFMT=F8.5 DFFMT=F4.0 T_BETAFMT=F8.2 DEFTFMT=F6.2
HRFMT=F7.2 LOWHRFMT=F6.2 UPHRFMT=F6.2 ADJWALDFFMT=F8.2;
run;
因為原始資料為多段機率不等抽樣，故使用Sudaan加權並校正抽樣標準誤!!
我想詢問的問題在於，
(1)觀察組的發生密度(新發個案數/總追蹤人年*1000)比參考組小(大)時，
為何用上述語法所分析出的觀察組HR會大於1(小於1)。
(2)y大所述的f(x)<g(x)時，f(x)/(1-F(x))不一定<g(x)/(1-G(x))又是什麼意思?
(3)若不用Sudaan，跑 proc phreg data=vo.dead93;
model follow_t*death(0)=plant_dietaM/risklimits;run;
中plant_dietaM數值小就會被直接當成參考組嗎? 還是可以改?
謝謝!!!
謝謝!!
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 210.60.122.8
→
10/30 19:51, , 1F
10/30 19:51, 1F
→
10/30 19:57, , 2F
10/30 19:57, 2F
→
10/30 20:03, , 3F
10/30 20:03, 3F
→
10/30 20:04, , 4F
10/30 20:04, 4F
→
10/30 20:06, , 5F
10/30 20:06, 5F
→
10/30 20:07, , 6F
10/30 20:07, 6F
→
10/30 20:09, , 7F
10/30 20:09, 7F
→
10/30 20:11, , 8F
10/30 20:11, 8F
討論串 (同標題文章)
完整討論串 (本文為第 2 之 2 篇)：