Page 330 Table 13.2 Kaplan-Meier estimates for data from Figure 13.2
data list list /id years dead. begin data. 1 2 1 2 6 1 3 2 0 4 3 0 5 1 1 end data. exe.
survival tables = years /status = dead(1) /interval=thru 1 by 1 thru 2 by 1 thru 6 by 4 thru 10 by 4.
This subfile contains: 5 observations
Life Table Survival Variable years
Number Number Number Number Cumul
Intrvl Entrng Wdrawn Exposd of Propn Propn Propn Proba-
Start this During to Termnl Termi- Sur- Surv bility Hazard
Time Intrvl Intrvl Risk Events nating viving at End Densty Rate
------ ------ ------ ------ ------ ------ ------ ------ ------ ------
.0 5.0 .0 5.0 .0 .0000 1.0000 1.0000 .0000 .0000
1.0 5.0 .0 5.0 1.0 .2000 .8000 .8000 .2000 .2222
2.0 4.0 2.0 3.0 1.0 .3333 .6667 .5333 .0667 .1000
6.0 1.0 .0 1.0 1.0 1.0000 .0000 .0000 .1333 .5000
The median survival time for these data is 6.25
SE of SE of
Intrvl Cumul Proba- SE of
Start Sur- bility Hazard
Time viving Densty Rate
------- ------ ------ ------
.0 .0000 .0000 .0000
1.0 .1789 .1789 .2208
2.0 .2483 .0564 .0980
6.0 .0000 .0621 .0000
Page 331 Figure 13.6 Kaplan-Meier estimates of the survival functions for lung cancer data
get file "c:pma5surv.sav". coxreg days /status = death(1) /method = enter staget /categorical = staget /plot = survival /pattern = by staget.
Page 336 Table 13.3 Percentage of deaths versus explanatory variables
get file 'c:pma5surv.sav'. crosstabs /tables = staget by death /tables = perfbl by death /tables = treat by death /tables = poinf by death /cells = row.
Page 337 Table 13.4 Log-linear model for lung cancer data: results
NOTE: We were unable to reproduce this table.
Page 339 Table 13.5 Cox’s model for lung cancer data: results
coxreg days /status = death(1) /method = enter staget perfbl poinf treat.
<some output omitted>
We have skipped the examples on pages 340-344 because they use a simulated data set.
Page 347 Figure 13.9 Computer-generated graph of log(-logS(t)) versus t for lung cancer data (A = large tumor, B = small tumor)
coxreg days /status = death(1) /method = enter staget perfbl poinf treat /categorical = staget perfbl poinf treat /plot = lml /pattern = perfbl(.156) poinf(.05) treat(.513) by staget.
