R version 2.12.1 (2010-12-16)
Copyright (C) 2010 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(36.68,36.77,36.78,36.78,37,37.12,37.3,37.34,37.4,37.4,37.34,37.29,37.39,37.42,37.42,43.49,44.3,44.36,44.52,44.66,44.77,44.77,44.82,44.97,45.28,45.5,45.52,45.52,45.17,45.25,45.32,45.41,45.44,45.44,45.46,45.53,45.61,45.8,45.83,45.83,45.96,46.01,46.18,46.32,46.51,46.51,46.56,46.54,46.62,46.76,46.82,46.82,46.7,46.72,46.47,46.74,46.89,46.89,46.96,47.6,47.67,47.58,47.58,47.58,47.57,47.53,47.68,47.56,47.81,47.81,47.81,47.81,48.12,48.13,48.01)
> 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] 75
> (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] 7 7 7 6 6 6 6 6 6 6 6 6
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]
 [1,] 36.68 37.39 45.28 45.61 46.62 47.67 48.12
 [2,] 36.77 37.42 45.50 45.80 46.76 47.58 48.13
 [3,] 36.78 37.42 45.52 45.83 46.82 47.58 48.01
 [4,] 36.78 43.49 45.52 45.83 46.82 47.58    NA
 [5,] 37.00 44.30 45.17 45.96 46.70 47.57    NA
 [6,] 37.12 44.36 45.25 46.01 46.72 47.53    NA
 [7,] 37.30 44.52 45.32 46.18 46.47 47.68    NA
 [8,] 37.34 44.66 45.41 46.32 46.74 47.56    NA
 [9,] 37.40 44.77 45.44 46.51 46.89 47.81    NA
[10,] 37.40 44.77 45.44 46.51 46.89 47.81    NA
[11,] 37.34 44.82 45.46 46.56 46.96 47.81    NA
[12,] 37.29 44.97 45.53 46.54 47.60 47.81    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/rcomp/tmp/181eo1323901745.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/rcomp/tmp/2whv71323901745.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/rcomp/tmp/33jys1323901745.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/rcomp/tmp/4a7pf1323901745.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,] 36.680 36.77 36.78 43.490 44.300 44.36 44.52 44.660 44.770 44.770 44.82
[2,] 41.335 41.46 41.47 43.490 44.300 44.36 44.52 44.660 44.770 44.770 44.82
[3,] 45.610 45.80 45.83 45.675 45.565 45.63 45.75 45.865 45.975 45.975 46.01
[4,] 47.145 47.17 47.20 46.820 46.700 46.72 46.47 46.740 46.890 46.890 46.96
[5,] 48.120 48.13 48.01 47.580 47.570 47.53 47.68 47.560 47.810 47.810 47.81
      [,12]
[1,] 44.970
[2,] 44.970
[3,] 46.035
[4,] 47.600
[5,] 47.810

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 42.14036 42.39008 42.40814 43.52704 44.01692 44.10772 44.49219 44.52333
[2,] 49.07964 49.20992 49.25186 47.82296 47.11308 47.15228 47.00781 47.20667
         [,9]    [,10]    [,11]    [,12]
[1,] 44.60753 44.60753 44.62963 44.33857
[2,] 47.34247 47.34247 47.39037 47.73143

$out
[1] 36.78 37.00 37.12 37.30 37.34 37.40 37.40 37.34 37.29

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

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

Warning message:
In bxp(list(stats = c(36.68, 41.335, 45.61, 47.145, 48.12, 36.77,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5kkut1323901745.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,] 36.680 37.390 45.17 45.610 46.47 47.530 48.010
[2,] 36.780 40.455 45.30 45.830 46.71 47.575 48.065
[3,] 37.205 44.440 45.44 46.095 46.79 47.625 48.120
[4,] 37.340 44.770 45.51 46.510 46.89 47.810 48.125
[5,] 37.400 44.970 45.53 46.560 46.96 47.810 48.130

$n
[1] 12 12 12 12 12 12  3

$conf
         [,1]    [,2]     [,3]     [,4]    [,5]     [,6]     [,7]
[1,] 36.94958 42.4719 45.34422 45.78485 46.7079 47.51781 48.06527
[2,] 37.46042 46.4081 45.53578 46.40515 46.8721 47.73219 48.17473

$out
[1] 47.6

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(36.68, 36.78, 37.205, 37.34, 37.4, 37.39,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6lyd51323901745.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,] 3.722679 10.220 1.67750
[2,] 3.780987 10.410 1.85750
[3,] 3.821354 10.495 1.97875
[4,] 4.373901 11.015 4.14250
[5,] 4.808950 11.440 5.81000

$n
[1] 12 12 12

$conf
         [,1]     [,2]      [,3]
[1,] 3.550922 10.21906 0.9365462
[2,] 4.091786 10.77094 3.0209538

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(3.72267869506175, 3.78098668968724, 3.82135410745374,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/181eo1323901745.ps tmp/181eo1323901745.png",intern=TRUE))
character(0)
> try(system("convert tmp/2whv71323901745.ps tmp/2whv71323901745.png",intern=TRUE))
character(0)
> try(system("convert tmp/33jys1323901745.ps tmp/33jys1323901745.png",intern=TRUE))
character(0)
> try(system("convert tmp/4a7pf1323901745.ps tmp/4a7pf1323901745.png",intern=TRUE))
character(0)
> try(system("convert tmp/5kkut1323901745.ps tmp/5kkut1323901745.png",intern=TRUE))
character(0)
> try(system("convert tmp/6lyd51323901745.ps tmp/6lyd51323901745.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.516   0.408   1.951