R version 2.8.0 (2008-10-20) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0 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. > x <- c(0.084299365852926 + ,-15.5692124815993 + ,-6.9953619930182 + ,-2.07575919037298 + ,6.06472185178652 + ,3.10063183929146 + ,-7.27872604143293 + ,11.9305592423578 + ,-15.1800926797913 + ,-4.96965957564265 + ,-3.06184332688462 + ,5.17876413221274 + ,-2.29638814463019 + ,-1.85853257792512 + ,-8.74637895628288 + ,-8.02346934029272 + ,4.80319999595127 + ,3.90173526693944 + ,-11.3069080265042 + ,4.26225040152088 + ,-6.88628010104842 + ,4.12797896412461 + ,4.24824928932737 + ,-5.91811547446539 + ,-2.62127541719435 + ,-5.11947651478644 + ,-6.19169292552164 + ,1.57682830128841 + ,6.93381888926023 + ,1.30435330636300 + ,0.701451547474772 + ,18.1487693268603 + ,-13.6251962127511 + ,-8.9684198944174 + ,-3.48024416476534 + ,0.133268755660886 + ,0.156896176204026 + ,-1.28412130463800 + ,-8.4586972345603 + ,-2.91518278408049 + ,0.150465578819962 + ,3.86877987666889 + ,7.70993797103159 + ,-4.11650976301151 + ,-2.74475852656699 + ,3.19286810860467 + ,-6.2212027713738 + ,1.25489355791219 + ,4.23585013027417 + ,0.153440583363572 + ,-2.05719003274709 + ,2.14422869581576 + ,-1.68156215590428 + ,5.0318450949173 + ,2.24823311061623 + ,-6.58372573676435 + ,-2.85420176875137 + ,5.0379752614482 + ,4.43711208741108 + ,-0.675511401722624 + ,-2.50684702275299 + ,2.35909448398834 + ,-6.31849429304567 + ,-3.81949862816358 + ,2.17307118170725 + ,4.13663986179168 + ,-6.96229654560136 + ,5.28008601863276 + ,4.2774558989149) > par2 = '0' > par1 = '0.95' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: Wessa P., (2007), Testing Population Mean with unknown Variance (confidence intervals) (v1.0.1) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_hypothesismeanu.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description > par1 <- as.numeric(par1) > par2 <- as.numeric(par2) > len <- length(x) > df <- len - 1 > sd <- sd(x) > mx <- mean(x) > > #Note: the /var/www/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab > load(file="/var/www/html/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/www/html/rcomp/tmp/1pa3o1230382690.tab") > > > > proc.time() user system elapsed 0.207 0.035 0.246