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

ns = 1:100
t_critical_values = qt( 1-0.05, ns ) # assume a one sided test 

#postscript("../../WriteUp/Graphics/Chapter4/prob_2_plot.eps", onefile=FALSE, horizontal=FALSE)
#
#
par( mfrow=c(1,2) )
plot( ns, t_critical_values, xlab="degrees of freedom", ylab="t critical value" )
plot( log(ns), log(t_critical_values), xlab="log(degrees of freedom)", ylab="log(t critical value)" )
par( mfrow=c(1,1) )
#
#
#dev.off()