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

df <- 46
pprime <- 5
n <- df + pprime 
sigmaHat <- 64.891 
eHatI <- c(-163.145, -137.599,-102.409,183.499,-49.452)
hII <- c(0.256,0.162,0.206,0.084,0.415)

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

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

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