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(-1575.00887827752 + ,-2263.82750005582 + ,7045.75768165492 + ,-2035.83149854133 + ,3886.02985519469 + ,158.539675832947 + ,960.457190377358 + ,-337.841379165606 + ,-7109.61280621946 + ,-7043.92461948593 + ,4685.51755601771 + ,977.408177407792 + ,918.328545859417 + ,589.375354065341 + ,-6891.67878957285 + ,1533.64184150137 + ,2210.55324487903 + ,-9280.29800578057 + ,8034.36112956057 + ,3659.22886246375 + ,7117.41595284467 + ,-2633.18383716224 + ,-4523.01908028135 + ,-10249.2788035485 + ,1517.80356687876 + ,1484.15119995578 + ,-1214.46781890457 + ,2169.24845396691 + ,-2253.92155265332 + ,-2631.74888936717 + ,-3849.26121145593 + ,821.620255775573 + ,-9408.87260464867 + ,1241.88804485160 + ,1179.93724525252 + ,-1962.37832237001 + ,516.653311386494 + ,-2767.97210236296 + ,5629.81911313198 + ,5803.53000967108 + ,-1873.10106680445 + ,-7023.12305076338 + ,-1722.94778568371 + ,-3011.92640473933 + ,-17439.1442995944 + ,-1770.63103827571 + ,-6148.68941133731 + ,7691.62364275427 + ,-4994.18050152824 + ,-4899.69881288656 + ,5524.45031429093 + ,-6451.65580766626 + ,-9155.30558471285 + ,9268.49118895307 + ,4067.36942971564 + ,-15460.4752115729 + ,9214.1160150789 + ,4009.46524142609 + ,8598.3415209042 + ,1171.18094385286 + ,-1585.96191465526 + ,-2329.27119138406 + ,5674.43930504972 + ,-9680.56803374957 + ,13269.3191881056 + ,-4706.51689196404 + ,-3552.85443447458 + ,-3123.75770966460 + ,2912.26032459412 + ,12338.9131048695 + ,10340.0965258702 + ,6543.69991954822 + ,7116.01012281237) > 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] 73 > (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 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -1575.0089 918.3285 1517.8036 516.6533 -4994.181 -1585.962 [2,] -2263.8275 589.3754 1484.1512 -2767.9721 -4899.699 -2329.271 [3,] 7045.7577 -6891.6788 -1214.4678 5629.8191 5524.450 5674.439 [4,] -2035.8315 1533.6418 2169.2485 5803.5300 -6451.656 -9680.568 [5,] 3886.0299 2210.5532 -2253.9216 -1873.1011 -9155.306 13269.319 [6,] 158.5397 -9280.2980 -2631.7489 -7023.1231 9268.491 -4706.517 [7,] 960.4572 8034.3611 -3849.2612 -1722.9478 4067.369 -3552.854 [8,] -337.8414 3659.2289 821.6203 -3011.9264 -15460.475 -3123.758 [9,] -7109.6128 7117.4160 -9408.8726 -17439.1443 9214.116 2912.260 [10,] -7043.9246 -2633.1838 1241.8880 -1770.6310 4009.465 12338.913 [11,] 4685.5176 -4523.0191 1179.9372 -6148.6894 8598.342 10340.097 [12,] 977.4082 -10249.2788 -1962.3783 7691.6236 1171.181 6543.700 [,7] [1,] 7116.01 [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [9,] NA [10,] NA [11,] NA [12,] NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [1,] -688.81862 -328.9532 -33.65237 -3284.6254 94.48169 -743.3093 [2,] 9309.58518 -7481.0541 -2698.61902 8397.7912 10424.14913 8003.7105 [3,] -9081.58918 8425.3206 3383.71627 173.7109 -11976.10612 -15355.0073 [4,] 5921.86135 676.9114 -4423.17001 -7676.6311 -2703.64978 22949.8872 [5,] -3727.49018 -11490.8513 -377.82734 -5150.0220 18423.79677 -17975.8361 [6,] 801.91751 17314.6591 -1217.51232 5300.1753 -5201.12176 1153.6625 [7,] -1298.29857 -4375.1323 4670.88147 -1288.9786 -19527.84464 429.0967 [8,] -6771.77143 3458.1871 -10230.49286 -14427.2179 24674.59123 6036.0180 [9,] 65.68819 -9750.5998 10650.76065 15668.5133 -5204.65077 9426.6528 [10,] 11729.44218 -1889.8352 -61.95080 -4378.0584 4588.87628 -1998.8166 [11,] -3708.10938 -5726.2597 -3142.31557 13840.3131 -7427.16058 -3796.3966 [12,] -59.07963 11767.0824 2479.03163 -12685.8041 -2757.14286 572.3102 [,7] [1,] NA [2,] NA [3,] NA [4,] NA [5,] NA [6,] NA [7,] NA [8,] NA [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/html/rcomp/tmp/1s46s1259947993.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/2bbia1259947993.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/371lm1259947993.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/4ddlr1259947993.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,] -4994.1805 -4899.6988 -6891.679 -9680.5680 -9155.3056 -9280.2980 [2,] -1580.4854 -2767.9721 -1214.468 -6451.6558 -2253.9216 -7023.1231 [3,] 516.6533 -2296.5493 5577.135 -251.0948 168.7261 -3669.1329 [4,] 1218.0661 589.3754 5674.439 2169.2485 3886.0299 158.5397 [5,] 1517.8036 1484.1512 7045.758 5803.5300 3886.0299 9268.4912 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -3849.2612 -3123.7577 -17439.144 -7043.9246 -6148.689 -10249.279 [2,] -3552.8544 -3123.7577 -9408.873 -2633.1838 -4523.019 -1962.378 [3,] -381.2453 -1674.8839 -2098.676 -264.3715 2932.727 1074.295 [4,] 4067.3694 821.6203 7117.416 4009.4652 8598.342 6543.700 [5,] 8034.3611 3659.2289 9214.116 12338.9131 10340.097 7691.624 $n [1] 7 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -1154.596 -4462.1469 1133.567 -5811.857 -3791.741 -8301.5373 -5296.536 [2,] 2187.903 -130.9518 10020.702 5309.667 4129.193 963.2715 4534.045 [,8] [,9] [,10] [,11] [,12] [1,] -4219.7801 -12758.666 -4549.095 -5530.974 -4412.401 [2,] 870.0123 8561.314 4020.352 11396.429 6560.990 $out [1] 7116.01 13269.32 -15460.48 $group [1] 1 5 8 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-4994.18050152824, -1580.48539646639, 516.653311386494, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5190v1259947993.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] [1,] -743.30928 -7481.054 -15355.007 -7676.631 -17975.8361 -5201.122 -4375.1323 [2,] -743.30928 -2698.619 -11976.106 -4423.170 -11490.8513 -1217.512 -4375.1323 [3,] -508.88591 8200.751 -4453.939 -1013.369 -4438.7561 977.790 -1293.6386 [4,] -33.65237 9309.585 3383.716 5921.861 -377.8273 5300.175 429.0967 [5,] 94.48169 10424.149 8425.321 5921.861 -377.8273 5300.175 4670.8815 [,8] [,9] [,10] [,11] [,12] [1,] -14427.218 -9750.600 -4378.058 -7427.161 -2757.1429 [2,] -10230.493 -5204.651 -1998.817 -5726.260 -2757.1429 [3,] -1656.792 4746.170 -975.893 -3752.253 256.6153 [4,] 6036.018 10650.761 4588.876 -3142.316 2479.0316 [5,] 24674.591 15668.513 11729.442 -3142.316 2479.0316 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -966.63753 455.0713 -14361.521 -7686.249 -11607.015 -3226.329 -4392.521 [2,] -51.13428 15946.4304 5453.642 5659.510 2729.503 5181.909 1805.244 [,8] [,9] [,10] [,11] [,12] [1,] -12149.217 -5481.082 -5225.168 -5418.980 -3120.886 [2,] 8835.633 14973.423 3273.382 -2085.526 3634.117 $out [1] -3284.625 22949.887 18423.797 17314.659 -19527.845 13840.313 11767.082 [8] -12685.804 $group [1] 1 4 5 6 7 11 12 12 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-743.3092767288, -743.3092767288, -508.885906786188, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/67cv11259947993.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,] -7109.61281 -10249.279 -3849.2612 -7023.123 -15460.475 -9680.5680 7116.01 [2,] -2149.82950 -5707.349 -2442.8352 -4580.308 -5722.918 -3338.3061 7116.01 [3,] -89.65085 753.852 -196.4238 -1821.866 2590.323 663.1492 7116.01 [4,] 2431.71902 2934.891 1363.0196 3073.236 7061.396 8441.8982 7116.01 [5,] 7045.75768 8034.361 2169.2485 7691.624 9268.491 13269.3192 7116.01 $n [1] 12 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2179.326 -3187.932 -1932.300 -5312.699 -3240.688 -4709.881 7116.01 [2,] 2000.024 4695.636 1539.452 1668.967 8421.335 6036.179 7116.01 $out [1] -9408.873 -17439.144 $group [1] 3 4 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(-7109.61280621946, -2149.82949929857, -89.6508516663295, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/77mfo1259947993.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,] -2908.8586 -3669.1329 -3491.50199 [2,] -2033.4916 -1886.7801 -1387.53821 [3,] 464.7826 -257.7332 -41.78894 [4,] 1018.8417 795.4739 1320.39569 [5,] 2628.0533 2932.7274 3066.77299 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -927.4072 -1481.1272 -1276.896 [2,] 1856.9723 965.6609 1193.318 $out [1] 5577.135 $group [1] 2 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-2908.85859781719, -2033.49158547979, 464.782568346804, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1s46s1259947993.ps tmp/1s46s1259947993.png") > system("convert tmp/2bbia1259947993.ps tmp/2bbia1259947993.png") > system("convert tmp/371lm1259947993.ps tmp/371lm1259947993.png") > system("convert tmp/4ddlr1259947993.ps tmp/4ddlr1259947993.png") > system("convert tmp/5190v1259947993.ps tmp/5190v1259947993.png") > system("convert tmp/67cv11259947993.ps tmp/67cv11259947993.png") > system("convert tmp/77mfo1259947993.ps tmp/77mfo1259947993.png") > > > proc.time() user system elapsed 1.204 0.943 1.446