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(-706.026542296181 + ,-59.4405579520267 + ,-969.833720383066 + ,1994.56903063601 + ,102.064814516947 + ,1282.18411025837 + ,-307.754606959943 + ,225.069461418321 + ,-862.471322666482 + ,-2850.16172979440 + ,-3875.76341994283 + ,2996.98129088558 + ,1931.51753614196 + ,488.179241950267 + ,-601.581627830989 + ,-3495.10356420977 + ,265.163825175723 + ,640.313274555833 + ,-3418.45097379814 + ,2072.03398243007 + ,2662.9984145817 + ,-1064.81719671349 + ,-9.57987333730985 + ,-2208.8867746199 + ,-4275.21498681124 + ,381.837227886395 + ,1914.99171248840 + ,1019.08609621188 + ,884.745394297871 + ,-1702.08983953302 + ,-2064.71644375347 + ,-2452.89205796152 + ,1313.39939683182 + ,-569.737683944086 + ,-187.9171774845 + ,-128.334928612711 + ,-2601.14827891639 + ,581.657297679053 + ,-299.792761503996 + ,2868.57371158154 + ,2723.0124492194 + ,-1846.18418996859 + ,-6526.92686969168 + ,-3652.78139505582 + ,2790.5729463211 + ,-7703.28365973859 + ,-1373.42544129638 + ,-2976.48819823291 + ,5585.1945937738 + ,-2568.07437379205 + ,-4750.48266905558 + ,1324.51875012340 + ,-1689.96832952026 + ,-7068.09981170815 + ,1855.55513995599 + ,2269.56153308754 + ,-4661.24493442205 + ,3260.97170059017 + ,2160.39821896397 + ,4793.18303923311 + ,-2209.06420471032 + ,-2213.50510931699 + ,-3114.14394163762 + ,3411.46075067066 + ,-4693.81993673823 + ,5920.52037235391 + ,-2077.91680939984 + ,-3642.83868232441 + ,332.648590721953 + ,2734.67402718683 + ,7903.1305482291 + ,5491.06198558059 + ,4354.83029956606 + ,1619.51905250655 + ,5543.22938817069 + ,-366.008933224093 + ,-3459.03179080274 + ,1134.72876550400 + ,-2525.56367971036 + ,-926.76612281528 + ,-2956.34221718809 + ,-885.143314774538) > 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] 82 > (np <- floor(n / par1)) [1] 6 > 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] 7 7 7 7 7 7 7 7 7 7 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -706.02654 1931.517536 -4275.2150 -2601.1483 5585.195 -2209.0642 [2,] -59.44056 488.179242 381.8372 581.6573 -2568.074 -2213.5051 [3,] -969.83372 -601.581628 1914.9917 -299.7928 -4750.483 -3114.1439 [4,] 1994.56903 -3495.103564 1019.0861 2868.5737 1324.519 3411.4608 [5,] 102.06481 265.163825 884.7454 2723.0124 -1689.968 -4693.8199 [6,] 1282.18411 640.313275 -1702.0898 -1846.1842 -7068.100 5920.5204 [7,] -307.75461 -3418.450974 -2064.7164 -6526.9269 1855.555 -2077.9168 [8,] 225.06946 2072.033982 -2452.8921 -3652.7814 2269.562 -3642.8387 [9,] -862.47132 2662.998415 1313.3994 2790.5729 -4661.245 332.6486 [10,] -2850.16173 -1064.817197 -569.7377 -7703.2837 3260.972 2734.6740 [11,] -3875.76342 -9.579873 -187.9172 -1373.4254 2160.398 7903.1305 [12,] 2996.98129 -2208.886775 -128.3349 -2976.4882 4793.183 5491.0620 [,7] [1,] 4354.8303 [2,] 1619.5191 [3,] 5543.2294 [4,] -366.0089 [5,] -3459.0318 [6,] 1134.7288 [7,] -2525.5637 [8,] -926.7661 [9,] -2956.3422 [10,] -885.1433 [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 646.5860 -1443.3383 4657.05221 3182.8056 -8153.2690 -4.440905 [2,] -910.3932 -1089.7609 1533.15448 -881.4501 -2182.4083 -900.638832 [3,] 2964.4028 -2893.5219 -895.90562 3168.3665 6075.0014 6525.604692 [4,] -1892.5042 3760.2674 -134.34070 -145.5613 -3014.4871 -8105.280687 [5,] 1180.1193 375.1494 -2586.83523 -4569.1966 -5378.1315 10614.340309 [6,] -1589.9387 -4058.7642 -362.62660 -4680.7427 8923.6550 -7998.437182 [7,] 532.8241 5490.4850 -388.17561 2874.1455 414.0064 -1564.921873 [8,] -1087.5408 590.9644 3766.29145 6443.3543 -6930.8065 3975.487273 [9,] -1987.6904 -3727.8156 -1883.13708 -10493.8566 7922.2166 2402.025436 [10,] -1025.6017 1055.2373 381.82051 6329.8582 -1100.5735 5168.456521 [11,] 6872.7447 -2199.3069 59.58225 -1603.0628 2632.7848 -2412.068563 [12,] -1065.4638 -2066.3282 -2472.81335 8561.6828 -7002.2472 -1136.231686 [,7] [1,] -2735.311 [2,] 3923.710 [3,] -5909.238 [4,] -3093.023 [5,] 4593.761 [6,] -3660.292 [7,] 1598.798 [8,] -2029.576 [9,] 2071.199 [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/html/rcomp/tmp/1m2ab1260307039.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/2yx951260307039.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/3sng41260307039.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/4ofxn1260307039.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] [1,] -4275.2150 -2568.0744 -4750.4827 -366.0089 -4693.8199 -1846.1842 [2,] -2405.1062 -1136.4728 -2041.9888 326.5386 -2574.5001 -1774.1370 [3,] -706.0265 381.8372 -601.5816 1324.5188 102.0648 640.3133 [4,] 3143.1739 534.9183 807.5995 2431.5714 574.9546 1208.4564 [5,] 5585.1946 1619.5191 1914.9917 3411.4608 2723.0124 1282.1841 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -3418.4510 -3652.7814 -4661.2449 -2850.1617 -3875.76342 -2976.488 [2,] -2972.0073 -3047.8654 -1909.4068 -1957.4895 -1373.42544 -2208.887 [3,] -2077.9168 -926.7661 332.6486 -885.1433 -98.74853 1434.323 [4,] -1186.2355 1148.5517 1988.1989 1082.4682 2160.39822 4793.183 [5,] -307.7546 2269.5615 2790.5729 3260.9717 2160.39822 5491.062 $n [1] 7 7 7 7 7 7 7 7 7 7 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -4019.370 -616.2906 -2303.310 67.42713 -1778.739 -1140.843 -3144.351 [2,] 2607.317 1379.9650 1100.147 2581.61037 1982.868 2421.470 -1011.483 [,8] [,9] [,10] [,11] [,12] [1,] -3432.799 -1994.939 -2700.5570 -2378.179 -3082.238 [2,] 1579.266 2660.236 930.2703 2180.682 5950.884 $out [1] 5543.229 -3495.104 -7068.100 5920.520 -6526.927 1855.555 -7703.284 [8] 7903.131 $group [1] 3 4 6 6 7 7 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-4275.21498681124, -2405.10624181336, -706.026542296181, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5k7fp1260307039.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,] -2735.311247 -2182.4083 -5909.238 -3093.0229 -5378.1315 -7998.4372 [2,] -2089.324771 -1000.0770 -1894.714 -3053.7550 -3578.0159 -4369.7535 [3,] -4.440905 -900.6388 2964.403 -1892.5042 375.1494 -3660.2924 [4,] 1914.695780 325.8522 4621.684 -139.9510 2886.9399 -976.2827 [5,] 4657.052215 1533.1545 6525.605 3760.2674 10614.3403 -362.6266 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -1564.92187 -6930.8065 -10493.857 -1100.5735 -2412.0686 -2472.813 [2,] 12.91539 -1558.5584 -2857.753 -1025.6017 -2199.3069 -2472.813 [3,] 532.82407 590.9644 -1883.137 718.5289 -771.7403 -1601.280 [4,] 2236.47152 3870.8894 2236.612 5168.4565 2632.7848 -1065.464 [5,] 5490.48496 6443.3543 7922.217 6329.8582 6872.7447 -1065.464 $n [1] 7 7 7 7 7 7 7 7 7 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2395.577 -1692.4624 -927.0848 -3632.5809 -3485.618 -5686.818 -795.0478 [2,] 2386.696 -108.8153 6855.8903 -152.4275 4235.917 -1633.766 1860.6959 [,8] [,9] [,10] [,11] [,12] [1,] -2651.414 -4925.410 -3276.839 -3888.595 -2509.066 [2,] 3833.343 1159.136 4713.897 2345.115 -693.494 $out [1] -8153.269 3923.710 -8105.281 8923.655 8561.683 -7002.247 $group [1] 1 2 4 6 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-2735.31124705951, -2089.3247706256, -4.44090460666985, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ifcb1260307039.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] [1,] -2850.1617 -3495.1036 -4275.2150 -7703.284 -7068.100 -4693.8199 -3459.0318 [2,] -916.1525 -1636.8520 -1883.4031 -3314.635 -3614.660 -2663.8245 -2525.5637 [3,] -183.5976 127.7920 -158.1261 -1609.805 1590.037 -872.6341 -625.5761 [4,] 753.6268 1285.9154 951.9157 1652.335 2765.267 4451.2614 1619.5191 [5,] 2996.9813 2662.9984 1914.9917 2868.574 5585.195 7903.1305 5543.2294 $n [1] 12 12 12 12 12 12 10 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -945.1951 -1205.302 -1451.334 -3875.2730 -1319.890 -4117.873 -2696.625 [2,] 578.0000 1460.886 1135.082 655.6634 4499.964 2372.604 1445.472 $out [1] -3875.763 $group [1] 1 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(-2850.1617297944, -916.152521524774, -183.597582455985, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7sshe1260307039.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,] -2152.2535 -2077.916809 -2079.1214 [2,] -855.4605 -795.584929 -783.4258 [3,] -243.4610 1.658145 -291.8088 [4,] 533.3146 511.075251 319.7307 [5,] 1327.9194 1434.323181 1379.0550 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -876.8907 -594.3183 -794.9659 [2,] 389.9686 597.6346 211.3483 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-2152.25346333678, -855.460489648028, -243.461038541175, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1m2ab1260307039.ps tmp/1m2ab1260307039.png") > system("convert tmp/2yx951260307039.ps tmp/2yx951260307039.png") > system("convert tmp/3sng41260307039.ps tmp/3sng41260307039.png") > system("convert tmp/4ofxn1260307039.ps tmp/4ofxn1260307039.png") > system("convert tmp/5k7fp1260307039.ps tmp/5k7fp1260307039.png") > system("convert tmp/6ifcb1260307039.ps tmp/6ifcb1260307039.png") > system("convert tmp/7sshe1260307039.ps tmp/7sshe1260307039.png") > > > proc.time() user system elapsed 1.190 0.942 1.451