R version 2.10.1 (2009-12-14) 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(101397,97994,94044,91159,87239,89235,118647,125620,125154,117529,109459,108483,107137,104699,100804,96066,91971,93228,120144,127233,127166,118194,109940,106683,102834,99882,96666,92540,88744,89321,115870,122401,122030,113802,105791,103076,98658,96945,92497,90687,88796,90015,113228,118711,117460,106556,97347,92657,93118,89037,83570,81693,75956,73993,97088,102394,96549,89727,82336,82653,82303,79596,74472,73562,66618,69029,89899,93774,90305,83799,80320,82497,84420,84646,84186,83269,77793,81145,101691,107357,104253,95963,91432,94324,93855,92183,87600,83641,78195,79604,100846,105293,102518,93132,87479,85476) > 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] 96 > (np <- floor(n / par1)) [1] 8 > 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] 8 8 8 8 8 8 8 8 8 8 8 8 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [1,] 101397 107137 102834 98658 93118 82303 84420 93855 NA [2,] 97994 104699 99882 96945 89037 79596 84646 92183 NA [3,] 94044 100804 96666 92497 83570 74472 84186 87600 NA [4,] 91159 96066 92540 90687 81693 73562 83269 83641 NA [5,] 87239 91971 88744 88796 75956 66618 77793 78195 NA [6,] 89235 93228 89321 90015 73993 69029 81145 79604 NA [7,] 118647 120144 115870 113228 97088 89899 101691 100846 NA [8,] 125620 127233 122401 118711 102394 93774 107357 105293 NA [9,] 125154 127166 122030 117460 96549 90305 104253 102518 NA [10,] 117529 118194 113802 106556 89727 83799 95963 93132 NA [11,] 109459 109940 105791 97347 82336 80320 91432 87479 NA [12,] 108483 106683 103076 92657 82653 82497 94324 85476 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/yougetitorg/rcomp/tmp/1x1w91304771735.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/yougetitorg/rcomp/tmp/2d71r1304771735.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/yougetitorg/rcomp/tmp/3hkuh1304771735.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/yougetitorg/rcomp/tmp/4ncai1304771735.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] [1,] 82303.0 79596.0 74472.0 73562.0 66618.0 69029.0 89899.0 93774.0 [2,] 88769.0 86841.5 83878.0 82481.0 76874.5 76798.5 98967.0 103843.5 [3,] 96256.5 94564.0 90048.5 87164.0 82717.0 85190.0 107459.5 113034.0 [4,] 102115.5 98938.0 95355.0 91849.5 88770.0 89668.0 117258.5 124010.5 [5,] 107137.0 104699.0 100804.0 96066.0 91971.0 93228.0 120144.0 127233.0 [,9] [,10] [,11] [,12] [1,] 90305.0 83799.0 80320.0 82497.0 [2,] 99533.5 91429.5 84907.5 84064.5 [3,] 110856.5 101259.5 94389.5 93490.5 [4,] 123592.0 115665.5 107625.0 104879.5 [5,] 127166.0 118194.0 109940.0 108483.0 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 88800.95 87806.72 83637.28 81930.62 76072 78000.91 97241.6 101768.4 [2,] 103712.05 101321.28 96459.72 92397.38 89362 92379.09 117677.4 124299.6 [,9] [,10] [,11] [,12] [1,] 97417.08 87720.92 81699.18 81862.94 [2,] 124295.92 114798.08 107079.82 105118.06 $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(82303, 88769, 96256.5, 102115.5, 107137, 79596, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/5j9m01304771735.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] [,8] [,9] [1,] 87239.0 91971 88744 88796 73993.0 66618.0 77793.0 78195.0 NA [2,] 92601.5 98435 94603 91592 82014.5 74017.0 83727.5 84558.5 NA [3,] 104940.0 106910 102955 97146 86303.5 81311.5 88039.0 89891.5 NA [4,] 118088.0 119169 114836 109892 94833.5 86849.0 98827.0 97350.5 NA [5,] 125620.0 127233 122401 118711 102394.0 93774.0 107357.0 105293.0 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] 93315.44 97453.08 93726.6 88799.25 80456.67 75458.74 81152.02 84056.98 [2,] 116564.56 116366.92 112183.4 105492.75 92150.33 87164.26 94925.98 95726.02 [,9] [1,] NA [2,] NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(87239, 92601.5, 104940, 118088, 125620, 91971, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/yougetitorg/rcomp/tmp/6whwv1304771735.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,] 7331.848 22504.0 8629.25 [2,] 8554.618 24968.5 10740.12 [3,] 9689.380 26159.0 14040.50 [4,] 12109.293 31852.0 19861.12 [5,] 14034.784 36861.0 22453.00 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 8068.068 23019.39 9880.35 [2,] 11310.691 29298.61 18200.65 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(7331.84789487421, 8554.61826785611, 9689.37961698622, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1x1w91304771735.ps tmp/1x1w91304771735.png",intern=TRUE)) character(0) > try(system("convert tmp/2d71r1304771735.ps tmp/2d71r1304771735.png",intern=TRUE)) character(0) > try(system("convert tmp/3hkuh1304771735.ps tmp/3hkuh1304771735.png",intern=TRUE)) character(0) > try(system("convert tmp/4ncai1304771735.ps tmp/4ncai1304771735.png",intern=TRUE)) character(0) > try(system("convert tmp/5j9m01304771735.ps tmp/5j9m01304771735.png",intern=TRUE)) character(0) > try(system("convert tmp/6whwv1304771735.ps tmp/6whwv1304771735.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 1.330 1.000 1.836