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(19246,17549,16428,16209,15235,16186,24971,30776,26416,23157,20155,19790,18849,17573,16597,16158,15507,16433,26325,31144,30535,27596,24064,23854,22407,21125,20226,19547,18933,20372,34331,37329,36761,32737,29321,28883,27436,25101,23776,23782,23027,25606,41328,44751,42855,37628,33544,33275,32009,30813,29143,28121,27007,29112,44067,48481,46581,41166,36824,35936,33633,31630,30434,28546,27660,29830,45599,49303,44417,40386,35544,35019,30400,29602,27701,27937,27283,29372,42821,45386,40170,34371,30077,29251,27202,25714,23784,22968,22243,24255,37282,38794,31828,27949,24605,25695,23338,21941,22034,20637,19418,22454,33261,34995,29132,26171,23828,25743,25204,25679,25281,25136,24794,28278,40062,42590,37885,34061,32412,34647,31750,31288,29331,28768,27780,30113,41240,43271,38108,34382,31551) > 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] 131 > (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 11 11 10 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] 19246 18849 22407 27436 32009 33633 30400 27202 23338 25204 31750 [2,] 17549 17573 21125 25101 30813 31630 29602 25714 21941 25679 31288 [3,] 16428 16597 20226 23776 29143 30434 27701 23784 22034 25281 29331 [4,] 16209 16158 19547 23782 28121 28546 27937 22968 20637 25136 28768 [5,] 15235 15507 18933 23027 27007 27660 27283 22243 19418 24794 27780 [6,] 16186 16433 20372 25606 29112 29830 29372 24255 22454 28278 30113 [7,] 24971 26325 34331 41328 44067 45599 42821 37282 33261 40062 41240 [8,] 30776 31144 37329 44751 48481 49303 45386 38794 34995 42590 43271 [9,] 26416 30535 36761 42855 46581 44417 40170 31828 29132 37885 38108 [10,] 23157 27596 32737 37628 41166 40386 34371 27949 26171 34061 34382 [11,] 20155 24064 29321 33544 36824 35544 30077 24605 23828 32412 31551 [12,] 19790 23854 28883 33275 35936 35019 29251 25695 25743 34647 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [1,] -1697 -1276 -1282 -2335 -1196 -2003 -798 -1488 -1397 475 -462 [2,] -1121 -976 -899 -1325 -1670 -1196 -1901 -1930 93 -398 -1957 [3,] -219 -439 -679 6 -1022 -1888 236 -816 -1397 -145 -563 [4,] -974 -651 -614 -755 -1114 -886 -654 -725 -1219 -342 -988 [5,] 951 926 1439 2579 2105 2170 2089 2012 3036 3484 2333 [6,] 8785 9892 13959 15722 14955 15769 13449 13027 10807 11784 11127 [7,] 5805 4819 2998 3423 4414 3704 2565 1512 1734 2528 2031 [8,] -4360 -609 -568 -1896 -1900 -4886 -5216 -6966 -5863 -4705 -5163 [9,] -3259 -2939 -4024 -5227 -5415 -4031 -5799 -3879 -2961 -3824 -3726 [10,] -3002 -3532 -3416 -4084 -4342 -4842 -4294 -3344 -2343 -1649 -2831 [11,] -365 -210 -438 -269 -888 -525 -826 1090 1915 2235 NA [12,] -941 -1447 -1447 -1266 -2303 -4619 -2049 -2357 -539 -2897 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/1s3z61293184054.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/html/rcomp/tmp/2s3z61293184054.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/html/rcomp/tmp/3kczr1293184054.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/html/rcomp/tmp/4kczr1293184054.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] [1,] 18849.0 17549.0 16428 16158 15235.0 16186 24971.0 30776.0 26416.0 23157.0 [2,] 22872.5 21533.0 21130 20092 19175.5 21413 33796.0 36162.0 31181.5 27772.5 [3,] 27202.0 25679.0 23784 23782 23027.0 25606 40062.0 42590.0 37885.0 34061.0 [4,] 31075.0 30207.5 28422 28029 27145.0 29242 42074.5 45068.5 41512.5 36005.0 [5,] 33633.0 31630.0 30434 28768 27780.0 30113 45599.0 49303.0 46581.0 41166.0 [,11] [,12] [1,] 20155.0 19790 [2,] 24334.5 25695 [3,] 30077.0 29067 [4,] 32978.0 34647 [5,] 36824.0 35936 $n [1] 11 11 11 11 11 11 11 11 11 11 11 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 23294.43 21546.57 20310.18 20000.91 19230.43 21876.36 36118.22 38347.05 [2,] 31109.57 29811.43 27257.82 27563.09 26823.57 29335.64 44005.78 46832.95 [,9] [,10] [,11] [,12] [1,] 32963.44 30139.14 25959.34 24594.22 [2,] 42806.56 37982.86 34194.66 33539.78 $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(18849, 22872.5, 27202, 31075, 33633, 17549, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5kczr1293184054.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,] -2335.0 -1957.0 -1888 -1219.0 926.0 8785 1512.0 -6966.0 -5799.0 -4842.0 [2,] -1592.5 -1785.5 -919 -981.0 1725.5 10967 2279.5 -5189.5 -4629.0 -4189.0 [3,] -1282.0 -1196.0 -563 -755.0 2105.0 13027 2998.0 -4705.0 -3879.0 -3416.0 [4,] -997.0 -937.5 -182 -652.5 2456.0 14457 4059.0 -1898.0 -3492.5 -2916.5 [5,] -462.0 93.0 236 -342.0 3484.0 15769 5805.0 -568.0 -2939.0 -1649.0 [,11] [,12] [1,] -888 -2897 [2,] -525 -2357 [3,] -317 -1748 [4,] 1090 -1266 [5,] 2235 -539 $n [1] 11 11 11 11 11 11 11 11 11 11 10 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1565.689 -1599.977 -914.0979 -911.4934 1756.999 11364.41 2150.268 [2,] -998.311 -792.023 -211.9021 -598.5066 2453.001 14689.59 3845.732 [,8] [,9] [,10] [,11] [,12] [1,] -6273.031 -4420.415 -4022.204 -1123.9184 -2293.107 [2,] -3136.969 -3337.585 -2809.796 489.9184 -1202.893 $out [1] 475 -4619 $group [1] 1 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-2335, -1592.5, -1282, -997, -462, -1957, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6vlyu1293184054.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,] 15235.0 15507.0 18933 23027.0 27007.0 27660.0 27283.0 22243.0 19418.0 [2,] 16318.5 16515.0 20299 24441.5 29127.5 30132.0 28594.0 24019.5 21987.5 [3,] 19518.0 21351.5 25645 30355.5 33972.5 34326.0 29839.5 25704.5 23583.0 [4,] 24064.0 26960.5 33534 39478.0 42616.5 42401.5 37270.5 29888.5 27651.5 [5,] 30776.0 31144.0 37329 44751.0 48481.0 49303.0 45386.0 37282.0 34995.0 [,10] [,11] [1,] 24794.0 27780 [2,] 25242.5 29722 [3,] 30345.0 31551 [4,] 36266.0 36245 [5,] 42590.0 43271 $n [1] 12 12 12 12 12 12 12 12 12 12 11 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 15985.23 16587.24 19608.43 23497.25 27820.08 28729.8 25882.09 23027.61 [2,] 23050.77 26115.76 31681.57 37213.75 40124.92 39922.2 33796.91 28381.39 [,9] [,10] [,11] [1,] 20999.61 25317.11 28443.52 [2,] 26166.39 35372.89 34658.48 $out [1] 38794 $group [1] 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" NA Warning message: In bxp(list(stats = c(15235, 16318.5, 19518, 24064, 30776, 15507, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7vlyu1293184054.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,] 22626.09 23027.0 23160.25 [2,] 24397.55 24695.0 25051.75 [3,] 27853.47 28134.5 27815.00 [4,] 34740.55 35973.0 34117.88 [5,] 40620.00 42590.0 40615.25 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 23135.96 22990.53 23679.88 [2,] 32570.98 33278.47 31950.12 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(22626.0909090909, 24397.5454545455, 27853.4681818182, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1s3z61293184054.ps tmp/1s3z61293184054.png",intern=TRUE)) character(0) > try(system("convert tmp/2s3z61293184054.ps tmp/2s3z61293184054.png",intern=TRUE)) character(0) > try(system("convert tmp/3kczr1293184054.ps tmp/3kczr1293184054.png",intern=TRUE)) character(0) > try(system("convert tmp/4kczr1293184054.ps tmp/4kczr1293184054.png",intern=TRUE)) character(0) > try(system("convert tmp/5kczr1293184054.ps tmp/5kczr1293184054.png",intern=TRUE)) character(0) > try(system("convert tmp/6vlyu1293184054.ps tmp/6vlyu1293184054.png",intern=TRUE)) character(0) > try(system("convert tmp/7vlyu1293184054.ps tmp/7vlyu1293184054.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.217 0.973 2.964