R version 2.8.1 (2008-12-22)
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(46,50,49,48,50,47,50,49,51,52,48,55,56,43,44,50,49,47,46,50,49,53,54,56,56,58,53,51,52,53,56,54,54,56,59,62,62,73,76,80,77,81,80,80,81,80,77,71,71,64,64,47,41,35,34,33,23,16,16,8,9)
> 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] 61
> (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] 6 5 5 5 5 5 5 5 5 5 5 5
> arr
      [,1] [,2] [,3] [,4] [,5] [,6]
 [1,]   46   56   56   62   71    9
 [2,]   50   43   58   73   64   NA
 [3,]   49   44   53   76   64   NA
 [4,]   48   50   51   80   47   NA
 [5,]   50   49   52   77   41   NA
 [6,]   47   47   53   81   35   NA
 [7,]   50   46   56   80   34   NA
 [8,]   49   50   54   80   33   NA
 [9,]   51   49   54   81   23   NA
[10,]   52   53   56   80   16   NA
[11,]   48   54   59   77   16   NA
[12,]   55   56   62   71    8   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/1wotw1258052742.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/rcomp/tmp/2ewt71258052742.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/rcomp/tmp/39ayo1258052742.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/rcomp/tmp/4f1wh1258052742.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] [,11] [,12]
[1,]   46   43   44   47   49   47   34   49   49    52    48    55
[2,]   46   50   49   48   49   47   46   49   49    52    48    55
[3,]   56   58   53   50   50   47   50   50   51    53    54    56
[4,]   62   64   64   51   52   53   56   54   54    56    59    62
[5,]   71   73   76   51   52   53   56   54   54    56    59    71

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 45.67948 48.10764 42.40104 47.88021 47.88021 42.76042 42.93403 46.46701
[2,] 66.32052 67.89236 63.59896 52.11979 52.11979 51.23958 57.06597 53.53299
         [,9]    [,10]    [,11]    [,12]
[1,] 47.46701 50.17361 46.22743 51.05382
[2,] 54.53299 55.82639 61.77257 60.94618

$out
 [1]  9 80 77 41 81 35 80 80 33 81 23 80 16 77 16  8

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

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

Warning message:
In bxp(list(stats = c(46, 46, 56, 62, 71, 43, 50, 58, 64, 73, 44,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/544e71258052742.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]
[1,] 46.0 43.0   51 71.0  8.0    9
[2,] 48.0 46.5   53 74.5 19.5    9
[3,] 49.5 49.5   55 78.5 34.5    9
[4,] 50.5 53.5   57 80.0 55.5    9
[5,] 52.0 56.0   62 81.0 71.0    9

$n
[1] 12 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 48.35973 46.30725 53.17557 75.99141 18.08016    9
[2,] 50.64027 52.69275 56.82443 81.00859 50.91984    9

$out
[1] 55 62

$group
[1] 1 4

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

Warning message:
In bxp(list(stats = c(46, 48, 49.5, 50.5, 52, 43, 46.5, 49.5, 53.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6pjx91258052742.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,] 11.71751 30.0  3.0
[2,] 13.79034 34.5  4.5
[3,] 17.08781 46.5  6.5
[4,] 21.98112 61.5 11.5
[5,] 24.54180 64.0 15.0

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 13.35194 34.18512 3.307253
[2,] 20.82368 58.81488 9.692747

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(11.7175082675456, 13.7903364554659, 17.0878110803264,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1wotw1258052742.ps tmp/1wotw1258052742.png")
> system("convert tmp/2ewt71258052742.ps tmp/2ewt71258052742.png")
> system("convert tmp/39ayo1258052742.ps tmp/39ayo1258052742.png")
> system("convert tmp/4f1wh1258052742.ps tmp/4f1wh1258052742.png")
> system("convert tmp/544e71258052742.ps tmp/544e71258052742.png")
> system("convert tmp/6pjx91258052742.ps tmp/6pjx91258052742.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.300   1.140   2.258