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(1639,1296,1063,1282,1365,1268,1532,1455,1393,1515,1510,1225,1577,1417,1224,1693,1633,1639,1914,1586,1552,2081,1500,1437,1470,1849,1387,1592,1590,1798,1935,1887,2027,2080,1556,1682,1785,1869,1781,2082,2571,1862,1938,1505,1767,1607,1578,1495,1615,1700,1337,1531,1623,1543,1638,1520,1416,1820,1596,1358,1267,1742,1402,1388,1646,1670,1531,1730,1407,1795,1504,1371,1734)
> 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] 73
> (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 6 6 6 6 6 6 6 6 6 6 6
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,] 1639 1577 1470 1785 1615 1267 1734
 [2,] 1296 1417 1849 1869 1700 1742   NA
 [3,] 1063 1224 1387 1781 1337 1402   NA
 [4,] 1282 1693 1592 2082 1531 1388   NA
 [5,] 1365 1633 1590 2571 1623 1646   NA
 [6,] 1268 1639 1798 1862 1543 1670   NA
 [7,] 1532 1914 1935 1938 1638 1531   NA
 [8,] 1455 1586 1887 1505 1520 1730   NA
 [9,] 1393 1552 2027 1767 1416 1407   NA
[10,] 1515 2081 2080 1607 1820 1795   NA
[11,] 1510 1500 1556 1578 1596 1504   NA
[12,] 1225 1437 1682 1495 1358 1371   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/1l10v1243580461.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/2xmjw1243580461.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/344l31243580461.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/4b5sw1243580461.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] [,12]
[1,] 1470.0 1296 1063 1282.0 1590 1268.0 1531 1455 1393 1515.0  1500  1225
[2,] 1523.5 1417 1224 1388.0 1590 1543.0 1532 1505 1407 1607.0  1504  1358
[3,] 1615.0 1721 1362 1561.5 1628 1654.5 1776 1553 1484 1807.5  1533  1404
[4,] 1686.5 1849 1402 1693.0 1646 1798.0 1935 1730 1767 2080.0  1578  1495
[5,] 1785.0 1869 1402 2082.0 1646 1862.0 1938 1887 2027 2081.0  1596  1682

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 1517.659 1442.346 1247.184 1364.765 1591.878 1490.017 1516.052 1407.868
[2,] 1712.341 1999.654 1476.816 1758.235 1664.122 1818.983 2035.948 1698.132
         [,9]    [,10]    [,11]    [,12]
[1,] 1251.788 1502.400 1485.268 1315.631
[2,] 1716.212 2112.600 1580.732 1492.369

$out
[1] 1267 1781 1365 2571

$group
[1] 1 3 5 5

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

Warning message:
In bxp(list(stats = c(1470, 1523.5, 1615, 1686.5, 1785, 1296, 1417,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5o1z51243580461.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,] 1063.0 1224.0 1387 1495.0 1337.0 1267.0 1734
[2,] 1275.0 1468.5 1573 1592.5 1468.0 1395.0 1734
[3,] 1379.0 1581.5 1740 1783.0 1569.5 1517.5 1734
[4,] 1512.5 1666.0 1911 1903.5 1630.5 1700.0 1734
[5,] 1639.0 1914.0 2080 2082.0 1820.0 1795.0 1734

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 1270.675 1491.419 1585.836 1641.151 1495.383 1378.387 1734
[2,] 1487.325 1671.581 1894.164 1924.849 1643.617 1656.613 1734

$out
[1] 2081 2571

$group
[1] 2 4

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

Warning message:
In bxp(list(stats = c(1063, 1275, 1379, 1512.5, 1639, 1224, 1468.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6s8by1243580461.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,] 153.9117 407.0  44.500
[2,] 168.7911 444.5 132.625
[3,] 222.6254 569.5 192.125
[4,] 247.4621 676.0 319.250
[5,] 279.8912 800.0 371.250

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 186.7430 463.9113 107.0041
[2,] 258.5078 675.0887 277.2459

$out
[1]  421.33977   41.54355 1206.00000   96.00000

$group
[1] 1 1 2 2

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

Warning message:
In bxp(list(stats = c(153.911662975877, 168.791126368919, 222.625400417204,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1l10v1243580461.ps tmp/1l10v1243580461.png")
> system("convert tmp/2xmjw1243580461.ps tmp/2xmjw1243580461.png")
> system("convert tmp/344l31243580461.ps tmp/344l31243580461.png")
> system("convert tmp/4b5sw1243580461.ps tmp/4b5sw1243580461.png")
> system("convert tmp/5o1z51243580461.ps tmp/5o1z51243580461.png")
> system("convert tmp/6s8by1243580461.ps tmp/6s8by1243580461.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.022   0.776   1.435