R version 2.10.1 (2009-12-14)
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(19097.1,19304.6,19601.7,16006.9,17681.2,19790.4,17014.2,17424.5,18908.9,15692.1,15160,15794.3,16032.1,16065,16236.8,12521,14762.1,15446.9,13635,14212.6,15021.7,14134.3,13721.4,14384.5,15638.6,19711.6,20359.8,16141.4,20056.9,20605.5,19325.8,20547.7,19211.2,19009.5,18746.8,16471.5,18957.2,20515.2,18374.4,16192.9,18147.5,19301.4,18344.7,17183.6,19630,17167.2,17428.5,16016.5,18466.5,18406.6,18174.1,14851.9,16260.7,18329.6,18003.8,15903.8,19554.2,16554.2,16198.9,16571.8,17535.2,16198.1,17487.5,13768,14915.8,17160.9,15607.4,16181.5,17413.2,15116.3,14544.5,15050.6,15535.4,15919.3,15853.1,12336.4,14355.5,16040.8,13867.7,14656.6)
> 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] 80
> (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 7 7 6 6 6 6
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]
 [1,] 19097.1 16032.1 15638.6 18957.2 18466.5 17535.2 15535.4
 [2,] 19304.6 16065.0 19711.6 20515.2 18406.6 16198.1 15919.3
 [3,] 19601.7 16236.8 20359.8 18374.4 18174.1 17487.5 15853.1
 [4,] 16006.9 12521.0 16141.4 16192.9 14851.9 13768.0 12336.4
 [5,] 17681.2 14762.1 20056.9 18147.5 16260.7 14915.8 14355.5
 [6,] 19790.4 15446.9 20605.5 19301.4 18329.6 17160.9 16040.8
 [7,] 17014.2 13635.0 19325.8 18344.7 18003.8 15607.4 13867.7
 [8,] 17424.5 14212.6 20547.7 17183.6 15903.8 16181.5 14656.6
 [9,] 18908.9 15021.7 19211.2 19630.0 19554.2 17413.2      NA
[10,] 15692.1 14134.3 19009.5 17167.2 16554.2 15116.3      NA
[11,] 15160.0 13721.4 18746.8 17428.5 16198.9 14544.5      NA
[12,] 15794.3 14384.5 16471.5 16016.5 16571.8 15050.6      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/yougetitorg/rcomp/tmp/1pjvx1305874874.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/yougetitorg/rcomp/tmp/2zfg11305874874.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/yougetitorg/rcomp/tmp/33j881305874874.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/yougetitorg/rcomp/tmp/4k3pb1305874874.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]
[1,] 15535.40 15919.30 15853.10 12336.40 14355.50 15446.90 13635.00 14212.60
[2,] 15835.35 16131.55 16862.15 13144.50 14838.95 16600.85 14737.55 15280.20
[3,] 17535.20 18406.60 18174.10 14851.90 16260.70 18329.60 17014.20 16181.50
[4,] 18711.85 19508.10 18988.05 16074.15 17914.35 19545.90 18174.25 17304.05
[5,] 19097.10 20515.20 20359.80 16192.90 20056.90 20605.50 19325.80 17424.50
         [,9]    [,10]    [,11]   [,12]
[1,] 15021.70 14134.30 13721.40 14384.5
[2,] 17413.20 15116.30 14544.50 15050.6
[3,] 19060.05 16123.15 15679.45 15905.4
[4,] 19554.20 17167.20 17428.50 16471.5
[5,] 19630.00 19009.50 18746.80 16571.8

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 15817.4 16390.18 16904.55 13102.36 14424.12 16570.86 14961.86 14972.89
[2,] 19253.0 20423.02 19443.65 16601.44 18097.28 20088.34 19066.54 17390.11
         [,9]    [,10]    [,11]    [,12]
[1,] 17679.04 14800.25 13819.18 14988.87
[2,] 20441.06 17446.05 17539.72 16821.93

$out
[1] 20547.7

$group
[1] 8

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

Warning message:
In bxp(list(stats = c(15535.4, 15835.35, 17535.2, 18711.85, 19097.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/5msg51305874874.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,] 15160.00 12521.00 15638.60 16016.5 14851.9 13768.00 12336.4
[2,] 15900.60 13927.85 17609.15 17175.4 16229.8 14983.20 14111.6
[3,] 17552.85 14573.30 19268.50 18246.1 17287.8 15894.45 15096.0
[4,] 19200.85 15739.50 20208.35 19129.3 18368.1 17287.05 15886.2
[5,] 19790.40 16236.80 20605.50 20515.2 19554.2 17535.20 16040.8

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]
[1,] 16047.58 13746.99 18082.99 17354.91 16312.51 14843.65 14104.68
[2,] 19058.12 15399.61 20454.01 19137.29 18263.09 16945.25 16087.32

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(15160, 15900.6, 17552.85, 19200.85, 19790.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/60bbs1305874874.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,] 1574.262 3561.70 1121.225
[2,] 1667.162 4181.60 1888.775
[3,] 1835.552 4741.75 2649.612
[4,] 2030.050 5424.70 3010.225
[5,] 2229.914 6335.10 3436.700

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1670.036 4174.764 2138.112
[2,] 2001.068 5308.736 3161.113

$out
[1]  850.323 2187.300

$group
[1] 1 2

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

Warning message:
In bxp(list(stats = c(1574.26180039740, 1667.16217589222, 1835.55196766708,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1pjvx1305874874.ps tmp/1pjvx1305874874.png",intern=TRUE))
character(0)
> try(system("convert tmp/2zfg11305874874.ps tmp/2zfg11305874874.png",intern=TRUE))
character(0)
> try(system("convert tmp/33j881305874874.ps tmp/33j881305874874.png",intern=TRUE))
character(0)
> try(system("convert tmp/4k3pb1305874874.ps tmp/4k3pb1305874874.png",intern=TRUE))
character(0)
> try(system("convert tmp/5msg51305874874.ps tmp/5msg51305874874.png",intern=TRUE))
character(0)
> try(system("convert tmp/60bbs1305874874.ps tmp/60bbs1305874874.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.310   1.000   1.772