par1<-as.numeric(par1) par2<-as.numeric(par2) par3<-as.numeric(par3) df <- par1 - 1 halfalpha <- (1 - par3) / 2 ub <- par2 * df / qchisq(halfalpha,df) ua <- par2 * df / qchisq(1-halfalpha,df) print(ua) print(ub) ur <- par2 * df / qchisq(1-par3,df) print(ur) ul <- par2 * df / qchisq(par3,df) print(ul) load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Testing Variance - Confidence Intervals for Population Variance',3,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Sample size',header=TRUE) a<-table.element(a,par1,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Sample variance',header=TRUE) a<-table.element(a,par2,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Confidence interval',header=TRUE) a<-table.element(a,par3,2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Type of Interval',header=TRUE) a<-table.element(a,'Left tail',header=TRUE) a<-table.element(a,'Right tail',header=TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('Two-sided confidence interval at ',par3), header=TRUE) a<-table.element(a,ua) a<-table.element(a,ub) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('Left one-sided confidence interval at ',par3), header=TRUE) a<-table.element(a,ul) a<-table.element(a,'+inf') a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,paste('Right one-sided confidence interval at ',par3), header=TRUE) a<-table.element(a,'0') a<-table.element(a,ur) a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|