R version 2.13.0 (2011-04-13)
Copyright (C) 2011 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(7.08,7.08,7.09,7.07,7.06,6.99,6.99,6.99,6.98,6.96,6.95,6.91,6.91,6.87,6.91,6.89,6.88,6.9,6.91,6.85,6.86,6.82,6.8,6.83,6.84,6.89,7.14,7.21,7.25,7.31,7.3,7.48,7.49,7.4,7.44,7.42,7.14,7.24,7.33,7.61,7.66,7.69,7.7,7.68,7.71,7.71,7.72,7.68,7.72,7.74,7.76,7.9,7.97,7.96,7.95,7.97,7.93,7.99,7.96,7.92,7.97,7.98,8,8.04,8.17,8.29,8.26,8.3,8.32,8.28,8.27,8.32,8.31,8.34,8.32,8.36,8.33,8.35,8.34,8.37,8.31,8.33,8.34,8.25,8.27,8.31,8.25,8.3,8.3,8.35,8.78,8.9,8.9,8.9,9,9.05)
> 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] 96
> (np <- floor(n / par1))
[1] 8
> 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 8 8
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9]
 [1,] 7.08 6.91 6.84 7.14 7.72 7.97 8.31 8.27   NA
 [2,] 7.08 6.87 6.89 7.24 7.74 7.98 8.34 8.31   NA
 [3,] 7.09 6.91 7.14 7.33 7.76 8.00 8.32 8.25   NA
 [4,] 7.07 6.89 7.21 7.61 7.90 8.04 8.36 8.30   NA
 [5,] 7.06 6.88 7.25 7.66 7.97 8.17 8.33 8.30   NA
 [6,] 6.99 6.90 7.31 7.69 7.96 8.29 8.35 8.35   NA
 [7,] 6.99 6.91 7.30 7.70 7.95 8.26 8.34 8.78   NA
 [8,] 6.99 6.85 7.48 7.68 7.97 8.30 8.37 8.90   NA
 [9,] 6.98 6.86 7.49 7.71 7.93 8.32 8.31 8.90   NA
[10,] 6.96 6.82 7.40 7.71 7.99 8.28 8.33 8.90   NA
[11,] 6.95 6.80 7.44 7.72 7.96 8.27 8.34 9.00   NA
[12,] 6.91 6.83 7.42 7.68 7.92 8.32 8.25 9.05   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/wessaorg/rcomp/tmp/148fc1322849938.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/wessaorg/rcomp/tmp/23p8u1322849938.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/wessaorg/rcomp/tmp/3q89q1322849938.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/wessaorg/rcomp/tmp/4c2dw1322849938.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,] 6.840 6.870 6.910 6.890 6.880 6.900 6.910 6.850 6.860 6.820 6.800 6.830
[2,] 6.995 6.985 7.115 7.140 7.155 7.150 7.145 7.235 7.235 7.180 7.195 7.165
[3,] 7.430 7.490 7.545 7.755 7.815 7.825 7.825 7.825 7.820 7.850 7.840 7.800
[4,] 8.120 8.145 8.125 8.170 8.235 8.320 8.300 8.335 8.315 8.305 8.305 8.285
[5,] 8.310 8.340 8.320 8.360 8.330 8.350 8.780 8.900 8.900 8.900 9.000 9.050

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]   [,7]     [,8]
[1,] 6.801559 6.842007 6.980799 7.179627 7.211696 7.171421 7.1798 7.210524
[2,] 8.058441 8.137993 8.109201 8.330373 8.418304 8.478579 8.4702 8.439476
         [,9]    [,10]    [,11]    [,12]
[1,] 7.216696 7.221559 7.219938 7.174352
[2,] 8.423304 8.478441 8.460062 8.425648

$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(6.84, 6.995, 7.43, 8.12, 8.31, 6.87, 6.985,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5fwqe1322849938.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] [,9]
[1,] 6.910 6.800 6.840 7.140 7.720 7.970 8.310 8.250   NA
[2,] 6.970 6.840 7.175 7.470 7.830 8.020 8.315 8.300   NA
[3,] 6.990 6.875 7.305 7.680 7.940 8.265 8.335 8.565   NA
[4,] 7.075 6.905 7.430 7.705 7.965 8.295 8.345 8.900   NA
[5,] 7.090 6.910 7.490 7.720 7.990 8.320 8.370 9.050   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 6.942109 6.845353 7.188693 7.572815 7.878426 8.139571 8.321317 8.291336
[2,] 7.037891 6.904647 7.421307 7.787185 8.001574 8.390429 8.348683 8.838664
     [,9]
[1,]   NA
[2,]   NA

$out
[1] 8.25

$group
[1] 7

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

Warning message:
In bxp(list(stats = c(6.91, 6.97, 6.99, 7.075, 7.09, 6.8, 6.84,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/66jlz1322849938.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.5538437 1.410 0.93000
[2,] 0.5906550 1.460 0.95500
[3,] 0.6444773 1.670 0.98750
[4,] 0.7110279 2.065 1.01875
[5,] 0.7503856 2.220 1.07500

$n
[1] 12 12 12

$conf
          [,1]     [,2]      [,3]
[1,] 0.5895744 1.394055 0.9584232
[2,] 0.6993802 1.945945 1.0165768

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.553843711838328, 0.590655005021854, 0.644477258502327,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/148fc1322849938.ps tmp/148fc1322849938.png",intern=TRUE))
character(0)
> try(system("convert tmp/23p8u1322849938.ps tmp/23p8u1322849938.png",intern=TRUE))
character(0)
> try(system("convert tmp/3q89q1322849938.ps tmp/3q89q1322849938.png",intern=TRUE))
character(0)
> try(system("convert tmp/4c2dw1322849938.ps tmp/4c2dw1322849938.png",intern=TRUE))
character(0)
> try(system("convert tmp/5fwqe1322849938.ps tmp/5fwqe1322849938.png",intern=TRUE))
character(0)
> try(system("convert tmp/66jlz1322849938.ps tmp/66jlz1322849938.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.354   0.274   1.638