R version 2.7.1 (2008-06-23) 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(100.58,118.48,79.58,81.97,127.13,120.76,120.26,74.9,67.59,87.73,102.87,144.94,110.48,96.34,100.43,90.88,128.28,101.21,73.76,73.64,66.4,57.34,113.59,123.53,102.87,102.99,95.8,98.43,102.65,129.55,100.37,101.93,101.94,93.87,100.91,92.64,101.67,88.67,129.86,98.07,166.45,176.52,82.07,92.18,95.02,84.69,103.01,107.9,204.13,101.99,119.23,95.65,160.95,111.06,150.41,94.79,160.34,104.08,101.07,111.5,136.9,141.71,153.98,134.27,124.71,72.89,101.2,73.28,174.05,111.9,97.06,105.23,109.13,84.04,118.82,90.84,144.28,110.16,86.09,59.87,108.97,94.93,87.36,143.52,108.7,121.13,210.25,110.2,161.46,99.41,132.72,174.29,69.93,83.43,127.53,187.58) > 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,] 100.58 110.48 102.87 101.67 204.13 136.90 109.13 108.70 NA [2,] 118.48 96.34 102.99 88.67 101.99 141.71 84.04 121.13 NA [3,] 79.58 100.43 95.80 129.86 119.23 153.98 118.82 210.25 NA [4,] 81.97 90.88 98.43 98.07 95.65 134.27 90.84 110.20 NA [5,] 127.13 128.28 102.65 166.45 160.95 124.71 144.28 161.46 NA [6,] 120.76 101.21 129.55 176.52 111.06 72.89 110.16 99.41 NA [7,] 120.26 73.76 100.37 82.07 150.41 101.20 86.09 132.72 NA [8,] 74.90 73.64 101.93 92.18 94.79 73.28 59.87 174.29 NA [9,] 67.59 66.40 101.94 95.02 160.34 174.05 108.97 69.93 NA [10,] 87.73 57.34 93.87 84.69 104.08 111.90 94.93 83.43 NA [11,] 102.87 113.59 100.91 103.01 101.07 97.06 87.36 127.53 NA [12,] 144.94 123.53 92.64 107.90 111.50 105.23 143.52 187.58 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/189mi1219064082.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/2g65b1219064082.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/3ttpq1219064082.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/45hqa1219064082.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] [1,] 100.580 84.040 79.580 81.970 102.650 72.890 73.760 59.87 66.400 [2,] 102.270 92.505 98.115 90.860 125.920 100.310 84.080 73.46 68.760 [3,] 108.915 102.490 119.025 96.860 136.280 110.610 100.785 83.54 98.480 [4,] 123.690 119.805 141.920 104.315 161.205 125.155 126.490 98.36 134.655 [5,] 136.900 141.710 153.980 110.200 166.450 129.550 150.410 101.93 174.050 [,10] [,11] [,12] [1,] 83.430 87.360 92.640 [2,] 84.060 98.985 106.565 [3,] 90.800 101.970 117.515 [4,] 99.505 108.300 144.230 [5,] 111.900 113.590 187.580 $n [1] 8 8 8 8 8 8 8 8 8 8 8 8 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 96.94948 87.23983 94.5549 89.34384 116.5693 96.73123 77.09417 [2,] 120.88052 117.74017 143.4951 104.37616 155.9907 124.48877 124.47583 [,8] [,9] [,10] [,11] [,12] [1,] 69.6305 61.67011 82.1722 96.7665 96.47479 [2,] 97.4495 135.28989 99.4278 107.1735 138.55521 $out [1] 204.13 210.25 134.27 176.52 174.29 57.34 127.53 $group [1] 1 3 4 6 8 10 11 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(100.58, 102.27, 108.915, 123.69, 136.9, 84.04, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/59xph1219064082.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] [,9] [1,] 67.590 57.340 92.640 82.070 94.790 72.890 59.870 69.930 NA [2,] 80.775 73.700 97.115 90.425 101.530 99.130 86.725 104.055 NA [3,] 101.725 98.385 101.420 99.870 111.280 118.305 101.950 124.330 NA [4,] 120.510 112.035 102.760 118.880 155.375 139.305 114.490 167.875 NA [5,] 144.940 128.280 102.990 129.860 204.130 174.050 144.280 210.250 NA $n [1] 12 12 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 83.6016 80.90015 98.84528 86.89148 86.72093 99.98091 89.2862 [2,] 119.8484 115.86985 103.99472 112.84852 135.83907 136.62909 114.6138 [,8] [,9] [1,] 95.22127 NA [2,] 153.43873 NA $out [1] 129.55 166.45 176.52 $group [1] 3 4 4 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" NA Warning message: In bxp(list(stats = c(67.59, 80.775, 101.725, 120.51, 144.94, 57.34, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/www/html/rcomp/tmp/6ig1j1219064082.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.89734 40.170 5.70750 [2,] 17.70085 56.115 13.67875 [3,] 28.26943 85.795 23.87250 [4,] 35.41618 105.640 36.63000 [5,] 41.52842 130.670 52.46750 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] 20.18935 63.20632 13.40428 [2,] 36.34951 108.38368 34.34072 $out numeric(0) $group numeric(0) $names [1] "S.D." "Range" "IQR" Warning message: In bxp(list(stats = c(11.8973418399718, 17.7008497789307, 28.2694272067760, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > system("convert tmp/189mi1219064082.ps tmp/189mi1219064082.png") > system("convert tmp/2g65b1219064082.ps tmp/2g65b1219064082.png") > system("convert tmp/3ttpq1219064082.ps tmp/3ttpq1219064082.png") > system("convert tmp/45hqa1219064082.ps tmp/45hqa1219064082.png") > system("convert tmp/59xph1219064082.ps tmp/59xph1219064082.png") > system("convert tmp/6ig1j1219064082.ps tmp/6ig1j1219064082.png") > > > proc.time() user system elapsed 1.268 0.801 1.516