survan.doc

完整的生存分析Fortran语言源代码。有5个大的功能模块。

***** LOG RANK TEST *****
Enter the set of groups to be included in the analysis:
(Terminate list with "END")
0 1 end

                     ===== Begin Commentary =====

    The specifications follow   the  same rule as  for  Kaplan   Meier
estimation above.

                      ===== End Commentary =====

***** Data Definitions *****

Survival variable: SURV    ,
Event variable:    STATUS  ,
Group variable:    GROUP   ,
   GROUP      0,
   GROUP      1,



          ***** Log rank test for equality of survival *****


                              Num       Num         Num
Group                       At risk  Obs. dead   Exp. dead      O/E


                                 19  17.000      12.238      1.3892
                                 21  19.000      23.762     0.79958







    Statistic      =  3.123
    Deg of freedom =  1.000
    Significance   = 0.7721E-01


    ***** Covariance Matrix *****

  7.2633     -7.2633

 -7.2633      7.2633



*** Log rank analysis complete ***

                     ===== Begin Commentary =====

      Option 10 performs a log rank analysis of the data.

                      ===== End Commentary =====

Press Return / Enter to continue:
?
Enter option number:
    (1) Print data
    (2) Define new survival variable
    (3) Define event for logistic regression model
    (4) Save current status to a file
    (5) Set/review/change group, strata, and compulsory variables
    (6) Define new print level for regressions
    (7) Compute frequencies
    (8) Compute condescriptive statistics
    (9) Run Kaplan-Meier estimation
   (10) Run log rank test
   (11) Run Gehan Breslow test
   (12) Run proportional hazards model estimation
   (13) Run logistic regression model estimation
   (14) End program
? 11
***** GEHAN BRESLOW TEST *****
Enter the set of groups to be included in the analysis:
(Terminate list with "END")
0 1 end

                     ===== Begin Commentary =====

    The  specifications follow    the same rule  as for   Kaplan Meier
estimation above.

                      ===== End Commentary =====






***** Data Definitions *****

Survival variable: SURV    ,
Event variable:    STATUS  ,
Group variable:    GROUP   ,
   GROUP      0,
   GROUP      1,



       ***** Gehan Breslow test for equality of survival *****


                              NUM       NUM         W
GROUP                       AT RISK  OBS. DEAD   SCORES


                                 19  17.000     -114.00
                                 21  19.000      114.00


    STATISTIC      =  2.651
    DEG OF FREEDOM =  1.000
    SIGNIFICANCE   = 0.1035


    ***** Covariance Matrix *****

  4902.2     -4902.2

 -4902.2      4902.2



*** Gehan Breslow analysis complete ***

                     ===== Begin Commentary =====

      Option 11 performs a Gehan Breslow analysis of the data.

                      ===== End Commentary =====






Press Return / Enter to continue:
?
Enter option number:
    (1) Print data
    (2) Define new survival variable
    (3) Define event for logistic regression model
    (4) Save current status to a file
    (5) Set/review/change group, strata, and compulsory variables
    (6) Define new print level for regressions
    (7) Compute frequencies
    (8) Compute condescriptive statistics
    (9) Run Kaplan-Meier estimation
   (10) Run log rank test
   (11) Run Gehan Breslow test
   (12) Run proportional hazards model estimation
   (13) Run logistic regression model estimation
   (14) End program
? 12
***** PROPORTIONAL HAZARDS MODEL ESTIMATION *****

Enter the variable name(s) of the independent variables to be
included in the regression. (Terminate list with 'END')
(A '?' causes the list of names to be displayed.)
group




          ***** Proportional Hazards Regression Model *****



DEPENDENT VARIABLE:      SURV    ,
Event variable:          STATUS  ,
Independent variable(s): GROUP   ,


Negative of the log likelihood with all
Regression coefficients equal to zero =  102.1583312986

Negative of the log likelihood =  100.7191290008

Test for the overall significance of the regression:

   -2.0 * LN(Likelihood ratio) =  2.878
   Degrees of freedom for approximate chi-square =     1
   Significance = 0.8978E-01


               Beta         Standard     Normal dev.      P-Values
Variable   Coefficient       Error       Beta/Std err     (1-Sided)


GROUP      -.595902       0.348404       -1.71038       0.435982E-01








***** Variance-Correlation Matrix *****


 0.12139



Want to compute the underlying survivor function? (Y/N)
? y

                     ===== Begin Commentary =====

    These estimates are for all covariate values held at zero.  Hence,
