R version 2.12.1 (2010-12-16)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i686-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.

  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.32,12.34,12.36,12.54,12.77,12.79,12.96,12.96,13,13.19,13.25,13.61,13.8,13.83,14.04,14.16,14.2,14.27,14.31,14.69,14.9,14.92,15.01,15.09,15.14,15.24,15.33,15.36,15.44,15.5,15.58,15.65,15.72,15.82,15.87,16.07,16.18,16.19,16.39,16.54,16.61,16.62,16.66,16.71,16.72,16.79,16.82,16.83,16.91,16.97,17.02,17.03,17.04,17.07,17.11,17.12,17.14,17.18,17.24,17.26,17.26,17.29,17.36,17.44,17.48,17.48,17.52,17.54,17.58,17.64,17.69,17.69,17.76,17.79,17.82,17.89,17.95,18,18.03,18.06,18.08,18.13,18.16,18.18,18.18,18.27,18.31,18.35,18.45,18.5)
> 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] 90
> (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 7 7 7 7 7 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,] 12.32 13.80 15.14 16.18 16.91 17.26 17.76 18.18
 [2,] 12.34 13.83 15.24 16.19 16.97 17.29 17.79 18.27
 [3,] 12.36 14.04 15.33 16.39 17.02 17.36 17.82 18.31
 [4,] 12.54 14.16 15.36 16.54 17.03 17.44 17.89 18.35
 [5,] 12.77 14.20 15.44 16.61 17.04 17.48 17.95 18.45
 [6,] 12.79 14.27 15.50 16.62 17.07 17.48 18.00 18.50
 [7,] 12.96 14.31 15.58 16.66 17.11 17.52 18.03    NA
 [8,] 12.96 14.69 15.65 16.71 17.12 17.54 18.06    NA
 [9,] 13.00 14.90 15.72 16.72 17.14 17.58 18.08    NA
[10,] 13.19 14.92 15.82 16.79 17.18 17.64 18.13    NA
[11,] 13.25 15.01 15.87 16.82 17.24 17.69 18.16    NA
[12,] 13.61 15.09 16.07 16.83 17.26 17.69 18.18    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/freestat/rcomp/tmp/10kql1305062748.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/freestat/rcomp/tmp/2uv461305062748.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/freestat/rcomp/tmp/3jaas1305062748.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/freestat/rcomp/tmp/4jlqo1305062748.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,] 12.320 12.340 12.360 12.540 12.770 12.790 12.960 12.96 13.00 13.19 13.250
[2,] 14.470 14.535 14.685 14.760 14.820 14.885 14.945 15.17 15.31 15.37 15.440
[3,] 16.545 16.580 16.705 16.785 16.825 16.845 16.660 16.71 16.72 16.79 16.820
[4,] 17.510 17.540 17.590 17.665 17.715 17.740 17.315 17.33 17.36 17.41 17.465
[5,] 18.180 18.270 18.310 18.350 18.450 18.500 18.030 18.06 18.08 18.13 18.160
      [,12]
[1,] 13.610
[2,] 15.580
[3,] 16.830
[4,] 17.475
[5,] 18.180

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 14.84681 14.90136 15.08223 15.16223 15.20781 15.25016 15.24467 15.42008
[2,] 18.24319 18.25864 18.32777 18.40777 18.44219 18.43984 18.07533 17.99992
         [,9]    [,10]    [,11]    [,12]
[1,] 15.49577 15.57174 15.61070 15.69834
[2,] 17.94423 18.00826 18.02930 17.96166

$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(12.32, 14.47, 16.545, 17.51, 18.18, 12.34,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/59qcz1305062748.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.320 13.80 15.140 16.180 16.910 17.26 17.760 18.18
[2,] 12.450 14.10 15.345 16.465 17.025 17.40 17.855 18.27
[3,] 12.875 14.29 15.540 16.640 17.090 17.50 18.015 18.33
[4,] 13.095 14.91 15.770 16.755 17.160 17.61 18.105 18.45
[5,] 13.610 15.09 16.070 16.830 17.260 17.69 18.180 18.50

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 12.58081 13.92055 15.34615 16.50773 17.02843 17.40422 17.90097 18.21389
[2,] 13.16919 14.65945 15.73385 16.77227 17.15157 17.59578 18.12903 18.44611

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(12.32, 12.45, 12.875, 13.095, 13.61, 13.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6r9yz1305062748.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,] 1.600281 4.570 1.89500
[2,] 1.748300 5.005 2.04500
[3,] 1.901324 5.390 2.39375
[4,] 2.022176 5.835 2.48000
[5,] 2.059105 5.950 2.58000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1.776407 5.011431 2.195344
[2,] 2.026240 5.768569 2.592156

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.60028122528511, 1.74829969779872, 1.9013238624931,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/10kql1305062748.ps tmp/10kql1305062748.png",intern=TRUE))
character(0)
> try(system("convert tmp/2uv461305062748.ps tmp/2uv461305062748.png",intern=TRUE))
character(0)
> try(system("convert tmp/3jaas1305062748.ps tmp/3jaas1305062748.png",intern=TRUE))
character(0)
> try(system("convert tmp/4jlqo1305062748.ps tmp/4jlqo1305062748.png",intern=TRUE))
character(0)
> try(system("convert tmp/59qcz1305062748.ps tmp/59qcz1305062748.png",intern=TRUE))
character(0)
> try(system("convert tmp/6r9yz1305062748.ps tmp/6r9yz1305062748.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.628   0.192   1.927