/*
** FILE BPM.E3
**
** Illustrates the use of BPROBRS and BPRS.
*/

library a:\bpm,pgraph;

bpmset;

r=3;
s=4;
x=ones(20,1)~rndn(20,2);
y=round(rndu(20,r*s));
let indvar=const covar1 covar2;
let depvar=dep1 dep2;

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

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

dataset="a:\\bpm1";
let indvar=x1 x2;
let depvar=y1 y2;

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

cls;
print "See the files e31.out and e32.out";