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(10284.5,12792,12823.61538,13845.66667,15335.63636,11188.5,13633.25,12298.46667,15353.63636,12696.15385,12213.93333,13683.72727,11214.14286,13950.23077,11179.13333,11801.875,11188.82353,16456.27273,11110.0625,16530.69231,10038.41176,11681.25,11148.88235,8631,9386.444444,9764.736842,12043.75,12948.06667,10987.125,11648.3125,10633.35294,10219.3,9037.6,10296.31579,11705.41176,10681.94444,9362.947368,11306.35294,10984.45,10062.61905,8118.583333,8867.48,8346.72,8529.307692,10697.18182,8591.84,8695.607143,8125.571429,7009.758621,7883.466667,7527.645161,6763.758621,6682.333333,7855.681818,6738.88,7895.434783,6361.884615,6935.956522,8344.454545,9107.944444) > 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] 60 > (np <- floor(n / par1)) [1] 5 > arr <- array(NA,dim=c(par1,np+1)) > ari <- array(0,dim=par1) > j <- 0 > for (i in 1:n) + { + j = j + 1 + ari[j] = ari[j] + 1 + arr[j,ari[j]] <- x[i] + if (j == par1) j = 0 + } > ari [1] 5 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 10284.50 11214.14 9386.444 9362.947 7009.759 NA [2,] 12792.00 13950.23 9764.737 11306.353 7883.467 NA [3,] 12823.62 11179.13 12043.750 10984.450 7527.645 NA [4,] 13845.67 11801.88 12948.067 10062.619 6763.759 NA [5,] 15335.64 11188.82 10987.125 8118.583 6682.333 NA [6,] 11188.50 16456.27 11648.312 8867.480 7855.682 NA [7,] 13633.25 11110.06 10633.353 8346.720 6738.880 NA [8,] 12298.47 16530.69 10219.300 8529.308 7895.435 NA [9,] 15353.64 10038.41 9037.600 10697.182 6361.885 NA [10,] 12696.15 11681.25 10296.316 8591.840 6935.957 NA [11,] 12213.93 11148.88 11705.412 8695.607 8344.455 NA [12,] 13683.73 8631.00 10681.944 8125.571 9107.944 NA > arr.sd <- array(NA,dim=par1) > arr.range <- array(NA,dim=par1) > arr.iqr <- array(NA,dim=par1) > for (j in 1:par1) + { + arr.sd[j] <- sqrt(var(arr[j,],na.rm=TRUE)) + arr.range[j] <- max(arr[j,],na.rm=TRUE) - min(arr[j,],na.rm=TRUE) + arr.iqr[j] <- quantile(arr[j,],0.75,na.rm=TRUE) - quantile(arr[j,],0.25,na.rm=TRUE) + } > overall.sd <- sqrt(var(x)) > overall.range <- max(x) - min(x) > overall.iqr <- quantile(x,0.75) - quantile(x,0.25) > postscript(file="/var/www/html/rcomp/tmp/1bgr61258117081.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.sd,type='b',ylab='S.D.',main='Standard Deviation Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.sd,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/28foy1258117081.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.range,type='b',ylab='range',main='Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.range,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/3os5g1258117081.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > plot(arr.iqr,type='b',ylab='IQR',main='Interquartile Range Plot',xlab='Periodic Index') > mtext(paste('# blocks = ',np)) > abline(overall.iqr,0) > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/4epoc1258117081.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,] 9362.947 7883.467 10984.45 6763.759 6682.333 7855.682 6738.88 [2,] 9362.947 9764.737 10984.45 10062.619 8118.583 8867.480 8346.72 [3,] 9386.444 11306.353 11179.13 11801.875 10987.125 11188.500 10633.35 [4,] 10284.500 12792.000 12043.75 12948.067 11188.824 11648.312 11110.06 [5,] 11214.143 13950.231 12823.62 13845.667 15335.636 11648.312 13633.25 [,8] [,9] [,10] [,11] [,12] [1,] 7895.435 9037.60 6935.957 8344.455 8125.571 [2,] 8529.308 9037.60 8591.840 8695.607 8631.000 [3,] 10219.300 10038.41 10296.316 11148.882 9107.944 [4,] 12298.467 10697.18 11681.250 11705.412 10681.944 [5,] 16530.692 10697.18 12696.154 12213.933 13683.727 $n [1] 5 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 8735.278 9167.296 10430.63 9763.025 8817.701 9223.57 8680.782 [2,] 10037.611 13445.409 11927.63 13840.725 13156.549 13153.43 12585.924 [,8] [,9] [,10] [,11] [,12] [1,] 7556.022 8865.755 8113.346 9022.162 7658.752 [2,] 12882.578 11211.068 12479.285 13275.603 10557.137 $out [1] 7009.759 7527.645 16456.273 15353.636 6361.885 $group [1] 1 3 6 9 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(9362.947368, 9362.947368, 9386.444444, 10284.5, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5hjy91258117081.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] [1,] 10284.50 8631.00 9037.600 8118.583 6361.885 NA [2,] 12256.20 11129.47 9992.018 8438.014 6751.319 NA [3,] 12807.81 11201.48 10657.649 8781.544 7268.702 NA [4,] 13764.70 12876.05 11676.862 10379.900 7889.451 NA [5,] 15353.64 13950.23 12948.067 11306.353 9107.944 NA $n [1] 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 12119.77 10404.86 9889.18 7895.836 6749.593 NA [2,] 13495.84 11998.11 11426.12 9667.251 7787.811 NA $out [1] 16456.27 16530.69 $group [1] 2 2 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(10284.5, 12256.2, 12807.80769, 13764.69697, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6gkur1258117081.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) > z <- data.frame(cbind(arr.sd,arr.range,arr.iqr)) > names(z) <- list('S.D.','Range','IQR') > (boxplot(z,notch=TRUE,col='grey',ylab='Overall Variability',main='Notched Box Plots')) $stats [,1] [,2] [,3] [1,] 1563.248 3869.479 921.5526 [2,] 2138.270 5427.063 1855.2631 [3,] 2529.510 6480.567 2833.1401 [4,] 3301.429 8617.924 3048.7517 [5,] 3485.037 8991.752 3769.1590 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1998.986 5025.194 2288.782 [2,] 3060.035 7935.940 3377.498 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(1563.24770834841, 2138.27009823341, 2529.51028049812, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1bgr61258117081.ps tmp/1bgr61258117081.png") > system("convert tmp/28foy1258117081.ps tmp/28foy1258117081.png") > system("convert tmp/3os5g1258117081.ps tmp/3os5g1258117081.png") > system("convert tmp/4epoc1258117081.ps tmp/4epoc1258117081.png") > system("convert tmp/5hjy91258117081.ps tmp/5hjy91258117081.png") > system("convert tmp/6gkur1258117081.ps tmp/6gkur1258117081.png") > > > proc.time() user system elapsed 1.071 0.825 1.355