read;file=a:ass1283.dat;nobs=1283;nvar=6;format=(7x,3f4.0,3f12.0);blanks;
names=size,ind,dy,ls85,ls80,lage$
dstats;rhs=size,ind,dy,ls85,ls80,lage$
reject;ls80<0;ls85<0$
crmodel;lhs=ls85;rhs=one,ls80$
reject;new$
reject;ls80<0;ls85<0&dy=1$
dstats;rhs=size,ind,dy,ls85,ls80,lage$
create;ls80sq=ls80*ls80;lagesq=lage*lage;lsage=ls80*lage$
probit;lhs=dy;rhs=one,ls80,ls80sq,lage,lagesq,lsage;hold$
select;lhs=ls85;rhs=one,ls80,lage;mle$