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(211868,229527,229139,198563,195722,202196,205816,212588,214320,220375,204442,206903,214126,226899,223532,195309,186005,188906,191563,189226,186413,178037,166827,169362,174330,187069,186530,158114,151001,159612,161914,164182,169701,171297,166444,173476,182516,202388,202300,168053,167302,172608,178106,185686,194581,194596,197922,208795,230580,240636,240048,211457,211142,214771,212610,219313,219277,231805,229245,241114,248624,265845,256446,219452,217142,221678,227184,230354,235243,237217,233575,244460,243324,260307,241476,203666,200237,204045,209465,213586,216234,213188,208679,217859,227247,243477,232571,191531,186029,189733,190420,194163,198770,195198,193111,195411,202108,215706,206348,166972,166070,169292,175041,177876,181140,179566,175335,184128,189917,194690,179612,150605,150569,153745,155511,159044,163095,159585,158644,166618,176512,200765,182698,153730,156145,161570,165688,173666,180144) > 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] [1,] 211868 214126 174330 182516 230580 248624 243324 227247 202108 189917 [2,] 229527 226899 187069 202388 240636 265845 260307 243477 215706 194690 [3,] 229139 223532 186530 202300 240048 256446 241476 232571 206348 179612 [4,] 198563 195309 158114 168053 211457 219452 203666 191531 166972 150605 [5,] 195722 186005 151001 167302 211142 217142 200237 186029 166070 150569 [6,] 202196 188906 159612 172608 214771 221678 204045 189733 169292 153745 [7,] 205816 191563 161914 178106 212610 227184 209465 190420 175041 155511 [8,] 212588 189226 164182 185686 219313 230354 213586 194163 177876 159044 [9,] 214320 186413 169701 194581 219277 235243 216234 198770 181140 163095 [10,] 220375 178037 171297 194596 231805 237217 213188 195198 179566 159585 [11,] 204442 166827 166444 197922 229245 233575 208679 193111 175335 158644 [12,] 206903 169362 173476 208795 241114 244460 217859 195411 184128 166618 [,11] [1,] 176512 [2,] 200765 [3,] 182698 [4,] 153730 [5,] 156145 [6,] 161570 [7,] 165688 [8,] 173666 [9,] 180144 [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [1,] 17659 12773 12739 19872 10056 17221 16983 16230 13598 4773 [2,] -388 -3367 -539 -88 -588 -9399 -18831 -10906 -9358 -15078 [3,] -30576 -28223 -28416 -34247 -28591 -36994 -37810 -41040 -39376 -29007 [4,] -2841 -9304 -7113 -751 -315 -2310 -3429 -5502 -902 -36 [5,] 6474 2901 8611 5306 3629 4536 3808 3704 3222 3176 [6,] 3620 2657 2302 5498 -2161 5506 5420 687 5749 1766 [7,] 6772 -2337 2268 7580 6703 3170 4121 3743 2835 3533 [8,] 1732 -2813 5519 8895 -36 4889 2648 4607 3264 4051 [9,] 6055 -8376 1596 15 12528 1974 -3046 -3572 -1574 -3510 [10,] -15933 -11210 -4853 3326 -2560 -3642 -4509 -2087 -4231 -941 [11,] 2461 2535 7032 10873 11869 10885 9180 2300 8793 7974 [12,] 7223 4968 9040 21785 7510 -1136 9388 6697 5789 9894 [,11] [1,] 24253 [2,] -18067 [3,] -28968 [4,] 2415 [5,] 5425 [6,] 4118 [7,] 7978 [8,] 6478 [9,] NA [10,] NA [11,] NA [12,] 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/1vgg21292593475.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/2vgg21292593475.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/357xn1292593475.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/457xn1292593475.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] [1,] 174330.0 187069.0 179612.0 150605.0 150569.0 153745.0 155511.0 159044 [2,] 186216.5 201576.5 194415.0 162543.0 161107.5 165431.0 170364.5 175771 [3,] 211868.0 226899.0 223532.0 191531.0 186005.0 188906.0 190420.0 189226 [4,] 228913.5 242056.5 236309.5 201114.5 197979.5 203120.5 207640.5 213087 [5,] 248624.0 265845.0 256446.0 219452.0 217142.0 221678.0 227184.0 230354 [,9] [,10] [,11] [,12] [1,] 163095 159585 158644.0 166618 [2,] 180642 178037 166827.0 173476 [3,] 194581 194897 195516.5 201157 [4,] 215277 220375 208679.0 217859 [5,] 235243 237217 233575.0 244460 $n [1] 11 11 11 11 11 11 11 11 11 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 191527.7 207614.8 203574.0 173156 168439.6 170951.2 172662.2 171449.1 [2,] 232208.3 246183.2 243490.0 209906 203570.4 206860.8 208177.8 207002.9 [,9] [,10] [,11] [,12] [1,] 178081.3 173743.2 174605.6 178981.5 [2,] 211080.7 216050.8 216427.4 223332.5 $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(174330, 186216.5, 211868, 228913.5, 248624, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/5ghe81292593475.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] [,10] [1,] 10056 -18831.0 -41040.0 -9304.0 2901.0 -2161 -2337.0 -36 -8376.0 -4853.0 [2,] 12756 -12992.0 -37402.0 -4465.5 3425.5 2034 3002.5 2190 -3510.0 -4853.0 [3,] 16230 -9358.0 -30576.0 -2310.0 3808.0 3620 3743.0 4051 -779.5 -3936.5 [4,] 17440 -563.5 -28779.5 -533.0 5365.5 5459 6737.5 5204 1974.0 -2087.0 [5,] 24253 -88.0 -28223.0 2415.0 6474.0 5749 7978.0 8895 6055.0 -941.0 [,11] [,12] [1,] 2300.0 4968.0 [2,] 2535.0 5789.0 [3,] 8383.5 7366.5 [4,] 10873.0 9388.0 [5,] 11869.0 9894.0 $n [1] 11 11 11 11 11 11 11 11 10 10 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 13998.6 -15278.787 -34683.65 -4183.3955 2883.807 1988.371 1963.691 [2,] 18461.4 -3437.213 -26468.35 -436.6045 4732.193 5251.629 5522.309 [,8] [,9] [,10] [,11] [,12] [1,] 2615.167 -3519.525 -5318.504 4217.503 5568.296 [2,] 5486.833 1960.525 -2554.496 12549.497 9164.704 $out [1] 4773 8611 -2813 12528 -15933 -11210 3326 21785 -1136 $group [1] 1 5 8 9 10 10 10 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(10056, 12756, 16230, 17440, 24253, -18831, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/6ghe81292593475.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] [1,] 195722.0 166827.0 151001.0 167302.0 211142.0 217142 200237.0 186029.0 [2,] 203319.0 182021.0 160763.0 175357.0 213690.5 224431 206362.0 190975.5 [3,] 209385.5 189066.0 168072.5 190133.5 224279.0 234409 213387.0 194680.5 [4,] 217347.5 204717.5 173903.0 200111.0 235926.5 246542 229667.5 213008.5 [5,] 229527.0 226899.0 187069.0 208795.0 241114.0 265845 260307.0 243477.0 [,9] [,10] [,11] [1,] 166070.0 150569.0 153730 [2,] 172166.5 154628.0 161570 [3,] 178721.0 159314.5 173666 [4,] 193118.0 173115.0 180144 [5,] 215706.0 194690.0 200765 $n [1] 12 12 12 12 12 12 12 12 12 12 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 202987 178714.0 162079.3 178843.0 214137 224324.0 202757.2 184631.1 [2,] 215784 199418.0 174065.7 201424.0 234421 244494.0 224016.8 204729.9 [,9] [,10] [,11] [1,] 169164.9 150882.5 163883.7 [2,] 188277.1 167746.5 183448.3 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(195722, 203319, 209385.5, 217347.5, 229527, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/rcomp/tmp/7ghe81292593475.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,] 180669.5 186005.0 179543.5 [2,] 186885.2 189823.0 186281.3 [3,] 194843.8 194739.0 195147.5 [4,] 205004.1 206512.5 203031.9 [5,] 224300.8 226899.0 221816.5 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 186579.7 187126.8 187507.5 [2,] 203108.0 202351.2 202787.5 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(180669.454545455, 186885.181818182, 194843.836363636, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1vgg21292593475.ps tmp/1vgg21292593475.png",intern=TRUE)) character(0) > try(system("convert tmp/2vgg21292593475.ps tmp/2vgg21292593475.png",intern=TRUE)) character(0) > try(system("convert tmp/357xn1292593475.ps tmp/357xn1292593475.png",intern=TRUE)) character(0) > try(system("convert tmp/457xn1292593475.ps tmp/457xn1292593475.png",intern=TRUE)) character(0) > try(system("convert tmp/5ghe81292593475.ps tmp/5ghe81292593475.png",intern=TRUE)) character(0) > try(system("convert tmp/6ghe81292593475.ps tmp/6ghe81292593475.png",intern=TRUE)) character(0) > try(system("convert tmp/7ghe81292593475.ps tmp/7ghe81292593475.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.63 1.09 2.71