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. 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(42.3,50.8,54.1,38.2,48.4,61.1,54.1,61.4,64.3,57.4,71.7,55.3,55.1,66.8,59.4,64.9,59.2,77.4,75.8,38.3,54,61.8,61.3,104.3,39.7,62.6,50.2,90.9,56.2,50.2,52.8,45.6,69,81.9,73.9,54.9,55.4,64.6,49.6,55.8,44.6,61.5,40.5,48.3,50.9,65.3,56.5,53.2,56.9,79.5,94,68.4,65.9,85.5,77.5,114.8,87.4,107.5,151.7,94.4,67.5,95.2,96.2,70.6,80.1,83.4,115.4,61.5,80.6,94.3,82.6,107.7,79.1,102.8,125.2,106.4,62.3,107.4,67.9,88,76.5,130.5,100.9,85.6) > 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,] 42.3 55.1 39.7 55.4 56.9 67.5 79.1 NA [2,] 50.8 66.8 62.6 64.6 79.5 95.2 102.8 NA [3,] 54.1 59.4 50.2 49.6 94.0 96.2 125.2 NA [4,] 38.2 64.9 90.9 55.8 68.4 70.6 106.4 NA [5,] 48.4 59.2 56.2 44.6 65.9 80.1 62.3 NA [6,] 61.1 77.4 50.2 61.5 85.5 83.4 107.4 NA [7,] 54.1 75.8 52.8 40.5 77.5 115.4 67.9 NA [8,] 61.4 38.3 45.6 48.3 114.8 61.5 88.0 NA [9,] 64.3 54.0 69.0 50.9 87.4 80.6 76.5 NA [10,] 57.4 61.8 81.9 65.3 107.5 94.3 130.5 NA [11,] 71.7 61.3 73.9 56.5 151.7 82.6 100.9 NA [12,] 55.3 104.3 54.9 53.2 94.4 107.7 85.6 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/1kk6y1211751531.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/20vtu1211751531.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/33juj1211751531.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/4sbvs1211751531.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] [,11] [1,] 39.7 50.80 49.60 38.20 44.6 50.20 40.50 38.30 50.90 57.40 56.50 [2,] 48.7 63.60 52.15 60.35 52.3 61.30 53.45 46.95 59.15 63.55 66.50 [3,] 55.4 66.80 59.40 68.40 59.2 77.40 67.90 61.40 69.00 81.90 73.90 [4,] 62.2 87.35 95.10 80.75 64.1 84.45 76.65 74.75 78.55 100.90 91.75 [5,] 79.1 102.80 125.20 106.40 80.1 107.40 77.50 114.80 87.40 130.50 100.90 [,12] [1,] 53.20 [2,] 55.10 [3,] 85.60 [4,] 99.35 [5,] 107.70 $n [1] 7 7 7 7 7 7 7 7 7 7 7 7 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 47.33802 52.61688 33.75095 56.21745 52.15323 63.57519 54.04533 44.79829 [2,] 63.46198 80.98312 85.04905 80.58255 66.24677 91.22481 81.75467 78.00171 [,9] [,10] [,11] [,12] [1,] 57.41463 59.59518 58.82111 59.17461 [2,] 80.58537 104.20482 88.97889 112.02539 $out [1] 115.4 151.7 $group [1] 7 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(39.7, 48.7, 55.4, 62.2, 79.1, 50.8, 63.6, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5ggsv1211751531.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,] 38.20 38.30 39.70 40.50 56.90 61.50 62.30 NA [2,] 49.60 57.15 50.20 48.95 72.95 75.35 77.80 NA [3,] 54.70 61.55 55.55 54.30 86.45 83.00 94.45 NA [4,] 61.25 71.30 71.45 59.00 100.95 95.70 106.90 NA [5,] 71.70 77.40 90.90 65.30 114.80 115.40 130.50 NA $n [1] 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 49.38636 55.09609 45.85773 49.71613 73.67901 73.71823 81.1773 NA [2,] 60.01364 68.00391 65.24227 58.88387 99.22099 92.28177 107.7227 NA $out [1] 104.3 151.7 $group [1] 2 5 $names [1] "1" "2" "3" "4" "5" "6" "7" NA Warning message: In bxp(list(stats = c(38.2, 49.6, 54.7, 61.25, 71.7, 38.3, 57.15, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ijrc1211751531.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,] 11.75694 35.50 11.800 [2,] 16.21529 45.70 19.900 [3,] 23.34502 62.70 23.475 [4,] 27.01831 75.25 32.575 [5,] 32.60789 95.20 44.250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 18.41769 49.22205 17.69385 [2,] 28.27235 76.17795 29.25615 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(11.7569432373264, 16.2152940263060, 23.3450247175464, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/1kk6y1211751531.ps tmp/1kk6y1211751531.png") > system("convert tmp/20vtu1211751531.ps tmp/20vtu1211751531.png") > system("convert tmp/33juj1211751531.ps tmp/33juj1211751531.png") > system("convert tmp/4sbvs1211751531.ps tmp/4sbvs1211751531.png") > system("convert tmp/5ggsv1211751531.ps tmp/5ggsv1211751531.png") > system("convert tmp/6ijrc1211751531.ps tmp/6ijrc1211751531.png") > > > proc.time() user system elapsed 1.318 0.844 5.900