#
# Written by:
# -- 
# John L. Weatherwax                2009-04-21
# 
# email: wax@alum.mit.edu
# 
# Please send comments and especially bug reports to the
# above email address.
# 
#-----

source('chap_4_prob_3.R')

n_sims = 50
X = matrix( 0, nrow=1, ncol=n_sims )
for( ii in 1:n_sims ){
  X[ii] = mean_zero_p_value()
}

#postscript("../../WriteUp/Graphics/Chapter4/prob_4_plot.eps", onefile=FALSE, horizontal=FALSE)
qqplot( X, qunif(ppoints(n_sims)) )
#dev.off()

print( min( X ) )
print( max( X ) )