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(-1065.28076055322,-2278.75083745204,-13750.7320730183,768.487940174182,480.329421548304,2441.84534097982,-2017.12122801278,-3677.24413007058,2376.07844077739,1991.57956111441,8177.8663887305,-9931.36901798728,1054.78346037674,-6334.48912633556,-10318.8992855633,-2985.24764365473,-4927.33781231944,2075.02534260383,-3870.63118024471,-2524.15412506803,725.810856004154,-6465.06532185037,-11625.1356478526,11621.1659938228,1767.57466477775,-10664.0704054107,6189.98081853049,2343.11555395316,3876.36005581636,1472.93269355467,118.383083333785,820.674136304355,1542.41090621852,-3650.11001150461,4906.01253028907,40.3638959194068,-1585.17421898316,-124.437694132677,2791.13055865956,5387.15230614261,5445.44402054302,2593.37824514089,6516.19939614798,6292.59289700375,2088.48664589605,2475.88763066489,3427.36476553255,-1086.08268469240,1510.63701057434,1414.23900224041) > 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] 50 > (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 4 4 4 4 4 4 4 4 4 4 > arr [,1] [,2] [,3] [,4] [,5] [1,] -1065.2808 1054.7835 1767.57466 -1585.1742 1510.637 [2,] -2278.7508 -6334.4891 -10664.07041 -124.4377 1414.239 [3,] -13750.7321 -10318.8993 6189.98082 2791.1306 NA [4,] 768.4879 -2985.2476 2343.11555 5387.1523 NA [5,] 480.3294 -4927.3378 3876.36006 5445.4440 NA [6,] 2441.8453 2075.0253 1472.93269 2593.3782 NA [7,] -2017.1212 -3870.6312 118.38308 6516.1994 NA [8,] -3677.2441 -2524.1541 820.67414 6292.5929 NA [9,] 2376.0784 725.8109 1542.41091 2088.4866 NA [10,] 1991.5796 -6465.0653 -3650.11001 2475.8876 NA [11,] 8177.8664 -11625.1356 4906.01253 3427.3648 NA [12,] -9931.3690 11621.1660 40.36390 -1086.0827 NA > darr [,1] [,2] [,3] [,4] [,5] [1,] -1213.4701 -7389.273 -12431.6451 1460.73652 -96.398 [2,] -11471.9812 -3984.410 16854.0512 2915.56825 NA [3,] 14519.2200 7333.652 -3846.8653 2596.02175 NA [4,] -288.1585 -1942.090 1533.2445 58.29171 NA [5,] 1961.5159 7002.363 -2403.4274 -2852.06578 NA [6,] -4458.9666 -5945.657 -1354.5496 3922.82115 NA [7,] -1660.1229 1346.477 702.2911 -223.60650 NA [8,] 6053.3226 3249.965 721.7368 -4204.10625 NA [9,] -384.4989 -7190.876 -5192.5209 387.40098 NA [10,] 6186.2868 -5160.070 8556.1225 951.47713 NA [11,] -18109.2354 23246.302 -4865.6486 -4513.44745 NA [12,] 10986.1525 -9853.591 -1625.5381 2596.71970 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/1nejo1260223579.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/2ld8x1260223579.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/345ec1260223579.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/47nrh1260223579.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] [1,] -1585.174 -10664.0704 -13750.732 -2985.248 -4927.338 1472.933 -3870.631 [2,] -1065.281 -6334.4891 -12034.816 -1108.380 -2223.504 1773.979 -2943.876 [3,] 1054.783 -2278.7508 -3763.884 1555.802 2178.345 2258.435 -949.369 [4,] 1510.637 -124.4377 4490.556 3865.134 4660.902 2517.612 3317.291 [5,] 1767.575 1414.2390 6189.981 5387.152 5445.444 2593.378 6516.199 [,8] [,9] [,10] [,11] [,12] [1,] -3677.244 725.8109 -6465.0653 -11625.136 -9931.3690 [2,] -3100.699 1134.1109 -5057.5877 -4098.885 -5508.7259 [3,] -851.740 1815.4488 -829.2652 4166.689 -522.8594 [4,] 3556.634 2232.2825 2233.7336 6541.939 5830.7649 [5,] 6292.593 2376.0784 2475.8876 8177.866 11621.1660 $n [1] 5 5 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -765.3535 -6666.758 -16818.928 -2373.274 -3260.336 1670.965 -5895.691 [2,] 2874.9205 2109.256 9291.159 5484.878 7617.026 2845.905 3996.953 [,8] [,9] [,10] [,11] [,12] [1,] -6111.033 947.8932 -6589.409 -4239.563 -9481.057 [2,] 4407.553 2683.0044 4930.879 12572.940 8435.338 $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(-1585.17421898316, -1065.28076055322, 1054.78346037674, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5tosc1260223579.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,] -12431.645 -11471.981 -3846.8653 -1942.0902 -2852.0658 -5945.657 [2,] -7389.273 -7728.196 -625.4218 -1115.1243 -2627.7466 -5202.312 [3,] -1213.470 -534.421 4964.8367 -114.9334 -220.9557 -2906.758 [4,] -96.398 9884.810 10926.4358 795.7681 4481.9395 1284.136 [5,] 1460.737 16854.051 14519.2200 1533.2445 7002.3632 3922.821 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -1660.1229 -4204.106 -7190.876178 -5160.070 -18109.235 -9853.5913 [2,] -941.8647 -1741.185 -6191.698548 -2104.297 -11487.442 -5739.5647 [3,] 239.3423 1985.851 -2788.509899 3568.882 -4689.548 485.5908 [4,] 1024.3841 4651.644 1.451053 7371.205 9366.427 6791.4361 [5,] 1346.4771 6053.323 387.400985 8556.123 23246.302 10986.1525 $n [1] 5 4 4 4 4 4 4 4 4 4 4 4 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -6366.597 -14448.70 -4161.131 -1624.538 -5837.608 -8031.051 -1313.994 [2,] 3939.657 13379.85 14090.804 1394.672 5395.696 2217.535 1792.679 [,8] [,9] [,10] [,11] [,12] [1,] -3064.484 -7681.098 -3916.764 -21164.10 -9413.90 [2,] 7036.185 2104.078 11054.528 11785.01 10385.08 $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(-12431.6450701884, -7389.2725867123, -1213.47007689882, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6gp9u1260223579.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,] -9931.3690 -11625.1356 -3650.11001 -1585.1742 1414.239 [2,] -2977.9975 -6399.7772 79.37349 982.0245 1414.239 [3,] -292.4757 -3427.9394 1507.67180 2692.2544 1462.438 [4,] 2183.8290 890.2972 3109.73780 5416.2982 1510.637 [5,] 8177.8664 11621.1660 6189.98082 6516.1994 1510.637 $n [1] 12 12 12 12 2 $conf [,1] [,2] [,3] [,4] [,5] [1,] -2646.819 -6752.9913 125.5023 669.7524 1354.739 [2,] 2061.868 -102.8876 2889.8413 4714.7564 1570.137 $out [1] -13750.73 -10664.07 $group [1] 1 3 $names [1] "1" "2" "3" "4" NA Warning message: In bxp(list(stats = c(-9931.36901798728, -2977.99748376131, -292.475669502458, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/7xo3z1260223579.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,] -1411.9270 -3763.8844 -3768.0072 [2,] -625.4537 -900.5545 -750.9635 [3,] 282.2376 265.9620 20.8791 [4,] 1299.9520 1996.8968 1723.9223 [5,] 2145.7954 4166.6886 2694.1078 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -595.9529 -1055.585 -1107.933 [2,] 1160.4281 1587.509 1149.691 $out [1] -3597.502 -3772.130 $group [1] 1 1 $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-1411.92703539392, -625.453744314144, 282.237612890432, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1nejo1260223579.ps tmp/1nejo1260223579.png") > system("convert tmp/2ld8x1260223579.ps tmp/2ld8x1260223579.png") > system("convert tmp/345ec1260223579.ps tmp/345ec1260223579.png") > system("convert tmp/47nrh1260223579.ps tmp/47nrh1260223579.png") > system("convert tmp/5tosc1260223579.ps tmp/5tosc1260223579.png") > system("convert tmp/6gp9u1260223579.ps tmp/6gp9u1260223579.png") > system("convert tmp/7xo3z1260223579.ps tmp/7xo3z1260223579.png") > > > proc.time() user system elapsed 1.193 0.947 1.665