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(3.11,3.1,3.13,3.14,3.16,3.17,3.19,3.18,3.17,3.17,3.2,3.21,3.25,3.22,3.25,3.31,3.35,3.37,3.38,3.38,3.39,3.44,3.56,3.65,3.69,3.71,3.71,3.74,3.75,3.77,3.75,3.78,3.8,3.78,3.79,3.8,3.82,3.82,3.84,3.86,3.8,3.85,3.78,3.79,3.77,3.78,3.77,3.76,3.78,3.76,3.75,3.71,3.72,3.7,3.69,3.7,3.72,3.73,3.73,3.69)
> 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] 60
> (np <- floor(n / par1))
[1] 5
> 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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,] 3.11 3.25 3.69 3.82 3.78   NA
 [2,] 3.10 3.22 3.71 3.82 3.76   NA
 [3,] 3.13 3.25 3.71 3.84 3.75   NA
 [4,] 3.14 3.31 3.74 3.86 3.71   NA
 [5,] 3.16 3.35 3.75 3.80 3.72   NA
 [6,] 3.17 3.37 3.77 3.85 3.70   NA
 [7,] 3.19 3.38 3.75 3.78 3.69   NA
 [8,] 3.18 3.38 3.78 3.79 3.70   NA
 [9,] 3.17 3.39 3.80 3.77 3.72   NA
[10,] 3.17 3.44 3.78 3.78 3.73   NA
[11,] 3.20 3.56 3.79 3.77 3.73   NA
[12,] 3.21 3.65 3.80 3.76 3.69   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/1vjyt1322592923.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/2zu851322592923.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/3wszw1322592923.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/4h4us1322592923.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,] 3.11 3.10 3.13 3.14 3.16 3.17 3.19 3.18 3.17  3.17  3.56  3.65
[2,] 3.25 3.22 3.25 3.31 3.35 3.37 3.38 3.38 3.39  3.44  3.56  3.65
[3,] 3.69 3.71 3.71 3.71 3.72 3.70 3.69 3.70 3.72  3.73  3.73  3.69
[4,] 3.78 3.76 3.75 3.74 3.75 3.77 3.75 3.78 3.77  3.78  3.77  3.76
[5,] 3.82 3.82 3.84 3.86 3.80 3.85 3.78 3.79 3.80  3.78  3.79  3.80

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 5

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 3.315503 3.328437 3.356701 3.406163 3.437361 3.417361 3.428559 3.417361
[2,] 4.064497 4.091563 4.063299 4.013837 4.002639 3.982639 3.951441 3.982639
         [,9]    [,10]    [,11]    [,12]
[1,] 3.451493 3.489757 3.581615 3.612274
[2,] 3.988507 3.970243 3.878385 3.767726

$out
[1] 3.20 3.21

$group
[1] 11 12

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

Warning message:
In bxp(list(stats = c(3.11, 3.25, 3.69, 3.78, 3.82, 3.1, 3.22, 3.71,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5t6op1322592923.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]
[1,] 3.100 3.220 3.690 3.760 3.69   NA
[2,] 3.135 3.280 3.725 3.775 3.70   NA
[3,] 3.170 3.375 3.760 3.795 3.72   NA
[4,] 3.185 3.415 3.785 3.830 3.74   NA
[5,] 3.210 3.560 3.800 3.860 3.78   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 3.147195 3.313426 3.732634 3.769914 3.701756   NA
[2,] 3.192805 3.436574 3.787366 3.820086 3.738244   NA

$out
[1] 3.65

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(3.1, 3.135, 3.17, 3.185, 3.21, 3.22, 3.28,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6lsby1322592923.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,] 0.2376342 0.590 0.210
[2,] 0.2647428 0.600 0.355
[3,] 0.2808524 0.635 0.400
[4,] 0.3158889 0.710 0.465
[5,] 0.3354400 0.720 0.540

$n
[1] 12 12 12

$conf
          [,1]      [,2]      [,3]
[1,] 0.2575243 0.5848283 0.3498283
[2,] 0.3041805 0.6851717 0.4501717

$out
[1] 0.11

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(0.23763417262675, 0.264742826288551, 0.280852401561609,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1vjyt1322592923.ps tmp/1vjyt1322592923.png",intern=TRUE))
character(0)
> try(system("convert tmp/2zu851322592923.ps tmp/2zu851322592923.png",intern=TRUE))
character(0)
> try(system("convert tmp/3wszw1322592923.ps tmp/3wszw1322592923.png",intern=TRUE))
character(0)
> try(system("convert tmp/4h4us1322592923.ps tmp/4h4us1322592923.png",intern=TRUE))
character(0)
> try(system("convert tmp/5t6op1322592923.ps tmp/5t6op1322592923.png",intern=TRUE))
character(0)
> try(system("convert tmp/6lsby1322592923.ps tmp/6lsby1322592923.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.426   0.324   1.758