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(102.1,102.86,102.99,103.73,105.02,104.43,104.63,104.93,105.87,105.66,106.76,106,107.22,107.33,107.11,108.86,107.72,107.88,108.38,107.72,108.41,109.9,111.45,112.18,113.34,113.46,114.06,115.54,116.39,115.94,116.97,115.94,115.91,116.43,116.26,116.35,117.9,117.7,117.53,117.86,117.65,116.51,115.93,115.31,115)
> 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] 45
> (np <- floor(n / par1))
[1] 3
> 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] 4 4 4 4 4 4 4 4 4 3 3 3
> arr
        [,1]   [,2]   [,3]   [,4]
 [1,] 102.10 107.22 113.34 117.90
 [2,] 102.86 107.33 113.46 117.70
 [3,] 102.99 107.11 114.06 117.53
 [4,] 103.73 108.86 115.54 117.86
 [5,] 105.02 107.72 116.39 117.65
 [6,] 104.43 107.88 115.94 116.51
 [7,] 104.63 108.38 116.97 115.93
 [8,] 104.93 107.72 115.94 115.31
 [9,] 105.87 108.41 115.91 115.00
[10,] 105.66 109.90 116.43     NA
[11,] 106.76 111.45 116.26     NA
[12,] 106.00 112.18 116.35     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/17zcy1257432801.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/2e0pb1257432801.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/3dutd1257432801.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/4qgce1257432801.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]
[1,] 102.10 102.860 102.990 103.730 105.020 104.430 104.630 104.930 105.870
[2,] 104.66 105.095 105.050 106.295 106.370 106.155 106.505 106.325 107.140
[3,] 110.28 110.395 110.585 112.200 112.055 111.910 112.155 111.515 111.705
[4,] 115.62 115.580 115.795 116.700 117.020 116.225 116.450 115.625 115.455
[5,] 117.90 117.700 117.530 117.860 117.650 116.510 116.970 115.940 115.910
       [,10]   [,11]   [,12]
[1,] 105.660 106.760 106.000
[2,] 107.780 109.105 109.090
[3,] 109.900 111.450 112.180
[4,] 113.165 113.855 114.265
[5,] 116.430 116.260 116.350

$n
 [1] 4 4 4 4 4 4 4 4 4 3 3 3

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 101.6216 102.1119 102.0965 103.9801 103.6415 103.9547 104.2984 104.168
[2,] 118.9384 118.6781 119.0736 120.4200 120.4685 119.8653 120.0116 118.862
         [,9]    [,10]    [,11]    [,12]
[1,] 105.1362 104.9877 107.1170 107.4593
[2,] 118.2738 114.8123 115.7830 116.9007

$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(102.1, 104.66, 110.28, 115.62, 117.9, 102.86,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5ce5m1257432801.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]
[1,] 102.100 107.110 113.34 115.00
[2,] 103.360 107.525 114.80 115.93
[3,] 104.780 108.130 115.94 117.53
[4,] 105.765 109.380 116.37 117.70
[5,] 106.760 111.450 116.97 117.90

$n
[1] 12 12 12  9

$conf
         [,1]     [,2]     [,3]     [,4]
[1,] 103.6831 107.2839 115.2239 116.5978
[2,] 105.8769 108.9761 116.6561 118.4622

$out
[1] 112.18

$group
[1] 2

$names
[1] "1" "2" "3" NA 

Warning message:
In bxp(list(stats = c(102.1, 103.36, 104.78, 105.765, 106.76, 107.11,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6jmrz1257432801.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,] 4.750126  9.500 4.75000
[2,] 5.316426 10.560 6.41875
[3,] 5.972052 12.210 8.49250
[4,] 6.493656 14.335 8.79750
[5,] 6.919191 15.800 9.66000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 5.435109 10.48820 7.407536
[2,] 6.508994 13.93180 9.577464

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(4.75012631410997, 5.31642605121679, 5.97205161318814,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/17zcy1257432801.ps tmp/17zcy1257432801.png")
> system("convert tmp/2e0pb1257432801.ps tmp/2e0pb1257432801.png")
> system("convert tmp/3dutd1257432801.ps tmp/3dutd1257432801.png")
> system("convert tmp/4qgce1257432801.ps tmp/4qgce1257432801.png")
> system("convert tmp/5ce5m1257432801.ps tmp/5ce5m1257432801.png")
> system("convert tmp/6jmrz1257432801.ps tmp/6jmrz1257432801.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.634   1.255   1.815