###################################################
### chunk number 1: 
###################################################
whitepine = c(61,17,38,32,30,38,25,46,38,27,43,31,34,41,27,22,40,22)


###################################################
### chunk number 2: 
###################################################
t.test(whitepine,mu=40)


###################################################
### chunk number 3: 
###################################################
t.test(whitepine,mu=40,alternative="less")


###################################################
### chunk number 4: 
###################################################
ybar =  mean( whitepine )
ybar
s = sd( whitepine )
s
n = 18
mu = 40
t.value = (ybar - mu ) / ( s / sqrt( n ))
t.value


###################################################
### chunk number 5: 
###################################################
2 * pt(abs(t.value), n-1, lower.tail=F)


###################################################
### chunk number 6: 
###################################################
alpha = 0.10
qt(alpha, n-1, lower.tail=F)


###################################################
### chunk number 7: 
###################################################
pop.var = 100
pnorm(ybar,mu,sqrt(pop.var/n))


###################################################
### chunk number 8: 
###################################################
n.draw = 200
alpha = 0.05
draws = matrix( rnorm( n.draw * n, mu, sqrt(pop.var) ), n )
get.p.value = function( x ) t.test( x, mu = mu )$p.value
pvalues = apply( draws, 2, get.p.value )
sum( pvalues <= alpha )


###################################################
### chunk number 9: 
###################################################
y = 83
n = 100
p = 0.75
prop.test(y,n,p)


###################################################
### chunk number 10: 
###################################################
n = 18
sample.var = var( whitepine)
sample.var
pop.var = 100
pchisq((n-1)*sample.var/pop.var,n-1,lower.tail=FALSE)