R version 2.8.0 (2008-10-20)
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(464,675,703,887,1139,1077,1318,1260,1120,963,996,960,530,883,894,1045,1199,1287,1565,1577,1076,918,1008,1063,544,635,804,980,1018,1064,1404,1286,1104,999,996,1015,615,722,832,977,1270,1437,1520,1708,1151,934,1159,1209,699,830,996,1124,1458,1270,1753,2258,1208,1241,1265,1828,809,997,1164,1205,1538,1513,1378,2083,1357,1536,1526,1376)
> 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,]  464  530  544  615  699  809   NA
 [2,]  675  883  635  722  830  997   NA
 [3,]  703  894  804  832  996 1164   NA
 [4,]  887 1045  980  977 1124 1205   NA
 [5,] 1139 1199 1018 1270 1458 1538   NA
 [6,] 1077 1287 1064 1437 1270 1513   NA
 [7,] 1318 1565 1404 1520 1753 1378   NA
 [8,] 1260 1577 1286 1708 2258 2083   NA
 [9,] 1120 1076 1104 1151 1208 1357   NA
[10,]  963  918  999  934 1241 1536   NA
[11,]  996 1008  996 1159 1265 1526   NA
[12,]  960 1063 1015 1209 1828 1376   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/18rv41242847113.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/2ah6w1242847113.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/346y21242847113.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/4i7cs1242847113.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,] 464.0  635  703  887.0 1018.0 1064.0 1318 1260.0 1076.0   918  996.0   960
[2,] 530.0  675  804  977.0 1139.0 1077.0 1378 1286.0 1104.0   934  996.0  1015
[3,] 579.5  776  863 1012.5 1234.5 1278.5 1462 1642.5 1135.5   981 1083.5  1136
[4,] 699.0  883  996 1124.0 1458.0 1437.0 1565 2083.0 1208.0  1241 1265.0  1376
[5,] 809.0  997 1164 1205.0 1538.0 1513.0 1753 2258.0 1357.0  1536 1526.0  1828

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

$conf
         [,1]     [,2]     [,3]      [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 470.4895 641.8333 739.1538  917.6803 1028.735 1046.288 1341.379 1128.409
[2,] 688.5105 910.1667 986.8462 1107.3197 1440.265 1510.712 1582.621 2156.591
         [,9]     [,10]     [,11]     [,12]
[1,] 1068.417  782.9751  909.9863  903.1433
[2,] 1202.583 1179.0249 1257.0137 1368.8567

$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(464, 530, 579.5, 699, 809, 635, 675, 776,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5xv1j1242847113.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,]  464.0  530  635  615.0  699.0  809.0   NA
[2,]  795.0  906  892  883.0 1060.0 1184.5   NA
[3,]  979.5 1054 1007 1155.0 1253.0 1377.0   NA
[4,] 1129.5 1243 1084 1353.5 1605.5 1531.0   NA
[5,] 1318.0 1577 1286 1708.0 2258.0 1538.0   NA

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

$conf
          [,1]     [,2]      [,3]      [,4]     [,5]     [,6] [,7]
[1,]  826.9323  900.292  919.4275  940.4018 1004.194 1218.959   NA
[2,] 1132.0677 1207.708 1094.5725 1369.5982 1501.806 1535.041   NA

$out
[1]  544 1404 2083

$group
[1] 3 3 6

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

Warning message:
In bxp(list(stats = c(464, 795, 979.5, 1129.5, 1318, 530, 906, 1054,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6ropy1242847113.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,] 102.4337 281.0  85.750
[2,] 132.0045 353.5 152.000
[3,] 172.4322 455.0 211.125
[4,] 227.7192 574.0 265.625
[5,] 324.2020 868.0 307.250

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 128.7762 354.4285 159.2999
[2,] 216.0883 555.5715 262.9501

$out
[1] 409.417 998.000 630.500

$group
[1] 1 2 3

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

Warning message:
In bxp(list(stats = c(102.433718406913, 132.004548433985, 172.432247845250,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/18rv41242847113.ps tmp/18rv41242847113.png")
> system("convert tmp/2ah6w1242847113.ps tmp/2ah6w1242847113.png")
> system("convert tmp/346y21242847113.ps tmp/346y21242847113.png")
> system("convert tmp/4i7cs1242847113.ps tmp/4i7cs1242847113.png")
> system("convert tmp/5xv1j1242847113.ps tmp/5xv1j1242847113.png")
> system("convert tmp/6ropy1242847113.ps tmp/6ropy1242847113.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.637   1.257   1.859