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(2834,4683,4120,3849,8435,12854,15883,10520,12562,5060,4520,2150,2905,4820,3950,4053,8700,13520,15400,11100,11950,4900,4633,2300,2945,3960,3900,3767,8820,11980,14085,11600,9814,4930,4360,2640,3050,5485,4366,4790,10100,14830,17930,13580,12490,6400,4980,4930,5856,5120,5100,5623,12035,19846,17030,15860,14890,8053,6080,5987,5682,4980,5450,6035,13240,18400,17689,16490,14062,9556,7555,4328) > 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,] 2834 2905 2945 3050 5856 5682 NA [2,] 4683 4820 3960 5485 5120 4980 NA [3,] 4120 3950 3900 4366 5100 5450 NA [4,] 3849 4053 3767 4790 5623 6035 NA [5,] 8435 8700 8820 10100 12035 13240 NA [6,] 12854 13520 11980 14830 19846 18400 NA [7,] 15883 15400 14085 17930 17030 17689 NA [8,] 10520 11100 11600 13580 15860 16490 NA [9,] 12562 11950 9814 12490 14890 14062 NA [10,] 5060 4900 4930 6400 8053 9556 NA [11,] 4520 4633 4360 4980 6080 7555 NA [12,] 2150 2300 2640 4930 5987 4328 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/18xzc1291719280.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/html/rcomp/tmp/28xzc1291719280.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/html/rcomp/tmp/38xzc1291719280.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/html/rcomp/tmp/406yx1291719280.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] [,11] [,12] [1,] 2834.0 4683 3900 3767.0 8435 11980 14085.0 10520 9814 4900 4360.0 2150 [2,] 2905.0 4683 3950 3849.0 8700 12854 15400.0 11100 11950 4930 4520.0 2300 [3,] 2997.5 4900 4243 4421.5 9460 14175 16456.5 12590 12526 5730 4806.5 3484 [4,] 5682.0 5120 5100 5623.0 12035 18400 17689.0 15860 14062 8053 6080.0 4930 [5,] 5856.0 5485 5450 6035.0 13240 19846 17930.0 16490 14890 9556 7555.0 5987 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 1206.245 4618.121 3501.213 3277.213 7308.817 10597.65 14980.02 9519.646 [2,] 4788.755 5181.879 4984.787 5565.787 11611.183 17752.35 17932.98 15660.354 [,9] [,10] [,11] [,12] [1,] 11163.69 3715.564 3800.250 1787.565 [2,] 13888.31 7744.436 5812.750 5180.435 $out [1] 3960 $group [1] 2 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(2834, 2905, 2997.5, 5682, 5856, 4683, 4683, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/506yx1291719280.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,] 2150.0 2300.0 2640.0 3050.0 5100.0 4328.0 NA [2,] 3984.5 4001.5 3833.5 4860.0 5739.5 5566.0 NA [3,] 4871.5 4860.0 4645.0 5942.5 7066.5 8555.5 NA [4,] 11541.0 11525.0 10707.0 13035.0 15375.0 15276.0 NA [5,] 15883.0 15400.0 14085.0 17930.0 19846.0 18400.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 1424.930 1428.481 1509.951 2213.828 2671.684 4126.704 NA [2,] 8318.070 8291.519 7780.049 9671.172 11461.316 12984.296 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(2150, 3984.5, 4871.5, 11541, 15883, 2300, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/606yx1291719280.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,] 512.5355 1525.0 367.75 [2,] 1103.4899 2645.0 1385.75 [3,] 1536.0442 3841.0 2056.25 [4,] 1968.1605 4940.5 2749.25 [5,] 3182.7006 7866.0 4487.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 1141.662 2794.007 1434.348 [2,] 1930.426 4887.993 2678.152 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" > dev.off() null device 1 > > try(system("convert tmp/18xzc1291719280.ps tmp/18xzc1291719280.png",intern=TRUE)) character(0) > try(system("convert tmp/28xzc1291719280.ps tmp/28xzc1291719280.png",intern=TRUE)) character(0) > try(system("convert tmp/38xzc1291719280.ps tmp/38xzc1291719280.png",intern=TRUE)) character(0) > try(system("convert tmp/406yx1291719280.ps tmp/406yx1291719280.png",intern=TRUE)) character(0) > try(system("convert tmp/506yx1291719280.ps tmp/506yx1291719280.png",intern=TRUE)) character(0) > try(system("convert tmp/606yx1291719280.ps tmp/606yx1291719280.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.070 0.852 18.651