x <- c(19,25,17,22,21,26,20,14,23,20,22,15,20,22,20,28,25,26,17,23,13,24,14,22,23,22,24,21,23,22,21,26,15,25,17,25,27,25,19,26,20,20,18,18,19,23,17,23,23,11,18,24,16,24,24,21,25,22,21,24,24,21,18)
par2 = '20'
par1 = '0.95'
par1 <- as.numeric(par1)
par2 <- as.numeric(par2)
len <- length(x)
df <- len - 1
sd <- sd(x)
mx <- mean(x)

#Note: the /var/fisher/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/fisher/rcomp/createtable")

delta2 <- abs(qt((1-par1)/2,df)) * sd / sqrt(len)
delta1 <- abs(qt((1-par1),df)) * sd / sqrt(len)
a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Hypothesis Test about the Mean - Confidence Interval',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Sample size',header=TRUE)
a<-table.element(a,len)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Sample standard deviation',header=TRUE)
a<-table.element(a,sd)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Confidence',header=TRUE)
a<-table.element(a,par1)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Null hypothesis',header=TRUE)
a<-table.element(a,par2)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Sample Mean',header=TRUE)
a<-table.element(a,mx)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'2-sided Confidence Interval',header=TRUE)
dum <- paste('[',mx-delta2)
dum <- paste(dum,',')
dum <- paste(dum,mx+delta2)
dum <- paste(dum,']')
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Left-sided Confidence Interval',header=TRUE)
dum <- paste('[',mx-delta1)
dum <- paste(dum,', +inf ]')
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'Right-sided Confidence Interval',header=TRUE)
dum <- paste('[ -inf, ',mx+delta1)
dum <- paste(dum,']')
a<-table.element(a,dum)
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file="/var/fisher/rcomp/tmp/1uq5z1383172407.tab")