R version 2.9.0 (2009-04-17)
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(3.65,3.59,3.31,3.89,4.31,4.35,4.11,3.90,3.75,3.75,3.88,3.93,3.97,3.97,4.33,4.16,4.93,3.86,4.06,4.18,4.08,4.38,4.48,4.41,4.37,4.56,4.71,4.94,5.03,5.08,5.05,4.83,4.68,4.69,4.58,4.54,4.75,4.71,4.50,4.62,4.69,5.05,4.93,4.53,4.33,4.33,3.87,3.74,3.31,3.21,2.93,3.19,3.46,3.73,3.60,3.46,3.25,3.19,2.82,1.89,1.98,2.30,2.42,2.47,2.81,3.37,3.14,3.21,3.02,2.96,2.92,3.07)
> 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,] 3.65 3.97 4.37 4.75 3.31 1.98   NA
 [2,] 3.59 3.97 4.56 4.71 3.21 2.30   NA
 [3,] 3.31 4.33 4.71 4.50 2.93 2.42   NA
 [4,] 3.89 4.16 4.94 4.62 3.19 2.47   NA
 [5,] 4.31 4.93 5.03 4.69 3.46 2.81   NA
 [6,] 4.35 3.86 5.08 5.05 3.73 3.37   NA
 [7,] 4.11 4.06 5.05 4.93 3.60 3.14   NA
 [8,] 3.90 4.18 4.83 4.53 3.46 3.21   NA
 [9,] 3.75 4.08 4.68 4.33 3.25 3.02   NA
[10,] 3.75 4.38 4.69 4.33 3.19 2.96   NA
[11,] 3.88 4.48 4.58 3.87 2.82 2.92   NA
[12,] 3.93 4.41 4.54 3.74 1.89 3.07   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/11o6q1291552269.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/www/html/rcomp/tmp/2bf5b1291552269.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/www/html/rcomp/tmp/3bf5b1291552269.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/www/html/rcomp/tmp/4bf5b1291552269.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,] 1.98 2.30 2.42 2.470 2.81 3.370 3.140 3.21 3.020  2.96 2.820 1.890
[2,] 3.31 3.21 2.93 3.190 3.46 3.730 3.600 3.46 3.250  3.19 2.920 3.070
[3,] 3.81 3.78 3.82 4.025 4.50 4.105 4.085 4.04 3.915  4.04 3.875 3.835
[4,] 4.37 4.56 4.50 4.620 4.93 5.050 4.930 4.53 4.330  4.38 4.480 4.410
[5,] 4.75 4.71 4.71 4.940 5.03 5.080 5.050 4.83 4.680  4.69 4.580 4.540

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 3.126266 2.909206 2.807299 3.102604 3.551803 3.253557 3.227107 3.349815
[2,] 4.493734 4.650794 4.832701 4.947396 5.448197 4.956443 4.942893 4.730185
         [,9]    [,10]    [,11]    [,12]
[1,] 3.218365 3.272412 2.868750 2.970657
[2,] 4.611635 4.807588 4.881250 4.699343

$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(1.98, 3.31, 3.81, 4.37, 4.75, 2.3, 3.21, 3.78,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/54p4w1291552269.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,] 3.310 3.860 4.370 3.740 2.82 1.980   NA
[2,] 3.700 4.015 4.570 4.330 3.06 2.445   NA
[3,] 3.885 4.170 4.700 4.575 3.23 2.940   NA
[4,] 4.020 4.395 4.985 4.730 3.46 3.105   NA
[5,] 4.350 4.930 5.080 5.050 3.73 3.370   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 3.739046 3.996679 4.510716 4.392557 3.047557 2.638970   NA
[2,] 4.030954 4.343321 4.889284 4.757443 3.412443 3.241030   NA

$out
[1] 1.89

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(3.31, 3.7, 3.885, 4.02, 4.35, 3.86, 4.015,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/64p4w1291552269.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.6206260 1.620 0.87250
[2,] 0.7058861 1.720 0.95125
[3,] 0.8192081 2.065 1.08000
[4,] 0.9302043 2.440 1.15625
[5,] 0.9872318 2.770 1.43250

$n
[1] 12 12 12

$conf
          [,1]     [,2]      [,3]
[1,] 0.7168950 1.736603 0.9864981
[2,] 0.9215211 2.393397 1.1735019

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.620626028028689, 0.705886056078201, 0.819208053257305,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/11o6q1291552269.ps tmp/11o6q1291552269.png",intern=TRUE))
character(0)
> try(system("convert tmp/2bf5b1291552269.ps tmp/2bf5b1291552269.png",intern=TRUE))
character(0)
> try(system("convert tmp/3bf5b1291552269.ps tmp/3bf5b1291552269.png",intern=TRUE))
character(0)
> try(system("convert tmp/4bf5b1291552269.ps tmp/4bf5b1291552269.png",intern=TRUE))
character(0)
> try(system("convert tmp/54p4w1291552269.ps tmp/54p4w1291552269.png",intern=TRUE))
character(0)
> try(system("convert tmp/64p4w1291552269.ps tmp/64p4w1291552269.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.065   0.807   2.537