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(15136 + ,16733 + ,20016 + ,17708 + ,18019 + ,19227 + ,22893 + ,23739 + ,21133 + ,22591 + ,26786 + ,29740 + ,15028 + ,17977 + ,20008 + ,21354 + ,19498 + ,22125 + ,25817 + ,28779 + ,20960 + ,22254 + ,27392 + ,29945 + ,16933 + ,17892 + ,20533 + ,23569 + ,22417 + ,22084 + ,26580 + ,27454 + ,24081 + ,23451 + ,28991 + ,31386 + ,16896 + ,20045 + ,23471 + ,21747 + ,25621 + ,23859 + ,25500 + ,30998 + ,24475 + ,23145 + ,29701 + ,34365 + ,17556 + ,22077 + ,25702 + ,22214 + ,26886 + ,23191 + ,27831 + ,35406 + ,23195 + ,25110 + ,30009 + ,36242 + ,18450 + ,21845 + ,26488 + ,22394 + ,28057 + ,25451 + ,24872 + ,33424 + ,24052 + ,28449 + ,33533 + ,37351 + ,19969 + ,21701 + ,26249 + ,24493 + ,24603 + ,26485 + ,30723 + ,34569 + ,26689 + ,26157 + ,32064 + ,38870 + ,21337 + ,19419 + ,23166 + ,28286 + ,24570 + ,24001 + ,33151 + ,24878 + ,26804 + ,28967 + ,33311 + ,40226 + ,20504 + ,23060 + ,23562 + ,27562 + ,23940 + ,24584 + ,34303 + ,25517 + ,23494 + ,29095 + ,32903 + ,34379 + ,16991 + ,21109 + ,23740 + ,25552 + ,21752 + ,20294 + ,29009 + ,25500 + ,24166 + ,26960 + ,31222 + ,38641 + ,14672 + ,17543 + ,25453 + ,32683 + ,22449 + ,22316 + ,27595 + ,25451 + ,25421 + ,25288 + ,32568 + ,35110 + ,16052 + ,22146 + ,21198 + ,19543 + ,22084 + ,23816 + ,29961 + ,26773 + ,26635 + ,26972 + ,30207 + ,38687 + ,16974 + ,21697 + ,24179 + ,23757 + ,25013 + ,24019 + ,30345 + ,24488 + ,25156 + ,25650 + ,30923 + ,37240 + ,17466 + ,19463 + ,24352 + ,26805 + ,25236 + ,24735 + ,29356 + ,31234 + ,22724 + ,28496 + ,32857 + ,37198 + ,13652 + ,22784 + ,23565 + ,26323 + ,23779 + ,27549 + ,29660 + ,23356) > 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] 176 > (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 15 14 14 14 14 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [1,] 15136 15028 16933 16896 17556 18450 19969 21337 20504 16991 14672 16052 [2,] 16733 17977 17892 20045 22077 21845 21701 19419 23060 21109 17543 22146 [3,] 20016 20008 20533 23471 25702 26488 26249 23166 23562 23740 25453 21198 [4,] 17708 21354 23569 21747 22214 22394 24493 28286 27562 25552 32683 19543 [5,] 18019 19498 22417 25621 26886 28057 24603 24570 23940 21752 22449 22084 [6,] 19227 22125 22084 23859 23191 25451 26485 24001 24584 20294 22316 23816 [7,] 22893 25817 26580 25500 27831 24872 30723 33151 34303 29009 27595 29961 [8,] 23739 28779 27454 30998 35406 33424 34569 24878 25517 25500 25451 26773 [9,] 21133 20960 24081 24475 23195 24052 26689 26804 23494 24166 25421 26635 [10,] 22591 22254 23451 23145 25110 28449 26157 28967 29095 26960 25288 26972 [11,] 26786 27392 28991 29701 30009 33533 32064 33311 32903 31222 32568 30207 [12,] 29740 29945 31386 34365 36242 37351 38870 40226 34379 38641 35110 38687 [,13] [,14] [,15] [1,] 16974 17466 13652 [2,] 21697 19463 22784 [3,] 24179 24352 23565 [4,] 23757 26805 26323 [5,] 25013 25236 23779 [6,] 24019 24735 27549 [7,] 30345 29356 29660 [8,] 24488 31234 23356 [9,] 25156 22724 NA [10,] 25650 28496 NA [11,] 30923 32857 NA [12,] 37240 37198 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 1597 2949 959 3149 4521 3395 1732 -1918 2556 4118 [2,] 3283 2031 2641 3426 3625 4643 4548 3747 502 2631 [3,] -2308 1346 3036 -1724 -3488 -4094 -1756 5120 4000 1812 [4,] 311 -1856 -1152 3874 4672 5663 110 -3716 -3622 -3800 [5,] 1208 2627 -333 -1762 -3695 -2606 1882 -569 644 -1458 [6,] 3666 3692 4496 1641 4640 -579 4238 9150 9719 8715 [7,] 846 2962 874 5498 7575 8552 3846 -8273 -8786 -3509 [8,] -2606 -7819 -3373 -6523 -12211 -9372 -7880 1926 -2023 -1334 [9,] 1458 1294 -630 -1330 1915 4397 -532 2163 5601 2794 [10,] 4195 5138 5540 6556 4899 5084 5907 4344 3808 4262 [11,] 2954 2553 2395 4664 6233 3818 6806 6915 1476 7419 [12,] -14712 -13012 -14490 -16809 -17792 -17382 -17533 -19722 -17388 -23969 [,11] [,12] [,13] [,14] [,15] [1,] 2871 6094 4723 1997 9132 [2,] 7910 -948 2482 4889 781 [3,] 7230 -1655 -422 2453 2758 [4,] -10234 2541 1256 -1569 -2544 [5,] -133 1732 -994 -501 3770 [6,] 5279 6145 6326 4621 2111 [7,] -2144 -3188 -5857 1878 -6304 [8,] -30 -138 668 -8510 NA [9,] -133 337 494 5772 NA [10,] 7280 3235 5273 4361 NA [11,] 2542 8480 6317 4341 NA [12,] -19058 -21713 -19774 -23546 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/1ndg11272358769.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/2ndg11272358769.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/3ndg11272358769.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/4385a1272358769.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] [1,] 13652 16733 20008.0 17708.0 18019.0 19227.0 22893.0 23356.0 20960.0 [2,] 15594 18698 22182.0 21980.5 22250.5 22220.5 26198.5 25164.5 23195.0 [3,] 16974 21109 23565.0 23757.0 23940.0 23859.0 29009.0 26773.0 24123.5 [4,] 18003 21961 24902.5 26564.0 25124.5 24659.5 30153.0 31116.0 25421.0 [5,] 21337 23060 26488.0 32683.0 28057.0 27549.0 34303.0 35406.0 26804.0 [,10] [,11] [,12] [1,] 22254.0 26786.0 29740 [2,] 23451.0 29701.0 34365 [3,] 25903.5 31072.5 36720 [4,] 28449.0 32857.0 38641 [5,] 29095.0 33533.0 40226 $n [1] 15 15 15 15 15 15 15 15 14 14 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 15991.24 19777.85 22455.16 21887.14 22767.54 22864 27395.74 24345.06 [2,] 17956.76 22440.15 24674.84 25626.86 25112.46 24854 30622.26 29200.94 [,9] [,10] [,11] [,12] [1,] 23183.52 23792.98 29739.81 34914.36 [2,] 25063.48 28014.02 32405.19 38525.64 $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(13652, 15594, 16974, 18003, 21337, 16733, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5385a1272358769.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] [1,] 959.0 502.0 -4094 -10234.0 -3695 1641.0 -8786 -12211.0 -1330 3235.0 [2,] 1864.5 2256.5 -1740 -3083.0 -1226 3679.0 -4683 -7880.0 -133 4262.0 [3,] 2949.0 3283.0 1346 -1152.0 -333 4621.0 846 -2989.5 1376 4991.5 [4,] 4319.5 4147.5 2897 1898.5 1470 6235.5 3404 -138.0 2794 5540.0 [5,] 6094.0 4889.0 7230 5663.0 3770 9719.0 8552 1926.0 5772 7280.0 [,11] [,12] [1,] 1476.0 -23969.0 [2,] 2553.0 -19774.0 [3,] 4502.5 -17662.5 [4,] 6806.0 -16809.0 [5,] 8480.0 -13012.0 $n [1] 15 15 15 15 15 15 15 14 14 14 14 14 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1947.472 2511.559 -545.6838 -3184.224 -1432.8446 3578.065 -2453.126 [2,] 3950.528 4054.441 3237.6838 880.224 766.8446 5663.935 4145.126 [,8] [,9] [,10] [,11] [,12] [1,] -6258.7357 140.0076 4451.835 2706.574 -18914.54 [2,] 279.7357 2611.9924 5531.165 6298.426 -16410.46 $out [1] -1918 9132 7910 -948 -579 $group [1] 1 1 2 2 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(959, 1864.5, 2949, 4319.5, 6094, 502, 2256.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6385a1272358769.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] [1,] 15136.0 15028.0 16933.0 16896 17556.0 18450.0 19969.0 19419.0 20504.0 [2,] 17863.5 19753.0 21308.5 22446 22702.5 23223.0 24548.0 23583.5 23528.0 [3,] 20574.5 21739.5 23510.0 24167 25406.0 25969.5 26367.0 25841.0 25050.5 [4,] 23316.0 26604.5 27017.0 27661 28920.0 30936.5 31393.5 31059.0 30999.0 [5,] 29740.0 29945.0 31386.0 34365 36242.0 37351.0 38870.0 40226.0 34379.0 [,10] [,11] [,12] [,13] [,14] [,15] [1,] 16991.0 14672.0 16052.0 21697.0 17466.0 22784 [2,] 21430.5 22382.5 21641.0 23888.0 23538.0 23070 [3,] 24833.0 25436.0 25225.5 24750.5 26020.5 23672 [4,] 27984.5 30081.5 28466.5 27997.5 30295.0 26936 [5,] 31222.0 35110.0 38687.0 30923.0 37198.0 29660 $n [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 18087.58 18614.48 20906.31 21788.40 22570.16 22451.32 23244.72 22431.37 [2,] 23061.42 24864.52 26113.69 26545.60 28241.84 29487.68 29489.28 29250.63 [,9] [,10] [,11] [,12] [,13] [,14] [,15] [1,] 21642.93 21843.68 21924.43 22112.34 22876.13 22938.59 21512.40 [2,] 28458.07 27822.32 28947.57 28338.66 26624.87 29102.41 25831.60 $out [1] 38641 16974 37240 13652 $group [1] 10 13 13 15 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" "13" "14" NA Warning message: In bxp(list(stats = c(15136, 17863.5, 20574.5, 23316, 29740, 15028, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7vz4d1272358769.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,] 17174.40 21109.00 16798.50 [2,] 23513.93 23661.00 23491.12 [3,] 24239.61 24031.75 24292.25 [4,] 28305.40 27891.00 28158.00 [5,] 30890.50 31072.50 31281.38 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 22054.19 22102.42 22163.66 [2,] 26425.03 25961.08 26420.84 $out [1] 35670.0 16974.0 36720.0 36343.5 $group [1] 1 2 2 3 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(17174.4, 23513.9333333333, 24239.6071428571, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1ndg11272358769.ps tmp/1ndg11272358769.png",intern=TRUE)) character(0) > try(system("convert tmp/2ndg11272358769.ps tmp/2ndg11272358769.png",intern=TRUE)) character(0) > try(system("convert tmp/3ndg11272358769.ps tmp/3ndg11272358769.png",intern=TRUE)) character(0) > try(system("convert tmp/4385a1272358769.ps tmp/4385a1272358769.png",intern=TRUE)) character(0) > try(system("convert tmp/5385a1272358769.ps tmp/5385a1272358769.png",intern=TRUE)) character(0) > try(system("convert tmp/6385a1272358769.ps tmp/6385a1272358769.png",intern=TRUE)) character(0) > try(system("convert tmp/7vz4d1272358769.ps tmp/7vz4d1272358769.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.216 0.942 7.302