# # Data from Feigl P, Zelen M (1965) Estimation of exponential probability # with concomitant information. Biometrics 1965; 21: 826-838. # (example from Aitkin et al (1989) Statistical Modelling in GLIM) # # Survival times in weeks of AML patients versus # WBC (x1000) and Ag factor (+/-) # No censoring # set loc surv qua set loc WBC qua set loc AG aff read pedigree inline 1 1 x x x 65 2.3 n 2 2 x x x 156 0.75 n 3 3 x x x 100 4.3 n 4 4 x x x 134 2.6 n 5 5 x x x 16 6.0 n 6 6 x x x 108 10.5 n 7 7 x x x 121 10.0 n 8 8 x x x 4 17.0 n 9 9 x x x 39 5.4 n 10 10 x x x 143 7.0 n 11 11 x x x 56 9.4 n 12 12 x x x 26 32.0 n 13 13 x x x 22 35.0 n 14 14 x x x 1 100.0 n 15 15 x x x 1 100.0 n 16 16 x x x 5 52.0 n 17 17 x x x 65 100.0 n 18 18 x x x 56 4.4 y 19 19 x x x 65 3.0 y 20 20 x x x 17 4.0 y 21 21 x x x 7 1.5 y 22 22 x x x 16 9.0 y 23 23 x x x 22 5.3 y 24 24 x x x 3 10.0 y 25 25 x x x 4 19.0 y 26 26 x x x 2 27.0 y 27 27 x x x 3 28.0 y 28 28 x x x 8 31.0 y 29 29 x x x 4 26.0 y 30 30 x x x 3 21.0 y 31 31 x x x 30 79.0 y 32 32 x x x 4 100.0 y 33 33 x x x 43 100.0 y ;;;; run set loc logt qua logt=log(surv) WBC = log(WBC) set loc WBC.AG qua WBC.AG=WBC*AG reg logt = WBC AG WBC.AG set ple 2 # fpm surv nqtl 0 g fix g gval 0 lik weibull cov WBC + AG + WBC.AG reg surv = WBC AG WBC.AG exponential reg surv = WBC AG WBC.AG weibull