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(-913.737775198605,23.1096757663404,456.385648731541,889.128571736973,836.385206042615,-2533.95352697662,-719.35167446997,-3745.06373199603,613.859653253457,659.119496405595,1059.93335663275,165.153227385284,-1259.25875793442,352.916648547945,-952.712549002952,3949.62792687456,2706.88419408659,-1983.44023600318,-5671.75150480296,-4491.18740262067,1977.74056046666,-7402.19548345163,-804.84956566237,-2561.90617124258,9086.02687097775,-2880.53187045752,-5510.78885909259,271.144630465973,-2568.10406759938,-6098.18357994379,5134.60266836802,4903.82157738512,-6164.11643961637,5797.99522829894,2868.48534028214,6908.93382436188,-2973.44012233353,-2043.11546028487,-2440.32738835400,1872.19672137889,-5813.85059764123,9609.8752579018,178.04755499601,-2432.7338999834,-214.086753766809,5446.18743629832,8510.17950026123,5524.61176329784,2640.19381267409,2753.13135132188,8434.31458855555,-3048.13516086747,-2363.04153284340,1655.57706698231,-2477.35519994173) > 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] 55 > (np <- floor(n / par1)) [1] 4 > 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] 5 5 5 5 5 5 5 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -913.73778 -1259.2588 9086.0269 -2973.4401 2640.194 [2,] 23.10968 352.9166 -2880.5319 -2043.1155 2753.131 [3,] 456.38565 -952.7125 -5510.7889 -2440.3274 8434.315 [4,] 889.12857 3949.6279 271.1446 1872.1967 -3048.135 [5,] 836.38521 2706.8842 -2568.1041 -5813.8506 -2363.042 [6,] -2533.95353 -1983.4402 -6098.1836 9609.8753 1655.577 [7,] -719.35167 -5671.7515 5134.6027 178.0476 -2477.355 [8,] -3745.06373 -4491.1874 4903.8216 -2432.7339 NA [9,] 613.85965 1977.7406 -6164.1164 -214.0868 NA [10,] 659.11950 -7402.1955 5797.9952 5446.1874 NA [11,] 1059.93336 -804.8496 2868.4853 8510.1795 NA [12,] 165.15323 -2561.9062 6908.9338 5524.6118 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] 936.84745 1612.175 -11966.5587 930.3247 112.9375 [2,] 433.27597 -1305.629 -2630.2570 -397.2119 5681.1832 [3,] 432.74292 4902.340 5781.9335 4312.5241 -11482.4497 [4,] -52.74337 -1242.744 -2839.2487 -7686.0473 685.0936 [5,] -3370.33873 -4690.324 -3530.0795 15423.7259 4018.6186 [6,] 1814.60185 -3688.311 11232.7862 -9431.8277 -4132.9323 [7,] -3025.71206 1180.564 -230.7811 -2610.7815 NA [8,] 4358.92339 6468.928 -11067.9380 2218.6471 NA [9,] 45.25984 -9379.936 11962.1117 5660.2742 NA [10,] 400.81386 6597.346 -2929.5099 3063.9921 NA [11,] -894.78013 -1757.057 4040.4485 -2985.5677 NA [12,] -1424.41199 11647.933 -9882.3739 -2884.4180 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/1lg2c1260015494.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/2zgsw1260015494.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/3jxxn1260015494.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/4jxfc1260015494.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,] -2973.4401 -2880.53187 -5510.7889 271.1446 -5813.8506 -6098.184 [2,] -1259.2588 -2043.11546 -2440.3274 271.1446 -2568.1041 -2533.954 [3,] -913.7378 23.10968 -952.7125 889.1286 -2363.0415 -1983.440 [4,] 2640.1938 352.91665 456.3856 1872.1967 836.3852 1655.577 [5,] 2640.1938 2753.13135 456.3856 3949.6279 2706.8842 1655.577 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -5671.7515 -4491.187 -6164.1164 -7402.195 -804.8496 -2561.906 [2,] -2477.3552 -4118.126 -3189.1016 -3371.538 127.5419 -1198.376 [3,] -719.3517 -3088.899 199.8864 3052.653 1964.2093 2844.882 [4,] 178.0476 1235.544 1295.8001 5622.091 5689.3324 6216.773 [5,] 178.0476 4903.822 1977.7406 5797.995 8510.1795 6908.934 $n [1] 5 5 5 5 5 5 5 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -3669.081 -1669.921 -2999.523 -242.1708 -4768.64508 -4943.7520 -2595.653 [2,] 1841.606 1716.140 1094.098 2020.4279 42.56201 976.8715 1156.949 [,8] [,9] [,10] [,11] [,12] [1,] -7318.298 -3343.186 -4052.314 -2429.605 -3013.085 [2,] 1140.500 3742.959 10157.621 6358.024 8702.850 $out [1] 9086.027 8434.315 -3048.135 9609.875 5134.603 $group [1] 1 3 4 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-2973.44012233353, -1259.25875793442, -913.737775198605, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5awq01260015494.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,] 112.9375 -2630.2570 432.7429 -2839.24870 -4690.324 -9431.828 -3025.7121 [2,] 112.9375 -1305.6292 432.7429 -2839.24870 -3530.080 -4132.932 -2818.2468 [3,] 930.3247 -397.2119 4312.5241 -1242.74373 -3370.339 -3688.311 -1420.7813 [4,] 936.8475 433.2760 4902.3405 -52.74337 4018.619 1814.602 474.8915 [5,] 1612.1754 433.2760 5781.9335 685.09363 4018.619 1814.602 1180.5641 [,8] [,9] [,10] [,11] [,12] [1,] -11067.938 -9379.936 -2929.510 -2985.568 -9882.374 [2,] -4424.645 -4667.338 -1264.348 -2371.312 -6383.396 [3,] 3288.785 2852.767 1732.403 -1325.918 -2154.415 [4,] 5413.926 8811.193 4830.669 1572.834 5111.761 [5,] 6468.928 11962.112 6597.346 4040.448 11647.933 $n [1] 5 5 5 5 5 5 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 348.152 -1625.9179 1154.318 -3211.6814 -8704.230 -7890.8239 -4022.360 [2,] 1512.497 831.4941 7470.730 726.1939 1963.552 514.2014 1180.798 [,8] [,9] [,10] [,11] [,12] [1,] -4483.686 -7795.272 -3082.660 -4441.794 -11235.589 [2,] 11061.256 13500.807 6547.466 1789.957 6926.759 $out [1] -11966.559 5681.183 -11482.450 -7686.047 15423.726 11232.786 $group [1] 1 2 3 4 5 6 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(112.93753864779, 112.93753864779, 930.32466204866, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/68fyg1260015494.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] [1,] -2533.9535 -7402.195 -6164.116 -5813.8506 -3048.135 [2,] -816.5447 -3526.547 -4195.660 -2436.5306 -2420.198 [3,] 310.7694 -1105.986 1569.815 -18.0196 1655.577 [4,] 747.7524 1165.329 5466.299 5485.3996 2696.663 [5,] 1059.9334 3949.628 9086.027 9609.8753 8434.315 $n [1] 12 12 12 12 7 $conf [,1] [,2] [,3] [,4] [,5] [1,] -402.7170 -3245.982 -2837.070 -3631.265 -1400.130 [2,] 1024.2558 1034.010 5976.699 3595.226 4711.284 $out [1] -3745.064 $group [1] 1 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-2533.95352697662, -816.544724834288, 310.769438058413, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7n8au1260015494.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,] -1441.29086 -3088.899 -2265.0948 [2,] -828.90619 -1468.076 -928.9152 [3,] 63.67464 -348.121 -406.2852 [4,] 1220.61674 1426.669 1580.3179 [5,] 2908.43716 3052.653 2677.0403 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -871.1265 -1668.4338 -1550.7632 [2,] 998.4758 972.1918 738.1928 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-1441.29086430374, -828.906188042946, 63.6746422798068, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1lg2c1260015494.ps tmp/1lg2c1260015494.png") > system("convert tmp/2zgsw1260015494.ps tmp/2zgsw1260015494.png") > system("convert tmp/3jxxn1260015494.ps tmp/3jxxn1260015494.png") > system("convert tmp/4jxfc1260015494.ps tmp/4jxfc1260015494.png") > system("convert tmp/5awq01260015494.ps tmp/5awq01260015494.png") > system("convert tmp/68fyg1260015494.ps tmp/68fyg1260015494.png") > system("convert tmp/7n8au1260015494.ps tmp/7n8au1260015494.png") > > > proc.time() user system elapsed 1.204 0.956 1.479