作者ofine (ofine)
看板Statistics
標題Re: [問題]Hazard ratio與發生密度
時間Tue Oct 30 11:41:24 2012
各位，以下是我分析所採用的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
1F：→ bmka:所以你是跟Cox model給的 hazard ratio 來比.. = = 10/30 19:51
2F：→ bmka:而且Cox model還做stratification 10/30 19:57
3F：→ bmka:你觀察到的現象有好幾個可能的解釋, 不過先釐清幾個問題 10/30 20:03
4F：→ bmka:(1)把data放一小部分上來,看看有沒有設定上的問題 10/30 20:04
5F：→ bmka:(2)每個人的follow-up time 怎麼設, 發生密度的分母怎麼算的 10/30 20:06
6F：→ bmka:(3)不做stratification的時候還是觀察到一樣的現象嗎? 10/30 20:07
7F：→ bmka:(4)HR跟發生密度雖然方向相反,但是數值大小差很多嗎? 10/30 20:09
8F：→ bmka:(5)兩個group的Kaplan-Meier survival curves 有交叉嗎 10/30 20:11