全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
4403 5
2017-05-25


一般cox,结局变量都是不利的事件,比如死亡时间,疾病复发时间,等等。如此作出的HR>1;

我做cox时,结局变量都是有益的事件(出院的时间),结果,cox出来的HR都<1,

这个讨论和解释起来本无可厚非,但不符合业界习惯!这样的表达也不被认可...>-< b (解释为该风险,比如吸烟,对出院的贡献为负)

请问各位老师,我该怎样修改代码,能把HR扭过来??使吸烟等风险>1  ??

谢谢老师!我的代码如下(因为是根据help自学的,(T-T)有错误一定要指出啊,谢谢老师!!)

. stset ex, failure( ex )                        //ex即结局变量:出院时间

     failure event:  ex != 0 & ex < .
obs. time interval:  (0, ex]
exit on or before:  failure

------------------------------------------------------------------------------
     1714  total observations
        0  exclusions
------------------------------------------------------------------------------
     1714  observations remaining, representing
     1714  failures in single-record/single-failure data
2845.103  total analysis time at risk and under observation
                                              at risk from t =         0
                                   earliest observed entry t =         0
                                        last observed exit t =        30

. stcox diameter3 blood2 age2 as


No. of subjects =         1114                     Number of obs   =      1114
No. of failures =         1714
Time at risk    =    2845.1029
                                                   LR chi2(11)     =    293.13
Log likelihood  =   -1153.947                     Prob > chi2     =    0.0000

------------------------------------------------------------------------------
          _t | Haz. Ratio   Std. Err.      z    P>|z|     [95% Conf. Interval]
-------------+----------------------------------------------------------------
   diameter3 |   .9261015   .0203681    -3.49   0.000     .8870289    .9668951
      blood2 |   .9861023   .0053189    -2.59   0.009     .9757323    .9965825
        age2 |   .9485863   .0184377    -2.72   0.007     .9131287    .9854207
         as  |   .5867303   .0992673    -3.15   0.002     .4211395    .8174309
------------------------------------------------------------------------------

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2017-5-25 12:51:32
help里的例子的结局变量也都是不利变量,比如肿瘤复发/死亡。。。我捯饬各种cox两天了。。。天天command栏一片红(T-T)
而且我也没找到stata里“删失”这一(0,1)变量该摆在代码的哪里。。。。?

谢谢各位前辈!!

(我ex取值1-30)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-5-25 13:03:39
另,我是想一起分析几个因素里,各个因素谁有影响,谁没影响。而不是help例子里分析某一个因素对结局的影响。

因为结局不止是0,1  ,还有时间轴,所以用cox  ,没用logistic。

(我一共10个变量,5个二分类变量,三个等级变量,2个连续变量)
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2017-5-25 16:43:03
贴错代码了,是. stset ex, failure( out ) ,但问题不变的,依旧HR<1..........
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2019-12-11 08:29:56
同问,遇到了同样的问题
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

2019-12-11 08:30:30
alasaa 发表于 2017-5-25 16:43
贴错代码了,是. stset ex, failure( out ) ,但问题不变的,依旧HR
请问解决了吗?
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群