R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-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.

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(11.8,11.87,11.86,12.11,11.92,12.61,13.34,13.31,13.47,13.24,13.18,13.3,13.15,13.47,13.65,13.52,14.13,14.84,15.29,15.51,15.43,15.42,15.56,15.43,15.36,15.18,15.41,15.15,15.21,15.09,15.09,15.5,15.41,15.42,15.47,15.23,15.59,15.22,15.45,15.02,15.5,15.59,15.98,15.76,15.43,15.45,15.32,15.4,15.42,15.54,15.6,15.67,15.61,16.01,16.06,16.15,15.87,15.89,15.73,15.78,16.07,16.2,16.42,16.61,16.89,17.62,17.83,17.94,18.07,17.85,17.86,17.85)
> 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,] 11.80 13.15 15.36 15.59 15.42 16.07   NA
 [2,] 11.87 13.47 15.18 15.22 15.54 16.20   NA
 [3,] 11.86 13.65 15.41 15.45 15.60 16.42   NA
 [4,] 12.11 13.52 15.15 15.02 15.67 16.61   NA
 [5,] 11.92 14.13 15.21 15.50 15.61 16.89   NA
 [6,] 12.61 14.84 15.09 15.59 16.01 17.62   NA
 [7,] 13.34 15.29 15.09 15.98 16.06 17.83   NA
 [8,] 13.31 15.51 15.50 15.76 16.15 17.94   NA
 [9,] 13.47 15.43 15.41 15.43 15.87 18.07   NA
[10,] 13.24 15.42 15.42 15.45 15.89 17.85   NA
[11,] 13.18 15.56 15.47 15.32 15.73 17.86   NA
[12,] 13.30 15.43 15.23 15.40 15.78 17.85   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/wessaorg/rcomp/tmp/1f9841326703528.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/wessaorg/rcomp/tmp/2lg7b1326703528.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/wessaorg/rcomp/tmp/3v8jh1326703528.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/wessaorg/rcomp/tmp/48ei71326703528.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]
[1,] 11.80 11.87 11.86 12.110 11.920 14.84 15.090 15.500 15.41 15.420 15.320
[2,] 13.15 13.47 13.65 13.520 14.130 14.84 15.090 15.500 15.41 15.420 15.320
[3,] 15.39 15.20 15.43 15.085 15.355 15.34 15.635 15.635 15.43 15.435 15.515
[4,] 15.59 15.54 15.60 15.670 15.610 16.01 16.060 16.150 15.87 15.890 15.730
[5,] 16.07 16.20 16.42 16.610 16.890 17.62 16.060 16.150 15.87 15.890 15.730
      [,12]
[1,] 15.230
[2,] 15.230
[3,] 15.415
[4,] 15.780
[5,] 15.780

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 13.81612 13.86478 14.17219 13.69818 14.40035 14.58531 15.00932 15.21573
[2,] 16.96388 16.53522 16.68781 16.47182 16.30965 16.09469 16.26068 16.05427
         [,9]    [,10]    [,11]    [,12]
[1,] 15.13329 15.13183 15.25054 15.06023
[2,] 15.72671 15.73817 15.77946 15.76977

$out
 [1] 12.61 13.34 17.83 13.31 17.94 13.47 18.07 13.24 17.85 13.18 17.86 13.30
[13] 17.85

$group
 [1]  6  7  7  8  8  9  9 10 10 11 11 12 12

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

Warning message:
In bxp(list(stats = c(11.8, 13.15, 15.39, 15.59, 16.07, 11.87, 13.47,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/51l291326703528.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,] 11.800 13.150 15.090 15.02 15.420 16.070   NA
[2,] 11.895 13.585 15.165 15.36 15.605 16.515   NA
[3,] 12.895 15.065 15.295 15.45 15.755 17.725   NA
[4,] 13.305 15.430 15.415 15.59 15.950 17.855   NA
[5,] 13.470 15.560 15.500 15.76 16.150 18.070   NA

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

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]     [,6] [,7]
[1,] 12.25189 14.22348 15.18097 15.3451 15.59764 17.11382   NA
[2,] 13.53811 15.90652 15.40903 15.5549 15.91236 18.33618   NA

$out
[1] 15.98

$group
[1] 4

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

Warning message:
In bxp(list(stats = c(11.8, 11.895, 12.895, 13.305, 13.47, 13.15,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6djfc1326703528.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,] 1.451116 4.270 0.33000
[2,] 1.469565 4.495 0.39000
[3,] 1.545810 4.580 0.95125
[4,] 1.657004 4.655 1.51750
[5,] 1.696487 4.680 1.84500

$n
[1] 12 12 12

$conf
         [,1]     [,2]      [,3]
[1,] 1.460318 4.507023 0.4369897
[2,] 1.631302 4.652977 1.4655103

$out
[1] 4.97 5.01

$group
[1] 2 2

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

Warning message:
In bxp(list(stats = c(1.45111566274597, 1.46956508205315, 1.54581046002781,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1f9841326703528.ps tmp/1f9841326703528.png",intern=TRUE))
character(0)
> try(system("convert tmp/2lg7b1326703528.ps tmp/2lg7b1326703528.png",intern=TRUE))
character(0)
> try(system("convert tmp/3v8jh1326703528.ps tmp/3v8jh1326703528.png",intern=TRUE))
character(0)
> try(system("convert tmp/48ei71326703528.ps tmp/48ei71326703528.png",intern=TRUE))
character(0)
> try(system("convert tmp/51l291326703528.ps tmp/51l291326703528.png",intern=TRUE))
character(0)
> try(system("convert tmp/6djfc1326703528.ps tmp/6djfc1326703528.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.386   0.313   1.705