作者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