R version 2.7.1 (2008-06-23)
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(115.6,118.5,118.1,117.7,117.9,118.5,117,115.6,113.7,112.1,113.5,113.5,112.2,112.7,111.9,111.5,110.8,111.2,110.9,114.1,121.3,121.5,122.5,122.4,122.5,122.8,122.9,123.6,123.4,122.9,121,119,119,114.3,113.4,113.7,113.5,114.4,114.6,111.6,111.8,110.8,112.6,115.3,115.4,115.3,115.5,115,116.3,117,116.9,117.2,117.1,116.7,117.4,126.7,128.6,129.2,131.1,131.1,131.5,133.2,133.7,135.7,134.9,135.7,131.2,127.4,127.5,124.1,124.1,124.1,124.5,123.7,125.9,126,125.5,125.3)
> 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] 78
> (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] 7 7 7 7 7 7 6 6 6 6 6 6
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]
 [1,] 115.6 112.2 122.5 113.5 116.3 131.5 124.5
 [2,] 118.5 112.7 122.8 114.4 117.0 133.2 123.7
 [3,] 118.1 111.9 122.9 114.6 116.9 133.7 125.9
 [4,] 117.7 111.5 123.6 111.6 117.2 135.7 126.0
 [5,] 117.9 110.8 123.4 111.8 117.1 134.9 125.5
 [6,] 118.5 111.2 122.9 110.8 116.7 135.7 125.3
 [7,] 117.0 110.9 121.0 112.6 117.4 131.2    NA
 [8,] 115.6 114.1 119.0 115.3 126.7 127.4    NA
 [9,] 113.7 121.3 119.0 115.4 128.6 127.5    NA
[10,] 112.1 121.5 114.3 115.3 129.2 124.1    NA
[11,] 113.5 122.5 113.4 115.5 131.1 124.1    NA
[12,] 113.5 122.4 113.7 115.0 131.1 124.1    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/freestat/rcomp/tmp/151h11219070876.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/freestat/rcomp/tmp/2e91i1219070876.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/freestat/rcomp/tmp/3v23p1219070876.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/freestat/rcomp/tmp/4eqoi1219070876.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] [,11]
[1,] 112.20 112.70 111.90 111.5 110.80 110.80 110.9 114.1 113.70 112.1 113.4
[2,] 114.55 115.70 115.75 114.4 114.45 113.95 112.6 115.3 115.40 114.3 113.5
[3,] 116.30 118.50 118.10 117.7 117.90 118.50 117.2 117.3 120.15 118.4 119.0
[4,] 123.50 123.25 124.40 124.8 124.45 124.10 121.0 126.7 127.50 124.1 124.1
[5,] 131.50 133.20 133.70 135.7 134.90 135.70 131.2 127.4 128.60 129.2 131.1
     [,12]
[1,] 113.5
[2,] 113.7
[3,] 118.7
[4,] 124.1
[5,] 131.1

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 110.9552 113.9913 112.9344 111.4893 111.9282 112.4386 111.7817 109.9466
[2,] 121.6448 123.0087 123.2656 123.9107 123.8718 124.5614 122.6183 124.6534
         [,9]    [,10]    [,11]    [,12]
[1,] 112.3451 112.0787 112.1627 111.9917
[2,] 127.9549 124.7213 125.8373 125.4083

$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(112.2, 114.55, 116.3, 123.5, 131.5, 112.7,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5hkgd1219070876.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,] 112.1 110.80 113.40 110.8 116.30 124.10 123.7
[2,] 113.6 111.35 116.65 112.2 116.95 125.75 124.5
[3,] 116.3 112.45 121.75 114.5 117.30 131.35 125.4
[4,] 118.0 121.40 122.90 115.3 128.90 134.30 125.9
[5,] 118.5 122.50 123.60 115.5 131.10 135.70 126.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]
[1,] 114.2931 107.8661 118.8993 113.0861 111.8495 127.4503 124.4970
[2,] 118.3069 117.0339 124.6007 115.9139 122.7505 135.2497 126.3030

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(112.1, 113.6, 116.3, 118, 118.5, 110.8, 111.35,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6h2ua1219070876.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,] 5.938153 13.30  7.550
[2,] 6.788005 17.35  8.775
[3,] 7.106628 19.80  9.525
[4,] 7.969213 22.95  9.850
[5,] 8.737249 24.90 10.400

$n
[1] 12 12 12

$conf
         [,1]     [,2]      [,3]
[1,] 6.567871 17.24580  9.034685
[2,] 7.645385 22.35420 10.015315

$out
[1] 6.4

$group
[1] 3

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

Warning message:
In bxp(list(stats = c(5.93815347281179, 6.78800472602348, 7.10662784514566,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/151h11219070876.ps tmp/151h11219070876.png")
> system("convert tmp/2e91i1219070876.ps tmp/2e91i1219070876.png")
> system("convert tmp/3v23p1219070876.ps tmp/3v23p1219070876.png")
> system("convert tmp/4eqoi1219070876.ps tmp/4eqoi1219070876.png")
> system("convert tmp/5hkgd1219070876.ps tmp/5hkgd1219070876.png")
> system("convert tmp/6h2ua1219070876.ps tmp/6h2ua1219070876.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.869   1.204   2.126