R version 2.6.2 (2008-02-08) 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. Natural language support but running in an English locale 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(117.99,118.13,117.06,118.03,118.57,118.44,118.02,118.18,116.79,116.91,117.19,117.3,116.35,117.04,116.26,115.91,115.75,115.21,115.07,114.17,113.96,113.59,114.47,113.52,113.56,113.59,113,112.07,111.11,109.61,109.51,109.21,108.86,107.99,108.74,108.74,107.74,107.8,106.96,106.73,106.27,105.69,105.39,105.09,104.94,104.75,104.38,103.65,103.49,103.39,103.48,103.33,103.61,103.21,102.85,103.25,102.88,102.51,101.99,101.6,101.21,100.74,100.36,100.26,100.37,100.23,100.02,99.8,99.58,99.69,99.53,99.37) > 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] 72 > (np <- floor(n / par1)) [1] 6 > 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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 117.99 116.35 113.56 107.74 103.49 101.21 NA [2,] 118.13 117.04 113.59 107.80 103.39 100.74 NA [3,] 117.06 116.26 113.00 106.96 103.48 100.36 NA [4,] 118.03 115.91 112.07 106.73 103.33 100.26 NA [5,] 118.57 115.75 111.11 106.27 103.61 100.37 NA [6,] 118.44 115.21 109.61 105.69 103.21 100.23 NA [7,] 118.02 115.07 109.51 105.39 102.85 100.02 NA [8,] 118.18 114.17 109.21 105.09 103.25 99.80 NA [9,] 116.79 113.96 108.86 104.94 102.88 99.58 NA [10,] 116.91 113.59 107.99 104.75 102.51 99.69 NA [11,] 117.19 114.47 108.74 104.38 101.99 99.53 NA [12,] 117.30 113.52 108.74 103.65 101.60 99.37 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/1w56c1211119949.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/2v2yw1211119949.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/3iuvf1211119949.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/430u31211119949.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] [,8] [,9] [,10] [1,] 101.21 100.740 100.36 100.26 100.37 100.23 100.02 99.80 99.58 99.69 [2,] 103.49 103.390 103.48 103.33 103.61 103.21 102.85 103.25 102.88 102.51 [3,] 110.65 110.695 109.98 109.40 108.69 107.65 107.45 107.15 106.90 106.37 [4,] 116.35 117.040 116.26 115.91 115.75 115.21 115.07 114.17 113.96 113.59 [5,] 117.99 118.130 117.06 118.03 118.57 118.44 118.02 118.18 116.79 116.91 [,11] [,12] [1,] 99.53 99.370 [2,] 101.99 101.600 [3,] 106.56 106.195 [4,] 114.47 113.520 [5,] 117.19 117.300 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 102.3549 101.8903 101.7365 101.2855 100.8593 99.90961 99.5677 100.1062 [2,] 118.9451 119.4997 118.2235 117.5145 116.5207 115.39039 115.3323 114.1938 [,9] [,10] [,11] [,12] [1,] 99.75304 99.22304 98.51 98.50621 [2,] 114.04696 113.51696 114.61 113.88379 $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(101.21, 103.49, 110.65, 116.35, 117.99, 100.74, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/59erb1211119949.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,] 116.790 113.520 107.990 103.650 101.600 99.370 NA [2,] 117.125 114.065 108.800 104.845 102.680 99.635 NA [3,] 118.005 115.140 109.560 105.540 103.230 100.125 NA [4,] 118.155 116.085 112.535 106.845 103.435 100.365 NA [5,] 118.570 117.040 113.590 107.800 103.610 101.210 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 117.5352 114.2187 107.8564 104.6278 102.8856 99.79204 NA [2,] 118.4748 116.0613 111.2636 106.4522 103.5744 100.45796 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(116.79, 117.125, 118.005, 118.155, 118.57, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/68yop1211119949.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,] 6.941964 16.700 9.12000 [2,] 6.944672 17.215 9.63500 [3,] 7.040348 17.715 10.26375 [4,] 7.077653 18.100 10.93250 [5,] 7.249096 18.380 11.68500 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 6.979694 17.31135 9.671952 [2,] 7.101001 18.11865 10.855548 $out [1] 6.630782 6.622730 $group [1] 1 1 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(6.94196417929873, 6.94467177066579, 7.04034795500649, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1w56c1211119949.ps tmp/1w56c1211119949.png") > system("convert tmp/2v2yw1211119949.ps tmp/2v2yw1211119949.png") > system("convert tmp/3iuvf1211119949.ps tmp/3iuvf1211119949.png") > system("convert tmp/430u31211119949.ps tmp/430u31211119949.png") > system("convert tmp/59erb1211119949.ps tmp/59erb1211119949.png") > system("convert tmp/68yop1211119949.ps tmp/68yop1211119949.png") > > > proc.time() user system elapsed 2.377 1.346 2.597