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(27.00,26.88,27.38,26.82,27.00,26.15,25.85,26.17,25.66,25.72,25.42,25.10,26.20,26.39,26.27,26.63,21.10,20.30,20.61,21.05,20.45,20.91,21.22,20.85,21.90,22.71,22.40,22.81,23.96,23.37,23.55,23.01,22.63,22.63,22.00,22.15,22.00,22.00,21.84,22.10,22.37,21.83,21.77,21.89,20.76,20.21,20.19,20.01,19.16,18.50,17.41,18.14,18.60,18.32,18.40,18.16,17.29,16.65,16.36,16.32,17.37,17.30,18.10,19.00,18.38,18.41,18.10,17.87,18.70,18.81,18.88,19.44,18.60,18.80,18.62,18.24,17.84,17.85,17.67,17.99,18.15,18.39,18.07,18.39)
> 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] 84
> (np <- floor(n / par1))
[1] 7
> 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 7 7 7 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7] [,8]
 [1,] 27.00 26.20 21.90 22.00 19.16 17.37 18.60   NA
 [2,] 26.88 26.39 22.71 22.00 18.50 17.30 18.80   NA
 [3,] 27.38 26.27 22.40 21.84 17.41 18.10 18.62   NA
 [4,] 26.82 26.63 22.81 22.10 18.14 19.00 18.24   NA
 [5,] 27.00 21.10 23.96 22.37 18.60 18.38 17.84   NA
 [6,] 26.15 20.30 23.37 21.83 18.32 18.41 17.85   NA
 [7,] 25.85 20.61 23.55 21.77 18.40 18.10 17.67   NA
 [8,] 26.17 21.05 23.01 21.89 18.16 17.87 17.99   NA
 [9,] 25.66 20.45 22.63 20.76 17.29 18.70 18.15   NA
[10,] 25.72 20.91 22.63 20.21 16.65 18.81 18.39   NA
[11,] 25.42 21.22 22.00 20.19 16.36 18.88 18.07   NA
[12,] 25.10 20.85 22.15 20.01 16.32 19.44 18.39   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/1yxat1243877432.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/21ih11243877432.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/3xk8l1243877432.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/4756v1243877432.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,] 17.37 17.30 17.410 18.14 17.840 17.850 17.67 17.870 17.290 16.65 16.360
[2,] 18.88 18.65 18.360 18.62 18.490 18.365 18.25 18.075 18.425 18.60 18.475
[3,] 21.90 22.00 21.840 22.10 21.100 20.300 20.61 21.050 20.450 20.21 20.190
[4,] 24.10 24.55 24.335 24.72 23.165 22.600 22.66 22.450 21.695 21.77 21.610
[5,] 27.00 26.88 27.380 26.82 27.000 26.150 25.85 26.170 25.660 25.72 25.420
      [,12]
[1,] 16.320
[2,] 18.915
[3,] 20.010
[4,] 21.500
[5,] 25.100

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 18.7827 18.47662 18.27183 18.45718 18.30817 17.77093 17.97642 18.43732
[2,] 25.0173 25.52338 25.40817 25.74282 23.89183 22.82907 23.24358 23.66268
         [,9]    [,10]    [,11]    [,12]
[1,] 18.49721 18.31693 18.31783 18.46628
[2,] 22.40279 22.10307 22.06217 21.55372

$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(17.37, 18.88, 21.9, 24.1, 27, 17.3, 18.65,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5dn711243877432.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] [,8]
[1,] 25.10 20.300 21.900 20.010 16.32 17.300 17.670   NA
[2,] 25.69 20.730 22.275 20.485 16.97 17.985 17.920   NA
[3,] 26.16 21.075 22.670 21.835 18.15 18.395 18.195   NA
[4,] 26.94 26.235 23.190 22.000 18.45 18.845 18.495   NA
[5,] 27.38 26.630 23.960 22.370 19.16 19.440 18.800   NA

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

$conf
         [,1]     [,2]     [,3]   [,4]     [,5]     [,6]     [,7] [,8]
[1,] 25.58987 18.56413 22.25266 21.144 17.47496 18.00275 17.93274   NA
[2,] 26.73013 23.58587 23.08734 22.526 18.82504 18.78725 18.45726   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(25.1, 25.69, 26.16, 26.94, 27.38, 20.3, 20.73,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6q4wz1243877432.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,] 2.800998 8.180 2.5850
[2,] 2.980946 8.335 3.2200
[3,] 3.101388 8.920 4.3925
[4,] 3.726346 9.370 5.5600
[5,] 3.968739 9.970 6.1000

$n
[1] 12 12 12

$conf
         [,1]    [,2]    [,3]
[1,] 2.761406 8.44793 3.32521
[2,] 3.441370 9.39207 5.45979

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2.80099812141859, 2.98094591482046, 3.10138797925708,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1yxat1243877432.ps tmp/1yxat1243877432.png")
> system("convert tmp/21ih11243877432.ps tmp/21ih11243877432.png")
> system("convert tmp/3xk8l1243877432.ps tmp/3xk8l1243877432.png")
> system("convert tmp/4756v1243877432.ps tmp/4756v1243877432.png")
> system("convert tmp/5dn711243877432.ps tmp/5dn711243877432.png")
> system("convert tmp/6q4wz1243877432.ps tmp/6q4wz1243877432.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.037   0.812   1.279