R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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.60773,0.58933,0.60039,0.61342,0.6348,0.634,0.62915,0.62168,0.61328,0.6089,0.60857,0.62672,0.62291,0.62393,0.61838,0.62012,0.61659,0.6116,0.61573,0.61407,0.62823,0.64405,0.6387,0.63633,0.63059,0.62994,0.63709,0.64217,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.66532,0.67157,0.66428,0.66576,0.66942,0.6813,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)
> 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] 88
> (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 7 7 7 7 7 7 7 7
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]
 [1,] 0.60773 0.62291 0.63059 0.69693 0.68130 0.67760 0.67511 0.68887
 [2,] 0.58933 0.62393 0.62994 0.69763 0.69144 0.68137 0.67254 0.69614
 [3,] 0.60039 0.61838 0.63709 0.69278 0.69862 0.67933 0.67397 0.70896
 [4,] 0.61342 0.62012 0.64217 0.70196 0.69500 0.67922 0.67286 0.72064
 [5,] 0.63480 0.61659 0.65711 0.69215 0.69867 0.68598 0.66341      NA
 [6,] 0.63400 0.61160 0.66977 0.67690 0.68968 0.68297 0.66800      NA
 [7,] 0.62915 0.61573 0.68255 0.67124 0.69233 0.68935 0.68021      NA
 [8,] 0.62168 0.61407 0.68902 0.66532 0.68293 0.69463 0.67934      NA
 [9,] 0.61328 0.62823 0.71322 0.67157 0.68399 0.68330 0.68136      NA
[10,] 0.60890 0.64405 0.70224 0.66428 0.66895 0.68666 0.67562      NA
[11,] 0.60857 0.63870 0.70045 0.66576 0.68756 0.68782 0.67440      NA
[12,] 0.62672 0.63633 0.69919 0.66942 0.68527 0.67669 0.67766      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/1gjg61230375544.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/261711230375544.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/3134r1230375544.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/4hhs81230375544.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.607730 0.589330 0.600390 0.613420 0.616590 0.611600 0.615730 0.614070
[2,] 0.626750 0.626935 0.627735 0.631145 0.645955 0.651000 0.650195 0.643500
[3,] 0.676355 0.676955 0.676650 0.676040 0.663410 0.669770 0.680210 0.679340
[4,] 0.685085 0.693790 0.695700 0.698480 0.689065 0.679935 0.685950 0.685975
[5,] 0.696930 0.697630 0.708960 0.720640 0.698670 0.689680 0.692330 0.694630
         [,9]    [,10]   [,11]    [,12]
[1,] 0.613280 0.644050 0.60857 0.626720
[2,] 0.649900 0.654165 0.65223 0.652875
[3,] 0.681360 0.668950 0.67440 0.676690
[4,] 0.683645 0.681140 0.68769 0.681465
[5,] 0.713220 0.702240 0.70045 0.699190

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.6437682 0.6396088 0.6386838 0.6384257 0.6376654 0.6524905 0.6588577
[2,] 0.7089418 0.7143012 0.7146162 0.7136543 0.6891546 0.6870495 0.7015623
          [,8]     [,9]    [,10]     [,11]     [,12]
[1,] 0.6539746 0.661208 0.652841 0.6532239 0.6596165
[2,] 0.7047054 0.701512 0.685059 0.6955761 0.6937635

$out
[1] 0.6089

$group
[1] 10

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

Warning message:
In bxp(list(stats = c(0.60773, 0.62675, 0.676355, 0.685085, 0.69693,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5s9ph1230375544.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]     [,8]
[1,] 0.589330 0.611600 0.62994 0.664280 0.668950 0.676690 0.668000 0.688870
[2,] 0.608150 0.616160 0.63963 0.667590 0.683460 0.679275 0.672700 0.692505
[3,] 0.613350 0.621515 0.67616 0.674235 0.688620 0.683135 0.674755 0.702550
[4,] 0.627935 0.632280 0.69982 0.694855 0.693665 0.687240 0.678500 0.714800
[5,] 0.634800 0.644050 0.71322 0.701960 0.698670 0.694630 0.681360 0.720640

$n
[1] 12 12 12 12 12 12 12  4

$conf
          [,1]      [,2]     [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.6043259 0.6141626 0.648707 0.6617993 0.6839654 0.6795021 0.6721096
[2,] 0.6223741 0.6288674 0.703613 0.6866707 0.6932746 0.6867679 0.6774004
         [,8]
[1,] 0.684937
[2,] 0.720163

$out
[1] 0.66341

$group
[1] 7

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

Warning message:
In bxp(list(stats = c(0.58933, 0.60815, 0.61335, 0.627935, 0.6348,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6e8i21230375544.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.02628612 0.07247 0.0269750
[2,] 0.03054222 0.07932 0.0313400
[3,] 0.03251189 0.09054 0.0391150
[4,] 0.03702837 0.10358 0.0573025
[5,] 0.04057223 0.10857 0.0641775

$n
[1] 12 12 12

$conf
           [,1]       [,2]       [,3]
[1,] 0.02955351 0.07947485 0.02727333
[2,] 0.03547026 0.10160515 0.05095667

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.0262861236251771, 0.0305422171393469, 0.0325118866451136,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1gjg61230375544.ps tmp/1gjg61230375544.png")
> system("convert tmp/261711230375544.ps tmp/261711230375544.png")
> system("convert tmp/3134r1230375544.ps tmp/3134r1230375544.png")
> system("convert tmp/4hhs81230375544.ps tmp/4hhs81230375544.png")
> system("convert tmp/5s9ph1230375544.ps tmp/5s9ph1230375544.png")
> system("convert tmp/6e8i21230375544.ps tmp/6e8i21230375544.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.089   0.826   1.420