R version 2.8.0 (2008-10-20) Copyright (C) 2008 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. Natural language support but running in an English locale 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(43880,43110,44496,44164,40399,36763,37903,35532,35533,32110,33374,35462,33508,36080,34560,38737,38144,37594,36424,36843,37246,38661,40454,44928,48441,48140,45998,47369,49554,47510,44873,45344,42413,36912,43452,42142,44382,43636,44167,44423,42868,43908,42013,38846,35087,33026,34646,37135,37985,43121,43722,43630,42234,39351,39327,35704,30466,28155,29257,29998,32529,34787,33855,34556,31348,30805,28353,24514,21106,21346,23335,24379,26290,30084,29429,30632,27349,27264,27474,24482,21453,18788,19282,19713,21917,23812,23785,24696,24562,23580,24939,23899,21454,19761,19815,20780,23462,25005,24725,26198,27543,26471,26558,25317,22896,22248,23406,25073,27691,30599,31948,32946,34012,32936,32974,30951,29812,29010,31068,32447,34844,35676,35387,36488,35652,33488,32914,29781,27951) > 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] 129 > (np <- floor(n / par1)) [1] 10 > 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] 11 11 11 11 11 11 11 11 11 10 10 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 43880 33508 48441 44382 37985 32529 26290 21917 23462 27691 34844 [2,] 43110 36080 48140 43636 43121 34787 30084 23812 25005 30599 35676 [3,] 44496 34560 45998 44167 43722 33855 29429 23785 24725 31948 35387 [4,] 44164 38737 47369 44423 43630 34556 30632 24696 26198 32946 36488 [5,] 40399 38144 49554 42868 42234 31348 27349 24562 27543 34012 35652 [6,] 36763 37594 47510 43908 39351 30805 27264 23580 26471 32936 33488 [7,] 37903 36424 44873 42013 39327 28353 27474 24939 26558 32974 32914 [8,] 35532 36843 45344 38846 35704 24514 24482 23899 25317 30951 29781 [9,] 35533 37246 42413 35087 30466 21106 21453 21454 22896 29812 27951 [10,] 32110 38661 36912 33026 28155 21346 18788 19761 22248 29010 NA [11,] 33374 40454 43452 34646 29257 23335 19282 19815 23406 31068 NA [12,] 35462 44928 42142 37135 29998 24379 19713 20780 25073 32447 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -770 2572 -301 -746 5136 2258 3794 1895 1543 2908 832 [2,] 1386 -1520 -2142 531 601 -932 -655 -27 -280 1349 -289 [3,] -332 4177 1371 256 -92 701 1203 911 1473 998 1101 [4,] -3765 -593 2185 -1555 -1396 -3208 -3283 -134 1345 1066 -836 [5,] -3636 -550 -2044 1040 -2883 -543 -85 -982 -1072 -1076 -2164 [6,] 1140 -1170 -2637 -1895 -24 -2452 210 1359 87 38 -574 [7,] -2371 419 471 -3167 -3623 -3839 -2992 -1040 -1241 -2023 -3133 [8,] 1 403 -2931 -3759 -5238 -3408 -3029 -2445 -2421 -1139 -1830 [9,] -3423 1415 -5501 -2061 -2311 240 -2665 -1693 -648 -802 NA [10,] 1264 1793 6540 1620 1102 1989 494 54 1158 2058 NA [11,] 2088 4474 -1310 2489 741 1044 431 965 1667 1379 NA [12,] -1954 3513 2240 850 2531 1911 2204 2682 2618 2397 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/freestat/rcomp/tmp/1c4su1289562258.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/freestat/rcomp/tmp/2c4su1289562258.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/freestat/rcomp/tmp/3c4su1289562258.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/freestat/rcomp/tmp/45d9x1289562258.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,] 21917.0 23812.0 23785.0 24696 24562.0 23580.0 24939.0 23899.0 21106 [2,] 26990.5 30341.5 30688.5 31789 29445.5 29034.5 27913.5 24915.5 22175 [3,] 33508.0 35676.0 34560.0 36488 35652.0 33488.0 32974.0 30951.0 29812 [4,] 40932.5 43115.5 43944.5 43897 41316.5 38472.5 38615.0 36273.5 35310 [5,] 48441.0 48140.0 45998.0 47369 49554.0 47510.0 44873.0 45344.0 42413 [,10] [,11] [,12] [1,] 18788.0 19282.0 19713.0 [2,] 21346.0 23335.0 24379.0 [3,] 28582.5 30162.5 31222.5 [4,] 33026.0 34646.0 37135.0 [5,] 38661.0 43452.0 44928.0 $n [1] 11 11 11 11 11 11 11 11 11 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 26866.2 29590.62 28245 30719.90 29996.8 28991.85 27875.93 25540.19 [2,] 40149.8 41761.38 40875 42256.10 41307.2 37984.15 38072.07 36361.81 [,9] [,10] [,11] [,12] [1,] 23554.64 22746.71 24511.07 24849.09 [2,] 36069.36 34418.29 35813.93 37595.91 $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(21917, 26990.5, 33508, 40932.5, 48441, 23812, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/55d9x1289562258.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,] -770.0 -2142.0 -332.0 -3765.0 -3636.0 -2637.0 -3839.0 -5238.0 -5501 54 [2,] 265.5 -793.5 478.5 -2381.5 -2104.0 -1532.5 -3150.0 -3218.5 -2665 1102 [3,] 1895.0 -280.0 998.0 -836.0 -1072.0 -24.0 -2371.0 -2445.0 -1877 1442 [4,] 2740.0 566.0 1287.0 466.0 -546.5 148.5 -1140.5 -1484.5 -648 1989 [5,] 5136.0 1386.0 1473.0 2185.0 1040.0 1359.0 471.0 403.0 1415 2058 [,11] [,12] [1,] 431.0 1911.0 [2,] 741.0 1911.0 [3,] 1211.5 2318.5 [4,] 2088.0 2618.0 [5,] 2489.0 3513.0 $n [1] 11 11 11 11 11 11 11 11 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 716.1781 -927.6494 612.8404 -2192.5146 -1813.9742 -824.8081 -3328.302 [2,] 3073.8219 367.6494 1383.1596 520.5146 -330.0258 776.8081 -1413.698 [,8] [,9] [,10] [,11] [,12] [1,] -3271.057 -2884.7736 998.8194 538.4851 1965.255 [2,] -1618.943 -869.2264 1885.1806 1884.5149 2671.745 $out [1] 4177 6540 4474 -1310 -1954 850 $group [1] 3 10 11 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-770, 265.5, 1895, 2740, 5136, -2142, -793.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/6xm901289562258.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,] 32110 33508 36912.0 33026.0 28155.0 21106 18788.0 19761.0 22248.0 27691.0 [2,] 35497 36252 42932.5 36111.0 30232.0 23857 20583.0 21117.0 23434.0 30205.5 [3,] 37333 37420 45671.0 42440.5 38656.0 29579 26777.0 23682.5 25039.0 31508.0 [4,] 43495 38699 47825.0 44037.5 42677.5 33192 28451.5 24230.5 26334.5 32941.0 [5,] 44496 40454 49554.0 44423.0 43722.0 34787 30632.0 24939.0 27543.0 34012.0 [,11] [1,] 29781 [2,] 32914 [3,] 34844 [4,] 35652 [5,] 36488 $n [1] 12 12 12 12 12 12 12 12 12 12 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 33685.06 36303.91 43439.5 38825.17 32979.52 25321.24 23188.12 22262.41 [2,] 40980.94 38536.09 47902.5 46055.83 44332.48 33836.76 30365.88 25102.59 [,9] [,10] [,11] [1,] 23716.06 30260.32 33401.99 [2,] 26361.94 32755.68 36286.01 $out [1] 44928 27951 $group [1] 2 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(32110, 35497, 37333, 43495, 44496, 33508, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/7xm901289562258.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,] 28001.70 28582.50 27184.25 [2,] 30507.30 30556.75 29717.44 [3,] 34030.95 33231.00 33508.88 [4,] 35715.32 35106.00 36054.75 [5,] 36712.64 36488.00 37843.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 31655.54 31156.06 30618.38 [2,] 36406.37 35305.94 36399.37 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(28001.7, 30507.3, 34030.9545454545, 35715.3181818182, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1c4su1289562258.ps tmp/1c4su1289562258.png",intern=TRUE)) character(0) > try(system("convert tmp/2c4su1289562258.ps tmp/2c4su1289562258.png",intern=TRUE)) character(0) > try(system("convert tmp/3c4su1289562258.ps tmp/3c4su1289562258.png",intern=TRUE)) character(0) > try(system("convert tmp/45d9x1289562258.ps tmp/45d9x1289562258.png",intern=TRUE)) character(0) > try(system("convert tmp/55d9x1289562258.ps tmp/55d9x1289562258.png",intern=TRUE)) character(0) > try(system("convert tmp/6xm901289562258.ps tmp/6xm901289562258.png",intern=TRUE)) character(0) > try(system("convert tmp/7xm901289562258.ps tmp/7xm901289562258.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.944 1.469 2.150