R version 4.2.2 Patched (2022-11-10 r83330) -- "Innocent and Trusting"
Copyright (C) 2022 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

R is free software and comes with ABSOLUTELY NO WARRANTY.
You are welcome to redistribute it under certain conditions.
Type 'license()' or 'licence()' for distribution details.

  Natural language support but running in an English locale

R is a collaborative project with many contributors.
Type 'contributors()' for more information and
'citation()' on how to cite R or R packages in publications.

Type 'demo()' for some demos, 'help()' for on-line help, or
'help.start()' for an HTML browser interface to help.
Type 'q()' to quit R.

> pdf("g01.pdf");options(width=64)
> #setwd("C:\\Users\\kolassa\\Class542")
> setwd("~/Taught1/960-542/Data")
> #*********************************************************/
> # The Stanford heart transplant data from                */
> # http://lib.stat.cmu.edu/datasets/stanford              */
> # represents a prototypical survival analysis data set.  */
> # Participants are waiting for a heart transplant.  Some */
> # of them get the transplant before death; others do not.*/
> # time is time in days to death or end of followup.      */
> # status is 1 if time is death, and 0 if end of followup.*/
> # wait is waiting time for transplant.                   */
> #*********************************************************/
> stanford<-read.table("stanford",skip=25,nrows=103,fill=TRUE,
+    col.names=c("id","start","age","status","time","prior",
+       "transp" ,"wait" ,"alleles" ,"antigen" ,"score"))
> #****************************************/
> # Exhibit complexities of survival data.*/
> #****************************************/
> #Print everything except the first column, with id
> print(stanford[,-1],row.names=FALSE)
 start age status time prior transp wait alleles antigen score
    68  53      1    1     0      0   NA      NA      NA    NA
    70  43      1    2     0      0   NA      NA      NA    NA
    71  52      1    2     0      0   NA      NA      NA    NA
    72  52      1    2     0      0   NA      NA      NA    NA
    68  54      1    3     0      0   NA      NA      NA    NA
    70  36      1    3     0      0   NA      NA      NA    NA
    71  47      1    3     0      0   NA      NA      NA    NA
    70  41      1    5     0      1    5       3       0  0.87
    73  47      1    5     0      0   NA      NA      NA    NA
    68  51      1    6     0      0   NA      NA      NA    NA
    67  39      1    6     0      0   NA      NA      NA    NA
    68  53      1    8     0      0   NA      NA      NA    NA
    71  56      1    9     0      0   NA      NA      NA    NA
    74  40      0   11     0      0   NA      NA      NA    NA
    70  43      1   12     0      0   NA      NA      NA    NA
    73  40      1   16     0      1    2       0       0  0.00
    69  54      1   16     0      0   NA      NA      NA    NA
    68  54      1   16     0      1    1       2       0  1.10
    72  29      1   17     0      1    5       1       0  0.69
    68  20      1   18     0      0   NA      NA      NA    NA
    72  41      1   21     0      0   NA      NA      NA    NA
    73  49      1   21     0      0   NA      NA      NA    NA
    69  55      1   28     0      1    1       3       1  2.76
    72  52      1   30     0      1    5       3       1  1.68
    74  49      0   31     0      0   NA      NA      NA    NA
    72  53      1   32     0      0   NA      NA      NA    NA
    69  50      1   35     0      0   NA      NA      NA    NA
    68  20      1   36     0      0   NA      NA      NA    NA
    68  59      1   37     0      0   NA      NA      NA    NA
    68  40      1   39     0      1   36       3       0    NA
    74  35      0   39     1      1   38       3       0    NA
    68  45      1   40     0      0   NA      NA      NA    NA
    70  42      1   40     0      0   NA      NA      NA    NA
    68  56      1   43     0      1   20       3       0    NA
    71  36      1   45     0      1    1       1       0    NA
    67  30      1   50     0      0   NA      NA      NA    NA
    69  42      1   51     0      1   12       3       0    NA
    70  50      1   53     0      1    2       0       0    NA
    68  42      1   58     0      1   12       2       1    NA
    71  52      1   61     0      1   10       2       1    NA
    70  61      1   66     0      1   19       3       0    NA
    72  45      1   68     0      1    3       3       0    NA
    71  49      1   68     0      1    3       3       0    NA
    71  39      1   69     0      0   NA      NA      NA    NA
    69  53      1   72     0      1   71       2       0    NA
    71  47      1   72     0      1   21       3       0    NA
    69  64      1   77     0      1   17       4       0    NA
    72  51      1   78     0      1   12       2       0    NA
    73  53      1   80     0      1   32       4       0    NA
    68  54      1   81     0      1   17       3       0    NA
    68  47      1   85     0      0   NA      NA      NA    NA
    72  56      1   90     0      1   27       3       1    NA
    72  53      1   96     0      1   67       2       0    NA
    70  48      1  100     0      1   46       2       0    NA
    71  41      1  102     0      0   NA      NA      NA    NA
    73  28      0  109     0      1   96       4       1    NA
    73  46      1  110     0      1   60       2       0    NA
    73  23      0  131     0      1   21       3       1    NA
    71  41      1  149     0      0   NA      NA      NA    NA
    68  47      1  153     0      1   26       1       0    NA
    73  43      1  165     1      1    4       3       0    NA
    73  26      0  180     0      1   13       2       0    NA
    73  52      1  186     1      1  160       3       1  0.82
    71  47      1  188     0      1   41       0       0  0.00
    73  51      1  207     0      1  139       4       1  1.33
    69  51      1  219     0      1   83       3       1  1.62
    69   8      1  263     0      0   NA      NA      NA    NA
    73  47      0  265     0      1   28       2       0  0.33
    71  48      1  285     0      1   32       4       1  1.08
    73  19      1  285     0      1   57       3       0  1.02
    68  49      1  308     0      1   28       2       0  1.12
    73  42      1  334     0      1   37       4       0  0.60
    73  47      1  340     0      0   NA      NA      NA    NA
    73  44      0  340     0      1  310       1       0  0.16
    71  47      1  342     1      1   21       2       1  1.82
    73  54      0  370     0      1   31       3       0  0.68
    73  48      0  397     0      1    8       3       1  1.44
    71  29      0  427     0      0   NA      NA      NA    NA
    73  52      0  445     0      1    6       4       1  1.94
    72  46      0  482     1      1   26       3       0  1.41
    72  48      0  515     0      1  210       3       0  0.81
    72  52      0  545     1      1   46       3       1  1.70
    72  48      1  583     1      1   32       1       0  0.12
    72  26      0  596     0      1    4       3       1  1.46
    72  47      0  630     0      1   31       3       0  0.97
    72  47      0  670     0      1   10       2       0  1.20
    68  50      1  675     0      1   51       4       0  1.32
    69  58      1  733     0      1    3       3       0  0.96
    71  32      0  841     0      1   27       3       1  1.93
    69  44      1  852     0      1   16       4       0  1.58
    71  41      0  915     0      1   78       2       0  0.19
    71  38      0  941     0      1   67       4       0  0.98
    71  45      1  979     1      1   83       0       0  0.00
    71  48      1  995     1      1    2       2       0  0.81
    69  43      1 1032     0      1    8       2       0  1.13
    71  36      0 1141     1      1   36       4       0  1.35
    70  45      0 1321     1      1   58       2       0  0.98
    68  53      1 1386     0      1   37       1       0  0.87
    69  30      0 1400     0      0   NA      NA      NA    NA
    70  48      0 1407     1      1   41       4       0  0.75
    69  40      0 1571     0      1   23       2       0  0.38
    69  48      0 1586     0      1   51       3       0  0.91
    69  33      0 1799     0      1   25       2       0  1.06
> plot(range(pretty(stanford$time)),c(0,.005),type="n",
+    xlab="Time (days)",ylab="Density",
+    main="Remaining lifetimes for Stanford heart patients")
> for(j in 0:1) hist(stanford$time[stanford$transp==j],
+    freq=FALSE,lty=j+1,add=TRUE)
> legend(400,.003,lty=c(2,1),
+    legend=c("Received Transplant","Did not Receive Transplant"))
> plot(range(pretty(log(stanford$time))),c(0,.5),type="n",
+     xlab="Log Time (days)",ylab="Density",
+     main="Remaining lifetimes for Stanford heart patients")
> for(j in 0:1) hist(log(stanford$time[stanford$transp==j]),
+    freq=FALSE,lty=j+1,add=TRUE)
> legend(0,.4,lty=c(2,1),
+    legend=c("Received Transplant","Did not Receive Transplant"))
>