R version 2.15.2 (2012-10-26) -- "Trick or Treat" Copyright (C) 2012 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i686-pc-linux-gnu (32-bit) R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > par3 = '0.65' > par2 = '13' > par1 = '63' > par1<-as.numeric(par1) > par2<-as.numeric(par2) > par3<-as.numeric(par3) > df <- par1 - 1 > halfalpha <- (1 - par3) / 2 > ua <- qchisq(halfalpha,df) * par2 / df > ub <- qchisq(1-halfalpha,df) * par2 / df > ua [1] 10.81229 > ub [1] 15.15259 > ul <- qchisq(1-par3,df) * par2 / df > ul [1] 11.98708 > ur <- qchisq(par3,df) * par2 / df > ur [1] 13.77536 > > #Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/wessaorg/rcomp/createtable") > > a<-table.start() > a<-table.row.start(a) > a<-table.element(a,hyperlink('http://www.xycoon.com/ht_variance.htm','Testing Variance - Confidence Intervals for Sample Variance','learn more about Statistical Hypothesis Testing about the 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,'Null hypothesis (H0)',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.row.start(a) > a<-table.element(a, hyperlink('http://www.xycoon.com/ht_variance.htm#ex3', 'more information about confidence interval','example'),3,TRUE) > a<-table.row.end(a) > a<-table.end(a) > table.save(a,file="/var/wessaorg/rcomp/tmp/12yri1356100699.tab") > > > > proc.time() user system elapsed 0.315 0.095 0.395