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(8.5968,8.5114,8.3884,8.2671,8.2410,8.3177,8.4070,8.3917,8.4145,8.5245,8.6289,8.6622,8.9055,8.9770,9.1264,9.1120,9.0576,9.2106,9.2637,9.3107,9.6744,9.5780,9.4166,9.4359,9.2275,9.1828,9.0594,9.1358,9.2208,9.1137,9.2689,9.2489,9.1679,9.1051,9.0818,9.0961,9.1733,9.1455,9.2265,9.1541,9.1559,9.1182,9.1856,9.2378,9.0682,9.0105,8.9939,9.0228,9.1368,9.1763,9.2346,9.1653,9.1277,9.1430,9.1962,9.1861,9.0920,9.0620,8.9981,8.9819,9.0476,9.0852,9.0884,9.1670,9.1931,9.2628,9.4276,9.3398,9.3342,9.4223,9.5614,9.4316,9.3111,9.3414,9.4017,9.3346,9.3310,9.2349,9.2170,9.2098,9.2665,9.2533,9.1008,9.0377,9.0795,9.1896,9.2992,9.2372,9.2061,9.3290,9.1842,9.3231,9.2835,9.1735,9.2889,9.4319)
> 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,] 8.5968 8.9055 9.2275 9.1733 9.1368 9.0476 9.3111 9.0795   NA
 [2,] 8.5114 8.9770 9.1828 9.1455 9.1763 9.0852 9.3414 9.1896   NA
 [3,] 8.3884 9.1264 9.0594 9.2265 9.2346 9.0884 9.4017 9.2992   NA
 [4,] 8.2671 9.1120 9.1358 9.1541 9.1653 9.1670 9.3346 9.2372   NA
 [5,] 8.2410 9.0576 9.2208 9.1559 9.1277 9.1931 9.3310 9.2061   NA
 [6,] 8.3177 9.2106 9.1137 9.1182 9.1430 9.2628 9.2349 9.3290   NA
 [7,] 8.4070 9.2637 9.2689 9.1856 9.1962 9.4276 9.2170 9.1842   NA
 [8,] 8.3917 9.3107 9.2489 9.2378 9.1861 9.3398 9.2098 9.3231   NA
 [9,] 8.4145 9.6744 9.1679 9.0682 9.0920 9.3342 9.2665 9.2835   NA
[10,] 8.5245 9.5780 9.1051 9.0105 9.0620 9.4223 9.2533 9.1735   NA
[11,] 8.6289 9.4166 9.0818 8.9939 8.9981 9.5614 9.1008 9.2889   NA
[12,] 8.6622 9.4359 9.0961 9.0228 8.9819 9.4316 9.0377 9.4319   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/1ui7n1244367460.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/23a4v1244367460.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/3iper1244367460.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/42nww1244367460.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,] 8.90550 8.9770 9.05940 9.1120 9.05760 9.11370 9.1842 9.18610 9.06820
[2,] 8.97655 9.0311 9.07390 9.1239 9.09265 9.11595 9.1849 9.19795 9.08010
[3,] 9.10815 9.1609 9.17645 9.1597 9.17450 9.17680 9.2066 9.24335 9.21720
[4,] 9.20040 9.1862 9.26690 9.2021 9.21345 9.24885 9.2663 9.31690 9.30885
[5,] 9.31110 9.3414 9.40170 9.2372 9.33100 9.32900 9.2689 9.33980 9.33420
       [,10]   [,11]   [,12]
[1,] 9.01050 8.62890 8.66220
[2,] 9.03625 8.99600 9.00235
[3,] 9.13930 9.09130 9.06690
[4,] 9.33780 9.35275 9.43175
[5,] 9.57800 9.56140 9.43590

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]    [,6]     [,7]     [,8]
[1,] 8.983104 9.074259 9.068637 9.116016 9.10702 9.10256 9.161129 9.176903
[2,] 9.233196 9.247541 9.284263 9.203384 9.24198 9.25104 9.252071 9.309797
         [,9]   [,10]    [,11]    [,12]
[1,] 9.089417 8.97085 8.892014 8.827031
[2,] 9.344983 9.30775 9.290586 9.306769

$out
 [1] 8.5968 8.5114 8.3884 8.2671 9.3346 8.2410 8.3177 8.4070 9.4276 8.3917
[11] 8.4145 9.6744 8.5245

$group
 [1]  1  2  3  4  4  5  6  7  7  8  9  9 10

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

Warning message:
In bxp(list(stats = c(8.9055, 8.97655, 9.10815, 9.2004, 9.3111,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5soga1244367460.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] [,9]
[1,] 8.24100 8.90550 9.05940 8.99390 8.9819 9.04760 9.0377 9.07950   NA
[2,] 8.35305 9.08480 9.10060 9.04550 9.0770 9.12770 9.2134 9.18690   NA
[3,] 8.41075 9.23715 9.15185 9.14980 9.1399 9.29850 9.2599 9.26035   NA
[4,] 8.56065 9.42625 9.22415 9.17945 9.1812 9.42495 9.3328 9.31115   NA
[5,] 8.66220 9.67440 9.26890 9.23780 9.2346 9.56140 9.4017 9.43190   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]    [,7]     [,8]
[1,] 8.316062 9.081412 9.095498 9.088705 9.092374 9.162922 9.20544 9.203679
[2,] 8.505438 9.392888 9.208202 9.210895 9.187426 9.434078 9.31436 9.317021
     [,9]
[1,]   NA
[2,]   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(8.241, 8.35305, 8.41075, 8.56065, 8.6622,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/69e4z1244367460.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.2232233 0.71430 0.0547000
[2,] 0.2828751 0.88125 0.1047625
[3,] 0.3119020 1.01230 0.1479750
[4,] 0.3269927 1.06050 0.2282750
[5,] 0.3566461 1.25990 0.3237750

$n
[1] 12 12 12

$conf
          [,1]      [,2]       [,3]
[1,] 0.2917797 0.9305429 0.09164012
[2,] 0.3320243 1.0940571 0.20430988

$out
[1] 0.4191

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(0.223223347073220, 0.282875108527754, 0.311901988562782,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1ui7n1244367460.ps tmp/1ui7n1244367460.png")
> system("convert tmp/23a4v1244367460.ps tmp/23a4v1244367460.png")
> system("convert tmp/3iper1244367460.ps tmp/3iper1244367460.png")
> system("convert tmp/42nww1244367460.ps tmp/42nww1244367460.png")
> system("convert tmp/5soga1244367460.ps tmp/5soga1244367460.png")
> system("convert tmp/69e4z1244367460.ps tmp/69e4z1244367460.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.057   0.800   1.438