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(8357,7454,8076,7248,7339,7292,7359,7537,7441,8057,8037,8257,8692,8119,8236,7432,7669,7453,7566,7731,7657,8130,8401,8737,9009,7919,8228,7903,7912,7857,7965,8091,8024,8772,8656,8953,9014,8103,8876,8231,8173,8087,8296,8007,8382,9168,9137,9321,9234,8451,9101,8279,8284,8225,8597,8305,8620,9102,9258,9652,9522,8874,9415,8525,8862,8421,8626,8750,8852,9412,9570,9513,9986,8907,9663,8799,8931,8732,8936,9127,9070,9773,9670,9929,10095,9025,9659,8954,9022,8855,9034,9196,9038,9650,9715,10052,10436,9314,9717,8997,9062,8885,9058,9095,9149,9857,9848,10269,10341,9690,10125,9349,9224,9224,9454,9347,9430,9933,10148,10677,10735,9760,10567,9333,9409,9502,9348,9319,9594,10160,10182,10810,11105,9874,10958,9311,9610,9398,9784,9425,9557,10166,10337,10770,11265,10183,10941,9628,9709,9637,9579,9741,9754,10508,10749,11079,11608,10668,10933,9703,9799,9656,9648,9712,9766,10540,10564,10911,11218,10230,10410,9227,9378,9105,9128) > 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] 175 > (np <- floor(n / par1)) [1] 14 > 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] 15 15 15 15 15 15 15 14 14 14 14 14 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] 8357 8692 9009 9014 9234 9522 9986 10095 10436 10341 10735 11105 11265 [2,] 7454 8119 7919 8103 8451 8874 8907 9025 9314 9690 9760 9874 10183 [3,] 8076 8236 8228 8876 9101 9415 9663 9659 9717 10125 10567 10958 10941 [4,] 7248 7432 7903 8231 8279 8525 8799 8954 8997 9349 9333 9311 9628 [5,] 7339 7669 7912 8173 8284 8862 8931 9022 9062 9224 9409 9610 9709 [6,] 7292 7453 7857 8087 8225 8421 8732 8855 8885 9224 9502 9398 9637 [7,] 7359 7566 7965 8296 8597 8626 8936 9034 9058 9454 9348 9784 9579 [8,] 7537 7731 8091 8007 8305 8750 9127 9196 9095 9347 9319 9425 9741 [9,] 7441 7657 8024 8382 8620 8852 9070 9038 9149 9430 9594 9557 9754 [10,] 8057 8130 8772 9168 9102 9412 9773 9650 9857 9933 10160 10166 10508 [11,] 8037 8401 8656 9137 9258 9570 9670 9715 9848 10148 10182 10337 10749 [12,] 8257 8737 8953 9321 9652 9513 9929 10052 10269 10677 10810 10770 11079 [,14] [,15] [1,] 11608 11218 [2,] 10668 10230 [3,] 10933 10410 [4,] 9703 9227 [5,] 9799 9378 [6,] 9656 9105 [7,] 9648 9128 [8,] 9712 NA [9,] 9766 NA [10,] 10540 NA [11,] 10564 NA [12,] 10911 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [1,] -903 -573 -1090 -911 -783 -648 -1079 -1070 -1122 -651 -975 -1231 -1082 [2,] 622 117 309 773 650 541 756 634 403 435 807 1084 758 [3,] -828 -804 -325 -645 -822 -890 -864 -705 -720 -776 -1234 -1647 -1313 [4,] 91 237 9 -58 5 337 132 68 65 -125 76 299 81 [5,] -47 -216 -55 -86 -59 -441 -199 -167 -177 0 93 -212 -72 [6,] 67 113 108 209 372 205 204 179 173 230 -154 386 -58 [7,] 178 165 126 -289 -292 124 191 162 37 -107 -29 -359 162 [8,] -96 -74 -67 375 315 102 -57 -158 54 83 275 132 13 [9,] 616 473 748 786 482 560 703 612 708 503 566 609 754 [10,] -20 271 -116 -31 156 158 -103 65 -9 215 22 171 241 [11,] 220 336 297 184 394 -57 259 337 421 529 628 433 330 [12,] 435 272 61 -87 -130 473 166 384 72 58 295 495 529 [,14] [,15] [1,] -940 -988 [2,] 265 180 [3,] -1230 -1183 [4,] 96 151 [5,] -143 -273 [6,] -8 23 [7,] 64 NA [8,] 54 NA [9,] 774 NA [10,] 24 NA [11,] 347 NA [12,] 307 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/1bwdu1259589807.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/2kbz21259589807.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/3o50a1259589807.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/46hee1259589807.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] [,9] [,10] [,11] [1,] 8357 7454 8076.0 7248 7339.0 7292 7359.0 7537 7441 8057.0 8037.0 [2,] 9124 8285 8988.5 8255 8228.5 8156 8446.5 8091 8382 9102.0 9137.0 [3,] 10095 9025 9663.0 8954 9022.0 8855 9034.0 9111 9054 9711.5 9692.5 [4,] 10920 9817 10488.5 9322 9393.5 9311 9401.0 9347 9557 10160.0 10182.0 [5,] 11608 10668 10958.0 9703 9799.0 9656 9784.0 9741 9766 10540.0 10749.0 [,12] [1,] 8257.0 [2,] 9321.0 [3,] 9990.5 [4,] 10770.0 [5,] 11079.0 $n [1] 15 15 15 15 15 15 15 14 14 14 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9362.314 8400.014 9051.069 8518.713 8546.733 8383.813 8644.608 8580.625 [2,] 10827.686 9649.986 10274.931 9389.287 9497.267 9326.187 9423.392 9641.375 [,9] [,10] [,11] [,12] [1,] 8557.83 9264.735 9251.225 9378.627 [2,] 9550.17 10158.265 10133.775 10602.373 $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(8357, 9124, 10095, 10920, 11608, 7454, 8285, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/57cko1259589807.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] [,7] [,8] [,9] [,10] [,11] [,12] [1,] -1231.0 117 -1647.0 -58.0 -273.0 -154 -359 -158 473 -116.0 184.0 -130.0 [2,] -1080.5 356 -1206.5 37.0 -205.5 45 -107 -67 560 -20.0 259.0 61.0 [3,] -975.0 622 -828.0 81.0 -143.0 173 94 54 614 44.5 336.5 283.5 [4,] -843.0 757 -748.0 141.5 -57.0 207 162 132 748 171.0 421.0 435.0 [5,] -573.0 1084 -325.0 237.0 93.0 386 191 375 786 271.0 628.0 529.0 $n [1] 15 15 15 15 15 15 14 14 14 14 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1071.8891 458.4103 -1015.047 38.36878 -203.5812 106.9114 -19.59137 [2,] -878.1109 785.5897 -640.953 123.63122 -82.4188 239.0886 207.59137 [,8] [,9] [,10] [,11] [,12] [1,] -30.03228 534.6127 -36.1541 268.0918 125.5700 [2,] 138.03228 693.3873 125.1541 404.9082 441.4300 $out [1] 337 -125 299 -441 -57 $group [1] 4 4 4 5 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-1231, -1080.5, -975, -843, -573, 117, 356, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6i3y31259589807.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] [,9] [,10] [1,] 7248.0 7432.0 7857.0 8007.0 8225.0 8421 8732.0 8855.0 8885.0 9224.0 [2,] 7349.0 7611.5 7915.5 8138.0 8294.5 8688 8919.0 9023.5 9060.0 9348.0 [3,] 7495.5 7925.0 8057.5 8339.0 8608.5 8868 9098.5 9117.0 9231.5 9572.0 [4,] 8066.5 8318.5 8714.0 9075.5 9168.0 9464 9721.5 9687.0 9852.5 10136.5 [5,] 8357.0 8737.0 9009.0 9321.0 9652.0 9570 9986.0 10095.0 10436.0 10677.0 [,11] [,12] [,13] [,14] [,15] [1,] 9319.0 9311.0 9579.0 9648.0 9105.0 [2,] 9378.5 9491.0 9673.0 9707.5 9177.5 [3,] 9677.0 9829.0 9968.5 10169.5 9378.0 [4,] 10374.5 10553.5 10845.0 10789.5 10320.0 [5,] 10810.0 11105.0 11265.0 11608.0 11218.0 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 7168.243 7602.533 7693.299 7911.4 8210.09 8514.061 8732.474 8814.373 [2,] 7822.757 8247.467 8421.701 8766.6 9006.91 9221.939 9464.526 9419.627 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [1,] 8870.035 9212.36 9222.718 9344.387 9433.943 9675.993 8695.717 [2,] 9592.965 9931.64 10131.282 10313.613 10503.057 10663.007 10060.283 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" NA Warning message: In bxp(list(stats = c(7248, 7349, 7495.5, 8066.5, 8357, 7432, 7611.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7eief1259589807.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,] 8688.600 8855.0 8733.500 [2,] 8819.136 9023.5 8799.750 [3,] 8992.867 9082.5 9017.188 [4,] 9625.595 9702.0 9704.438 [5,] 10041.133 10095.0 10057.875 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 8625.035 8773.032 8604.553 [2,] 9360.698 9391.968 9429.822 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(8688.6, 8819.13571428572, 8992.86666666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1bwdu1259589807.ps tmp/1bwdu1259589807.png") > system("convert tmp/2kbz21259589807.ps tmp/2kbz21259589807.png") > system("convert tmp/3o50a1259589807.ps tmp/3o50a1259589807.png") > system("convert tmp/46hee1259589807.ps tmp/46hee1259589807.png") > system("convert tmp/57cko1259589807.ps tmp/57cko1259589807.png") > system("convert tmp/6i3y31259589807.ps tmp/6i3y31259589807.png") > system("convert tmp/7eief1259589807.ps tmp/7eief1259589807.png") > > > proc.time() user system elapsed 1.236 0.967 1.998