/*
Example program: calls mfem_cox.src
*/

new;
/* Gehan Data: Time on Remission for Matched Pairs of Leukemia Patients */
let dta[42,5] =
 1  0 0  6 1
 2  0 1  6 2
 3  0 1  6 3
 4  0 1  6 4
 5  0 1  7 5
 6  0 0  9 6
 7  0 0 10 7
 8  0 1 10 8
 9  0 0 11 9
10  0 1 13 10
11  0 1 16 11
12  0 0 17 12
13  0 0 19 13
14  0 0 20 14
15  0 1 22 15
16  0 1 23 16
17  0 0 25 17
18  0 0 32 18
19  0 0 32 19
20  0 0 34 20
21  0 0 35 21
22  1 1  1 1
23  1 1  1 2
24  1 1  2 3
25  1 1  2 4
26  1 1  3 5
27  1 1  4 6
28  1 1  4 7
29  1 1  5 8
30  1 1  5 9
31  1 1  8 10
32  1 1  8 11
33  1 1  8 12
34  1 1  8 13
35  1 1 11 14
36  1 1 11 15
37  1 1 12 16
38  1 1 12 17
39  1 1 15 18
40  1 1 17 19
41  1 1 22 20
42  1 1 23 21;

output file=surv_pair_cox.out reset;

dta=sortc(dta,5);

  m=rows(dta);
 id=dta[.,5];
  d=dta[.,3];
  t=dta[.,4];
  c=unique(id,1); /* id number for pair */
  n=rows(c);

ndta = dta[.,2];

let dlist =  TMP6; /* = 1 defines control group */
let model =  TMP6 ;

print_it=1;
write_frail=0;
ffname="frail_file.dat";
alpha_st= 1;
maxit=150;
maxit1=20;
tol = 1.e-5;


#include mfem_cox.src;