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(100,87.14054095,112.0054296,112.312101,109.474134,104.9746116,100.4926851,104.2154743,120.1768388,112.1028355,108.1481575,116.802197,102.1699512,95.15358705,120.6707808,111.5234277,119.9669448,113.3697401,110.0717661,111.5567342,132.424212,107.900558,122.1626615,124.3992258,110.4450505,101.5874013,122.3203962,125.2582826,125.4411543,108.9902468,118.9243879,116.7242723,134.1724901,116.8530994,124.5732995,130.9914031,123.4239103,111.4536725,124.5135991,139.2589613,129.8596099,112.3460359,131.381655,133.0004776,134.3220552,144.2379719,134.1278719,150.1891559,140.722563,114.8389975,143.1973003,140.2738676,112.1248303,102.8951536,100.5090242,103.3513901,111.4134533,104.5887587,101.7840983,114.7007441,108.7426474) > 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] 61 > (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] 6 5 5 5 5 5 5 5 5 5 5 5 > arr [,1] [,2] [,3] [,4] [,5] [,6] [1,] 100.00000 102.16995 110.4451 123.4239 140.7226 108.7426 [2,] 87.14054 95.15359 101.5874 111.4537 114.8390 NA [3,] 112.00543 120.67078 122.3204 124.5136 143.1973 NA [4,] 112.31210 111.52343 125.2583 139.2590 140.2739 NA [5,] 109.47413 119.96694 125.4412 129.8596 112.1248 NA [6,] 104.97461 113.36974 108.9902 112.3460 102.8952 NA [7,] 100.49269 110.07177 118.9244 131.3817 100.5090 NA [8,] 104.21547 111.55673 116.7243 133.0005 103.3514 NA [9,] 120.17684 132.42421 134.1725 134.3221 111.4135 NA [10,] 112.10284 107.90056 116.8531 144.2380 104.5888 NA [11,] 108.14816 122.16266 124.5733 134.1279 101.7841 NA [12,] 116.80220 124.39923 130.9914 150.1892 114.7007 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/1vrba1257955876.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/2g6gy1257955876.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/3v1fy1257955876.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/4wf0e1257955876.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] [1,] 100.0000 87.14054 120.6708 111.5234 109.4741 102.8952 100.4927 103.3514 [2,] 102.1700 95.15359 120.6708 112.3121 112.1248 104.9746 100.5090 104.2155 [3,] 109.5938 101.58740 122.3204 125.2583 119.9669 108.9902 110.0718 111.5567 [4,] 123.4239 111.45367 124.5136 139.2590 125.4412 112.3460 118.9244 116.7243 [5,] 140.7226 114.83900 124.5136 140.2739 129.8596 113.3697 131.3817 133.0005 [,9] [,10] [,11] [,12] [1,] 111.4135 104.5888 101.7841 114.7007 [2,] 120.1768 107.9006 108.1482 116.8022 [3,] 132.4242 112.1028 122.1627 124.3992 [4,] 134.1725 116.8531 124.5733 130.9914 [5,] 134.3221 116.8531 134.1279 150.1892 $n [1] 6 5 5 5 5 5 5 5 5 5 5 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 95.88436 90.0698 119.6051 106.2177 110.5577 103.7816 97.05952 102.7180 [2,] 123.30334 113.1050 125.0357 144.2989 129.3762 114.1989 123.08402 120.3954 [,9] [,10] [,11] [,12] [1,] 122.5349 105.7770 110.5567 114.3732 [2,] 142.3135 118.4287 133.7686 134.4253 $out [1] 112.0054 143.1973 144.2380 $group [1] 3 3 10 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(100, 102.1699512, 109.59384895, 123.4239103, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/5k1061257955876.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,] 100.0000 95.15359 101.5874 111.4537 100.5090 108.7426 [2,] 102.3541 108.98616 113.5847 123.9688 103.1233 108.7426 [3,] 108.8111 112.46324 120.6224 132.1911 111.7691 108.7426 [4,] 112.2075 121.41672 125.3497 136.7905 127.5564 108.7426 [5,] 120.1768 132.42421 134.1725 150.1892 143.1973 108.7426 $n [1] 12 12 12 12 12 1 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] 104.3169 106.7936 115.2563 126.3430 100.6250 108.7426 [2,] 113.3053 118.1329 125.9885 138.0392 122.9133 108.7426 $out [1] 87.14054 $group [1] 1 $names [1] "1" "2" "3" "4" "5" NA Warning message: In bxp(list(stats = c(100, 102.3540797, 108.81114575, 112.20746825, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6avfz1257955876.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,] 8.626525 20.38548 3.842818 [2,] 10.852244 25.30353 10.730670 [3,] 12.559517 30.26903 14.092429 [4,] 14.103406 33.91609 16.395606 [5,] 15.832374 40.72256 18.415364 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 11.07664 26.34078 11.50861 [2,] 14.04239 34.19728 16.67624 $out [1] 4.543903 10.474587 26.946860 $group [1] 1 2 3 $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(8.62652545574986, 10.8522440633783, 12.5595172468239, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > system("convert tmp/1vrba1257955876.ps tmp/1vrba1257955876.png") > system("convert tmp/2g6gy1257955876.ps tmp/2g6gy1257955876.png") > system("convert tmp/3v1fy1257955876.ps tmp/3v1fy1257955876.png") > system("convert tmp/4wf0e1257955876.ps tmp/4wf0e1257955876.png") > system("convert tmp/5k1061257955876.ps tmp/5k1061257955876.png") > system("convert tmp/6avfz1257955876.ps tmp/6avfz1257955876.png") > > > proc.time() user system elapsed 1.085 0.837 1.601