library(foreign)
library(tigerstats)
time<-read.dta(file="http://kanggc.iptime.org/book/data/chap11-2-1.dta")
time$var1
(n<-length(time$var1))
(m<-mean(time$var1))
(s<-sd(time$var1))
(tc<-(m-34.5)/(s/sqrt(n)))  
t.test(time$var1, mu=34.5, conf.level=0.9)
(t14<-qt(0.05, df=14, lower.tail=T))
(t14<-qt(0.05, df=14, lower.tail=F))
par(mfrow=c(2,1))
ptGC(c(-4.3136,4.3136),region="outside",df=14, graph=T)
ptGC(c(-1.76131,1.76131),region="outside",df=14, graph=T)