# Statistics Lab Math M6 - simulated teaching model
xbar <- 68; s <- 10; n <- 25
tcrit <- qt(0.975, df=n-1)
se <- s/sqrt(n)
interval <- xbar + c(-1,1)*tcrit*se
c(tcrit=tcrit, se=se, lower=interval[1], upper=interval[2])