x <- c(1,1,1,1,1,0,0,1,0,1,0,0,0,0,0,0,0,1,0,0,1,1,1,1,1,1,1,1,0,0,1,0,0,0,1)
par2 = '0'
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/1jpeg1350986052.tab") 


