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 + ,108.1560276 + ,114.0150276 + ,102.1880309 + ,110.3672031 + ,96.8602511 + ,94.1944583 + ,99.51621961 + ,94.06333487 + ,97.5541476 + ,78.15062422 + ,81.2434643 + ,92.36262465 + ,96.06324371 + ,114.0523777 + ,110.6616666 + ,104.9171949 + ,90.00187193 + ,95.7008067 + ,86.02741157 + ,84.85287668 + ,100.04328 + ,80.91713823 + ,74.06539709 + ,77.30281369 + ,97.23043249 + ,90.75515676 + ,100.5614455 + ,92.01293267 + ,99.24012138 + ,105.8672755 + ,90.9920463 + ,93.30624423 + ,91.17419413 + ,77.33295039 + ,91.1277721 + ,85.01249943 + ,83.90390242 + ,104.8626302 + ,110.9039108 + ,95.43714373 + ,111.6238727 + ,108.8925403 + ,96.17511682 + ,101.9740205 + ,99.11953031 + ,86.78158147 + ,118.4195003 + ,118.7441447 + ,106.5296192 + ,134.7772694 + ,104.6778714 + ,105.2954304 + ,139.4139849 + ,103.6060491 + ,99.78182974 + ,103.4610301 + ,120.0594945 + ,96.71377168 + ,107.1308929 + ,105.3608372 + ,111.6942359 + ,132.0519998 + ,126.8037879 + ,154.4824253 + ,141.5570984 + ,109.9506882 + ,127.904198 + ,133.0888617 + ,120.0796299 + ,117.5557142 + ,143.0362309 + ,159.982927 + ,128.5991124 + ,149.7373327 + ,126.8169313 + ,140.9639674 + ,137.6691981 + ,117.9402337 + ,122.3095247 + ,127.7804207 + ,136.1677176 + ,116.2405856 + ,123.1576893 + ,116.3400234 + ,108.6119282 + ,125.8982264 + ,112.8003105 + ,107.5182447 + ,135.0955413 + ,115.5096488 + ,115.8640759 + ,104.5883906 + ,163.7213386 + ,113.4482275 + ,98.0428844 + ,116.7868521 + ,126.5330444 + ,113.0336597 + ,124.3392163 + ,109.8298759 + ,124.4434777 + ,111.5039454 + ,102.0350019 + ,116.8726598 + ,112.2073122 + ,101.1513902 + ,124.4255108) > 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] [1,] 100.00000 92.36262 77.30281 85.01250 118.74414 105.3608 159.9829 [2,] 108.15603 96.06324 97.23043 83.90390 106.52962 111.6942 128.5991 [3,] 114.01503 114.05238 90.75516 104.86263 134.77727 132.0520 149.7373 [4,] 102.18803 110.66167 100.56145 110.90391 104.67787 126.8038 126.8169 [5,] 110.36720 104.91719 92.01293 95.43714 105.29543 154.4824 140.9640 [6,] 96.86025 90.00187 99.24012 111.62387 139.41398 141.5571 137.6692 [7,] 94.19446 95.70081 105.86728 108.89254 103.60605 109.9507 117.9402 [8,] 99.51622 86.02741 90.99205 96.17512 99.78183 127.9042 122.3095 [9,] 94.06333 84.85288 93.30624 101.97402 103.46103 133.0889 127.7804 [10,] 97.55415 100.04328 91.17419 99.11953 120.05949 120.0796 136.1677 [11,] 78.15062 80.91714 77.33295 86.78158 96.71377 117.5557 116.2406 [12,] 81.24346 74.06540 91.12777 118.41950 107.13089 143.0362 123.1577 [,8] [,9] [,10] [1,] 116.34002 116.7869 NA [2,] 108.61193 126.5330 NA [3,] 125.89823 113.0337 NA [4,] 112.80031 124.3392 NA [5,] 107.51824 109.8299 NA [6,] 135.09554 124.4435 NA [7,] 115.50965 111.5039 NA [8,] 115.86408 102.0350 NA [9,] 104.58839 116.8727 NA [10,] 163.72134 112.2073 NA [11,] 113.44823 101.1514 NA [12,] 98.04288 124.4255 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/180vx1257855511.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/2w9qh1257855511.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/3sx9y1257855511.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/494u41257855511.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] [1,] 77.30281 83.90390 90.75516 100.5614 104.9172 90.00187 94.19446 [2,] 92.36262 97.23043 113.03366 104.6779 104.9172 99.24012 103.60605 [3,] 105.36084 108.15603 114.05238 110.9039 107.5182 124.44348 108.89254 [4,] 116.78685 111.69424 132.05200 124.3392 110.3672 137.66920 111.50395 [5,] 118.74414 128.59911 149.73733 126.8169 110.3672 141.55710 117.94023 [,8] [,9] [,10] [,11] [,12] [1,] 86.02741 84.85288 91.17419 77.33295 74.06540 [2,] 96.17512 94.06333 99.11953 80.91714 91.12777 [3,] 99.78183 103.46103 112.20731 96.71377 107.13089 [4,] 115.86408 116.87266 120.07963 113.44823 123.15769 [5,] 127.90420 133.08886 136.16772 117.55571 143.03623 $n [1] 9 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 92.49741 100.5384 104.0361 100.5489 104.6479 104.2042 104.7330 89.41231 [2,] 118.22426 115.7736 124.0687 121.2589 110.3886 144.6828 113.0521 110.15135 [,9] [,10] [,11] [,12] [1,] 91.44812 101.1683 79.58073 90.2618 [2,] 115.47394 123.2463 113.84681 124.0000 $out [1] 159.98293 92.01293 95.43714 154.48243 140.96397 163.72134 $group [1] 1 5 5 5 5 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(77.30281369, 92.36262465, 105.3608372, 116.7868521, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/51gu31257855511.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] [1,] 78.15062 74.06540 90.75516 83.90390 96.71377 105.3608 116.2406 [2,] 94.12890 85.44014 90.87360 91.10936 103.53354 114.6250 122.7336 [3,] 98.53518 94.03172 91.59356 100.54678 105.91252 127.3540 128.1898 [4,] 105.17203 102.48024 98.23528 109.89823 119.40182 137.3230 139.3166 [5,] 114.01503 114.05238 105.86728 118.41950 139.41398 154.4824 159.9829 [,8] [,9] [,10] [1,] 98.04288 101.1514 NA [2,] 108.06509 110.6669 NA [3,] 114.47894 114.9103 NA [4,] 121.11912 124.3824 NA [5,] 135.09554 126.5330 NA $n [1] 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 93.49834 86.25961 88.23585 91.97705 98.6749 117.0013 120.6262 108.5249 [2,] 103.57203 101.80382 94.95127 109.11650 113.1502 137.7067 135.7534 120.4330 [,9] [,10] [1,] 108.6545 NA [2,] 121.1660 NA $out [1] 77.30281 77.33295 163.72134 $group [1] 3 3 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(78.15062422, 94.128896585, 98.535183605, 105.17202925, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6349s1257855511.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,] 8.132721 23.74578 5.450008 [2,] 14.276057 41.04978 16.741072 [3,] 17.075445 49.89561 20.324529 [4,] 21.623161 65.72016 28.227072 [5,] 24.368679 82.68011 38.429077 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 13.72438 38.64328 15.08569 [2,] 20.42651 61.14794 25.56337 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(8.13272107100048, 14.2760565917519, 17.0754452349667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > system("convert tmp/180vx1257855511.ps tmp/180vx1257855511.png") > system("convert tmp/2w9qh1257855511.ps tmp/2w9qh1257855511.png") > system("convert tmp/3sx9y1257855511.ps tmp/3sx9y1257855511.png") > system("convert tmp/494u41257855511.ps tmp/494u41257855511.png") > system("convert tmp/51gu31257855511.ps tmp/51gu31257855511.png") > system("convert tmp/6349s1257855511.ps tmp/6349s1257855511.png") > > > proc.time() user system elapsed 1.074 0.806 1.430