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(-0.0375732130599521 + ,-1.42842903073985 + ,0.344188230160782 + ,-0.163082724149866 + ,-0.0816939481183535 + ,0.252095350579819 + ,0.289469636284585 + ,-0.274908034110628 + ,-0.457959047417949 + ,-0.120161814446042 + ,0.115277151044634 + ,0.52594879230336 + ,-0.574167421833416 + ,-0.977009688083828 + ,0.28705772607985 + ,0.868376961764599 + ,-1.01003929038058 + ,1.06363782814538 + ,0.657427733182704 + ,0.566693320822928 + ,-1.64827620985414 + ,1.21153674122654 + ,-2.01659873949968 + ,1.08949523651952 + ,-1.10332128763400 + ,0.393057450103926 + ,-0.366241722302763 + ,0.0740797820367478 + ,0.899595501277396 + ,-0.559421942167576 + ,-1.79385744860224 + ,0.280148270173554 + ,0.603531479858657 + ,-0.196938655171879 + ,-0.0491476438699476 + ,1.27471298429356 + ,-0.0205523160493216 + ,0.638508206542742 + ,-0.906094943253174 + ,-0.711904745736901 + ,1.03912815618083 + ,0.797303418582777 + ,-0.57685992535319 + ,1.06628808049365 + ,-0.00187118162893192 + ,-0.993835979441977 + ,1.53717318164927 + ,-0.552769095442449 + ,0.835516821617911 + ,-0.767777920380284 + ,-0.464942973832369 + ,0.190103523651543 + ,-1.10193377119234 + ,-3.13783784070054 + ,-2.74591790271741 + ,-1.20624604013675 + ,0.661408396700863 + ,-0.0544913767182588 + ,-0.782616089336004 + ,0.122339478528019 + ,0.487830064783113) > 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] 61 > (np <- floor(n / par1)) [1] 5 > arr <- array(NA,dim=c(par1,np+1)) > darr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > dx <- diff(x) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + darr[j,ari[j]] <- dx[i] + if (j == par1) j = 0 + } > ari [1] 6 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.03757321 -0.5741674 -1.10332129 -0.020552316 0.83551682 0.4878301 [2,] -1.42842903 -0.9770097 0.39305745 0.638508207 -0.76777792 NA [3,] 0.34418823 0.2870577 -0.36624172 -0.906094943 -0.46494297 NA [4,] -0.16308272 0.8683770 0.07407978 -0.711904746 0.19010352 NA [5,] -0.08169395 -1.0100393 0.89959550 1.039128156 -1.10193377 NA [6,] 0.25209535 1.0636378 -0.55942194 0.797303419 -3.13783784 NA [7,] 0.28946964 0.6574277 -1.79385745 -0.576859925 -2.74591790 NA [8,] -0.27490803 0.5666933 0.28014827 1.066288080 -1.20624604 NA [9,] -0.45795905 -1.6482762 0.60353148 -0.001871182 0.66140840 NA [10,] -0.12016181 1.2115367 -0.19693866 -0.993835979 -0.05449138 NA [11,] 0.11527715 -2.0165987 -0.04914764 1.537173182 -0.78261609 NA [12,] 0.52594879 1.0894952 1.27471298 -0.552769095 0.12233948 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -1.39085582 -0.40284227 1.4963787 0.6590605 -1.6032947 NA [2,] 1.77261726 1.26406741 -0.7592992 -1.5446031 0.3028349 NA [3,] -0.50727095 0.58131924 0.4403215 0.1941902 0.6550465 NA [4,] 0.08138878 -1.87841625 0.8255157 1.7510329 -1.2920373 NA [5,] 0.33378930 2.07367712 -1.4590174 -0.2418247 -2.0359041 NA [6,] 0.03737429 -0.40621009 -1.2344355 -1.3741633 0.3919199 NA [7,] -0.56437767 -0.09073441 2.0740057 1.6431480 1.5396719 NA [8,] -0.18305101 -2.21496953 0.3233832 -1.0681593 1.8676544 NA [9,] 0.33779723 2.85981295 -0.8004701 -0.9919648 -0.7158998 NA [10,] 0.23543897 -3.22813548 0.1477910 2.5310092 -0.7281247 NA [11,] 0.41067164 3.10609398 1.3238606 -2.0899423 0.9049556 NA [12,] -1.10011621 -2.19281652 -1.2952653 1.3882859 0.3654906 NA > arr.mean <- array(NA,dim=par1) > arr.median <- array(NA,dim=par1) > arr.midrange <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.mean[j] <- mean(arr[j,],na.rm=TRUE) + arr.median[j] <- median(arr[j,],na.rm=TRUE) + arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2 + } > overall.mean <- mean(x) > overall.median <- median(x) > overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2 > postscript(file="/var/www/html/rcomp/tmp/1zpq11260460877.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.mean,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/2dpdb1260460877.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.median,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3r3qm1260460877.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index') > mtext(paste('#blocks = ',np)) > abline(overall.midrange,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4ljip1260460877.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] [1,] -1.10332129 -1.4284290 -0.9060949 -0.16308272 -1.10193377 -0.5594219 [2,] -0.57416742 -0.9770097 -0.4649430 -0.16308272 -1.01003929 -0.5594219 [3,] -0.02906276 -0.7677779 -0.3662417 0.07407978 -0.08169395 0.2520954 [4,] 0.48783006 0.3930575 0.2870577 0.19010352 0.89959550 0.7973034 [5,] 0.83551682 0.6385082 0.3441882 0.19010352 1.03912816 1.0636378 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -2.7459179 -1.2062460 -1.648276210 -0.19693866 -2.01659874 -0.5527691 [2,] -1.7938574 -0.2749080 -0.457959047 -0.19693866 -0.78261609 0.1223395 [3,] -0.5768599 0.2801483 -0.001871182 -0.12016181 -0.04914764 0.5259488 [4,] 0.2894696 0.5666933 0.603531480 -0.05449138 0.11527715 1.0894952 [5,] 0.6574277 1.0662881 0.661408397 -0.05449138 0.11527715 1.2747130 $n [1] 6 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.7140854 -1.7358639 -0.8976035 -0.1754807 -1.431037 -0.7065634 [2,] 0.6559599 0.2003081 0.1651201 0.3236403 1.267649 1.2107541 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -2.0489336 -0.3145251 -0.7519177 -0.22081470 -0.6835967 -0.1574410 [2,] 0.8952137 0.8748217 0.7481754 -0.01950893 0.5853015 1.2093386 $out [1] 0.8683770 -0.7119047 -3.1378378 1.2115367 -0.9938360 1.5371732 $group [1] 4 4 6 10 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1.103321287634, -0.574167421833416, -0.0290627645546368, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/56lmw1260460877.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(t(darr)) > names(z) <- c(1:par1) > (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries')) $stats [,1] [,2] [,3] [,4] [,5] [,6] [1,] -1.6032947 -1.5446031 0.1941902 -1.87841625 -2.0359041 -1.37416334 [2,] -1.3908558 -0.7592992 0.1941902 -1.29203729 -1.4590174 -1.23443551 [3,] -0.4028423 0.3028349 0.4403215 0.08138878 -0.2418247 -0.40621009 [4,] 0.6590605 1.2640674 0.5813192 0.82551572 0.3337893 0.03737429 [5,] 1.4963787 1.7726173 0.6550465 1.75103290 2.0736771 0.39191994 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -0.56437767 -2.2149695 -0.9919648 -0.7281247 0.4106716 -2.1928165 [2,] -0.09073441 -1.0681593 -0.8004701 -0.7281247 0.4106716 -1.2952653 [3,] 1.53967186 -0.1830510 -0.7158998 0.1477910 0.9049556 -1.1001162 [4,] 1.64314801 0.3233832 0.3377972 0.2354390 1.3238606 0.3654906 [5,] 2.07400572 1.8676544 0.3377972 0.2354390 1.3238606 1.3882859 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1.851308 -1.126871 0.1667771 -1.414869 -1.508617 -1.3048677 0.3145149 [2,] 1.045623 1.732541 0.7138659 1.577646 1.024968 0.4924475 2.7648288 [,8] [,9] [,10] [,11] [,12] [1,] -1.1663114 -1.52019663 -0.5330607 0.2596985 -2.27360214 [2,] 0.8002094 0.08839708 0.8286427 1.5502126 0.07336971 $out [1] -0.507271 2.859813 -3.228135 2.531009 3.106094 -2.089942 $group [1] 3 9 10 10 11 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1.60329474199819, -1.39085581767990, -0.402842266250412, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/617og1260460877.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] [1,] -0.45795905 -2.0165987 -1.79385745 -0.99383598 -2.7459179 0.4878301 [2,] -0.21899538 -0.9935245 -0.46283183 -0.64438234 -1.1540899 0.4878301 [3,] -0.05963358 0.4268755 0.01246607 -0.01121175 -0.6163604 0.4878301 [4,] 0.27078249 0.9660074 0.49829446 0.91821579 0.1562215 0.4878301 [5,] 0.52594879 1.2115367 1.27471298 1.53717318 0.8355168 0.4878301 $n [1] 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.2830246 -0.4668801 -0.4259101 -0.7239232 -1.21400228 0.4878301 [2,] 0.1637574 1.3206312 0.4508422 0.7014997 -0.01871862 0.4878301 $out [1] -1.428429 -3.137838 $group [1] 1 5 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(-0.457959047417949, -0.218995379130247, -0.0596335805891528, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7q1n21260460877.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.mean,arr.median,arr.midrange)) > names(z) <- list('mean','median','midrange') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] -0.42833020 -0.76777792 -0.33366947 [2,] -0.27801353 -0.24320177 -0.20884557 [3,] -0.11867227 -0.03910520 -0.04743205 [4,] 0.01036817 0.16308757 0.09586348 [5,] 0.08639512 0.52594879 0.14589264 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -0.25020510 -0.2244165 -0.1864119 [2,] 0.01286056 0.1462061 0.0915478 $out [1] -0.8339476 0.4919455 -0.7521939 0.6059174 $group [1] 1 1 3 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-0.428330196511459, -0.278013532557187, -0.118672268915456, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1zpq11260460877.ps tmp/1zpq11260460877.png") > system("convert tmp/2dpdb1260460877.ps tmp/2dpdb1260460877.png") > system("convert tmp/3r3qm1260460877.ps tmp/3r3qm1260460877.png") > system("convert tmp/4ljip1260460877.ps tmp/4ljip1260460877.png") > system("convert tmp/56lmw1260460877.ps tmp/56lmw1260460877.png") > system("convert tmp/617og1260460877.ps tmp/617og1260460877.png") > system("convert tmp/7q1n21260460877.ps tmp/7q1n21260460877.png") > > > proc.time() user system elapsed 1.199 0.955 1.535