/* program example to call coxfit */


new;
/* Namboodiri & Suchindran 21.1 pp. 212 */

let dta[12,4]=
1   1    1   1
2   1    1   3
3   1    0   3
4   1    1   4
5   1    1   5
6   1    1   5
7   1    0   9
8   0    1   2
9   0    0   3
10  0    1   4
11  0    1   6
12  0    1   6;


id=dta[.,1];
 x=dta[.,2];
 d=dta[.,3];
 t=dta[.,4];
offs=ones(rows(t),1);
 k=cols(x);
 b=zeros(k,1);

{m1, m2, h ,b}=coxfit(b,t,d,x,offs,.00001,10);

b~sqrt(diag(invpd(h)))~b./sqrt(diag(invpd(h)));

#include coxfit.src;