R version 2.6.2 (2008-02-08)
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.

  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(44.13,44.13,44.14,44.14,44.14,44.14,44.15,44.17,44.19,44.26,44.27,44.29,44.29,44.29,44.29,44.32,44.33,44.34,44.34,44.34,44.37,44.47,44.51,44.51,44.51,44.52,44.7,44.84,44.9,44.91,44.94,44.94,44.95,45.28,45.34,45.34,45.34,45.36,45.44,45.62,45.75,45.77,45.77,45.77,46.09,46.25,46.28,46.29,46.29,46.29,46.3,46.34,46.34,46.35,46.42,46.52,46.59,46.66,46.67,46.72,46.72,46.72,46.76,46.89,47.02,47.02,47.04,47.18,47.22,47.8,47.88,47.91)
> 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,] 44.13 44.29 44.51 45.34 46.29 46.72   NA
 [2,] 44.13 44.29 44.52 45.36 46.29 46.72   NA
 [3,] 44.14 44.29 44.70 45.44 46.30 46.76   NA
 [4,] 44.14 44.32 44.84 45.62 46.34 46.89   NA
 [5,] 44.14 44.33 44.90 45.75 46.34 47.02   NA
 [6,] 44.14 44.34 44.91 45.77 46.35 47.02   NA
 [7,] 44.15 44.34 44.94 45.77 46.42 47.04   NA
 [8,] 44.17 44.34 44.94 45.77 46.52 47.18   NA
 [9,] 44.19 44.37 44.95 46.09 46.59 47.22   NA
[10,] 44.26 44.47 45.28 46.25 46.66 47.80   NA
[11,] 44.27 44.51 45.34 46.28 46.67 47.88   NA
[12,] 44.29 44.51 45.34 46.29 46.72 47.91   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/1eq5s1211213798.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/2j8q01211213798.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/3ay5i1211213798.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/4blsn1211213798.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]
[1,] 44.130 44.13 44.14 44.14 44.140 44.14 44.150 44.170 44.19 44.260 44.27
[2,] 44.290 44.29 44.29 44.32 44.330 44.34 44.340 44.340 44.37 44.470 44.51
[3,] 44.925 44.94 45.07 45.23 45.325 45.34 45.355 45.355 45.52 45.765 45.81
[4,] 46.290 46.29 46.30 46.34 46.340 46.35 46.420 46.520 46.59 46.660 46.67
[5,] 46.720 46.72 46.76 46.89 47.020 47.02 47.040 47.180 47.22 47.800 47.88
      [,12]
[1,] 44.290
[2,] 44.510
[3,] 45.815
[4,] 46.720
[5,] 47.910

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 43.63494 43.64994 43.77349 43.92703 44.02849 44.04349 44.01333 43.94883
[2,] 46.21506 46.23006 46.36651 46.53297 46.62151 46.63651 46.69667 46.76117
         [,9]    [,10]    [,11]    [,12]
[1,] 44.08803 44.35238 44.41673 44.38948
[2,] 46.95197 47.17762 47.20327 47.24052

$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(44.13, 44.29, 44.925, 46.29, 46.72, 44.13,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5jyil1211213798.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,] 44.130 44.290 44.510 45.34 46.290 46.720   NA
[2,] 44.140 44.305 44.770 45.53 46.320 46.825   NA
[3,] 44.145 44.340 44.925 45.77 46.385 47.030   NA
[4,] 44.225 44.420 45.115 46.17 46.625 47.510   NA
[5,] 44.290 44.510 45.340 46.29 46.720 47.910   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 44.10623 44.28755 44.76764 45.47809 46.24589 46.71757   NA
[2,] 44.18377 44.39245 45.08236 46.06191 46.52411 47.34243   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(44.13, 44.14, 44.145, 44.225, 44.29, 44.29,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6aacg1211213798.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,] 1.083705 2.590 1.6925
[2,] 1.103368 2.685 1.7100
[3,] 1.156082 2.885 1.7450
[4,] 1.306790 3.285 1.8700
[5,] 1.391627 3.620 1.9500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1.063299 2.611336 1.672023
[2,] 1.248864 3.158664 1.817977

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.08370506442789, 1.10336768960458, 1.15608150768698,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1eq5s1211213798.ps tmp/1eq5s1211213798.png")
> system("convert tmp/2j8q01211213798.ps tmp/2j8q01211213798.png")
> system("convert tmp/3ay5i1211213798.ps tmp/3ay5i1211213798.png")
> system("convert tmp/4blsn1211213798.ps tmp/4blsn1211213798.png")
> system("convert tmp/5jyil1211213798.ps tmp/5jyil1211213798.png")
> system("convert tmp/6aacg1211213798.ps tmp/6aacg1211213798.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.366   1.346   2.588