/*
** FILE BPM.E1
**
** Illustrates the use of BIVPROB and PROFLIK.
*/

library a:\bpm,pgraph;

bpmset;

x={1 22.50,
   1 27.50,
   1 32.50,
   1 37.50,
   1 42.50,
   1 47.50,
   1 52.50,
   1 57.50,
   1 62.50};

y={1841   7  95   9,
   1654   9 105  23,
   1863  19 177  54,
   2357  48 257 121,
   1778  54 273 169,
   1712  88 324 269,
   1324 117 245 404,
    967 152 225 406,
    526 106 132 372};

let indvar=constant age;
let depvar=breathl wheeze;

_restric=1;
_rpar=0;
_rhonorm=1;
_initval=1;
_signif=1;
_tsiatis=1;

{teta,ll,tetra,llt,teta1,ll1,teta2,ll2}=bivprob(x,y,indvar,depvar);
dos copy bprobit.out e11.out;

_restric=1;
_rpar=1;
_rhonorm=0;
_initval=0;
_signif=0;
_tsiatis=0;

{teta,ll,tetra,llt,teta1,ll1,teta2,ll2}=bivprob(x,y,indvar,depvar);
dos copy bprobit.out e12.out;
dos del bprobit.out;

cls;
hardcopy("5");
_rpar=0;
_rhonorm=1;
call proflik(x,y,indvar,depvar,0.025,0.6,0.9);


cls;
print "See the files e11.out and e12.out";