R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(100 + ,99.94940551 + ,102.0743739 + ,102.0237794 + ,102.6309132 + ,102.8838857 + ,103.0103719 + ,104.3258285 + ,105.211232 + ,104.7305844 + ,104.174045 + ,103.6933974 + ,104.3258285 + ,105.211232 + ,105.3124209 + ,105.5400961 + ,106.7037693 + ,106.6025803 + ,105.008854 + ,104.3258285 + ,104.174045 + ,103.4151278 + ,102.3273463 + ,101.8214015 + ,103.5669112 + ,103.8957754 + ,104.5788009 + ,105.0341513 + ,105.6665823 + ,105.6918796 + ,105.7171768 + ,105.843663 + ,105.7930686 + ,105.3883127 + ,105.9701493 + ,106.5013913 + ,107.1338224 + ,109.3599798 + ,109.3599798 + ,108.4239818 + ,107.9433342 + ,108.0951176 + ,108.3986845 + ,110.5489502 + ,111.8138123 + ,112.5474323 + ,111.6620288 + ,111.3837592 + ,113.1039717 + ,115.1783456 + ,121.0726031 + ,123.0710853 + ,123.3999494 + ,122.9193018 + ,122.3880597 + ,123.5517329 + ,124.9683784 + ,124.8671895 + ,123.2734632 + ,121.9580066 + ,122.4892487 + ,125.6767012 + ,126.7644827 + ,126.4356185 + ,125.3478371 + ,126.0055654 + ,127.4475082 + ,130.5843663 + ,133.0887933 + ,133.3417657 + ,132.8358209 + ,133.7971161 + ,136.6304073 + ,138.8818619 + ,140.6526689 + ,143.9160132 + ,149.0260562 + ,149.3296231 + ,152.6182646 + ,161.5987857 + ,162.6612699 + ,166.1269922 + ,165.9752087 + ,165.4692639 + ,166.9365039 + ,167.4930433 + ,169.2891475 + ,170.2757399 + ,171.5911966 + ,171.3635214 + ,171.641791 + ,170.0986592 + ,168.302555 + ,168.9602833 + ,170.0227675 + ,167.5942322 + ,172.6030863 + ,175.082216 + ,177.7637238 + ,182.4943081 + ,185.4793827 + ,189.1474829 + ,186.5924614 + ,184.4927903 + ,178.62383 + ,170.5793069 + ,166.7088287 + ,163.3442955) > par1 = '12' > #'GNU S' R Code compiled by R2WASP v. 1.0.44 () > #Author: Prof. Dr. P. Wessa > #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/ > #Source of accompanying publication: Office for Research, Development, and Education > #Technical description: Write here your technical program description (don't use hard returns!) > par1 <- as.numeric(par1) > (n <- length(x)) [1] 108 > (np <- floor(n / par1)) [1] 9 > arr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + if (j == par1) j = 0 + } > ari [1] 9 9 9 9 9 9 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 100.0000 104.3258 103.5669 107.1338 113.1040 122.4892 136.6304 166.9365 [2,] 99.9494 105.2112 103.8958 109.3600 115.1783 125.6767 138.8819 167.4930 [3,] 102.0744 105.3124 104.5788 109.3600 121.0726 126.7645 140.6527 169.2891 [4,] 102.0238 105.5401 105.0342 108.4240 123.0711 126.4356 143.9160 170.2757 [5,] 102.6309 106.7038 105.6666 107.9433 123.3999 125.3478 149.0261 171.5912 [6,] 102.8839 106.6026 105.6919 108.0951 122.9193 126.0056 149.3296 171.3635 [7,] 103.0104 105.0089 105.7172 108.3987 122.3881 127.4475 152.6183 171.6418 [8,] 104.3258 104.3258 105.8437 110.5490 123.5517 130.5844 161.5988 170.0987 [9,] 105.2112 104.1740 105.7931 111.8138 124.9684 133.0888 162.6613 168.3026 [10,] 104.7306 103.4151 105.3883 112.5474 124.8672 133.3418 166.1270 168.9603 [11,] 104.1740 102.3273 105.9701 111.6620 123.2735 132.8358 165.9752 170.0228 [12,] 103.6934 101.8214 106.5014 111.3838 121.9580 133.7971 165.4693 167.5942 [,9] [,10] [1,] 172.6031 NA [2,] 175.0822 NA [3,] 177.7637 NA [4,] 182.4943 NA [5,] 185.4794 NA [6,] 189.1475 NA [7,] 186.5925 NA [8,] 184.4928 NA [9,] 178.6238 NA [10,] 170.5793 NA [11,] 166.7088 NA [12,] 163.3443 NA > arr.sd <- array(NA,dim=par1) > arr.range <- array(NA,dim=par1) > arr.iqr <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) + arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) + arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) + } > overall.sd <- sqrt(var(x)) > overall.range <- max(x) - min(x) > overall.iqr <- quantile(x,0.75) - quantile(x,0.25) > postscript(file="/var/www/html/rcomp/tmp/1n9eq1257976264.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.sd,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2iywb1257976264.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.range,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3tfng1257976264.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.iqr,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4zewz1257976264.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(arr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 100.0000 99.9494 102.0744 102.0238 102.6309 102.8839 103.0104 104.3258 [2,] 104.3258 105.2112 105.3124 105.5401 106.7038 106.6026 105.7172 105.8437 [3,] 113.1040 115.1783 121.0726 123.0711 123.3999 122.9193 122.3881 123.5517 [4,] 136.6304 138.8819 140.6527 143.9160 149.0261 149.3296 152.6183 161.5988 [5,] 172.6031 175.0822 177.7637 182.4943 185.4794 189.1475 186.5925 184.4928 [,9] [,10] [,11] [,12] [1,] 104.1740 103.4151 102.3273 101.8214 [2,] 105.7931 105.3883 105.9701 106.5014 [3,] 124.9684 124.8672 123.2735 121.9580 [4,] 162.6613 166.1270 165.9752 163.3443 [5,] 178.6238 170.5793 170.0228 167.5942 $n [1] 9 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 96.09023 97.44515 102.4601 102.8598 101.1102 100.4164 97.68682 [2,] 130.11772 132.91154 139.6851 143.2824 145.6897 145.4222 147.08930 [,8] [,9] [,10] [,11] [,12] [1,] 94.18737 95.0178 92.87815 91.6708 92.02074 [2,] 152.91610 154.9190 156.85623 154.8761 151.89527 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(100, 104.3258285, 113.1039717, 136.6304073, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5x3fv1257976264.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(arr) > names(z) <- c(1:np) > (boxplot(z,notch=TRUE,col='grey',xlab='Block Index',ylab='Value',main='Notched Box Plots - Sequential Blocks')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 99.9494 101.8214 103.5669 107.1338 121.0726 122.4892 136.6304 166.9365 [2,] 102.0491 103.7946 104.8065 108.2469 121.5153 125.8411 142.2843 167.9484 [3,] 102.9471 104.6673 105.6792 109.3600 122.9952 127.1060 150.9739 169.6560 [4,] 104.2499 105.4263 105.8184 111.5229 123.4758 132.9623 164.0653 170.8196 [5,] 105.2112 106.7038 106.5014 112.5474 124.9684 133.7971 166.1270 171.6418 [,9] [,10] [1,] 163.3443 NA [2,] 171.5912 NA [3,] 178.1938 NA [4,] 184.9861 NA [5,] 189.1475 NA $n [1] 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 101.9433 103.9231 105.2177 107.8658 122.1010 123.8580 141.0395 168.3464 [2,] 103.9510 105.4116 106.1408 110.8542 123.8894 130.3540 160.9084 170.9655 [,9] [,10] [1,] 172.0843 NA [2,] 184.3033 NA $out [1] 113.1040 115.1783 $group [1] 5 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(99.94940551, 102.04907665, 102.9471288, 104.24993675, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6l4ym1257976264.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) > names(z) <- list('S.D.','Range','IQR') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 27.70030 65.77283 32.30458 [2,] 28.23511 70.14925 36.85808 [3,] 29.34725 75.41108 44.81407 [4,] 31.04860 81.65950 56.85555 [5,] 31.59592 86.26360 60.73868 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 28.06400 70.16118 35.69308 [2,] 30.63051 80.66098 53.93505 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(27.7003000009886, 28.2351050509087, 29.3472541800025, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1n9eq1257976264.ps tmp/1n9eq1257976264.png") > system("convert tmp/2iywb1257976264.ps tmp/2iywb1257976264.png") > system("convert tmp/3tfng1257976264.ps tmp/3tfng1257976264.png") > system("convert tmp/4zewz1257976264.ps tmp/4zewz1257976264.png") > system("convert tmp/5x3fv1257976264.ps tmp/5x3fv1257976264.png") > system("convert tmp/6l4ym1257976264.ps tmp/6l4ym1257976264.png") > > > proc.time() user system elapsed 1.107 0.851 1.630