R version 2.6.2 (2008-02-08)
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.70291,0.6885,0.67127,0.66502,0.65825,0.65025,0.65779,0.66014,0.64683,0.64587,0.63702,0.62651,0.61834,0.61466,0.61063,0.59802,0.60151,0.62927,0.62304,0.6071,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,0.74725)
> 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] 109
> (np <- floor(n / par1))
[1] 9
> 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] 10  9  9  9  9  9  9  9  9  9  9  9
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]    [,9]
 [1,] 0.70291 0.61834 0.63480 0.61659 0.65711 0.69215 0.69867 0.68598 0.66341
 [2,] 0.68850 0.61466 0.63400 0.61160 0.66977 0.67690 0.68968 0.68297 0.66800
 [3,] 0.67127 0.61063 0.62915 0.61573 0.68255 0.67124 0.69233 0.68935 0.68021
 [4,] 0.66502 0.59802 0.62168 0.61407 0.68902 0.66532 0.68293 0.69463 0.67934
 [5,] 0.65825 0.60151 0.61328 0.62823 0.71322 0.67157 0.68399 0.68330 0.68136
 [6,] 0.65025 0.62927 0.60890 0.64405 0.70224 0.66428 0.66895 0.68666 0.67562
 [7,] 0.65779 0.62304 0.60857 0.63870 0.70045 0.66576 0.68756 0.68782 0.67440
 [8,] 0.66014 0.60710 0.62672 0.63633 0.69919 0.66942 0.68527 0.67669 0.67766
 [9,] 0.64683 0.60773 0.62291 0.63059 0.69693 0.68130 0.67760 0.67511 0.68887
[10,] 0.64587 0.58933 0.62393 0.62994 0.69763 0.69144 0.68137 0.67254 0.69614
[11,] 0.63702 0.60039 0.61838 0.63709 0.69278 0.69862 0.67933 0.67397 0.70896
[12,] 0.62651 0.61342 0.62012 0.64217 0.70196 0.69500 0.67922 0.67286 0.72064
        [,10]
 [1,] 0.74725
 [2,]      NA
 [3,]      NA
 [4,]      NA
 [5,]      NA
 [6,]      NA
 [7,]      NA
 [8,]      NA
 [9,]      NA
[10,]      NA
[11,]      NA
[12,]      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/1l8i11211218145.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/2xxow1211218145.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/3v2os1211218145.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/46vge1211218145.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]    [,9]
[1,] 0.616590 0.61160 0.61063 0.59802 0.60151 0.60890 0.60857 0.60710 0.60773
[2,] 0.634800 0.63400 0.62915 0.62168 0.62823 0.64405 0.63870 0.63633 0.63059
[3,] 0.674695 0.66977 0.67127 0.66532 0.67157 0.66428 0.66576 0.66942 0.67511
[4,] 0.698670 0.68297 0.68255 0.68293 0.68330 0.67562 0.68756 0.67766 0.68130
[5,] 0.747250 0.68968 0.69233 0.69463 0.71322 0.70224 0.70045 0.69919 0.69693
       [,10]   [,11]   [,12]
[1,] 0.58933 0.60039 0.61342
[2,] 0.62994 0.63702 0.62651
[3,] 0.67254 0.67397 0.67286
[4,] 0.69144 0.69278 0.69500
[5,] 0.69763 0.70896 0.72064

$n
 [1] 10  9  9  9  9  9  9  9  9  9  9  9

$conf
         [,1]      [,2]     [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 0.642783 0.6439791 0.643146 0.6330617 0.6425665 0.6476531 0.6400271
[2,] 0.706607 0.6955609 0.699394 0.6975783 0.7005735 0.6809069 0.6914929
          [,8]      [,9]   [,10]     [,11]     [,12]
[1,] 0.6476529 0.6484027 0.64015 0.6446031 0.6367886
[2,] 0.6911871 0.7018173 0.70493 0.7033369 0.7089314

$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(0.61659, 0.6348, 0.674695, 0.69867, 0.74725,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5whqx1211218145.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.626510 0.58933 0.608570 0.611600 0.669770 0.664280 0.668950 0.672540
[2,] 0.646350 0.60095 0.615830 0.616160 0.685785 0.667590 0.679275 0.674540
[3,] 0.658020 0.60918 0.622295 0.630265 0.697280 0.674235 0.683460 0.683135
[4,] 0.668145 0.61650 0.627935 0.637895 0.701205 0.691795 0.688620 0.687240
[5,] 0.688500 0.62927 0.634800 0.644050 0.713220 0.698620 0.698670 0.694630
         [,9]   [,10]
[1,] 0.663410 0.74725
[2,] 0.675010 0.74725
[3,] 0.679775 0.74725
[4,] 0.692505 0.74725
[5,] 0.708960 0.74725

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

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]     [,6]      [,7]
[1,] 0.6480792 0.6020875 0.6167738 0.6203515 0.6902468 0.663195 0.6791977
[2,] 0.6679608 0.6162725 0.6278162 0.6401785 0.7043132 0.685275 0.6877223
          [,8]      [,9]   [,10]
[1,] 0.6773424 0.6717954 0.74725
[2,] 0.6889276 0.6877546 0.74725

$out
[1] 0.70291 0.65711 0.72064

$group
[1] 1 5 9

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

Warning message:
In bxp(list(stats = c(0.62651, 0.64635, 0.65802, 0.668145, 0.6885,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6bm001211218145.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.02900810 0.078080 0.04133000
[2,] 0.03124114 0.090540 0.04891500
[3,] 0.03417960 0.094975 0.05423500
[4,] 0.03838376 0.108435 0.05895625
[5,] 0.04157853 0.130660 0.06849000

$n
[1] 12 12 12

$conf
           [,1]       [,2]       [,3]
[1,] 0.03092180 0.08681297 0.04965512
[2,] 0.03743740 0.10313703 0.05881488

$out
[1] 0.03157

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(0.0290080997309372, 0.0312411390057939, 0.0341796001463451,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1l8i11211218145.ps tmp/1l8i11211218145.png")
> system("convert tmp/2xxow1211218145.ps tmp/2xxow1211218145.png")
> system("convert tmp/3v2os1211218145.ps tmp/3v2os1211218145.png")
> system("convert tmp/46vge1211218145.ps tmp/46vge1211218145.png")
> system("convert tmp/5whqx1211218145.ps tmp/5whqx1211218145.png")
> system("convert tmp/6bm001211218145.ps tmp/6bm001211218145.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.343   0.871   1.601