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(1.2638,1.2640,1.2261,1.1989,1.2000,1.2146,1.2266,1.2191,1.2224,1.2507,1.2997,1.3406,1.3123,1.3013,1.3185,1.2943,1.2697,1.2155,1.2041,1.2295,1.2234,1.2022,1.1789,1.1861,1.2126,1.1940,1.2028,1.2273,1.2767,1.2661,1.2681,1.2810,1.2722,1.2617,1.2888,1.3205,1.2993,1.3080,1.3246,1.3513,1.3518,1.3421,1.3726,1.3626,1.3910,1.4233,1.4683,1.4559,1.4728,1.4759,1.5520,1.5754,1.5554,1.5562,1.5759,1.4955,1.4342,1.3266,1.2744,1.3511,1.3244,1.2797,1.3050,1.3199,1.3646,1.4014,1.4092,1.4266,1.4575,1.4821,1.4908,1.4579,1.4266,1.3680,1.3570,1.3417,1.2563,1.2223,1.2811,1.2903,1.3103,1.3901,1.3654,1.3221)
> 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] 84
> (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] 7 7 7 7 7 7 7 7 7 7 7 7
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7] [,8]
 [1,] 1.2638 1.3123 1.2126 1.2993 1.4728 1.3244 1.4266   NA
 [2,] 1.2640 1.3013 1.1940 1.3080 1.4759 1.2797 1.3680   NA
 [3,] 1.2261 1.3185 1.2028 1.3246 1.5520 1.3050 1.3570   NA
 [4,] 1.1989 1.2943 1.2273 1.3513 1.5754 1.3199 1.3417   NA
 [5,] 1.2000 1.2697 1.2767 1.3518 1.5554 1.3646 1.2563   NA
 [6,] 1.2146 1.2155 1.2661 1.3421 1.5562 1.4014 1.2223   NA
 [7,] 1.2266 1.2041 1.2681 1.3726 1.5759 1.4092 1.2811   NA
 [8,] 1.2191 1.2295 1.2810 1.3626 1.4955 1.4266 1.2903   NA
 [9,] 1.2224 1.2234 1.2722 1.3910 1.4342 1.4575 1.3103   NA
[10,] 1.2507 1.2022 1.2617 1.4233 1.3266 1.4821 1.3901   NA
[11,] 1.2997 1.1789 1.2888 1.4683 1.2744 1.4908 1.3654   NA
[12,] 1.3406 1.1861 1.3205 1.4559 1.3511 1.4579 1.3221   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/1ir2x1304872178.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/26g4a1304872178.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/3r8dp1304872178.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/48tys1304872178.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]
[1,] 1.21260 1.19400 1.20280 1.1989 1.2000 1.21460 1.20410 1.21910 1.2224
[2,] 1.28155 1.27185 1.26555 1.2608 1.2630 1.21890 1.24735 1.25525 1.2478
[3,] 1.31230 1.30130 1.31850 1.3199 1.2767 1.26610 1.28110 1.29030 1.3103
[4,] 1.37550 1.33800 1.34080 1.3465 1.3582 1.37175 1.39090 1.39460 1.4126
[5,] 1.47280 1.36800 1.35700 1.3513 1.3646 1.55620 1.57590 1.49550 1.4575
      [,10]   [,11]  [,12]
[1,] 1.2022 1.17890 1.3205
[2,] 1.2562 1.28160 1.3213
[3,] 1.3266 1.29970 1.3406
[4,] 1.4067 1.41685 1.4035
[5,] 1.4821 1.49080 1.4579

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 1.256195 1.261796 1.273562 1.268721 1.219848 1.174820 1.195374 1.207082
[2,] 1.368405 1.340804 1.363438 1.371079 1.333552 1.357380 1.366826 1.373518
         [,9]    [,10]    [,11]    [,12]
[1,] 1.211884 1.236724 1.218931 1.291511
[2,] 1.408716 1.416476 1.380469 1.389689

$out
[1] 1.4759 1.5520 1.5754 1.5554 1.1861

$group
[1]  2  3  4  5 12

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

Warning message:
In bxp(list(stats = c(1.2126, 1.28155, 1.3123, 1.3755, 1.4728, 1.194,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5bd2o1304872178.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,] 1.19890 1.17890 1.19400 1.29930 1.27440 1.27970 1.2223   NA
[2,] 1.21685 1.20315 1.21995 1.33335 1.39265 1.32215 1.2857   NA
[3,] 1.22635 1.22645 1.26710 1.35720 1.48570 1.40530 1.3319   NA
[4,] 1.26390 1.29780 1.27885 1.40715 1.55580 1.45770 1.3667   NA
[5,] 1.29970 1.31850 1.32050 1.46830 1.57590 1.49080 1.4266   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7] [,8]
[1,] 1.204890 1.183279 1.240235 1.323539 1.411286 1.343475 1.294955   NA
[2,] 1.247810 1.269621 1.293965 1.390861 1.560114 1.467125 1.368845   NA

$out
[1] 1.3406

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(1.1989, 1.21685, 1.22635, 1.2639, 1.2997,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6o8xp1304872178.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.08887650 0.2351 0.066150
[2,] 0.09523787 0.2741 0.083950
[3,] 0.10735134 0.2969 0.115225
[4,] 0.11934370 0.3523 0.147025
[5,] 0.12987852 0.3765 0.164800

$n
[1] 12 12 12

$conf
           [,1]      [,2]       [,3]
[1,] 0.09635651 0.2612325 0.08645607
[2,] 0.11834617 0.3325675 0.14399393

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.0888764953661621, 0.0952378736030116, 0.107351340526094,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1ir2x1304872178.ps tmp/1ir2x1304872178.png",intern=TRUE))
character(0)
> try(system("convert tmp/26g4a1304872178.ps tmp/26g4a1304872178.png",intern=TRUE))
character(0)
> try(system("convert tmp/3r8dp1304872178.ps tmp/3r8dp1304872178.png",intern=TRUE))
character(0)
> try(system("convert tmp/48tys1304872178.ps tmp/48tys1304872178.png",intern=TRUE))
character(0)
> try(system("convert tmp/5bd2o1304872178.ps tmp/5bd2o1304872178.png",intern=TRUE))
character(0)
> try(system("convert tmp/6o8xp1304872178.ps tmp/6o8xp1304872178.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.816   0.568   3.367