使用Stata进行生存分析实例

海蓝无魂

. stset time,f(outcome)

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

——————————————————————————
      33  total obs.
       0  exclusions
——————————————————————————
      33  obs. remaining, representing
      18  failures in single record/single failure data
   11370  total analysis time at risk, at risk from t =         0
                            earliest observed entry t =         0
                                 last observed exit t =      1045

. sts test treat,logrank

        failure _d:  outcome
  analysis time _t:  time
Log-rank test for equality of survivor functions

     |   Events         Events
treat |  observed       expected
——+————————-
1     |        14           8.57
2     |         4           9.43
——+————————-
Total |        18          18.00

           chi2(1) =       6.71
           Pr>chi2 =     0.0096

. sts graph,by(treat)

        failure _d:  outcome
  analysis time _t:  time