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(9.733,9.259,9.864,9.215,10.103,9.380,9.896,10.117,9.451,9.700,9.081,9.084,9.743,8.587,9.731,9.563,9.998,9.437,10.038,9.918,9.252,9.737,9.035,9.133,9.487,8.700,9.627,8.947,9.283,8.829,9.947,9.628,9.318,9.605,8.640,9.214,9.567,8.547,9.185,9.470,9.123,9.278,10.170,9.434,9.655,9.429,8.739,9.552,9.687,9.019,9.672,9.206,9.069,9.788,10.312,10.105,9.863,9.656,9.295,9.946,9.701,9.049,10.190,9.706,9.765,9.893,9.994,10.433,10.073,10.112,9.266,9.820,10.097,9.115,10.411,9.678,10.408,10.153,10.368,10.581,10.597,10.680,9.738,9.556) > 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,] 9.733 9.743 9.487 9.567 9.687 9.701 10.097 NA [2,] 9.259 8.587 8.700 8.547 9.019 9.049 9.115 NA [3,] 9.864 9.731 9.627 9.185 9.672 10.190 10.411 NA [4,] 9.215 9.563 8.947 9.470 9.206 9.706 9.678 NA [5,] 10.103 9.998 9.283 9.123 9.069 9.765 10.408 NA [6,] 9.380 9.437 8.829 9.278 9.788 9.893 10.153 NA [7,] 9.896 10.038 9.947 10.170 10.312 9.994 10.368 NA [8,] 10.117 9.918 9.628 9.434 10.105 10.433 10.581 NA [9,] 9.451 9.252 9.318 9.655 9.863 10.073 10.597 NA [10,] 9.700 9.737 9.605 9.429 9.656 10.112 10.680 NA [11,] 9.081 9.035 8.640 8.739 9.295 9.266 9.738 NA [12,] 9.084 9.133 9.214 9.552 9.946 9.820 9.556 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/19hzb1243510249.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/2vfyc1243510249.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/3ubq61243510249.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/48vam1243510249.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,] 9.487 8.5470 9.1850 8.9470 9.0690 8.8290 9.8960 9.434 9.2520 9.4290 [2,] 9.627 8.6435 9.6495 9.2105 9.2030 9.3290 9.9705 9.773 9.3845 9.6305 [3,] 9.701 9.0190 9.7310 9.4700 9.7650 9.4370 10.0380 10.105 9.6550 9.7000 [4,] 9.738 9.0820 10.0270 9.6205 10.0505 9.8405 10.2410 10.275 9.9680 9.9245 [5,] 9.743 9.2590 10.4110 9.7060 10.4080 10.1530 10.3680 10.581 10.5970 10.1120 [,11] [,12] [1,] 8.6400 9.0840 [2,] 8.8870 9.1735 [3,] 9.0810 9.5520 [4,] 9.2805 9.6880 [5,] 9.7380 9.9460 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9.634713 8.757135 9.505563 9.225155 9.258887 9.13154 9.876462 9.805214 [2,] 9.767287 9.280865 9.956437 9.714845 10.271113 9.74246 10.199538 10.404786 [,9] [,10] [,11] [,12] [1,] 9.306543 9.524428 8.846008 9.244749 [2,] 10.003457 9.875572 9.315992 9.859251 $out [1] 10.097 10.680 $group [1] 1 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(9.487, 9.627, 9.701, 9.738, 9.743, 8.547, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/56e1l1243510249.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,] 9.0810 8.5870 8.6400 8.5470 9.0190 9.2660 9.1150 NA [2,] 9.2370 9.1925 8.8880 9.1540 9.2505 9.7035 9.7080 NA [3,] 9.5755 9.6470 9.3005 9.4315 9.6795 9.8565 10.2605 NA [4,] 9.8800 9.8305 9.6160 9.5595 9.9045 10.0925 10.4960 NA [5,] 10.1170 10.0380 9.9470 9.6550 10.3120 10.4330 10.6800 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 9.282223 9.356004 8.968454 9.246549 9.381206 9.679074 9.901088 NA [2,] 9.868777 9.937996 9.632546 9.616451 9.977794 10.033926 10.619912 NA $out [1] 10.170 9.049 $group [1] 4 6 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(9.081, 9.237, 9.5755, 9.88, 10.117, 8.587, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6nn9v1243510249.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,] 0.1833411 0.4720 0.11100 [2,] 0.2818662 0.7355 0.33575 [3,] 0.3843833 1.1225 0.42425 [4,] 0.4319088 1.2875 0.51300 [5,] 0.5257413 1.3450 0.58350 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 0.3159479 0.8707291 0.3434051 [2,] 0.4528188 1.3742709 0.5050949 $out [1] 0.8475 $group [1] 3 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(0.183341082087328, 0.281866213108328, 0.384383326227625, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/19hzb1243510249.ps tmp/19hzb1243510249.png") > system("convert tmp/2vfyc1243510249.ps tmp/2vfyc1243510249.png") > system("convert tmp/3ubq61243510249.ps tmp/3ubq61243510249.png") > system("convert tmp/48vam1243510249.ps tmp/48vam1243510249.png") > system("convert tmp/56e1l1243510249.ps tmp/56e1l1243510249.png") > system("convert tmp/6nn9v1243510249.ps tmp/6nn9v1243510249.png") > > > proc.time() user system elapsed 1.095 0.850 1.578