R version 2.8.0 (2008-10-20) Copyright (C) 2008 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(1664.81,2397.53,2840.71,3547.29,3752.96,3714.74,4349.61,3566.34,5021.82,6423.48,7600.60,19756.21,2499.81,5198.24,7225.14,4806.03,5900.88,4951.34,6179.12,4752.15,5496.43,5835.10,12600.08,28541.72,4717.02,5702.63,9957.58,5304.78,6492.43,6630.80,7349.62,8176.62,8573.17,9690.50,15151.84,34061.01,5921.10,5814.58,12421.25,6369.77,7609.12,7224.75,8121.22,7979.25,8093.06,8476.70,17914.66,30114.41,4826.64,6470.23,9638.77,8821.17,8722.37,10209.48,11276.55,12552.22,11637.39,13606.89,21822.11,45060.69,7615.03,9849.69,14558.40,11587.33,9332.56,13082.09,16732.78,19888.61,23933.38,25391.35,36024.80,80721.71,10243.24,11266.88,21826.84,17357.33,15997.79,18601.53,26155.15,28586.52,30505.41,30821.33,46634.38,104660.67) > 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] 84 > (np <- floor(n / par1)) [1] 7 > 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] 7 7 7 7 7 7 7 7 7 7 7 7 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1664.81 2499.81 4717.02 5921.10 4826.64 7615.03 10243.24 NA [2,] 2397.53 5198.24 5702.63 5814.58 6470.23 9849.69 11266.88 NA [3,] 2840.71 7225.14 9957.58 12421.25 9638.77 14558.40 21826.84 NA [4,] 3547.29 4806.03 5304.78 6369.77 8821.17 11587.33 17357.33 NA [5,] 3752.96 5900.88 6492.43 7609.12 8722.37 9332.56 15997.79 NA [6,] 3714.74 4951.34 6630.80 7224.75 10209.48 13082.09 18601.53 NA [7,] 4349.61 6179.12 7349.62 8121.22 11276.55 16732.78 26155.15 NA [8,] 3566.34 4752.15 8176.62 7979.25 12552.22 19888.61 28586.52 NA [9,] 5021.82 5496.43 8573.17 8093.06 11637.39 23933.38 30505.41 NA [10,] 6423.48 5835.10 9690.50 8476.70 13606.89 25391.35 30821.33 NA [11,] 7600.60 12600.08 15151.84 17914.66 21822.11 36024.80 46634.38 NA [12,] 19756.21 28541.72 34061.01 30114.41 45060.69 80721.71 104660.67 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/freestat/rcomp/tmp/125zs1244216580.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/freestat/rcomp/tmp/2lip61244216580.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/freestat/rcomp/tmp/37fa01244216580.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/freestat/rcomp/tmp/4l5p71244216580.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,] 1664.810 2397.530 2840.710 3547.290 3752.960 3714.74 4349.61 [2,] 3608.415 5450.435 8431.955 5055.405 6196.655 5791.07 6764.37 [3,] 4826.640 5814.580 9957.580 6369.770 7609.120 7224.75 8121.22 [4,] 6768.065 8159.960 13489.825 10204.250 9027.465 11645.78 14004.66 [5,] 10243.240 11266.880 14558.400 17357.330 9332.560 18601.53 16732.78 [,8] [,9] [,10] [,11] [,12] [1,] 3566.34 5021.820 5835.10 7600.60 19756.21 [2,] 6365.70 6794.745 7450.09 13875.96 29328.07 [3,] 8176.62 8573.170 9690.50 17914.66 34061.01 [4,] 16220.42 17785.385 19499.12 28923.46 62891.20 [5,] 28586.52 30505.410 30821.33 46634.38 104660.67 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 2939.748 4196.495 6937.102 3294.963 5918.606 3728.409 3797.433 [2,] 6713.532 7432.665 12978.058 9444.577 9299.634 10721.091 12445.007 [,8] [,9] [,10] [,11] [,12] [1,] 2291.543 2009.737 2495.014 8928.539 14017.65 [2,] 14061.697 15136.603 16885.986 26900.781 54104.37 $out [1] 21826.84 15997.79 26155.15 $group [1] 3 5 7 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(1664.81, 3608.415, 4826.64, 6768.065, 10243.24, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/574671244216580.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] [,8] [1,] 1664.81 2499.810 4717.02 5814.580 4826.64 7615.03 10243.24 NA [2,] 3194.00 4878.685 6097.53 6797.260 8771.77 10718.51 16677.56 NA [3,] 3733.85 5665.765 7763.12 8036.155 10743.01 15645.59 23991.00 NA [4,] 5722.65 6702.130 9824.04 10448.975 13079.56 24662.36 30663.37 NA [5,] 7600.60 7225.140 15151.84 12421.250 13606.89 36024.80 46634.38 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 2580.516 4834.079 6063.434 6370.583 8778.205 9285.704 17611.97 NA [2,] 4887.184 6497.451 9462.806 9701.727 12707.825 22005.476 30370.02 NA $out [1] 19756.21 12600.08 28541.72 34061.01 17914.66 30114.41 21822.11 [8] 45060.69 80721.71 104660.67 $group [1] 1 2 2 3 4 4 5 5 6 7 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(1664.81, 3194, 3733.85, 5722.65, 7600.6, 2499.81, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/freestat/rcomp/tmp/660ea1244216580.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,] 2937.655 8578.43 2709.525 [2,] 4363.967 13027.44 4108.760 [3,] 6815.826 20395.83 6547.505 [4,] 9899.623 25251.89 11519.835 [5,] 13909.592 39033.78 15047.495 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 4290.976 14820.18 3167.264 [2,] 9340.677 25971.49 9927.746 $out [1] 31552.80 84904.46 33563.14 $group [1] 1 2 3 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(2937.65464975848, 4363.96659986192, 6815.8263548304, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/125zs1244216580.ps tmp/125zs1244216580.png") > system("convert tmp/2lip61244216580.ps tmp/2lip61244216580.png") > system("convert tmp/37fa01244216580.ps tmp/37fa01244216580.png") > system("convert tmp/4l5p71244216580.ps tmp/4l5p71244216580.png") > system("convert tmp/574671244216580.ps tmp/574671244216580.png") > system("convert tmp/660ea1244216580.ps tmp/660ea1244216580.png") > > > proc.time() user system elapsed 1.632 1.244 1.930