R version 2.12.2 (2011-02-25)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: x86_64-redhat-linux-gnu (64-bit)

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.

  Natural language support but running in an English locale

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(190.2,180.7,193.6,192.8,195.5,197.2,196.9,178.9,172.4,156.4,143.7,153.6,168.8,185.8,199.9,205.4,197.5,199.6,200.5,193.7,179.6,169.1,169.8,195.5,194.8,204.5,203.8,204.8,204.9,240,248.3,258.4,254.9,288.3,333.6,346.3,357.5,490.7,468.2,471.2,517.1,609.2,682,614,554.2,406.8,348.6,298.8,313.7,282.1,232.9,239.3,241.9,265.7,276,271.5,254.6,269.9,293.5,306.1,365.4,347.9,352.1,377.9,377.4,372.2,362.5,341.9,354.8,369.2,406.7,454.7)
> 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] 72
> (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] 6 6 6 6 6 6 6 6 6 6 6 6
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6] [,7]
 [1,] 190.2 168.8 194.8 357.5 313.7 365.4   NA
 [2,] 180.7 185.8 204.5 490.7 282.1 347.9   NA
 [3,] 193.6 199.9 203.8 468.2 232.9 352.1   NA
 [4,] 192.8 205.4 204.8 471.2 239.3 377.9   NA
 [5,] 195.5 197.5 204.9 517.1 241.9 377.4   NA
 [6,] 197.2 199.6 240.0 609.2 265.7 372.2   NA
 [7,] 196.9 200.5 248.3 682.0 276.0 362.5   NA
 [8,] 178.9 193.7 258.4 614.0 271.5 341.9   NA
 [9,] 172.4 179.6 254.9 554.2 254.6 354.8   NA
[10,] 156.4 169.1 288.3 406.8 269.9 369.2   NA
[11,] 143.7 169.8 333.6 348.6 293.5 406.7   NA
[12,] 153.6 195.5 346.3 298.8 306.1 454.7   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/wessaorg/rcomp/tmp/1jyg71304866463.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2eqt11304866463.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3xsh51304866463.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4srdc1304866463.ps",horizontal=F,onefile=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,] 168.80 180.7 193.60 192.80 195.5 197.20 196.90 178.90 172.40 156.4 143.70
[2,] 190.20 185.8 199.90 204.80 197.5 199.60 200.50 193.70 179.60 169.1 169.80
[3,] 254.25 243.3 218.35 222.35 223.4 252.85 262.15 264.95 254.75 279.1 313.55
[4,] 357.50 347.9 352.10 377.90 377.4 372.20 362.50 341.90 354.80 369.2 348.60
[5,] 365.40 490.7 468.20 471.20 517.1 609.20 362.50 341.90 554.20 406.8 406.70
      [,12]
[1,] 153.60
[2,] 195.50
[3,] 302.45
[4,] 346.30
[5,] 454.70

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 146.3361 138.7403 120.1761 110.6949 107.3587 141.5174 157.6548 169.3562
[2,] 362.1639 347.8597 316.5239 334.0051 339.4413 364.1826 366.6452 360.5438
         [,9]    [,10]    [,11]    [,12]
[1,] 141.7403 150.0290 198.2182 205.1791
[2,] 367.7597 408.1710 428.8818 399.7209

$out
[1] 682 614

$group
[1] 7 8

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

Warning message:
In bxp(list(stats = c(168.8, 190.2, 254.25, 357.5, 365.4, 180.7,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/55x001304866463.ps",horizontal=F,onefile=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,] 143.70 168.80 194.80 298.80 232.90 341.90   NA
[2,] 164.40 174.70 204.65 382.15 248.25 353.45   NA
[3,] 185.45 194.60 244.15 480.95 270.70 367.30   NA
[4,] 194.55 199.75 273.35 581.70 287.80 377.65   NA
[5,] 197.20 205.40 346.30 682.00 313.70 406.70   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6] [,7]
[1,] 171.6984 183.1745 212.8155 389.9339 252.661 356.2622   NA
[2,] 199.2016 206.0255 275.4845 571.9661 288.739 378.3378   NA

$out
[1] 454.7

$group
[1] 6

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

Warning message:
In bxp(list(stats = c(143.7, 164.4, 185.45, 194.55, 197.2, 168.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/wessaorg/rcomp/tmp/6u1sz1304866463.ps",horizontal=F,onefile=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,]  90.31677 196.60 114.4250
[2,] 106.11052 268.80 124.9250
[3,] 118.29560 305.55 137.0875
[4,] 150.61494 396.90 144.1500
[5,] 183.85062 485.10 155.2000

$n
[1] 12 12 12

$conf
          [,1]     [,2]     [,3]
[1,]  97.99684 247.1227 128.3188
[2,] 138.59437 363.9773 145.8562

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(90.3167684689098, 106.110523949565, 118.295602987503,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1jyg71304866463.ps tmp/1jyg71304866463.png",intern=TRUE))
character(0)
> try(system("convert tmp/2eqt11304866463.ps tmp/2eqt11304866463.png",intern=TRUE))
character(0)
> try(system("convert tmp/3xsh51304866463.ps tmp/3xsh51304866463.png",intern=TRUE))
character(0)
> try(system("convert tmp/4srdc1304866463.ps tmp/4srdc1304866463.png",intern=TRUE))
character(0)
> try(system("convert tmp/55x001304866463.ps tmp/55x001304866463.png",intern=TRUE))
character(0)
> try(system("convert tmp/6u1sz1304866463.ps tmp/6u1sz1304866463.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  4.030   0.420   4.564