x <- c(1,1,1,1,1,1,1,0,1,1,1,1,0,1,1,1,0,1,1,1,1,0,0,1,1,1,0,0,1,1,0,1,1)
par20 = ''
par19 = ''
par18 = ''
par17 = ''
par16 = ''
par15 = ''
par14 = ''
par13 = ''
par12 = ''
par11 = ''
par10 = ''
par9 = ''
par8 = ''
par7 = ''
par6 = ''
par5 = ''
par4 = ''
par3 = ''
par2 = '50'
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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/1w25k1350133692.tab")