they  are more useful for processing  by a subsequent program than for
viewing interactively.

                      ===== End Commentary =====

                ***** Estimated Survivor Function TABLE *****


                                                          Survivor
                                           Hazard         Function
              Time          Alpha          Values          F(Time)


          0.00000                                      1.0000000
         42.00000      0.96753114     0.32468855E-01  0.96753114
         43.00000      0.96669030     0.33309698E-01  0.93530297
         56.00000      0.96580940     0.34190595E-01  0.90332441
         63.00000      0.96515286     0.34847140E-01  0.87184614
         64.00000      0.96418232     0.35817683E-01  0.84061863
         88.00000      0.92570353     0.74296474E-01  0.77816363
         90.00000      0.95987028     0.40129721E-01  0.74693614
         92.00000      0.95819253     0.41807473E-01  0.71570863
         98.00000      0.95679641     0.43203592E-01  0.68478744
         105.0000      0.95463663     0.45363367E-01  0.65372318
         106.0000      0.95297652     0.47023475E-01  0.62298284
         109.0000      0.95065624     0.49343765E-01  0.59224252
         113.0000      0.94809502     0.51904976E-01  0.56150219
         116.0000      0.94525343     0.54746568E-01  0.53076187
         120.0000      0.93852198     0.61478019E-01  0.49813168
         127.0000      0.93449485     0.65505147E-01  0.46550149
         130.0000      0.92990315     0.70096850E-01  0.43287131
         132.0000      0.85436708     0.14563292      0.36983099
         133.0000      0.69571728     0.30428272      0.25729781
         134.0000      0.89959425     0.10040575      0.23146363
         139.0000      0.89119482     0.10880518      0.20627919
         140.0000      0.88426095     0.11573905      0.18240463
         146.0000      0.85417336     0.14582664      0.15580518
         161.0000      0.83585900     0.16414100      0.13023116
         165.0000      0.81154805     0.18845195      0.10568884
         180.0000      0.58482063     0.41517937      0.61809016E-01
         196.0000      0.46536219     0.53463781      0.28763579E-01
         204.0000      0.52429223     0.47570777      0.15080521E-01
         223.0000      0.28426993     0.71573007      0.42869387E-02








*** Proportional hazards analysis complete ***

                     ===== Begin Commentary =====

      Option 12 performs a proportional hazards analysis on the data.

                      ===== End Commentary =====

Press Return / Enter to continue:
?
Enter option number:
    (1) Print data
    (2) Define new survival variable
    (3) Define event for logistic regression model
    (4) Save current status to a file
    (5) Set/review/change group, strata, and compulsory variables
    (6) Define new print level for regressions
    (7) Compute frequencies
    (8) Compute condescriptive statistics
    (9) Run Kaplan-Meier estimation
   (10) Run log rank test
   (11) Run Gehan Breslow test
   (12) Run proportional hazards model estimation
   (13) Run logistic regression model estimation
   (14) End program
? 13
***** LOGISTIC REGRESSION MODEL ESTIMATION *****

Enter the variable name(s) of the independent variables to be
included in the regression. (Terminate list with 'END')
(A '?' causes the list of names to be displayed.)
group

Do you want an intercept? (Y/N)
? y



               ***** Logistic regression model *****



DEPENDENT VARIABLE:      STATUS  ,
Independent variable(s): GROUP   ,


Negative of the log likelihood with
constant probability of event =  13.00331893566
     Constant =  2.197224577336

Negative of the log likelihood =  12.99775562100






Test for the overall significance of the regression:

   -2.0 * LN(Likelihood ratio) = 0.1113E-01
   Degrees of freedom for approximate chi-square =     1
   Significance = 0.9160


               Beta         Standard     Normal dev.      P-values
Variable   Coefficient       Error       Beta/Std err     (1-Sided)


GROUP      0.111226        1.05426       0.105502       0.457989

INT         2.14007       0.747545        2.86279       0.209963E-02




***** Variance-correlation matrix *****


  1.1115     -.70907

 -.70907     0.55882



*** Logistic regression analysis complete ***

                     ===== Begin Commentary =====

      Option 13 performs a logistic regression analysis on the data.

                      ===== End Commentary =====






Press Return / Enter to continue:
?
Enter option number:
    (1) Print data
    (2) Define new survival variable
    (3) Define event for logistic regression model
    (4) Save current status to a file
    (5) Set/review/change group, strata, and compulsory variables
    (6) Define new print level for regressions
    (7) Compute frequencies
    (8) Compute condescriptive statistics
    (9) Run Kaplan-Meier estimation
   (10) Run log rank test
   (11) Run Gehan Breslow test
   (12) Run proportional hazards model estimation
   (13) Run logistic regression model estimation
   (14) End program
? 14
*** PROGRAM RUN ENDED ***
 FORTRAN STOP

                     ===== Begin Commentary =====

     Option 14 exits SURVAN.

                      ===== End Commentary =====