#
# 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