par1<-as.numeric(par1) par2<-as.numeric(par2) par3<-as.numeric(par3) par4<-as.numeric(par4) df <- par1 - 1 if (par2 > par3) { myc <- par3 / df * qchisq(1-par4,df) } else { myc <- par3 / df * qchisq(par4,df) } myc load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,hyperlink('ht_variance.htm','Testing Variance - Critical value (region)','learn more about Statistical Hypothesis Testing about the Variance'),2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Sample size',header=TRUE) a<-table.element(a,par1) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Sample variance',header=TRUE) a<-table.element(a,par2) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Null hypothesis (H0)',header=TRUE) a<-table.element(a,par3) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Type I error (alpha)',header=TRUE) a<-table.element(a,par4) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Critical value',header=TRUE) a<-table.element(a,myc) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Conclusion',header=TRUE) if (par2 > par3) { if (par2 < myc) a<-table.element(a,'There is no reason to reject the null hypothesis') else a<-table.element(a,'Reject the null hypothesis') } else { if (par2 > myc) a<-table.element(a,'There is no reason to reject the null hypothesis') else a<-table.element(a,'Reject the null hypothesis') } a<-table.row.end(a) a<-table.end(a) table.save(a,file='mytable.tab')
|