par1<-as.numeric(par1) par2<-as.numeric(par2) par3<-as.numeric(par3) par4<-as.numeric(par4) df <- par1 - 1 myc <- df * par2 / par3 myc if (par2 > par3) { myp <- 1 - pchisq(myc,df) } else { myp <- pchisq(myc,df) } myp load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,hyperlink('ht_variance.htm','Testing Variance - p-value (probability)','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,'p-value (probability)',header=TRUE) a<-table.element(a,myp) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'Conclusion',header=TRUE) if (myp > par4) 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')
|