#
# 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.
# 
#-----

#postscript("../../WriteUp/Graphics/Chapter4/prob_5_plot.eps", onefile=FALSE, horizontal=FALSE)

oldpar = par(mfrow=c(2,2))

# groups of 1000
tenfold1000 = rep( 1:10, rep(1000,10) )
boxplot( split( rnorm(1000*10), tenfold1000 ), ylab="normal - 1000" )
boxplot( split( rt(1000*10, 7), tenfold1000 ), ylab="t[7] - 1000" )

# groups of 100 
tenfold100 = rep( 1:10, rep(100,10) )
boxplot( split( rnorm(100*10), tenfold100 ), ylab="normal - 100" )
boxplot( split( rt(100*10, 7), tenfold100 ), ylab="t[7] - 100" )

par(oldpar)

#dev.off()