survan.doc

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

the strata.  Strata variables  should follow the  same rules  as group
variables.  Currently, SURVAN allows the log rank test to be  run over
strata  and Kaplan-Meier estimation  can  be performed separately  for
groups  by  strata.   Although  there   is a  stratified  version   of
proportional hazards estimation, it is not yet included in Survan.

      Choice 3  allows  the  user to specify the compulsory variables.
If any compulsory variable has a missing value in proportional hazards
or  logistic  regression,   the  case  will  not  be  included in  the
estimation.  This feature allows differing models to be tested  on the
same cases.

      Choice 4  allows the user  to review the  variable  name list as
well  as    the variable definitions  for    the group, event, strata,
compulsory, and logistic regression event variables.

      Choice 5 returns the user to the main menu.






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

SURVAN environment:
    Variable list:
     SURV       STATUS     GROUP
     Number of variables =     3

    Survival variable    = SURV
    Event variable       = STATUS
    Group variable       = NOT DEFINED
    Strata variable      = NOT DEFINED
    Compulsory variable  = NOT DEFINED
    Logit event variable = NOT DEFINED

Want to:
    (1) enter/change group variable
    (2) enter/change strata variable
    (3) enter/change compulsory variables
    (4) view group, strata, compulsory variables
    (5) return to option list
? 1
Enter group variable name.
(A '?' causes the list of names to be displayed.)
group
Want to:
    (1) enter/change group variable
    (2) enter/change strata variable
    (3) enter/change compulsory variables
    (4) view group, strata, compulsory variables
    (5) return to option list
? 5
*** group,strata,compulsory variables set/reviewed/changed ***
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
? 6
***** PRINT LEVEL DEFINITION *****
Enter desired print level:
    (1) Print only log likelihood
    (2) Print log likelihood and coefficient estimates
    (3) Print log likelihood, coefficient estimates and
        variance-correlation matrix
? 3
*** Print level defined ***





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

      Option   6 allows the  user to  specify  how much information is
printed after a regression is performed. If this option is not invoked
before a regression is performed, then option 3 is in effect.

                      ===== 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
? 7
***** COMPUTE FREQUENCIES *****
Enter the variable name on which to compute frequencies.
(A '?' causes the list of names to be displayed.)
group



          ***** Frequency Table *****



  Variable: GROUP


                      Relative    Cum
             Absolute   Freq      Freq
  Value        Freq     (Pct)     (Pct)


 0.00000         19      47.5      47.5
 1.00000         21      52.5     100.0


Number of observations =         40
Missing observations =          0



*** Frequency table computed ***






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

      Menu option   7  computes the  frequency  table for  a specified
variable.

                      ===== 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
? 8
***** COMPUTE CONDESCRIPTIVE STATISTICS *****
Enter the variable name(s) on which to compute the condescriptive statistics.
(A '?' causes the list of names to be displayed.)
all



               ***** Condescriptive Statistics *****



  Variable: SURV


  Mean       127.950
  Variance   2150.82      Std dev   46.3769      Std err   1.18915
  Minimum    42.0000      Maximum   244.000      Range     202.000


  Number of observations =    40
  Missing observations   =     0



  Variable: STATUS


  Mean      0.900000
  Variance  0.923077E-01  Std dev  0.303822      Std err  0.779030E-02
  Minimum    0.00000      Maximum   1.00000      Range     1.00000







  Number of observations =    40
  Missing observations   =     0



  Variable: GROUP


  Mean      0.525000
  Variance  0.255769      Std dev  0.505736      Std err  0.129676E-01
  Minimum    0.00000      Maximum   1.00000      Range     1.00000


  Number of observations =    40
  Missing observations   =     0



*** Condescriptive statistics computed ***

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

      Menu  option 8 prints   the condescriptive   statistics  for the
variables specified.

                      ===== 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
? 9
***** KAPLAN-MEIER ESTIMATION *****
Enter the set of groups to be included in the analysis:
(Terminate list with "END")
0 1 end






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

    A subset of groups can be specified by entering the  values of the
group  variable as is done here.   ALL can also be  specified and is a
frequent choice.

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

Want to:
    (1) Analyze each group in the set, or
    (2) Analyze the set of groups as a whole
? 2

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

    Kaplan Meier estimation can be performed  on each group separately
or on the subset of groups considered as a single data set.

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

***** CONFIDENCE INTERVAL SPECIFICATIONS *****
Want to compute the confidence interval for the proportion
surviving? (Y/N)
? y
Enter the desired confidence level:
? .95
Want to compute the confidence intervals over a range of quantiles,
such as 25%, 50% and 75% where 50% is the median? (Y/N)
? n
Do you want the life table to be printed? (Y/N)
? y



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

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







                 ***** KAPLAN MEIER LIFE TABLE *****



    The confidence interval for the proportion surviving
    computed, with the confidence level = 0.9500



                                      Lower             Upper
              Num      Num     Num    Limit     Cum     Limit      SD
   Time     Censored   Dead  At risk  of CI     Surv    of CI     Surv

 0.0000           0       0      40   1.0000   1.0000   1.0000   0.0000
 42.000           0       1      40   0.8732   0.9750   0.9955   0.0247
 43.000           0       1      39   0.8370   0.9500   0.9860   0.0345
 56.000           0       1      38   0.8035   0.9250   0.9738   0.0416
 63.000           0       1      37   0.7717   0.9000   0.9599   0.0474
 64.000           0       1      36   0.7412   0.8750   0.9448   0.0523
 88.000           0       2      35   0.6853   0.8250   0.9108   0.0601
 90.000           0       1      33   0.6549   0.8000   0.8940   0.0632
 92.000           0       1      32   0.6275   0.7750   0.8756   0.0660
 98.000           0       1      31   0.6007   0.7500   0.8568   0.0685
 104.00           1       0      30   0.5981   0.7500   0.8581   0.0685
 105.00           0       1      29   0.5708   0.7241   0.8382   0.0708
 106.00           0       1      28   0.5441   0.6983   0.8178   0.0729
 109.00           0       1      27   0.5179   0.6724   0.7969   0.0746
 113.00           0       1      26   0.4920   0.6466   0.7755   0.0761
 116.00           1       1      25   0.4666   0.6207   0.7537   0.0773
 120.00           0       1      23   0.4375   0.5937   0.7330   0.0785
 127.00           0       1      22   0.4119   0.5667   0.7095   0.0795
 130.00           0       1      21   0.3867   0.5397   0.6856   0.0801
 132.00           0       2      20   0.3410   0.4858   0.6330   0.0807
 133.00           0       4      18   0.2542   0.3778   0.5197   0.0788
 134.00           0       1      14   0.2219   0.3508   0.5059   0.0776
 139.00           0       1      13   0.2001   0.3238   0.4784   0.0762
 140.00           0       1      12   0.1787   0.2969   0.4502   0.0745
 144.00           1       0      11   0.1746   0.2969   0.4572   0.0745
 146.00           0       1      10   0.1521   0.2672   0.4256   0.0727
 161.00           0       1       9   0.1303   0.2375   0.3930   0.0704
 165.00           0       1       8   0.1093   0.2078   0.3594   0.0676
 180.00           0       2       7   0.0738   0.1484   0.2760   0.0599
 196.00           0       2       5   0.0386   0.0891   0.1921   0.0485
 204.00           0       1       3   0.0204   0.0594   0.1606   0.0404
 223.00           0       1       2   0.0082   0.0297   0.1019   0.0291
 244.00           1       0       1   0.0053   0.0297   0.1504   0.0291



*** Kaplan-Meier analysis complete ***






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

      Option 9 performs a Kaplan-Meier 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
? 10