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(47.87,47.87,47.89,47.88,47.91,47.92,47.92,47.91,47.93,48.05,48.03,48.04,48.04,48.06,48.04,48.09,48.12,48.16,48.16,48.16,48.08,48.13,48.16,48.15,48.15,48.15,48.27,48.47,48.51,48.53,48.53,48.53,48.68,48.64,48.67,48.66,48.66,48.67,48.71,48.96,49.01,49.04,49.04,49.04,49.06,49.13,49.19,49.26,49.26,49.26,49.29,49.43,49.43,49.45,49.45,49.46,49.57,49.68,49.71,49.7,49.7,49.8,49.84,50.09,50.2,50.16,50.16,50.29,50.36,51.02,51.03,51.04)
> 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,] 47.87 48.04 48.15 48.66 49.26 49.70   NA
 [2,] 47.87 48.06 48.15 48.67 49.26 49.80   NA
 [3,] 47.89 48.04 48.27 48.71 49.29 49.84   NA
 [4,] 47.88 48.09 48.47 48.96 49.43 50.09   NA
 [5,] 47.91 48.12 48.51 49.01 49.43 50.20   NA
 [6,] 47.92 48.16 48.53 49.04 49.45 50.16   NA
 [7,] 47.92 48.16 48.53 49.04 49.45 50.16   NA
 [8,] 47.91 48.16 48.53 49.04 49.46 50.29   NA
 [9,] 47.93 48.08 48.68 49.06 49.57 50.36   NA
[10,] 48.05 48.13 48.64 49.13 49.68 51.02   NA
[11,] 48.03 48.16 48.67 49.19 49.71 51.03   NA
[12,] 48.04 48.15 48.66 49.26 49.70 51.04   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/1kysr1210958947.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/24e9j1210958947.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/3zjty1210958947.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/4m75w1210958947.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,] 47.870 47.87 47.89 47.880 47.91 47.920 47.920 47.910 47.93 48.050 48.03
[2,] 48.040 48.06 48.04 48.090 48.12 48.160 48.160 48.160 48.08 48.130 48.16
[3,] 48.405 48.41 48.49 48.715 48.76 48.785 48.785 48.785 48.87 48.885 48.93
[4,] 49.260 49.26 49.29 49.430 49.43 49.450 49.450 49.460 49.57 49.680 49.71
[5,] 49.700 49.80 49.84 50.090 50.20 50.160 50.160 50.290 50.36 51.020 51.03
     [,12]
[1,] 48.04
[2,] 48.15
[3,] 48.96
[4,] 49.70
[5,] 51.04

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 47.61806 47.63596 47.68371 47.85066 47.91501 47.95291 47.95291 47.94646
[2,] 49.19194 49.18404 49.29629 49.57934 49.60499 49.61709 49.61709 49.62354
        [,9]   [,10]   [,11]   [,12]
[1,] 47.9089 47.8852 47.9302 47.9602
[2,] 49.8311 49.8848 49.9298 49.9598

$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(47.87, 48.04, 48.405, 49.26, 49.7, 47.87,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5nb7d1210958947.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,] 47.870 48.040 48.15 48.660 49.260 49.700   NA
[2,] 47.885 48.070 48.37 48.835 49.360 49.965   NA
[3,] 47.915 48.125 48.53 49.040 49.450 50.180   NA
[4,] 47.980 48.160 48.65 49.095 49.625 50.690   NA
[5,] 48.050 48.160 48.68 49.260 49.710 51.040   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6] [,7]
[1,] 47.87167 48.08395 48.40229 48.92141 49.32913 49.84932   NA
[2,] 47.95833 48.16605 48.65771 49.15859 49.57087 50.51068   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(47.87, 47.885, 47.915, 47.98, 48.05, 48.04,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/69f7e1210958947.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,] 0.7344567 1.830 1.03000
[2,] 0.8030357 2.080 1.07125
[3,] 0.8524048 2.265 1.10500
[4,] 1.0215143 2.700 1.24875
[5,] 1.1276953 3.000 1.31250

$n
[1] 12 12 12

$conf
          [,1]     [,2]     [,3]
[1,] 0.7527552 1.982214 1.024041
[2,] 0.9520543 2.547786 1.185959

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(0.734456715311848, 0.803035734333046, 0.85240477441256,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1kysr1210958947.ps tmp/1kysr1210958947.png")
> system("convert tmp/24e9j1210958947.ps tmp/24e9j1210958947.png")
> system("convert tmp/3zjty1210958947.ps tmp/3zjty1210958947.png")
> system("convert tmp/4m75w1210958947.ps tmp/4m75w1210958947.png")
> system("convert tmp/5nb7d1210958947.ps tmp/5nb7d1210958947.png")
> system("convert tmp/69f7e1210958947.ps tmp/69f7e1210958947.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.367   1.362   2.582