R version 2.12.2 (2011-02-25) Copyright (C) 2011 The R Foundation for Statistical Computing ISBN 3-900051-07-0 Platform: x86_64-redhat-linux-gnu (64-bit) 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(12544,12264,13783,11214,11453,10883,10381,10348,10024,10805,10796,11907,12261,11377,12689,11474,10992,10764,12164,10409,10398,10349,10865,11630,12221,10884,12019,11021,10799,10423,10484,10450,9906,11049,11281,12485,12849,11380,12079,11366,11328,10444,10854,10434,10137,10992,10906,12367,14371,11695,11546,10922,10670,10254,10573,10239,10253,11176,10719,11817,12503,11510,12012,10941,11252,10662,11114,10415,10626,11411,10936,12513) > 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,] 12544 12261 12221 12849 14371 12503 NA [2,] 12264 11377 10884 11380 11695 11510 NA [3,] 13783 12689 12019 12079 11546 12012 NA [4,] 11214 11474 11021 11366 10922 10941 NA [5,] 11453 10992 10799 11328 10670 11252 NA [6,] 10883 10764 10423 10444 10254 10662 NA [7,] 10381 12164 10484 10854 10573 11114 NA [8,] 10348 10409 10450 10434 10239 10415 NA [9,] 10024 10398 9906 10137 10253 10626 NA [10,] 10805 10349 11049 10992 11176 11411 NA [11,] 10796 10865 11281 10906 10719 10936 NA [12,] 11907 11630 12485 12367 11817 12513 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/wessaorg/rcomp/tmp/1ai3a1305024753.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2a9h01305024753.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3foea1305024753.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4wzo41305024753.ps",horizontal=F,onefile=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,] 12221.0 11377 11546 10922.0 10670 10254 10381.0 10239 9906 10349.0 [2,] 12261.0 11377 12012 10941.0 10799 10423 10484.0 10348 10024 10805.0 [3,] 12523.5 11445 12049 11117.5 11122 10553 10713.5 10412 10195 11020.5 [4,] 12849.0 11695 12689 11366.0 11328 10764 11114.0 10434 10398 11176.0 [5,] 12849.0 11695 12689 11474.0 11453 10883 11114.0 10450 10626 11411.0 [,11] [,12] [1,] 10719.0 11630 [2,] 10796.0 11817 [3,] 10885.5 12137 [4,] 10936.0 12485 [5,] 10936.0 12513 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 12144.22 11239.88 11612.31 10843.36 10780.78 10333.04 10307.13 10356.53 [2,] 12902.78 11650.12 12485.69 11391.64 11463.22 10772.96 11119.87 10467.47 [,9] [,10] [,11] [,12] [1,] 9953.758 10781.19 10795.20 11706.12 [2,] 10436.242 11259.81 10975.80 12567.88 $out [1] 14371 12264 10884 13783 12164 11281 $group [1] 1 2 2 3 7 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(12221, 12261, 12523.5, 12849, 12849, 11377, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/5u9r21305024753.ps",horizontal=F,onefile=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,] 10024.0 10349.0 9906.0 10137.0 10239.0 10415 NA [2,] 10588.5 10586.5 10467.0 10649.0 10413.5 10799 NA [3,] 11048.5 11184.5 10952.5 11160.0 10820.5 11183 NA [4,] 12085.5 11897.0 11650.0 11729.5 11620.5 11761 NA [5,] 13783.0 12689.0 12485.0 12849.0 11817.0 12513 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 10365.71 10586.77 10412.93 10667.18 10269.98 10744.23 NA [2,] 11731.29 11782.23 11492.07 11652.82 11371.02 11621.77 NA $out [1] 14371 $group [1] 5 $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(10024, 10588.5, 11048.5, 12085.5, 13783, 10349, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/wessaorg/rcomp/tmp/60qc51305024753.ps",horizontal=F,onefile=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,] 78.4213 211.0 66.000 [2,] 234.3294 595.5 281.750 [3,] 336.8670 833.0 338.625 [4,] 557.7209 1581.5 492.250 [5,] 806.1657 2237.0 616.000 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 189.3659 383.2788 242.6145 [2,] 484.3681 1282.7212 434.6355 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(78.4212981274857, 234.329352945183, 336.867004123272, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1ai3a1305024753.ps tmp/1ai3a1305024753.png",intern=TRUE)) character(0) > try(system("convert tmp/2a9h01305024753.ps tmp/2a9h01305024753.png",intern=TRUE)) character(0) > try(system("convert tmp/3foea1305024753.ps tmp/3foea1305024753.png",intern=TRUE)) character(0) > try(system("convert tmp/4wzo41305024753.ps tmp/4wzo41305024753.png",intern=TRUE)) character(0) > try(system("convert tmp/5u9r21305024753.ps tmp/5u9r21305024753.png",intern=TRUE)) character(0) > try(system("convert tmp/60qc51305024753.ps tmp/60qc51305024753.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 3.940 0.430 4.483