R version 2.12.0 (2010-10-15) Copyright (C) 2010 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: i486-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. > x <- c(7992,6114,5965,8460,8323,6333,5675,10090,9035,6976,6459,10896,9978,7466,7199,10977,9412,6341,7784,11911,10079,7721,8197,12038,11963,8033,8618,13625,11734,8895,8727,13974,12583,9525,9662,15490,13839,10047,9788,14978,13045,9489,8741,13149,14106,9998,10034,15081,13266,9997,9027,14324,13149,11209,10332,15354,13800,11786,10550,16114,13255,11403,10269,14009,15847,12967,11328,15814,18626,13219,13818,18062,15722,12111,11702,15589,14852,13612,12380,15501,16322,12157,11124,14621,14035,11159,10944,15824,14378,11816,12233,17344,16812,12181,13275,18458,17375,14609,13323,18327,16053,15070,13806,18245,17461,14999,16022,20564,16372,15854,15115,18207,19488,16644,18631,21093,22212,19762,19403,21227,23176,20823,20647,21336,23458,22003,21647,26416,25226,24723,19945,24040,25034,24885,21168,23541,26019,24657,20599,24534,28717,26138,22968,26577,28660,30430,27356,25454,30194) > par1 = '4' > #'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] 149 > (np <- floor(n / par1)) [1] 37 > 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] 38 37 37 37 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 7992 8323 9035 9978 9412 10079 11963 11734 12583 13839 13045 14106 [2,] 6114 6333 6976 7466 6341 7721 8033 8895 9525 10047 9489 9998 [3,] 5965 5675 6459 7199 7784 8197 8618 8727 9662 9788 8741 10034 [4,] 8460 10090 10896 10977 11911 12038 13625 13974 15490 14978 13149 15081 [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 13266 13149 13800 13255 15847 18626 15722 14852 16322 14035 14378 16812 [2,] 9997 11209 11786 11403 12967 13219 12111 13612 12157 11159 11816 12181 [3,] 9027 10332 10550 10269 11328 13818 11702 12380 11124 10944 12233 13275 [4,] 14324 15354 16114 14009 15814 18062 15589 15501 14621 15824 17344 18458 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [1,] 17375 16053 17461 16372 19488 22212 23176 23458 25226 25034 26019 28717 [2,] 14609 15070 14999 15854 16644 19762 20823 22003 24723 24885 24657 26138 [3,] 13323 13806 16022 15115 18631 19403 20647 21647 19945 21168 20599 22968 [4,] 18327 18245 20564 18207 21093 21227 21336 26416 24040 23541 24534 26577 [,37] [,38] [1,] 28660 30194 [2,] 30430 NA [3,] 27356 NA [4,] 25454 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/rcomp/tmp/18puf1304952262.ps",horizontal=F,onefile=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/rcomp/tmp/2l6kj1304952262.ps",horizontal=F,onefile=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/rcomp/tmp/38fpb1304952262.ps",horizontal=F,onefile=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/rcomp/tmp/4z34a1304952262.ps",horizontal=F,onefile=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] [1,] 7992 6114 5675 8460 [2,] 13045 9525 9027 14009 [3,] 15287 12111 11328 15814 [4,] 19488 15854 16022 20564 [5,] 28717 24885 22968 26577 $n [1] 38 37 37 37 $conf [,1] [,2] [,3] [,4] [1,] 13635.60 10467.04 9511.046 14111.34 [2,] 16938.40 13754.96 13144.954 17516.66 $out [1] 30194 26138 30430 27356 $group [1] 1 2 2 3 $names [1] "1" "2" "3" "4" > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5yh7d1304952262.ps",horizontal=F,onefile=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] [,9] [1,] 5965.0 5675.0 6459.0 7199.0 6341.0 7721.0 8033.0 8727.0 9525.0 [2,] 6039.5 6004.0 6717.5 7332.5 7062.5 7959.0 8325.5 8811.0 9593.5 [3,] 7053.0 7328.0 8005.5 8722.0 8598.0 9138.0 10290.5 10314.5 11122.5 [4,] 8226.0 9206.5 9965.5 10477.5 10661.5 11058.5 12794.0 12854.0 14036.5 [5,] 8460.0 10090.0 10896.0 10977.0 11911.0 12038.0 13625.0 13974.0 15490.0 [,10] [,11] [,12] [,13] [,14] [,15] [,16] [,17] [,18] [,19] [1,] 9788.0 8741 9998.0 9027.0 10332.0 10550 10269 11328.0 13219.0 11702.0 [2,] 9917.5 9115 10016.0 9512.0 10770.5 11168 10836 12147.5 13518.5 11906.5 [3,] 11943.0 11267 12070.0 11631.5 12179.0 12793 12329 14390.5 15940.0 13850.0 [4,] 14408.5 13097 14593.5 13795.0 14251.5 14957 13632 15830.5 18344.0 15655.5 [5,] 14978.0 13149 15081.0 14324.0 15354.0 16114 14009 15847.0 18626.0 15722.0 [,20] [,21] [,22] [,23] [,24] [,25] [,26] [,27] [,28] [1,] 12380.0 11124.0 10944.0 11816.0 12181.0 13323 13806.0 14999.0 15115.0 [2,] 12996.0 11640.5 11051.5 12024.5 12728.0 13966 14438.0 15510.5 15484.5 [3,] 14232.0 13389.0 12597.0 13305.5 15043.5 15992 15561.5 16741.5 16113.0 [4,] 15176.5 15471.5 14929.5 15861.0 17635.0 17851 17149.0 19012.5 17289.5 [5,] 15501.0 16322.0 15824.0 17344.0 18458.0 18327 18245.0 20564.0 18207.0 [,29] [,30] [,31] [,32] [,33] [,34] [,35] [,36] [,37] [1,] 16644.0 19403.0 20647.0 21647.0 19945.0 21168.0 20599.0 22968.0 25454 [2,] 17637.5 19582.5 20735.0 21825.0 21992.5 22354.5 22566.5 24553.0 26405 [3,] 19059.5 20494.5 21079.5 22730.5 24381.5 24213.0 24595.5 26357.5 28008 [4,] 20290.5 21719.5 22256.0 24937.0 24974.5 24959.5 25338.0 27647.0 29545 [5,] 21093.0 22212.0 23176.0 26416.0 25226.0 25034.0 26019.0 28717.0 30430 [,38] [1,] 30194 [2,] 30194 [3,] 30194 [4,] 30194 [5,] 30194 $n [1] 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 5325.665 4798.025 5439.58 6237.45 5754.79 6689.395 6760.385 7120.53 [2,] 8780.335 9857.975 10571.42 11206.55 11441.21 11586.605 13820.615 13508.47 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [,16] [1,] 7612.53 8395.11 8121.22 8453.775 8247.93 9429.01 9799.69 10120.16 [2,] 14632.47 15490.89 14412.78 15686.225 15015.07 14928.99 15786.31 14537.84 [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [1,] 11480.93 12127.85 10888.29 12509.41 10362.51 9533.38 10274.67 11166.97 [2,] 17300.07 19752.15 16811.71 15954.59 16415.49 15660.62 16336.33 18920.03 [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32] [1,] 12922.85 13419.81 13974.92 14687.05 16963.63 18806.27 19877.91 20272.02 [2,] 19061.15 17703.19 19508.08 17538.95 21155.37 22182.73 22281.09 25188.98 [,33] [,34] [,35] [,36] [,37] [,38] [1,] 22025.72 22155.05 22406.01 23913.24 25527.4 30194 [2,] 26737.28 26270.95 26784.99 28801.76 30488.6 30194 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" "15" [16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30" [31] "31" "32" "33" "34" "35" "36" "37" NA Warning message: In bxp(list(stats = c(5965, 6039.5, 7053, 8226, 8460, 5675, 6004, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6x9gg1304952262.ps",horizontal=F,onefile=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,] 4706.780 18117.0 6201.50 [2,] 5071.848 19899.0 6265.25 [3,] 5717.301 21941.5 6442.00 [4,] 6144.566 23259.0 6775.00 [5,] 6291.447 24316.0 6995.00 $n [1] 4 4 4 $conf [,1] [,2] [,3] [1,] 4869.854 19287.1 6039.297 [2,] 6564.748 24595.9 6844.703 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(4706.7797449602, 5071.8480690672, 5717.30099367094, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/18puf1304952262.ps tmp/18puf1304952262.png",intern=TRUE)) character(0) > try(system("convert tmp/2l6kj1304952262.ps tmp/2l6kj1304952262.png",intern=TRUE)) character(0) > try(system("convert tmp/38fpb1304952262.ps tmp/38fpb1304952262.png",intern=TRUE)) character(0) > try(system("convert tmp/4z34a1304952262.ps tmp/4z34a1304952262.png",intern=TRUE)) character(0) > try(system("convert tmp/5yh7d1304952262.ps tmp/5yh7d1304952262.png",intern=TRUE)) character(0) > try(system("convert tmp/6x9gg1304952262.ps tmp/6x9gg1304952262.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.340 1.000 2.332