#
# Survival in a randomised trial comparing two treatments for ovarian cancer.
# Dataset included with Therneau's survival5 library.
# futime (f-u time) fustat (died) age (years) resid.ds (residual disease)
# rx (treatment arm) ecog.ps (ECOG performance status)
#
# Using the O'Brien test for effects of covariates:
#           coef exp(coef) se(coef)      z      p
# age      0.551     1.735    0.170  3.243 0.0012
# rx2     -0.937     0.392    0.619 -1.514 0.1300
# ecog.ps  0.308     1.361    0.610  0.506 0.6100
#
# And the equivalent Weibull regression 
# survreg(formula = Surv(futime, fustat) ~ age + rx + ecog.ps, 
#         data = ovarian, dist = "weibull", scale = 1)
#               Value Std. Error      z        p
# (Intercept) 12.2157     2.4690  4.948 7.51e-07
# age         -0.1044     0.0322 -3.242 1.19e-03
# rx           0.6655     0.6154  1.081 2.80e-01
# ecog.ps     -0.0908     0.5916 -0.154 8.78e-01
#
set data ../pedigrees
set loc futime qua
set loc fustat aff
set loc age qua
set loc resid.ds aff
set loc rx aff
set loc ecog.ps qua
read ped ovarian.ped
run
kap futime fustat res
his futime
regress futime on age resid.ds rx ecog.ps
regress futime on age rx ecog.ps
regress futime on age
