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(78973,46146,46492,60656,21898,36555,74680,22807,61282,37981,41553,45081,38557,51641,30658,52924,79256,53462,68950,53639,67819,48333,28001,51665,39019,46221,65792,39858,19574,41829,78688,36781,44314,24874,56911,37048,48426,33388,26998,46502,41507,40001,33144,29501,43059,43249,29272,49821,98341,44372,42448,5950,64839,32551,30767,62046,71930,67328,67253,35373,85544,88087,30621,50580,49670,25456,69245,43787,53638,35683,38008,18801,44324,51408,53880,55708,63858,183643,35660,41664,29883,62047,33321,46553,56622,15430,49379,58215,38253,77786,21331,55292,30105,37651,59370,46216,73122,93927,55935,93308,74344,78094,25625,43750,28995,47336,57582,60875,165877,32984,61638,36367,1168,40530,21427,15024,39088,855,80455,14116,43915,76705,40112,41821,8773,52045,51491,53470,53211,63091,131634,41745,23656,51442,54574,35708,66627,39585,50029,25266,34860,62759,62307,37238,42452,59820,75075,97567,0,6023,0,0,0,0,42420,31116,0,0,1644,6179,3926,23238,0,38818) > 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] 164 > (np <- floor(n / par1)) [1] 13 > 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] 14 14 14 14 14 14 14 14 13 13 13 13 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 78973 38557 39019 48426 98341 85544 44324 56622 73122 165877 43915 [2,] 46146 51641 46221 33388 44372 88087 51408 15430 93927 32984 76705 [3,] 46492 30658 65792 26998 42448 30621 53880 49379 55935 61638 40112 [4,] 60656 52924 39858 46502 5950 50580 55708 58215 93308 36367 41821 [5,] 21898 79256 19574 41507 64839 49670 63858 38253 74344 1168 8773 [6,] 36555 53462 41829 40001 32551 25456 183643 77786 78094 40530 52045 [7,] 74680 68950 78688 33144 30767 69245 35660 21331 25625 21427 51491 [8,] 22807 53639 36781 29501 62046 43787 41664 55292 43750 15024 53470 [9,] 61282 67819 44314 43059 71930 53638 29883 30105 28995 39088 53211 [10,] 37981 48333 24874 43249 67328 35683 62047 37651 47336 855 63091 [11,] 41553 28001 56911 29272 67253 38008 33321 59370 57582 80455 131634 [12,] 45081 51665 37048 49821 35373 18801 46553 46216 60875 14116 41745 [,12] [,13] [,14] [1,] 23656 42452 0 [2,] 51442 59820 0 [3,] 54574 75075 1644 [4,] 35708 97567 6179 [5,] 66627 0 3926 [6,] 39585 6023 23238 [7,] 50029 0 0 [8,] 25266 0 38818 [9,] 34860 0 NA [10,] 62759 0 NA [11,] 62307 42420 NA [12,] 37238 31116 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] -32827 13084 7202 -15038 -53969 2543 7084 -41192 20805 -132893 [2,] 346 -20983 19571 -6390 -1924 -57466 2472 33949 -37992 28654 [3,] 14164 22266 -25934 19504 -36498 19959 1828 8836 37373 -25271 [4,] -38758 26332 -20284 -4995 58889 -910 8150 -19962 -18964 -35199 [5,] 14657 -25794 22255 -1506 -32288 -24214 119785 39533 3750 39362 [6,] 38125 15488 36859 -6857 -1784 43789 -147983 -56455 -52469 -19103 [7,] -51873 -15311 -41907 -3643 31279 -25458 6004 33961 18125 -6403 [8,] 38475 14180 7533 13558 9884 9851 -11781 -25187 -14755 24064 [9,] -23301 -19486 -19440 190 -4602 -17955 32164 7546 18341 -38233 [10,] 3572 -20332 32037 -13977 -75 2325 -28726 21719 10246 79600 [11,] 3528 23664 -19863 20549 -31880 -19207 13232 -13154 3293 -66339 [12,] -6524 -12646 11378 48520 50171 25523 10069 26906 105002 29799 [,11] [,12] [,13] [,14] [1,] 32790 27786 17368 0 [2,] -36593 3132 15255 1644 [3,] 1709 -18866 22492 4535 [4,] -33048 30919 -97567 -2253 [5,] 43272 -27042 6023 19312 [6,] -554 10444 -6023 -23238 [7,] 1979 -24763 0 38818 [8,] -259 9594 0 NA [9,] 9880 27899 0 NA [10,] 68543 -452 42420 NA [11,] -89889 -25069 -11304 NA [12,] -18089 5214 -31116 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/rcomp/tmp/1rb7m1321350306.ps",horizontal=F,onefile=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/rcomp/tmp/2yheq1321350306.ps",horizontal=F,onefile=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/rcomp/tmp/3xa4w1321350306.ps",horizontal=F,onefile=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/rcomp/tmp/4hfrr1321350306.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] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 0 0.0 1644.0 5950 0 6023.0 0 0 0 0 28001 [2,] 39019 33388.0 30658.0 36367 8773 32551.0 21427 25266 30105 35683 38008 [3,] 46375 48814.5 47935.5 48541 39880 40265.5 34402 40241 43059 43249 56911 [4,] 78973 59820.0 55935.0 58215 64839 53462.0 68950 53470 53638 62047 62307 [5,] 98341 93927.0 75075.0 60656 79256 78094.0 78688 62046 71930 67328 80455 [,12] [1,] 18801 [2,] 35373 [3,] 41745 [4,] 46553 [5,] 60875 $n [1] 14 14 14 14 14 14 14 14 13 13 13 13 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 29503.51 37652.99 37261.71 39315.18 16204.86 31435.35 14334.33 28331.22 [2,] 63246.49 59976.01 58609.29 57766.82 63555.14 49095.65 54469.67 52150.78 [,9] [,10] [,11] [,12] [1,] 32746.53 31695.95 46262.86 36845.78 [2,] 53371.47 54802.05 67559.14 46644.22 $out [1] 165877 93308 97567 183643 131634 14116 $group [1] 1 4 4 6 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(0, 39019, 46375, 78973, 98341, 0, 33388, 48814.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5sata1321350306.ps",horizontal=F,onefile=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] [,7] [,8] [,9] [1,] -53969.0 -57466 -36498.0 -38758.0 -32288 -56455.0 -51873.0 -14755 -38233 [2,] -32827.0 -20983 -18866.0 -33048.0 -24214 -23238.0 -24763.0 -259 -19440 [3,] 4813.5 995 6685.5 -11979.5 10340 -3903.5 -1821.5 9594 0 [4,] 17368.0 15255 19959.0 8150.0 39362 15488.0 18125.0 13558 9880 [5,] 32790.0 33949 37373.0 58889.0 119785 43789.0 38818.0 24064 32164 [,10] [,11] [,12] [1,] -28726 -66339 -31116 [2,] -452 -25069 -6524 [3,] 3572 -13154 11378 [4,] 32037 3528 29799 [5,] 79600 23664 50171 $n [1] 14 14 14 14 14 14 14 13 13 13 13 13 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -16382.48 -14307.32 -9709.24 -29376.293 -16506.41 -20256.43 -19931.93 [2,] 26009.48 16297.32 23080.24 5417.293 37186.41 12449.43 16288.93 [,8] [,9] [,10] [,11] [,12] [1,] 3539.209 -12848.41 -10665.11 -25685.5816 -4539.216 [2,] 15648.791 12848.41 17809.11 -622.4184 27295.216 $out [1] -132893 -97567 -147983 38475 -25187 -89889 105002 $group [1] 1 4 6 8 8 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-53969, -32827, 4813.5, 17368, 32790, -57466, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/658oq1321350306.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,] 21898.0 28001.0 19574.0 26998.0 5950.0 18801.0 29883.0 15430.0 25625.0 [2,] 37268.0 43445.0 36914.5 31322.5 33962.0 33152.0 38662.0 33878.0 45543.0 [3,] 45613.5 52294.5 40843.5 40754.0 53209.0 46728.5 48980.5 47797.5 59228.5 [4,] 60969.0 60729.0 51566.0 44875.5 67290.5 61441.5 58877.5 57418.5 76219.0 [5,] 78973.0 79256.0 65792.0 49821.0 98341.0 88087.0 63858.0 77786.0 93927.0 [,10] [,11] [,12] [,13] [,14] [1,] 855.0 40112.0 23656.0 0.0 0.0 [2,] 14570.0 41783.0 35284.0 0.0 0.0 [3,] 34675.5 51768.0 44807.0 18569.5 2785.0 [4,] 51084.0 58280.5 58440.5 51136.0 14708.5 [5,] 80455.0 76705.0 66627.0 97567.0 23238.0 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 34803.31 44411.15 34160.85 34572.39 38007.65 33825.47 39760.07 37060.52 [2,] 56423.69 60177.85 47526.15 46935.61 68410.35 59631.53 58200.93 58534.48 [,9] [,10] [,11] [,12] [,13] [,14] [1,] 45236.97 18021.22 44243.38 34245.16 -4753.973 -5431.379 [2,] 73220.03 51329.78 59292.62 55368.84 41892.973 11001.379 $out [1] 78688 183643 165877 8773 131634 38818 $group [1] 3 7 10 11 11 14 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" NA Warning message: In bxp(list(stats = c(21898, 37268, 45613.5, 60969, 78973, 28001, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/79eqp1321350306.ps",horizontal=F,onefile=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,] 37274.64 34402.00 38033.50 [2,] 39869.65 40253.25 41417.25 [3,] 44155.97 43154.00 43919.50 [4,] 50798.89 48238.25 48139.50 [5,] 59916.29 56911.00 50157.50 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 39171.07 39511.99 40853.44 [2,] 49140.87 46796.01 46985.56 $out [1] 58693.75 $group [1] 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(37274.6428571429, 39869.6510989011, 44155.9725274725, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1rb7m1321350306.ps tmp/1rb7m1321350306.png",intern=TRUE)) character(0) > try(system("convert tmp/2yheq1321350306.ps tmp/2yheq1321350306.png",intern=TRUE)) character(0) > try(system("convert tmp/3xa4w1321350306.ps tmp/3xa4w1321350306.png",intern=TRUE)) character(0) > try(system("convert tmp/4hfrr1321350306.ps tmp/4hfrr1321350306.png",intern=TRUE)) character(0) > try(system("convert tmp/5sata1321350306.ps tmp/5sata1321350306.png",intern=TRUE)) character(0) > try(system("convert tmp/658oq1321350306.ps tmp/658oq1321350306.png",intern=TRUE)) character(0) > try(system("convert tmp/79eqp1321350306.ps tmp/79eqp1321350306.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.43 0.30 2.71