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(132,131.4,132.7,130.9,126,109.7,68.3,70.6,75.3,74.1,74.9,74,74.2,76,76.2,74.9,74.1,76.5,57.8,59.2,57.3,57.5,60.4,59.9,59.9,60,60.2,65.4,62.4,78.8,65.6,64.4,67.4,65.3,66.7,66.8,69.4,71.7,77.1,81.1,82.1,92.1,77.1,78.2,77.7,77.3,78.5,78.8,78.7,79.8,82.2,84,81.7,77.6,64.3,72.6,73.8,73.8,70.1,70,72.3,72.1,73.3,79.1,77,76.1,66.4,72.7,73.2,70.7,73.6,74.2,72.6,73.6,79.1,79.6,78,85.4,82,91.9,89.4,92.1,93.8,93.6,95.6,99.9,103.7,99.2,93.7,93.5,80.7,91.8,105.8,111.3,110.3,109.4,111.4,111.6,111.8,106.6,104.3,105.5,98.5,108.5,106,101.8,101.3,92.4,88.9,84.9,86.4,90.7,86.8,90.6,88.3,95.4,93.6,91.3,91.3,89.3)
> 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] 120
> (np <- floor(n / par1))
[1] 10
> 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] 10 10 10 10 10 10 10 10 10 10 10 10
> arr
       [,1] [,2] [,3] [,4] [,5] [,6] [,7]  [,8]  [,9] [,10] [,11]
 [1,] 132.0 74.2 59.9 69.4 78.7 72.3 72.6  95.6 111.4  88.9    NA
 [2,] 131.4 76.0 60.0 71.7 79.8 72.1 73.6  99.9 111.6  84.9    NA
 [3,] 132.7 76.2 60.2 77.1 82.2 73.3 79.1 103.7 111.8  86.4    NA
 [4,] 130.9 74.9 65.4 81.1 84.0 79.1 79.6  99.2 106.6  90.7    NA
 [5,] 126.0 74.1 62.4 82.1 81.7 77.0 78.0  93.7 104.3  86.8    NA
 [6,] 109.7 76.5 78.8 92.1 77.6 76.1 85.4  93.5 105.5  90.6    NA
 [7,]  68.3 57.8 65.6 77.1 64.3 66.4 82.0  80.7  98.5  88.3    NA
 [8,]  70.6 59.2 64.4 78.2 72.6 72.7 91.9  91.8 108.5  95.4    NA
 [9,]  75.3 57.3 67.4 77.7 73.8 73.2 89.4 105.8 106.0  93.6    NA
[10,]  74.1 57.5 65.3 77.3 73.8 70.7 92.1 111.3 101.8  91.3    NA
[11,]  74.9 60.4 66.7 78.5 70.1 73.6 93.8 110.3 101.3  91.3    NA
[12,]  74.0 59.9 66.8 78.8 70.0 74.2 93.6 109.4  92.4  89.3    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/1iu681291793805.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/html/rcomp/tmp/2bl6b1291793805.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/html/rcomp/tmp/3bl6b1291793805.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/html/rcomp/tmp/4bl6b1291793805.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] [,12]
[1,]  59.90  60.0  60.20  65.40  62.4  76.1 57.8  59.20  57.3  57.5  60.4  59.9
[2,]  72.30  72.1  76.20  79.10  77.0  77.6 65.6  70.60  73.2  70.7  70.1  70.0
[3,]  76.45  77.9  80.65  82.55  81.9  88.0 72.7  75.45  76.5  75.7  76.7  76.5
[4,]  95.60  99.9 103.70  99.20  93.7  93.5 82.0  91.90  93.6  92.1  93.8  92.4
[5,] 111.40 131.4 132.70 106.60 104.3 109.7 98.5 108.50 106.0 111.3 110.3 109.4

$n
 [1] 10 10 10 10 10 10 10 10 10 10 10 10

$conf
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6]    [,7]     [,8]
[1,] 64.80839 64.01001 66.9099 72.50724 73.55601 80.05573 64.5059 64.80767
[2,] 88.09161 91.78999 94.3901 92.59276 90.24399 95.94427 80.8941 86.09233
         [,9]    [,10]    [,11]    [,12]
[1,] 66.30735 65.00771 64.85854 65.30807
[2,] 86.69265 86.39229 88.54146 87.69193

$out
[1] 132.0 130.9 126.0

$group
[1] 1 4 5

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

Warning message:
In bxp(list(stats = c(59.9, 72.3, 76.45, 95.6, 111.4, 60, 72.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5mu5e1291793805.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]   [,8]   [,9] [,10] [,11]
[1,]  68.30 57.30 59.90 77.10 64.30 70.70 72.60  80.70  92.40 84.90    NA
[2,]  74.05 58.50 61.30 77.10 71.35 72.20 78.55  93.60 101.55 87.55    NA
[3,]  92.50 67.25 65.35 77.95 75.70 73.25 83.70  99.55 105.75 89.95    NA
[4,] 131.15 75.45 66.75 79.95 80.75 75.15 92.00 107.60 109.95 91.30    NA
[5,] 132.70 76.50 67.40 82.10 84.00 79.10 93.80 111.30 111.80 95.40    NA

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

$conf
          [,1]     [,2]     [,3]    [,4]    [,5]     [,6]     [,7]     [,8]
[1,]  66.45631 59.51899 62.86422 76.6501 71.4126 71.90449 77.56536  93.1645
[2,] 118.54369 74.98101 67.83578 79.2499 79.9874 74.59551 89.83464 105.9355
         [,9]   [,10] [,11]
[1,] 101.9187 88.2396    NA
[2,] 109.5813 91.6604    NA

$out
[1] 78.8 69.4 71.7 92.1 66.4

$group
[1] 3 4 4 4 6

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

Warning message:
In bxp(list(stats = c(68.3, 74.05, 92.5, 131.15, 132.7, 57.3, 58.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6mu5e1291793805.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,] 12.02320 33.60 14.7250
[2,] 15.33553 49.00 16.8625
[3,] 16.61216 51.85 20.5250
[4,] 20.19653 68.45 21.8750
[5,] 22.04989 72.50 23.6750

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 14.39503 42.97872 18.23877
[2,] 18.82930 60.72128 22.81123

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(12.0231997959501, 15.335530504781, 16.6121636081030,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1iu681291793805.ps tmp/1iu681291793805.png",intern=TRUE))
character(0)
> try(system("convert tmp/2bl6b1291793805.ps tmp/2bl6b1291793805.png",intern=TRUE))
character(0)
> try(system("convert tmp/3bl6b1291793805.ps tmp/3bl6b1291793805.png",intern=TRUE))
character(0)
> try(system("convert tmp/4bl6b1291793805.ps tmp/4bl6b1291793805.png",intern=TRUE))
character(0)
> try(system("convert tmp/5mu5e1291793805.ps tmp/5mu5e1291793805.png",intern=TRUE))
character(0)
> try(system("convert tmp/6mu5e1291793805.ps tmp/6mu5e1291793805.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.065   0.829   6.464