R version 2.8.0 (2008-10-20)
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(113.9000,112.0000,113.8500,113.0800,111.7200,110.6900,113.5300,113.9900,112.7400,112.1500,115.8200,118.3800,118.8100,123.8500,117.9600,120.1600,118.7400,119.8400,124.8100,121.3300,120.2000,118.3200,129.5800,130.2000,127.1900,133.1000,129.1200,123.2800,123.3600,124.1300,126.9700,127.1400,123.7000,123.6700,130.1900,134.0100,124.9600,129.9600,128.3200,132.3800,126.2500,128.9100,131.4200,129.4400,126.8600,126.7100,131.6300,132.7800,126.6100,132.8400,123.1400,128.1300,125.4900,126.4800,130.8600,127.3200,126.5600,126.6400,129.2600,126.4700,135.4000,135.5000,132.2200,122.6200,125.1600,128.5000,133.8600,128.8700,125.0700,125.2500,132.1600,130.2400)
> 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,] 113.90 118.81 127.19 124.96 126.61 135.40   NA
 [2,] 112.00 123.85 133.10 129.96 132.84 135.50   NA
 [3,] 113.85 117.96 129.12 128.32 123.14 132.22   NA
 [4,] 113.08 120.16 123.28 132.38 128.13 122.62   NA
 [5,] 111.72 118.74 123.36 126.25 125.49 125.16   NA
 [6,] 110.69 119.84 124.13 128.91 126.48 128.50   NA
 [7,] 113.53 124.81 126.97 131.42 130.86 133.86   NA
 [8,] 113.99 121.33 127.14 129.44 127.32 128.87   NA
 [9,] 112.74 120.20 123.70 126.86 126.56 125.07   NA
[10,] 112.15 118.32 123.67 126.71 126.64 125.25   NA
[11,] 115.82 129.58 130.19 131.63 129.26 132.16   NA
[12,] 118.38 130.20 134.01 132.78 126.47 130.24   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/15q4g1229883759.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/2iyh41229883759.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/3aqnj1229883759.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/4kh5x1229883759.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]  [,10]
[1,] 113.900 112.00 113.85 113.08 111.72 110.690 124.810 113.99 112.740 112.15
[2,] 118.810 123.85 117.96 120.16 118.74 119.840 124.810 121.33 120.200 118.32
[3,] 125.785 131.40 125.73 122.95 124.26 125.305 128.915 127.23 124.385 124.46
[4,] 127.190 133.10 129.12 128.13 125.49 128.500 131.420 128.87 126.560 126.64
[5,] 135.400 135.50 132.22 132.38 126.25 128.910 133.860 129.44 126.860 126.71
       [,11]  [,12]
[1,] 129.260 118.38
[2,] 129.260 126.47
[3,] 129.885 130.22
[4,] 131.630 132.78
[5,] 132.160 134.01

$n
 [1] 6 6 6 6 6 6 6 6 6 6 6 6

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 120.3796 125.4335 118.5314 117.8091 119.9060 119.7190 124.6513 122.3665
[2,] 131.1904 137.3665 132.9286 128.0909 128.6140 130.8910 133.1787 132.0935
         [,9]    [,10]    [,11]    [,12]
[1,] 120.2826 119.0933 128.3563 126.1498
[2,] 128.4874 129.8267 131.4137 134.2902

$out
[1] 113.53 115.82

$group
[1]  7 11

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12"

Warning message:
In bxp(list(stats = c(113.9, 118.81, 125.785, 127.19, 135.4, 112,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5tfvg1229883759.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]
[1,] 110.690 117.960 123.280 124.960 125.490 122.620   NA
[2,] 112.075 118.775 123.685 126.785 126.475 125.205   NA
[3,] 113.305 120.180 127.055 129.175 126.625 129.555   NA
[4,] 113.945 124.330 129.655 131.525 128.695 133.040   NA
[5,] 115.820 130.200 134.010 132.780 130.860 135.500   NA

$n
[1] 12 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 112.4521 117.6463 124.3320 127.0131 125.6124 125.9814   NA
[2,] 114.1579 122.7137 129.7780 131.3369 127.6376 133.1286   NA

$out
[1] 118.38 132.84 123.14

$group
[1] 1 5 5

$names
[1] "1" "2" "3" "4" "5" "6" NA 

Warning message:
In bxp(list(stats = c(110.69, 112.075, 113.305, 113.945, 115.82,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6yrfi1229883759.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,] 5.368107 14.120 4.74250
[2,] 5.730991 15.005 5.31250
[3,] 6.383770 17.280 6.03625
[4,] 7.212496 19.815 6.89000
[5,] 8.748211 23.500 7.65750

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 5.708046 15.08613 5.316742
[2,] 7.059495 19.47387 6.755758

$out
[1] 9.665 1.930

$group
[1] 3 3

$names
[1] "S.D."  "Range" "IQR"  

Warning message:
In bxp(list(stats = c(5.36810736355624, 5.73099083434044, 6.38377040768932,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/15q4g1229883759.ps tmp/15q4g1229883759.png")
> system("convert tmp/2iyh41229883759.ps tmp/2iyh41229883759.png")
> system("convert tmp/3aqnj1229883759.ps tmp/3aqnj1229883759.png")
> system("convert tmp/4kh5x1229883759.ps tmp/4kh5x1229883759.png")
> system("convert tmp/5tfvg1229883759.ps tmp/5tfvg1229883759.png")
> system("convert tmp/6yrfi1229883759.ps tmp/6yrfi1229883759.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.061   0.821   1.329