R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(0.63709,0.64218,0.65711,0.66977,0.68255,0.68902,0.71322,0.70224,0.70045,0.69919,0.69693,0.69763,0.69278,0.70196,0.69215,0.6769,0.67124,0.66533,0.67157,0.66428,0.66576,0.66942,0.68130,0.69144,0.69862,0.695,0.69867,0.68968,0.69233,0.68293,0.68399,0.66895,0.68756,0.68527,0.6776,0.68137,0.67933,0.67922,0.68598,0.68297,0.68935,0.69463,0.6833,0.68666,0.68782,0.67669,0.67511,0.67254,0.67397,0.67286,0.66341,0.668,0.68021,0.67934,0.68136,0.67562,0.6744,0.67766,0.68887,0.69614,0.70896,0.72064,0.74725,0.75094,0.77494,0.79487,0.79209,0.79152,0.79308,0.79279,0.79924,0.78668)
> 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] 72
> (np <- floor(n / par1))
[1] 6
> 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 6 6 6 6 6 6 6 6 6 6 6
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6] [,7]
 [1,] 0.63709 0.69278 0.69862 0.67933 0.67397 0.70896   NA
 [2,] 0.64218 0.70196 0.69500 0.67922 0.67286 0.72064   NA
 [3,] 0.65711 0.69215 0.69867 0.68598 0.66341 0.74725   NA
 [4,] 0.66977 0.67690 0.68968 0.68297 0.66800 0.75094   NA
 [5,] 0.68255 0.67124 0.69233 0.68935 0.68021 0.77494   NA
 [6,] 0.68902 0.66533 0.68293 0.69463 0.67934 0.79487   NA
 [7,] 0.71322 0.67157 0.68399 0.68330 0.68136 0.79209   NA
 [8,] 0.70224 0.66428 0.66895 0.68666 0.67562 0.79152   NA
 [9,] 0.70045 0.66576 0.68756 0.68782 0.67440 0.79308   NA
[10,] 0.69919 0.66942 0.68527 0.67669 0.67766 0.79279   NA
[11,] 0.69693 0.68130 0.67760 0.67511 0.68887 0.79924   NA
[12,] 0.69763 0.69144 0.68137 0.67254 0.69614 0.78668   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/rcomp/tmp/1ibg81244128797.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/html/rcomp/tmp/2v9rz1244128797.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/html/rcomp/tmp/39ft51244128797.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/html/rcomp/tmp/4z1gm1244128797.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]
[1,] 0.637090 0.64218 0.657110 0.668000 0.67124 0.665330 0.671570 0.66428
[2,] 0.673970 0.67286 0.663410 0.669770 0.68021 0.679340 0.681360 0.66895
[3,] 0.686055 0.68711 0.689065 0.679935 0.68595 0.685975 0.683645 0.68114
[4,] 0.698620 0.70196 0.698670 0.689680 0.69233 0.694630 0.713220 0.70224
[5,] 0.708960 0.72064 0.747250 0.689680 0.69233 0.694630 0.713220 0.70224
        [,9]    [,10]    [,11]   [,12]
[1,] 0.66576 0.669420 0.675110 0.67254
[2,] 0.67440 0.676690 0.677600 0.68137
[3,] 0.68769 0.681465 0.685085 0.69379
[4,] 0.70045 0.699190 0.696930 0.69763
[5,] 0.70045 0.699190 0.696930 0.69763

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

$conf
         [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.670155 0.6683396 0.6663212 0.6670924 0.6781322 0.6761125 0.6630943
[2,] 0.701955 0.7058804 0.7118088 0.6927776 0.6937678 0.6958375 0.7041957
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,] 0.6596669 0.6708869 0.6669518 0.6726165 0.6833018
[2,] 0.7026131 0.7044931 0.6959782 0.6975535 0.7042782

$out
[1] 0.75094 0.77494 0.79487 0.79209 0.79152 0.79308 0.79279 0.79924 0.78668

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

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

Warning message:
In bxp(list(stats = c(0.63709, 0.67397, 0.686055, 0.69862, 0.70896,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5s2211244128797.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] [,7]
[1,] 0.637090 0.664280 0.668950 0.672540 0.663410 0.708960   NA
[2,] 0.663440 0.667590 0.682150 0.677955 0.673415 0.749095   NA
[3,] 0.692975 0.674235 0.686415 0.683135 0.676640 0.789100   NA
[4,] 0.699820 0.691795 0.693665 0.687240 0.680785 0.792935   NA
[5,] 0.713220 0.701960 0.698670 0.694630 0.688870 0.799240   NA

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

$conf
          [,1]     [,2]      [,3]    [,4]      [,5]      [,6] [,7]
[1,] 0.6763818 0.663195 0.6811629 0.67890 0.6732785 0.7691043   NA
[2,] 0.7095682 0.685275 0.6916671 0.68737 0.6800015 0.8090957   NA

$out
[1] 0.69614

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(0.63709, 0.66344, 0.692975, 0.69982, 0.71322,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6hzv31244128797.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,] 0.02532798 0.071870 0.01079000
[2,] 0.03159688 0.086540 0.01488000
[3,] 0.04335811 0.117330 0.01919000
[4,] 0.04676317 0.125685 0.02491875
[5,] 0.04773334 0.129540 0.02798750

$n
[1] 12 12 12

$conf
           [,1]      [,2]       [,3]
[1,] 0.03644066 0.0994757 0.01461126
[2,] 0.05027555 0.1351843 0.02376874

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.0253279761660237, 0.0315968847307013, 0.0433581071351971,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1ibg81244128797.ps tmp/1ibg81244128797.png")
> system("convert tmp/2v9rz1244128797.ps tmp/2v9rz1244128797.png")
> system("convert tmp/39ft51244128797.ps tmp/39ft51244128797.png")
> system("convert tmp/4z1gm1244128797.ps tmp/4z1gm1244128797.png")
> system("convert tmp/5s2211244128797.ps tmp/5s2211244128797.png")
> system("convert tmp/6hzv31244128797.ps tmp/6hzv31244128797.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.031   0.778   1.468