#
# Diabetic Retinopathy Study
# Huster et al. , and Therneau and Grambsch analysed these data
#
# id:   patient number
# trteye:   treated eye (1=right 2=left)
# ageonset:   age at diagnosis of diabetes
# typediab:  type of diabetes (1= juvenile age at onset < 20 years , 2=adult)
# status:   status for eye (0=censored, 1=blindness)
# time:   follow-up time in months
# trt:   treatment of eye (1=control eye, 2=treated eye) 
#
set loc eye      cat . 1=right 2=left
set loc trteye   cat . 1=right 2=left
set loc ageonset qua . age at diagnosis of diabetes
set loc typediab cat . 1=T1D 2=T2D
set loc status   qua . blindness
set loc time     qua . time to blindness
set loc trt      cat . laser treatment 1=control 2=laser
read pedigree inline
#   id trteye ageonset typediab status  time trt 
5    right x x x 1 2 28 2 0 46.23 1 
5    left  x x x 2 2 28 2 0 46.23 2 
14   right x x x 1 1 12 1 0 42.50 2 
14   left  x x x 2 1 12 1 1 31.30 1 
16   right x x x 1 1 9  1 0 42.27 2 
16   left  x x x 2 1 9  1 0 42.27 1 
25   right x x x 1 2 9  1 0 20.60 1 
25   left  x x x 2 2 9  1 0 20.60 2 
29   right x x x 1 2 13 1 1 0.30  1 
29   left  x x x 2 2 13 1 0 38.77 2 
46   right x x x 1 1 12 1 0 65.23 2 
46   left  x x x 2 1 12 1 1 54.27 1 
49   right x x x 1 1 8  1 0 63.50 2 
49   left  x x x 2 1 8  1 1 10.80 1 
56   right x x x 1 1 12 1 0 23.17 2 
56   left  x x x 2 1 12 1 0 23.17 1 
61   right x x x 1 1 16 1 0 1.47  2 
61   left  x x x 2 1 16 1 0 1.47  1 
71   right x x x 1 1 21 2 0 58.07 2 
71   left  x x x 2 1 21 2 1 13.83 1 
100  right x x x 1 2 23 2 0 48.53 1 
100  left  x x x 2 2 23 2 1 46.43 2 
112  right x x x 1 1 44 2 0 44.40 2 
112  left  x x x 2 1 44 2 1 7.90  1 
120  right x x x 1 2 47 2 0 39.57 1 
120  left  x x x 2 2 47 2 0 39.57 2 
127  right x x x 1 1 48 2 1 30.83 2 
127  left  x x x 2 1 48 2 1 38.57 1 
133  right x x x 1 1 26 2 0 66.27 2 
133  left  x x x 2 1 26 2 1 14.10 1 
150  right x x x 1 2 10 1 1 6.90  1 
150  left  x x x 2 2 10 1 1 20.17 2 
167  right x x x 1 2 23 2 1 41.40 1 
167  left  x x x 2 2 23 2 0 58.43 2 
176  right x x x 1 1 5  1 0 58.20 2 
176  left  x x x 2 1 5  1 0 58.20 1 
185  right x x x 1 1 46 2 0 57.43 2 
185  left  x x x 2 1 46 2 0 57.43 1 
190  right x x x 1 1 5  1 0 56.03 2 
190  left  x x x 2 1 5  1 0 56.03 1 
202  right x x x 1 2 13 1 0 67.53 1 
202  left  x x x 2 2 13 1 0 67.53 2 
214  right x x x 1 1 45 2 0 61.40 2 
214  left  x x x 2 1 45 2 1 0.60  1 
220  right x x x 1 1 11 1 1 10.27 2 
220  left  x x x 2 1 11 1 1 1.63  1 
243  right x x x 1 2 1  1 0 66.20 1 
243  left  x x x 2 2 1  1 0 66.20 2 
255  right x x x 1 2 10 1 1 13.83 1 
255  left  x x x 2 2 10 1 1 5.67  2 
264  right x x x 1 1 12 1 0 58.83 2 
264  left  x x x 2 1 12 1 1 29.97 1 
266  right x x x 1 1 36 2 0 60.27 2 
266  left  x x x 2 1 36 2 1 26.37 1 
284  right x x x 1 2 53 2 1 1.33  1 
284  left  x x x 2 2 53 2 1 5.77  2 
295  right x x x 1 2 10 1 1 35.53 1 
295  left  x x x 2 2 10 1 1 5.90  2 
300  right x x x 1 2 25 2 1 21.90 1 
300  left  x x x 2 2 25 2 1 25.63 2 
302  right x x x 1 2 14 1 1 14.80 1 
302  left  x x x 2 2 14 1 1 33.90 2 
315  right x x x 1 2 16 1 1 6.20  1 
315  left  x x x 2 2 16 1 1 1.73  2 
324  right x x x 1 2 38 2 1 22.00 1 
324  left  x x x 2 2 38 2 0 46.90 2 
328  right x x x 1 2 14 1 0 31.13 1 
328  left  x x x 2 2 14 1 0 31.13 2 
335  right x x x 1 2 10 1 1 22.00 1 
335  left  x x x 2 2 10 1 1 30.20 2 
342  right x x x 1 2 17 1 0 70.90 1 
342  left  x x x 2 2 17 1 0 70.90 2 
349  right x x x 1 1 44 2 1 25.80 2 
349  left  x x x 2 1 44 2 1 13.87 1 
357  right x x x 1 2 21 2 1 48.30 1 
357  left  x x x 2 2 21 2 1 5.73  2 
368  right x x x 1 2 19 1 0 53.43 1 
368  left  x x x 2 2 19 1 0 53.43 2 
385  right x x x 1 1 13 1 1 1.90  2 
385  left  x x x 2 1 13 1 0 51.10 1 
396  right x x x 1 2 40 2 1 9.90  1 
396  left  x x x 2 2 40 2 1 9.90  2 
405  right x x x 1 1 9  1 0 34.20 2 
405  left  x x x 2 1 9  1 0 34.20 1 
409  right x x x 1 2 48 2 1 2.67  1 
409  left  x x x 2 2 48 2 0 46.73 2 
419  right x x x 1 1 42 2 0 18.73 2 
419  left  x x x 2 1 42 2 1 13.83 1 
429  right x x x 1 1 24 2 0 32.03 2 
429  left  x x x 2 1 24 2 1 4.27  1 
433  right x x x 1 2 55 2 1 13.90 1 
433  left  x x x 2 2 55 2 0 69.87 2 
445  right x x x 1 2 17 1 0 66.80 1 
445  left  x x x 2 2 17 1 0 66.80 2 
454  right x x x 1 1 5  1 0 64.73 2 
454  left  x x x 2 1 5  1 0 64.73 1 
468  right x x x 1 1 6  1 1 1.70  2 
468  left  x x x 2 1 6  1 1 1.70  1 
480  right x x x 1 1 19 1 1 1.77  2 
480  left  x x x 2 1 19 1 1 43.03 1 
485  right x x x 1 1 12 1 0 29.03 2 
485  left  x x x 2 1 12 1 0 29.03 1 
491  right x x x 1 2 45 2 0 56.57 1 
491  left  x x x 2 2 45 2 0 56.57 2 
503  right x x x 1 2 27 2 1 8.30  1 
503  left  x x x 2 2 27 2 1 8.30  2 
515  right x x x 1 1 43 2 0 21.57 2 
515  left  x x x 2 1 43 2 1 18.43 1 
522  right x x x 1 2 4  1 0 31.57 1 
522  left  x x x 2 2 4  1 0 31.57 2 
538  right x x x 1 2 45 2 1 31.63 1 
538  left  x x x 2 2 45 2 0 31.63 2 
547  right x x x 1 1 32 2 0 39.77 2 
547  left  x x x 2 1 32 2 0 39.77 1 
550  right x x x 1 2 3  1 1 6.53  1 
550  left  x x x 2 2 3  1 1 18.70 2 
554  right x x x 1 1 14 1 0 18.90 2 
554  left  x x x 2 1 14 1 0 18.90 1 
557  right x x x 1 1 13 1 0 56.80 2 
557  left  x x x 2 1 13 1 1 22.23 1 
561  right x x x 1 1 15 1 0 55.60 2 
561  left  x x x 2 1 15 1 1 14.00 1 
568  right x x x 1 2 10 1 1 42.17 1 
568  left  x x x 2 2 10 1 1 42.17 2 
572  right x x x 1 2 6  1 1 5.33  1 
572  left  x x x 2 2 6  1 0 10.70 2 
576  right x x x 1 2 17 1 1 59.80 1 
576  left  x x x 2 2 17 1 0 66.33 2 
581  right x x x 1 2 37 2 1 5.83  1 
581  left  x x x 2 2 37 2 0 52.33 2 
606  right x x x 1 1 18 1 0 58.17 2 
606  left  x x x 2 1 18 1 1 2.17  1 
610  right x x x 1 2 13 1 1 48.43 1 
610  left  x x x 2 2 13 1 1 14.30 2 
615  right x x x 1 2 14 1 0 25.83 1 
615  left  x x x 2 2 14 1 0 25.83 2 
618  right x x x 1 1 12 1 0 45.40 2 
618  left  x x x 2 1 12 1 0 45.40 1 
624  right x x x 1 1 9  1 0 47.60 2 
624  left  x x x 2 1 9  1 0 47.60 1 
631  right x x x 1 2 11 1 1 9.60  1 
631  left  x x x 2 2 11 1 1 13.33 2 
636  right x x x 1 2 10 1 0 42.10 1 
636  left  x x x 2 2 10 1 0 42.10 2 
645  right x x x 1 1 5  1 0 39.93 2 
645  left  x x x 2 1 5  1 0 39.93 1 
653  right x x x 1 2 15 1 1 7.60  1 
653  left  x x x 2 2 15 1 1 14.27 2 
662  right x x x 1 2 7  1 1 1.80  1 
662  left  x x x 2 2 7  1 1 34.57 2 
664  right x x x 1 2 2  1 1 4.30  1 
664  left  x x x 2 2 2  1 0 65.80 2 
683  right x x x 1 2 22 2 1 12.20 1 
683  left  x x x 2 2 22 2 1 4.10  2 
687  right x x x 1 2 5  1 0 60.93 1 
687  left  x x x 2 2 5  1 0 60.93 2 
701  right x x x 1 2 4  1 0 57.20 1 
701  left  x x x 2 2 4  1 0 57.20 2 
706  right x x x 1 1 27 2 0 38.07 2 
706  left  x x x 2 1 27 2 1 12.73 1 
717  right x x x 1 1 53 2 0 54.10 2 
717  left  x x x 2 1 53 2 1 54.10 1 
722  right x x x 1 1 10 1 0 59.27 2 
722  left  x x x 2 1 10 1 1 9.40  1 
731  right x x x 1 2 13 1 1 9.90  1 
731  left  x x x 2 2 13 1 1 21.57 2 
740  right x x x 1 2 12 1 0 54.10 1 
740  left  x x x 2 2 12 1 0 54.10 2 
749  right x x x 1 2 24 2 0 50.47 1 
749  left  x x x 2 2 24 2 0 50.47 2 
757  right x x x 1 2 17 1 0 46.17 1 
757  left  x x x 2 2 17 1 0 46.17 2 
760  right x x x 1 2 8  1 0 46.30 1 
760  left  x x x 2 2 8  1 0 46.30 2 
766  right x x x 1 1 58 2 0 38.83 2 
766  left  x x x 2 1 58 2 0 38.83 1 
769  right x x x 1 2 17 1 0 44.60 1 
769  left  x x x 2 2 17 1 0 44.60 2 
772  right x x x 1 2 12 1 0 43.07 1 
772  left  x x x 2 2 12 1 0 43.07 2 
778  right x x x 1 2 25 2 0 40.03 1 
778  left  x x x 2 2 25 2 1 26.23 2 
780  right x x x 1 1 15 1 0 41.60 2 
780  left  x x x 2 1 15 1 1 18.03 1 
793  right x x x 1 2 21 2 0 38.07 1 
793  left  x x x 2 2 21 2 0 38.07 2 
800  right x x x 1 1 20 2 0 65.23 2 
800  left  x x x 2 1 20 2 0 65.23 1 
804  right x x x 1 1 23 2 1 7.07  2 
804  left  x x x 2 1 23 2 0 66.77 1 
810  right x x x 1 2 13 1 1 13.77 1 
810  left  x x x 2 2 13 1 1 13.77 2 
815  right x x x 1 2 45 2 1 9.63  1 
815  left  x x x 2 2 45 2 0 9.63  2 
832  right x x x 1 1 5  1 0 46.23 2 
832  left  x x x 2 1 5  1 0 46.23 1 
834  right x x x 1 2 8  1 1 1.50  1 
834  left  x x x 2 2 8  1 0 45.73 2 
838  right x x x 1 2 30 2 1 33.63 1 
838  left  x x x 2 2 30 2 1 33.63 2 
857  right x x x 1 2 7  1 0 40.17 1 
857  left  x x x 2 2 7  1 0 40.17 2 
866  right x x x 1 2 39 2 1 27.60 1 
866  left  x x x 2 2 39 2 1 63.33 2 
887  right x x x 1 1 26 2 1 38.47 2 
887  left  x x x 2 1 26 2 1 1.63  1 
903  right x x x 1 1 50 2 0 55.23 2 
903  left  x x x 2 1 50 2 0 55.23 1 
910  right x x x 1 2 34 2 1 25.30 1 
910  left  x x x 2 2 34 2 0 52.77 2 
920  right x x x 1 2 10 1 1 46.20 1 
920  left  x x x 2 2 10 1 0 57.17 2 
925  right x x x 1 1 40 2 0 9.87  2 
925  left  x x x 2 1 40 2 1 1.70  1 
931  right x x x 1 2 13 1 0 57.90 1 
931  left  x x x 2 2 13 1 0 57.90 2 
936  right x x x 1 2 7  1 0 5.90  1 
936  left  x x x 2 2 7  1 0 5.90  2 
945  right x x x 1 2 11 1 0 32.20 1 
945  left  x x x 2 2 11 1 0 32.20 2 
949  right x x x 1 1 13 1 1 10.33 2 
949  left  x x x 2 1 13 1 1 0.83  1 
952  right x x x 1 2 9  1 0 50.90 1 
952  left  x x x 2 2 9  1 1 6.13  2 
962  right x x x 1 2 5  1 1 25.93 1 
962  left  x x x 2 2 5  1 0 43.67 2 
964  right x x x 1 2 10 1 0 38.30 1 
964  left  x x x 2 2 10 1 0 38.30 2 
971  right x x x 1 1 23 2 0 38.77 2 
971  left  x x x 2 1 23 2 1 19.40 1 
978  right x x x 1 2 2  1 1 21.97 1 
978  left  x x x 2 2 2  1 0 38.07 2 
983  right x x x 1 1 12 1 0 38.30 2 
983  left  x x x 2 1 12 1 0 38.30 1 
987  right x x x 1 2 7  1 0 70.03 1 
987  left  x x x 2 2 7  1 1 26.20 2 
1002 right x x x 1 2 13 1 1 18.03 1 
1002 left  x x x 2 2 13 1 0 62.57 2 
1017 right x x x 1 2 50 2 1 1.57  1 
1017 left  x x x 2 2 50 2 1 13.83 2 
1029 right x x x 1 1 20 2 0 46.50 2 
1029 left  x x x 2 1 20 2 1 13.37 1 
1034 right x x x 1 2 15 1 1 1.97  1 
1034 left  x x x 2 2 15 1 1 11.07 2 
1037 right x x x 1 1 30 2 0 42.47 2 
1037 left  x x x 2 1 30 2 1 22.20 1 
1042 right x x x 1 2 32 2 0 38.73 1 
1042 left  x x x 2 2 32 2 0 38.73 2 
1069 right x x x 1 2 39 2 0 51.13 1 
1069 left  x x x 2 2 39 2 0 51.13 2 
1074 right x x x 1 2 4  1 0 46.50 1 
1074 left  x x x 2 2 4  1 1 6.10  2 
1098 right x x x 1 2 3  1 1 11.30 1 
1098 left  x x x 2 2 3  1 1 2.10  2 
1102 right x x x 1 1 10 1 1 17.73 2 
1102 left  x x x 2 1 10 1 0 42.30 1 
1112 right x x x 1 2 6  1 0 26.47 1 
1112 left  x x x 2 2 6  1 0 26.47 2 
1117 right x x x 1 1 15 1 0 10.77 2 
1117 left  x x x 2 1 15 1 0 10.77 1 
1126 right x x x 1 1 33 2 0 55.33 2 
1126 left  x x x 2 1 33 2 0 55.33 1 
1135 right x x x 1 1 15 1 0 58.67 2 
1135 left  x x x 2 1 15 1 0 58.67 1 
1145 right x x x 1 2 44 2 1 4.97  1 
1145 left  x x x 2 2 44 2 1 12.93 2 
1148 right x x x 1 2 48 2 1 26.47 1 
1148 left  x x x 2 2 48 2 0 54.20 2 
1167 right x x x 1 2 4  1 0 49.57 1 
1167 left  x x x 2 2 4  1 0 49.57 2 
1184 right x x x 1 2 46 2 1 9.87  1 
1184 left  x x x 2 2 46 2 1 24.43 2 
1191 right x x x 1 2 25 2 0 50.23 1 
1191 left  x x x 2 2 25 2 0 50.23 2 
1205 right x x x 1 2 12 1 1 30.40 1 
1205 left  x x x 2 2 12 1 1 13.97 2 
1213 right x x x 1 1 12 1 0 43.33 2 
1213 left  x x x 2 1 12 1 1 43.33 1 
1228 right x x x 1 1 26 2 0 42.23 2 
1228 left  x x x 2 1 26 2 0 42.23 1 
1247 right x x x 1 2 11 1 0 74.93 1 
1247 left  x x x 2 2 11 1 0 74.93 2 
1250 right x x x 1 1 36 2 0 66.93 2 
1250 left  x x x 2 1 36 2 0 66.93 1 
1253 right x x x 1 1 12 1 0 73.43 2 
1253 left  x x x 2 1 12 1 0 73.43 1 
1267 right x x x 1 1 50 2 0 67.47 2 
1267 left  x x x 2 1 50 2 1 38.57 1 
1281 right x x x 1 1 44 2 0 3.67  2 
1281 left  x x x 2 1 44 2 1 3.67  1 
1287 right x x x 1 2 8  1 0 67.03 1 
1287 left  x x x 2 2 8  1 1 48.87 2 
1293 right x x x 1 1 14 1 0 65.60 2 
1293 left  x x x 2 1 14 1 0 65.60 1 
1296 right x x x 1 2 18 1 1 15.83 1 
1296 left  x x x 2 2 18 1 0 15.83 2 
1309 right x x x 1 1 56 2 0 20.07 2 
1309 left  x x x 2 1 56 2 1 8.83  1 
1312 right x x x 1 2 9  1 0 67.43 1 
1312 left  x x x 2 2 9  1 0 67.43 2 
1317 right x x x 1 1 15 1 0 1.47  2 
1317 left  x x x 2 1 15 1 0 1.47  1 
1321 right x x x 1 1 5  1 0 62.93 2 
1321 left  x x x 2 1 5  1 1 22.13 1 
1333 right x x x 1 1 1  1 1 6.30  2 
1333 left  x x x 2 1 1  1 0 56.97 1 
1347 right x x x 1 1 1  1 0 59.70 2 
1347 left  x x x 2 1 1  1 1 18.93 1 
1361 right x x x 1 2 14 1 1 19.00 1 
1361 left  x x x 2 2 14 1 1 13.80 2 
1366 right x x x 1 1 57 2 0 55.13 2 
1366 left  x x x 2 1 57 2 0 55.13 1 
1373 right x x x 1 2 8  1 1 5.43  1 
1373 left  x x x 2 2 8  1 1 13.57 2 
1397 right x x x 1 1 33 2 0 42.20 2 
1397 left  x x x 2 1 33 2 0 42.20 1 
1410 right x x x 1 2 46 2 0 38.27 1 
1410 left  x x x 2 2 46 2 0 38.27 2 
1413 right x x x 1 1 3  1 0 7.10  2 
1413 left  x x x 2 1 3  1 1 7.10  1 
1425 right x x x 1 2 35 2 1 26.17 1 
1425 left  x x x 2 2 35 2 0 63.63 2 
1447 right x x x 1 2 8  1 1 24.73 1 
1447 left  x x x 2 2 8  1 0 59.00 2 
1461 right x x x 1 1 30 2 0 54.37 2 
1461 left  x x x 2 1 30 2 0 54.37 1 
1469 right x x x 1 1 51 2 0 54.60 2 
1469 left  x x x 2 1 51 2 1 10.97 1 
1480 right x x x 1 2 42 2 1 21.10 1 
1480 left  x x x 2 2 42 2 0 63.87 2 
1487 right x x x 1 1 20 2 0 62.37 2 
1487 left  x x x 2 1 20 2 1 43.70 1 
1491 right x x x 1 2 23 2 0 62.80 1 
1491 left  x x x 2 2 23 2 0 62.80 2 
1499 right x x x 1 1 22 2 0 63.33 2 
1499 left  x x x 2 1 22 2 1 14.37 1 
1503 right x x x 1 2 25 2 0 58.53 1 
1503 left  x x x 2 2 25 2 0 58.53 2 
1513 right x x x 1 2 45 2 0 58.07 1 
1513 left  x x x 2 2 45 2 0 58.07 2 
1524 right x x x 1 2 20 2 0 58.50 1 
1524 left  x x x 2 2 20 2 0 58.50 2 
1533 right x x x 1 2 41 2 0 14.37 1 
1533 left  x x x 2 2 41 2 1 1.50  2 
1537 right x x x 1 1 19 1 0 54.73 2 
1537 left  x x x 2 1 19 1 1 38.40 1 
1552 right x x x 1 1 4  1 0 50.63 2 
1552 left  x x x 2 1 4  1 1 2.83  1 
1554 right x x x 1 2 36 2 0 51.10 1 
1554 left  x x x 2 2 36 2 0 51.10 2 
1562 right x x x 1 1 20 2 0 49.93 2 
1562 left  x x x 2 1 20 2 1 6.57  1 
1572 right x x x 1 1 24 2 0 46.27 2 
1572 left  x x x 2 1 24 2 1 46.27 1 
1581 right x x x 1 1 28 2 0 10.60 2 
1581 left  x x x 2 1 28 2 0 10.60 1 
1585 right x x x 1 1 51 2 0 42.77 2 
1585 left  x x x 2 1 51 2 0 42.77 1 
1596 right x x x 1 1 16 1 1 34.37 2 
1596 left  x x x 2 1 16 1 0 42.27 1 
1600 right x x x 1 1 16 1 0 42.07 2 
1600 left  x x x 2 1 16 1 0 42.07 1 
1603 right x x x 1 2 10 1 0 38.77 1 
1603 left  x x x 2 2 10 1 0 38.77 2 
1619 right x x x 1 2 20 2 1 61.83 1 
1619 left  x x x 2 2 20 2 0 74.97 2 
1627 right x x x 1 2 10 1 0 66.97 1 
1627 left  x x x 2 2 10 1 1 6.57  2 
1636 right x x x 1 1 16 1 1 38.87 2 
1636 left  x x x 2 1 16 1 0 68.30 1 
1640 right x x x 1 2 10 1 1 46.63 1 
1640 left  x x x 2 2 10 1 1 42.43 2 
1643 right x x x 1 1 11 1 0 67.07 2 
1643 left  x x x 2 1 11 1 0 67.07 1 
1649 right x x x 1 1 1  1 1 2.70  2 
1649 left  x x x 2 1 1  1 0 2.70  1 
1666 right x x x 1 2 17 1 0 63.80 1 
1666 left  x x x 2 2 17 1 0 63.80 2 
1672 right x x x 1 2 7  1 0 32.63 1 
1672 left  x x x 2 2 7  1 0 32.63 2 
1683 right x x x 1 1 29 2 0 62.00 2 
1683 left  x x x 2 1 29 2 0 62.00 1 
1688 right x x x 1 2 5  1 0 54.80 1 
1688 left  x x x 2 2 5  1 1 13.10 2 
1705 right x x x 1 2 1  1 0 8.00  1 
1705 left  x x x 2 2 1  1 0 8.00  2 
1717 right x x x 1 2 22 2 1 42.33 1 
1717 left  x x x 2 2 22 2 0 51.60 2 
1727 right x x x 1 1 33 2 0 49.97 2 
1727 left  x x x 2 1 33 2 1 2.90  1 
1746 right x x x 1 1 3  1 0 45.90 2 
1746 left  x x x 2 1 3  1 1 1.43  1 
1749 right x x x 1 1 32 2 0 41.93 2 
1749 left  x x x 2 1 32 2 0 41.93 1 
;;;;
run
recast status affection
tab trt typediab
set loc inter aff
if (trt == 2 and typediab == 2) then inter=y else inter=n
reg time = trt + typediab + inter weibull status
set chain 4
set burn 1000
set iter 5000
fpm time nqtl 0 c lik weibull cens status cov trt + typediab + inter
