/*
** FILE BPM.E2
**
** Illustrates the use of BPROBIT.
*/

library a:\bpm,pgraph;

bpmset;

x=ones(50,1)~rndn(50,2);
y=round(rndu(50,2));

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

{teta,ll,tetra,llt,teta1,ll1,teta2,ll2}=bprobit("",x,y);
dos copy bprobit.out e21.out;

dataset="a:\\bpm";
let indvar=x1 x2 x3 x5;
let depvar=y1 y2;

{teta,ll,tetra,llt,teta1,ll1,teta2,ll2}=bprobit(dataset,indvar,depvar);
dos copy bprobit.out e22.out;
dos del bprobit.out;



_rpar=2;

hardcopy("5");
{rho,par,pl,tvec,pvec,step,min,max}=bprobit("",x,y);


cls;
print "See the files e21.out and e22.out";