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(100 + ,99.94940551 + ,102.0743739 + ,102.0237794 + ,102.6309132 + ,102.8838857 + ,103.0103719 + ,104.3258285 + ,105.211232 + ,104.7305844 + ,104.174045 + ,103.6933974 + ,104.3258285 + ,105.211232 + ,105.3124209 + ,105.5400961 + ,106.7037693 + ,106.6025803 + ,105.008854 + ,104.3258285 + ,104.174045 + ,103.4151278 + ,102.3273463 + ,101.8214015 + ,103.5669112 + ,103.8957754 + ,104.5788009 + ,105.0341513 + ,105.6665823 + ,105.6918796 + ,105.7171768 + ,105.843663 + ,105.7930686 + ,105.3883127 + ,105.9701493 + ,106.5013913 + ,107.1338224 + ,109.3599798 + ,109.3599798 + ,108.4239818 + ,107.9433342 + ,108.0951176 + ,108.3986845 + ,110.5489502 + ,111.8138123 + ,112.5474323 + ,111.6620288 + ,111.3837592 + ,113.1039717 + ,115.1783456 + ,121.0726031 + ,123.0710853 + ,123.3999494 + ,122.9193018 + ,122.3880597 + ,123.5517329 + ,124.9683784 + ,124.8671895 + ,123.2734632 + ,121.9580066 + ,122.4892487 + ,125.6767012 + ,126.7644827 + ,126.4356185 + ,125.3478371 + ,126.0055654 + ,127.4475082 + ,130.5843663 + ,133.0887933 + ,133.3417657 + ,132.8358209 + ,133.7971161 + ,136.6304073 + ,138.8818619 + ,140.6526689 + ,143.9160132 + ,149.0260562 + ,149.3296231 + ,152.6182646 + ,161.5987857 + ,162.6612699 + ,166.1269922 + ,165.9752087 + ,165.4692639 + ,166.9365039 + ,167.4930433 + ,169.2891475 + ,170.2757399 + ,171.5911966 + ,171.3635214 + ,171.641791 + ,170.0986592 + ,168.302555 + ,168.9602833 + ,170.0227675 + ,167.5942322 + ,172.6030863 + ,175.082216 + ,177.7637238 + ,182.4943081 + ,185.4793827 + ,189.1474829 + ,186.5924614 + ,184.4927903 + ,178.62383 + ,170.5793069 + ,166.7088287 + ,163.3442955) > 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] 108 > (np <- floor(n / par1)) [1] 9 > 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] 9 9 9 9 9 9 9 9 9 9 9 9 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 100.0000 104.3258 103.5669 107.1338 113.1040 122.4892 136.6304 166.9365 [2,] 99.9494 105.2112 103.8958 109.3600 115.1783 125.6767 138.8819 167.4930 [3,] 102.0744 105.3124 104.5788 109.3600 121.0726 126.7645 140.6527 169.2891 [4,] 102.0238 105.5401 105.0342 108.4240 123.0711 126.4356 143.9160 170.2757 [5,] 102.6309 106.7038 105.6666 107.9433 123.3999 125.3478 149.0261 171.5912 [6,] 102.8839 106.6026 105.6919 108.0951 122.9193 126.0056 149.3296 171.3635 [7,] 103.0104 105.0089 105.7172 108.3987 122.3881 127.4475 152.6183 171.6418 [8,] 104.3258 104.3258 105.8437 110.5490 123.5517 130.5844 161.5988 170.0987 [9,] 105.2112 104.1740 105.7931 111.8138 124.9684 133.0888 162.6613 168.3026 [10,] 104.7306 103.4151 105.3883 112.5474 124.8672 133.3418 166.1270 168.9603 [11,] 104.1740 102.3273 105.9701 111.6620 123.2735 132.8358 165.9752 170.0228 [12,] 103.6934 101.8214 106.5014 111.3838 121.9580 133.7971 165.4693 167.5942 [,9] [,10] [1,] 172.6031 NA [2,] 175.0822 NA [3,] 177.7637 NA [4,] 182.4943 NA [5,] 185.4794 NA [6,] 189.1475 NA [7,] 186.5925 NA [8,] 184.4928 NA [9,] 178.6238 NA [10,] 170.5793 NA [11,] 166.7088 NA [12,] 163.3443 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -0.05059449 0.8854035 0.3288642 2.2261574 2.0743739 3.1874525 [2,] 2.12496839 0.1011889 0.6830255 0.0000000 5.8942575 1.0877815 [3,] -0.05059450 0.2276752 0.4553504 -0.9359980 1.9984822 -0.3288642 [4,] 0.60713380 1.1636732 0.6324310 -0.4806476 0.3288641 -1.0877814 [5,] 0.25297250 -0.1011890 0.0252973 0.1517834 -0.4806476 0.6577283 [6,] 0.12648620 -1.5937263 0.0252972 0.3035669 -0.5312421 1.4419428 [7,] 1.31545660 -0.6830255 0.1264862 2.1502657 1.1636732 3.1368581 [8,] 0.88540350 -0.1517835 -0.0505944 1.2648621 1.4166455 2.5044270 [9,] -0.48064760 -0.7589172 -0.4047559 0.7336200 -0.1011889 0.2529724 [10,] -0.55653940 -1.0877815 0.5818366 -0.8854035 -1.5937263 -0.5059448 [11,] -0.48064760 -0.5059448 0.5312420 -0.2782696 -1.3154566 0.9612952 [12,] 0.63243110 1.7455097 0.6324311 1.7202125 0.5312421 2.8332912 [,7] [,8] [,9] [,10] [1,] 2.2514546 0.5565394 2.479130 NA [2,] 1.7708070 1.7961042 2.681508 NA [3,] 3.2633443 0.9865924 4.730584 NA [4,] 5.1100430 1.3154567 2.985075 NA [5,] 0.3035669 -0.2276752 3.668100 NA [6,] 3.2886415 0.2782696 -2.555022 NA [7,] 8.9805211 -1.5431318 -2.099671 NA [8,] 1.0624842 -1.7961042 -5.868960 NA [9,] 3.4657223 0.6577283 -8.044523 NA [10,] -0.1517835 1.0624842 -3.870478 NA [11,] -0.5059448 -2.4285353 -3.364533 NA [12,] 1.4672400 5.0088541 NA 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/19le51257975842.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/2euc11257975842.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/3xqir1257975842.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/4zpoo1257975842.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] [1,] 100.0000 99.9494 102.0744 102.0238 102.6309 102.8839 103.0104 104.3258 [2,] 104.3258 105.2112 105.3124 105.5401 106.7038 106.6026 105.7172 105.8437 [3,] 113.1040 115.1783 121.0726 123.0711 123.3999 122.9193 122.3881 123.5517 [4,] 136.6304 138.8819 140.6527 143.9160 149.0261 149.3296 152.6183 161.5988 [5,] 172.6031 175.0822 177.7637 182.4943 185.4794 189.1475 186.5925 184.4928 [,9] [,10] [,11] [,12] [1,] 104.1740 103.4151 102.3273 101.8214 [2,] 105.7931 105.3883 105.9701 106.5014 [3,] 124.9684 124.8672 123.2735 121.9580 [4,] 162.6613 166.1270 165.9752 163.3443 [5,] 178.6238 170.5793 170.0228 167.5942 $n [1] 9 9 9 9 9 9 9 9 9 9 9 9 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 96.09023 97.44515 102.4601 102.8598 101.1102 100.4164 97.68682 [2,] 130.11772 132.91154 139.6851 143.2824 145.6897 145.4222 147.08930 [,8] [,9] [,10] [,11] [,12] [1,] 94.18737 95.0178 92.87815 91.6708 92.02074 [2,] 152.91610 154.9190 156.85623 154.8761 151.89527 $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(100, 104.3258285, 113.1039717, 136.6304073, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/54g2l1257975842.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] [1,] -0.05059449 0.0000000 -0.9359980 -1.0877814 -0.4806476 -1.5937263 [2,] 0.55653940 0.6830255 -0.0505945 0.3288641 -0.1011890 -0.5312421 [3,] 2.07437390 1.7708070 0.4553504 0.6324310 0.1517834 0.1264862 [4,] 2.25145460 2.1249684 1.9984822 1.3154567 0.3035669 0.3035669 [5,] 3.18745250 2.6815078 4.7305843 1.3154567 0.6577283 1.4419428 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -2.0996711 -1.7961042 -0.7589172 -1.5937263 -2.4285353 0.5312421 [2,] -0.6830255 -0.1517835 -0.4806476 -1.0877815 -1.3154566 0.6324311 [3,] 1.1636732 0.8854035 -0.1011889 -0.5565394 -0.5059448 1.5937263 [4,] 2.1502657 1.2648621 0.6577283 -0.1517835 -0.2782696 2.2894005 [5,] 3.1368581 2.5044270 0.7336200 1.0624842 0.9612952 2.8332912 $n [1] 9 9 9 9 9 9 9 9 9 9 9 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1.181719 1.011384 -0.623830 0.1128256 -0.06138804 -0.3131799 -0.3285268 [2,] 2.967029 2.530230 1.534531 1.1520364 0.36495484 0.5661523 2.6558732 [,8] [,9] [,10] [,11] [,12] [1,] 0.1393035 -0.7007335 -1.04949835 -1.05219662 0.6681194 [2,] 1.6315035 0.4983557 -0.06358045 0.04030702 2.5193331 $out [1] 5.894257 5.110043 2.985075 3.668100 3.288642 -2.555022 8.980521 [8] -5.868960 3.465722 -8.044523 -3.870478 -3.364533 5.008854 $group [1] 2 4 4 5 6 6 7 8 9 9 10 11 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-0.0505944899999946, 0.55653940000002, 2.07437390000000, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6e40t1257975842.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] [1,] 99.9494 101.8214 103.5669 107.1338 121.0726 122.4892 136.6304 166.9365 [2,] 102.0491 103.7946 104.8065 108.2469 121.5153 125.8411 142.2843 167.9484 [3,] 102.9471 104.6673 105.6792 109.3600 122.9952 127.1060 150.9739 169.6560 [4,] 104.2499 105.4263 105.8184 111.5229 123.4758 132.9623 164.0653 170.8196 [5,] 105.2112 106.7038 106.5014 112.5474 124.9684 133.7971 166.1270 171.6418 [,9] [,10] [1,] 163.3443 NA [2,] 171.5912 NA [3,] 178.1938 NA [4,] 184.9861 NA [5,] 189.1475 NA $n [1] 12 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 101.9433 103.9231 105.2177 107.8658 122.1010 123.8580 141.0395 168.3464 [2,] 103.9510 105.4116 106.1408 110.8542 123.8894 130.3540 160.9084 170.9655 [,9] [,10] [1,] 172.0843 NA [2,] 184.3033 NA $out [1] 113.1040 115.1783 $group [1] 5 5 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" NA Warning message: In bxp(list(stats = c(99.94940551, 102.04907665, 102.9471288, 104.24993675, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7gt841257975842.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,] 125.1989 121.0726 120.4781 [2,] 129.1157 121.5153 123.8553 [3,] 131.1016 122.9952 128.5669 [4,] 131.8281 123.4758 134.5750 [5,] 132.8190 124.9684 135.9727 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 129.8644 122.1010 123.6776 [2,] 132.3387 123.8894 133.4562 $out [1] 113.1040 115.1783 $group [1] 2 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(125.198864444444, 129.115720838889, 131.101554377778, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/19le51257975842.ps tmp/19le51257975842.png") > system("convert tmp/2euc11257975842.ps tmp/2euc11257975842.png") > system("convert tmp/3xqir1257975842.ps tmp/3xqir1257975842.png") > system("convert tmp/4zpoo1257975842.ps tmp/4zpoo1257975842.png") > system("convert tmp/54g2l1257975842.ps tmp/54g2l1257975842.png") > system("convert tmp/6e40t1257975842.ps tmp/6e40t1257975842.png") > system("convert tmp/7gt841257975842.ps tmp/7gt841257975842.png") > > > proc.time() user system elapsed 1.214 0.991 2.182