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(108.01,101.21,119.93,94.76,95.26,117.96,115.86,111.44,108.16,108.77,109.45,124.83,115.31,109.49,124.24,92.85,98.42,120.88,111.72,116.1,109.37,111.65,114.29,133.68,114.27,126.49,131,104,108.88,128.48,132.44,128.04,116.35,120.93,118.59,133.1,121.05,127.62,135.44,114.88,114.34,128.85,138.9,129.44,114.96,127.98,127.03,128.75,137.91,128.37,135.9,122.19,113.08,136.2,138,115.24,110.95,99.23,102.39,112.67)
> 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] 60
> (np <- floor(n / par1))
[1] 5
> 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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
        [,1]   [,2]   [,3]   [,4]   [,5] [,6]
 [1,] 108.01 115.31 114.27 121.05 137.91   NA
 [2,] 101.21 109.49 126.49 127.62 128.37   NA
 [3,] 119.93 124.24 131.00 135.44 135.90   NA
 [4,]  94.76  92.85 104.00 114.88 122.19   NA
 [5,]  95.26  98.42 108.88 114.34 113.08   NA
 [6,] 117.96 120.88 128.48 128.85 136.20   NA
 [7,] 115.86 111.72 132.44 138.90 138.00   NA
 [8,] 111.44 116.10 128.04 129.44 115.24   NA
 [9,] 108.16 109.37 116.35 114.96 110.95   NA
[10,] 108.77 111.65 120.93 127.98  99.23   NA
[11,] 109.45 114.29 118.59 127.03 102.39   NA
[12,] 124.83 133.68 133.10 128.75 112.67   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/1wxca1257763013.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/2vzek1257763013.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/3dptl1257763013.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/4svv91257763013.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]  [,10]
[1,] 108.01 101.21 119.93  92.85  95.26 117.96 111.72 111.44 108.16  99.23
[2,] 114.27 109.49 124.24  94.76  98.42 120.88 115.86 115.24 109.37 108.77
[3,] 115.31 126.49 131.00 104.00 108.88 128.48 132.44 116.10 110.95 111.65
[4,] 121.05 127.62 135.44 114.88 113.08 128.85 138.00 128.04 114.96 120.93
[5,] 121.05 128.37 135.90 122.19 114.34 136.20 138.90 129.44 116.35 127.98
      [,11]  [,12]
[1,] 102.39 112.67
[2,] 109.45 124.83
[3,] 114.29 128.75
[4,] 118.59 133.10
[5,] 127.03 133.68

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 5

$conf
         [,1]     [,2]     [,3]      [,4]      [,5]     [,6]     [,7]     [,8]
[1,] 110.5193 113.6794 123.0861  89.78326  98.52128 122.8484 116.7959 107.0556
[2,] 120.1007 139.3006 138.9139 118.21674 119.23872 134.1116 148.0841 125.1444
         [,9]    [,10]    [,11]    [,12]
[1,] 107.0001 103.0578 107.8317 122.9064
[2,] 114.8999 120.2422 120.7483 134.5936

$out
[1] 137.91

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(108.01, 114.27, 115.31, 121.05, 121.05, 101.21,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5ts891257763013.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]
[1,]  94.76  98.420 104.00 114.340  99.230   NA
[2,] 104.61 109.430 115.31 118.005 111.810   NA
[3,] 109.11 113.005 123.71 127.800 118.715   NA
[4,] 116.91 118.490 129.74 129.145 136.050   NA
[5,] 124.83 124.240 133.10 138.900 138.000   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5] [,6]
[1,] 103.4999 108.8727 117.1284 122.7190 107.6590   NA
[2,] 114.7201 117.1373 130.2916 132.8810 129.7710   NA

$out
[1]  92.85 133.68

$group
[1] 2 2

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

Warning message:
In bxp(list(stats = c(94.76, 104.61, 109.11, 116.91, 124.83, 98.42,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6nvh41257763013.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,]  3.551066  8.190  5.590
[2,]  7.674752 18.120  8.120
[3,]  8.984782 22.825 11.680
[4,] 11.937653 27.965 16.395
[5,] 12.738174 29.900 22.140

$n
[1] 12 12 12

$conf
          [,1]     [,2]      [,3]
[1,]  7.040443 18.33463  7.905717
[2,] 10.929120 27.31537 15.454283

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(3.55106603712181, 7.6747515059802, 8.98478161611802,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1wxca1257763013.ps tmp/1wxca1257763013.png")
> system("convert tmp/2vzek1257763013.ps tmp/2vzek1257763013.png")
> system("convert tmp/3dptl1257763013.ps tmp/3dptl1257763013.png")
> system("convert tmp/4svv91257763013.ps tmp/4svv91257763013.png")
> system("convert tmp/5ts891257763013.ps tmp/5ts891257763013.png")
> system("convert tmp/6nvh41257763013.ps tmp/6nvh41257763013.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.044   0.790   1.771