* -------------------------------------
* Exercise 2.16 in Stock and Watson, 2e
* -------------------------------------

drawnorm z1 z2 z3 z4 z5 z6 z7 z8 z9 z10, n(100) means(5 5 5 5 5 5 5 5 5 5) sds(10 10 10 10 10 10 10 10 10 10) clear
summarize z1 if z1 < 3.6
scalar p1 = r(N)/100
summarize z2 if z2 < 3.6
scalar p2 = r(N)/100
summarize z3 if z3 < 3.6
scalar p3 = r(N)/100
summarize z4 if z4 < 3.6
scalar p4 = r(N)/100
summarize z5 if z5 < 3.6
scalar p5 = r(N)/100
summarize z6 if z6 < 3.6
scalar p6 = r(N)/100
summarize z7 if z7 < 3.6
scalar p7 = r(N)/100
summarize z8 if z8 < 3.6
scalar p8 = r(N)/100
summarize z9 if z9 < 3.6
scalar p9 = r(N)/100
summarize z10 if z10 < 3.6
scalar p10 = r(N)/100

scalar list p1 p2 p3 p4 p5 p6 p7 p8 p9 p10 
scalar prop = (p1+p2+p3+p4+p5+p6+p7+p8+p9+p10)/10
scalar list prop

scalar z = normal((3.6-5)/10)
scalar list z