set more off
set scheme sj

sjlog using ../../mvprobit1, replace
use http://www.stata-press.com/data/r7/school, clear
biprobit (private=years logptax loginc) (vote=years logptax loginc), nolog
sjlog close, replace

sjlog using ../../mvprobit2, replace
predict p11, p11
predict p00, p00
predict xbb1, xb1
predict xbb2, xb2
predict stdpb1, stdp1
predict stdpb2, stdp2
predict pmargb1, pmarg1
predict pmargb2, pmarg2
sjlog close, replace

sjlog using ../../mvprobit3, replace
mvprobit (private = years logptax loginc) (vote=years logptax loginc), nolog
sjlog close, replace

sjlog using ../../mvprobit3b, replace
mvprobit (private = years logptax loginc) (vote=years logptax loginc), nolog draws(100)
sjlog close, replace

sjlog using ../../mvprobit4, replace
mvppred pall, pall
mvppred xbm, xb
mvppred stdpm, stdp
mvppred pmargm, pmarg
sjlog close, replace

sjlog using ../../mvprobit5, replace
summarize pall1s p11 pall0s p00 xbm1 xbb1 xbm2 xbb2 /*
        */ stdpm1 stdpb1 stdpm2 stdpb2 pmargm1 pmargb1 pmargm2 pmargb2
sjlog close, replace

sjlog using ../../mvprobit6, replace
set seed 12309
set obs 5000
matrix R = (1, .25, .5, .75 \ .25, 1, .75, .5 \ .5, .75, 1, .75 \ .75, .5, .75, 1)
drawnorm u1 u2 u3 u4, corr(R)
correlate u*
generate x1 = uniform()-.5
generate x2 = uniform() + 1/3
generate x3 = 2*uniform() + .5
generate x4 = .5*uniform() - 1/3
*Equations
ge y1s = .5 + 4*x1 + u1
ge y2s = 3 + .5*x1 - 3*x2 + u2
ge y3s = 1 - 2*x1 + .4*x2 - .75*x3 + u3
ge y4s = -6 + 1*x1 - .3*x2 + 3*x3 - .4*x4 + u4
ge y1 = y1s>0
ge y2 = y2s>0
ge y3 = y3s>0
ge y4 = y4s>0
sjlog close, replace

sjlog using ../../mvprobit7, replace
mvprobit (y1=x1) (y2=x1 x2) (y3 = x1 x2 x3) (y4=x1 x2 x3 x4), dr(75)
sjlog close, replace