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(109.87,95.74,123.06,123.39,120.28,115.33,110.4,114.49,132.03,123.16,118.82,128.32,112.24,104.53,132.57,122.52,131.8,124.55,120.96,122.6,145.52,118.57,134.25,136.7,121.37,111.63,134.42,137.65,137.86,119.77,130.69,128.28,147.45,128.42,136.9,143.95,135.64,122.48,136.83,153.04,142.71,123.46,144.37,146.15,147.61,158.51,147.4,165.05,154.64,126.2,157.36,154.15,123.21,113.07,110.45,113.57,122.44,114.93,111.85,126.04,121.34)
> 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] 61
> (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] 6 5 5 5 5 5 5 5 5 5 5 5
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]
 [1,] 109.87 112.24 121.37 135.64 154.64 121.34
 [2,]  95.74 104.53 111.63 122.48 126.20     NA
 [3,] 123.06 132.57 134.42 136.83 157.36     NA
 [4,] 123.39 122.52 137.65 153.04 154.15     NA
 [5,] 120.28 131.80 137.86 142.71 123.21     NA
 [6,] 115.33 124.55 119.77 123.46 113.07     NA
 [7,] 110.40 120.96 130.69 144.37 110.45     NA
 [8,] 114.49 122.60 128.28 146.15 113.57     NA
 [9,] 132.03 145.52 147.45 147.61 122.44     NA
[10,] 123.16 118.57 128.42 158.51 114.93     NA
[11,] 118.82 134.25 136.90 147.40 111.85     NA
[12,] 128.32 136.70 143.95 165.05 126.04     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/1z30a1258147813.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=8.3333333333333) 
> 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/2kttn1258147813.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=8.3333333333333) 
> 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/3v5jy1258147813.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=8.3333333333333) 
> 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/4seed1258147813.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=8.3333333333333) 
> 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,] 109.870  95.74 132.57 122.52 120.28 113.07 110.40 113.57 122.44 114.93
[2,] 112.240 104.53 132.57 123.39 123.21 115.33 110.45 114.49 132.03 118.57
[3,] 121.355 111.63 134.42 137.65 131.80 119.77 120.96 122.60 145.52 123.16
[4,] 135.640 122.48 136.83 153.04 137.86 123.46 130.69 128.28 147.45 128.42
[5,] 154.640 126.20 136.83 154.15 142.71 124.55 144.37 146.15 147.61 128.42
      [,11]  [,12]
[1,] 111.85 126.04
[2,] 118.82 128.32
[3,] 134.25 136.70
[4,] 136.90 143.95
[5,] 147.40 165.05

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

$conf
         [,1]      [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 106.2612  98.94658 131.4099 116.6994 121.4483 114.0254 106.6585 112.8560
[2,] 136.4488 124.31342 137.4301 158.6006 142.1517 125.5146 135.2615 132.3440
         [,9]    [,10]    [,11]    [,12]
[1,] 134.6243 116.2000 121.4747 125.6559
[2,] 156.4157 130.1200 147.0253 147.7441

$out
[1] 123.06 157.36 158.51

$group
[1]  3  3 10

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

Warning message:
In bxp(list(stats = c(109.87, 112.24, 121.355, 135.64, 154.64, 95.74,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/500sp1258147813.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=8.3333333333333) 
> 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,] 109.870 104.530 111.630 122.480 110.450 121.34
[2,] 112.445 119.765 124.825 136.235 113.320 121.34
[3,] 119.550 123.575 132.555 145.260 122.825 121.34
[4,] 123.275 133.410 137.755 150.325 140.175 121.34
[5,] 132.030 145.520 147.450 165.050 157.360 121.34

$n
[1] 12 12 12 12 12  1

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]   [,6]
[1,] 114.6104 117.3514 126.6575 138.8335 110.5763 121.34
[2,] 124.4896 129.7986 138.4525 151.6865 135.0737 121.34

$out
[1] 95.74

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(109.87, 112.445, 119.55, 123.275, 132.03,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6wlsd1258147813.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=8.3333333333333) 
> 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,]  9.490499 22.430  4.260
[2,] 11.934100 27.815 11.820
[3,] 13.808496 33.275 15.525
[4,] 15.512696 37.280 18.015
[5,] 17.404269 44.770 20.240

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 12.17627 28.95795 12.69942
[2,] 15.44072 37.59205 18.35058

$out
[1]  4.99105 11.48000 29.65000

$group
[1] 1 2 3

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

Warning message:
In bxp(list(stats = c(9.49049893314362, 11.9341000780153, 13.8084956457466,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1z30a1258147813.ps tmp/1z30a1258147813.png")
> system("convert tmp/2kttn1258147813.ps tmp/2kttn1258147813.png")
> system("convert tmp/3v5jy1258147813.ps tmp/3v5jy1258147813.png")
> system("convert tmp/4seed1258147813.ps tmp/4seed1258147813.png")
> system("convert tmp/500sp1258147813.ps tmp/500sp1258147813.png")
> system("convert tmp/6wlsd1258147813.ps tmp/6wlsd1258147813.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.040   0.810   1.388