R version 2.12.2 (2011-02-25)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-redhat-linux-gnu (64-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.

  Natural language support but running in an English locale

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(12.94,12.79,12.82,12.85,12.85,12.72,12.62,12.67,12.6,12.54,12.64,12.67,12.51,12.59,12.52,12.5,12.58,12.51,12.47,12.44,12.51,12.27,12.51,12.41,12.35,12.39,12.31,12.31,12.21,12.1,12.01,11.85,12.12,11.96,11.99,11.93,11.91,11.83,11.92,11.86,11.94,11.87,11.86,11.92,11.82,11.85,11.77,11.82,11.61,11.56,11.45,11.4,11.38,11.33,11.19,11.15,10.98,10.92,10.99,11,10.9,10.99,11.04,11.03,10.99,11,10.87,10.88,10.91,10.92,10.83,10.9,10.82,10.79,10.77,10.72,10.71,10.63,10.61,10.57,10.65,10.57,10.57,10.57,10.52,10.43,10.35,10.2,10.2,10.17,10.14,10.05,10.12,10.12)
> 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] 94
> (np <- floor(n / par1))
[1] 7
> 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] 8 8 8 8 8 8 8 8 8 8 7 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,] 12.94 12.51 12.35 11.91 11.61 10.90 10.82 10.52
 [2,] 12.79 12.59 12.39 11.83 11.56 10.99 10.79 10.43
 [3,] 12.82 12.52 12.31 11.92 11.45 11.04 10.77 10.35
 [4,] 12.85 12.50 12.31 11.86 11.40 11.03 10.72 10.20
 [5,] 12.85 12.58 12.21 11.94 11.38 10.99 10.71 10.20
 [6,] 12.72 12.51 12.10 11.87 11.33 11.00 10.63 10.17
 [7,] 12.62 12.47 12.01 11.86 11.19 10.87 10.61 10.14
 [8,] 12.67 12.44 11.85 11.92 11.15 10.88 10.57 10.05
 [9,] 12.60 12.51 12.12 11.82 10.98 10.91 10.65 10.12
[10,] 12.54 12.27 11.96 11.85 10.92 10.92 10.57 10.12
[11,] 12.64 12.51 11.99 11.77 10.99 10.83 10.57    NA
[12,] 12.67 12.41 11.93 11.82 11.00 10.90 10.57    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/wessaorg/rcomp/tmp/199g11305021409.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/wessaorg/rcomp/tmp/2fmsz1305021409.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/wessaorg/rcomp/tmp/3l01z1305021409.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/wessaorg/rcomp/tmp/47hcw1305021409.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,] 10.52 10.430 10.350 10.200 10.200 10.170 10.140 10.050 10.120 10.120 10.57
[2,] 10.86 10.890 10.905 10.875 10.850 10.815 10.740 10.725 10.780 10.745 10.91
[3,] 11.76 11.695 11.685 11.630 11.660 11.600 11.525 11.500 11.400 11.385 11.77
[4,] 12.43 12.490 12.415 12.405 12.395 12.305 12.240 12.180 12.315 12.115 12.25
[5,] 12.94 12.790 12.820 12.850 12.850 12.720 12.620 12.670 12.600 12.540 12.64
     [,12]
[1,] 10.57
[2,] 10.95
[3,] 11.82
[4,] 12.17
[5,] 12.67

$n
 [1] 8 8 8 8 8 8 8 8 8 8 7 7

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 10.88298 10.80122 10.84149 10.77532 10.79694 10.76766 10.68708 10.68722
[2,] 12.63702 12.58878 12.52851 12.48468 12.52306 12.43234 12.36292 12.31278
         [,9]    [,10]    [,11]    [,12]
[1,] 10.54253 10.61970 10.96977 11.09144
[2,] 12.25747 12.15030 12.57023 12.54856

$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(10.52, 10.86, 11.76, 12.43, 12.94, 10.43,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/5csb41305021409.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]   [,8]
[1,] 12.540 12.410 11.850 11.770 10.920 10.830 10.570 10.050
[2,] 12.630 12.455 11.975 11.825 10.995 10.890 10.570 10.120
[3,] 12.695 12.510 12.110 11.860 11.260 10.915 10.640 10.185
[4,] 12.835 12.515 12.310 11.915 11.425 10.995 10.745 10.350
[5,] 12.940 12.590 12.390 11.940 11.610 11.040 10.820 10.520

$n
[1] 12 12 12 12 12 12 12 10

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 12.60150 12.48263 11.95720 11.81895 11.06387 10.86711 10.56018 10.07008
[2,] 12.78850 12.53737 12.26280 11.90105 11.45613 10.96289 10.71982 10.29992

$out
[1] 12.27

$group
[1] 2

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

Warning message:
In bxp(list(stats = c(12.54, 12.63, 12.695, 12.835, 12.94, 12.41,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/6fplw1305021409.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.8281678 2.100 1.20500
[2,] 0.8779667 2.390 1.27125
[3,] 0.8953421 2.475 1.35625
[4,] 0.9238827 2.585 1.39750
[5,] 0.9398138 2.650 1.51000

$n
[1] 12 12 12

$conf
          [,1]     [,2]     [,3]
[1,] 0.8743995 2.386059 1.298667
[2,] 0.9162847 2.563941 1.413833

$out
[1] 0.8028876 2.0700000

$group
[1] 1 2

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

Warning message:
In bxp(list(stats = c(0.828167799882309, 0.877966747972528, 0.895342134011842,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/199g11305021409.ps tmp/199g11305021409.png",intern=TRUE))
character(0)
> try(system("convert tmp/2fmsz1305021409.ps tmp/2fmsz1305021409.png",intern=TRUE))
character(0)
> try(system("convert tmp/3l01z1305021409.ps tmp/3l01z1305021409.png",intern=TRUE))
character(0)
> try(system("convert tmp/47hcw1305021409.ps tmp/47hcw1305021409.png",intern=TRUE))
character(0)
> try(system("convert tmp/5csb41305021409.ps tmp/5csb41305021409.png",intern=TRUE))
character(0)
> try(system("convert tmp/6fplw1305021409.ps tmp/6fplw1305021409.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  3.990   0.380   4.473