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(49.32,58.53,49.09,24.68,16.71,19.86,38.23,36.11,19.59,14.91,15.74,15.40,13.06,19.07,15.28,15.82,12.77,12.05,11.69,13.85,13.85,10.07,9.17,10.79,13.44,21.17,18.64,13.21,15.54,21.94,23.11,18.64,14.94,16.90,15.46,11.15,13.13,12.48,12.95,12.59,10.58,10.58,12.39,15.53,13.06,10.22,16.33,19.72,21.31,18.84,24.84,15.67,15.57,12.73,13.56,15.54,17.22,12.14,11.07,12.02,11.55,6.92,10.33,8.38,12.11,11.46,12.75,13.32,13.00,11.90,11.79,12.55,11.84,11.25,11.15,10.99,11.70,14.01,17.51,17.27,16.90,15.79,15.45,16.24,16.71,16.77,16.64,17.80,16.87,16.13,15.76,15.66,15.54,15.30,15.05,14.69,14.39,14.18,13.70,13.66,13.27,13.56,13.14,14.19,22.57,23.09,23.31,22.91,22.36,43.06,64.67,64.68,56.90,48.79,45.21,41.40,22.17,25.52,20.28,22.87,27.63,22.95,21.35,18.38,17.15,18.27,19.40,20.52)
> 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] 128
> (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] 11 11 11 11 11 11 11 11 10 10 10 10
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11]
 [1,] 49.32 13.06 13.44 13.13 21.31 11.55 11.84 16.71 14.39 22.36 27.63
 [2,] 58.53 19.07 21.17 12.48 18.84  6.92 11.25 16.77 14.18 43.06 22.95
 [3,] 49.09 15.28 18.64 12.95 24.84 10.33 11.15 16.64 13.70 64.67 21.35
 [4,] 24.68 15.82 13.21 12.59 15.67  8.38 10.99 17.80 13.66 64.68 18.38
 [5,] 16.71 12.77 15.54 10.58 15.57 12.11 11.70 16.87 13.27 56.90 17.15
 [6,] 19.86 12.05 21.94 10.58 12.73 11.46 14.01 16.13 13.56 48.79 18.27
 [7,] 38.23 11.69 23.11 12.39 13.56 12.75 17.51 15.76 13.14 45.21 19.40
 [8,] 36.11 13.85 18.64 15.53 15.54 13.32 17.27 15.66 14.19 41.40 20.52
 [9,] 19.59 13.85 14.94 13.06 17.22 13.00 16.90 15.54 22.57 22.17    NA
[10,] 14.91 10.07 16.90 10.22 12.14 11.90 15.79 15.30 23.09 25.52    NA
[11,] 15.74  9.17 15.46 16.33 11.07 11.79 15.45 15.05 23.31 20.28    NA
[12,] 15.40 10.79 11.15 19.72 12.02 12.55 16.24 14.69 22.91 22.87    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/1u60j1244446214.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/29g181244446214.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/3ckz11244446214.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/447711244446214.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,] 11.550  6.92 10.330  8.38 10.58 10.580 11.690 13.32 13.00 10.070  9.170
[2,] 13.095 13.33 13.325 12.90 12.44 12.390 12.945 14.86 13.85 11.900 11.790
[3,] 14.390 18.84 16.640 15.67 15.54 14.010 15.760 15.66 16.22 15.105 15.455
[4,] 21.835 22.06 23.095 18.09 16.79 19.065 21.255 19.58 19.59 16.900 16.330
[5,] 27.630 22.95 24.840 24.68 17.15 21.940 23.110 20.52 22.57 23.090 20.280
      [,12]
[1,] 10.790
[2,] 12.020
[3,] 15.045
[4,] 19.720
[5,] 22.910

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 10.22637 14.68113 11.98569 13.19755 13.46771 10.83011 11.80122 13.41145
[2,] 18.55363 22.99887 21.29431 18.14245 17.61229 17.18989 19.71878 17.90855
         [,9]   [,10]    [,11]    [,12]
[1,] 13.35207 12.6068 13.18663 11.19777
[2,] 19.08793 17.6032 17.72337 18.89223

$out
 [1] 49.32 58.53 43.06 49.09 64.67 64.68 56.90 48.79 38.23 45.21 36.11 41.40
[13] 25.52 23.31

$group
 [1]  1  2  2  3  3  4  5  6  7  7  8  8 10 11

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

Warning message:
In bxp(list(stats = c(11.55, 13.095, 14.39, 21.835, 27.63, 6.92,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/53l251244446214.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]   [,9]  [,10]
[1,] 14.910  9.170 11.150 10.220 11.070  8.380 10.990 14.690 13.140 20.280
[2,] 16.225 11.240 14.190 11.485 12.435 10.895 11.475 15.420 13.610 22.615
[3,] 22.270 12.915 16.220 12.770 15.555 11.845 14.730 15.945 14.185 42.230
[4,] 43.660 14.565 19.905 14.330 18.030 12.650 16.570 16.740 22.740 52.845
[5,] 58.530 19.070 23.110 16.330 24.840 13.320 17.510 17.800 23.310 64.680
      [,11]
[1,] 17.150
[2,] 18.325
[3,] 19.960
[4,] 22.150
[5,] 27.630

$n
 [1] 12 12 12 12 12 12 12 12 12 12  8

$conf
          [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,]  9.756712 11.39845 13.61335 11.47238 13.00308 11.04453 12.40614 15.34294
[2,] 34.783288 14.43155 18.82665 14.06762 18.10692 12.64547 17.05386 16.54706
         [,9]    [,10]   [,11]
[1,] 10.02075 28.44189 17.8233
[2,] 18.34925 56.01811 22.0967

$out
[1] 19.72  6.92

$group
[1] 4 6

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

Warning message:
In bxp(list(stats = c(14.91, 16.225, 22.27, 43.66, 58.53, 9.17,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6w4u21244446214.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,]  3.530927  9.570 3.57750
[2,]  4.884231 14.795 4.69125
[3,] 10.972229 35.645 5.93250
[4,] 14.148880 48.965 8.52000
[5,] 17.400560 56.300 9.77000

$n
[1] 12 12 12

$conf
          [,1]     [,2]     [,3]
[1,]  6.746561 20.05983 4.186181
[2,] 15.197898 51.23017 7.678819

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(3.53092748281367, 4.88423122131986, 10.9722291409897,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1u60j1244446214.ps tmp/1u60j1244446214.png")
> system("convert tmp/29g181244446214.ps tmp/29g181244446214.png")
> system("convert tmp/3ckz11244446214.ps tmp/3ckz11244446214.png")
> system("convert tmp/447711244446214.ps tmp/447711244446214.png")
> system("convert tmp/53l251244446214.ps tmp/53l251244446214.png")
> system("convert tmp/6w4u21244446214.ps tmp/6w4u21244446214.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.026   0.771   1.361