read; nobs=122; nvar=14; file=a:top9388.wks; format=wks ;names  $
dstat; rhs=usd,asl1,asl2,asd1,asd2,tsl1,tsl2,tsd1,tsd2 $
dstat; rhs=p1,p2,e1,e2 $
 
create;las1=log(asd1);las2=log(asd2); lts1=log(tsd1); lts2=log(tsd2)$
create;le1=log(e1); le2=log(e2) $
create;shd1=asd1/tsd1;shd2=asd2/tsd2 $
dstat; rhs=las1,las2,lts1,lts2,le1,le2,shd1,shd2 $
 
sample; all$
create; lx1=las1;lx2=las2 $
 
reject; lx1<0; lx2<0 $
dstat; rhs=lx1,lx2 $
crmodel;lhs=lx1;rhs=one,lx2; het $
crmodel;lhs=lx1;rhs=one,lx2,usd; het $
crmodel;lhs=lx1;rhs=one,lx2,usd,shd2; het $
plot; lhs=lx2; rhs=lx1 $

sample; all$
sample; 1-122 $
create; dum=1 $
reject; lx1<0 $
create; dum=0 $
dstat; rhs=dum $
 
sample; all$
sample; 1-122 $
create; lx1sq=lx1*lx1; lx2sq=lx2*lx2 $
 
reject;lx2<0 $
dstat;rhs=dum,lx1,lx2 $
list; dum,lx2 $
logit;lhs=dum;rhs=one,lx2 ; hold $
? probit;lhs=dum;rhs=one,lx2; list; alg=dfp; hold $
? probit;lhs=dum;rhs=one,lx2;start=0.0; hold $
select; lhs=lx1;rhs=one,lx2 $

sample; all$
sample; 1-122 $
dstat; rhs=dum,lx1,lx2 $
probit; lhs=dum; rhs=one,lx2,lx2sq;hold $
select; lhs=lx1; rhs=one,lx2 $