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

n <- 54
pprime <- 5
sigmaHat <- 4.0 
eHatI <- c(1.,1.732,9.,10.295)
hII <- c(0.9,0.75,0.25,0.185)

rI <- eHatI / ( sigmaHat * sqrt( 1 - hII ) )

DI <- (rI^2) * hII / ( pprime * ( 1 - hII ) ) 

tI <- rI * sqrt( ( n-pprime-1 ) / ( n - pprime - rI^2 ) )

pValues <- 2* ( 1 - pt( tI, n-pprime-1 ) )