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(26.663
+ ,23.598
+ ,26.931
+ ,24.740
+ ,25.806
+ ,24.364
+ ,24.477
+ ,23.901
+ ,23.175
+ ,23.227
+ ,21.672
+ ,21.870
+ ,21.439
+ ,21.089
+ ,23.709
+ ,21.669
+ ,21.752
+ ,20.761
+ ,23.479
+ ,23.824
+ ,23.105
+ ,23.110
+ ,21.759
+ ,22.073
+ ,21.937
+ ,20.035
+ ,23.590
+ ,21.672
+ ,22.222
+ ,22.123
+ ,23.950
+ ,23.504
+ ,22.238
+ ,23.142
+ ,21.059
+ ,21.573
+ ,21.548
+ ,20.000
+ ,22.424
+ ,20.615
+ ,21.761
+ ,22.874
+ ,24.104
+ ,23.748
+ ,23.262
+ ,22.907
+ ,21.519
+ ,22.025
+ ,22.604
+ ,20.894
+ ,24.677
+ ,23.673
+ ,25.320
+ ,23.583
+ ,24.671
+ ,24.454
+ ,24.122
+ ,24.252
+ ,22.084
+ ,22.991
+ ,23.287
+ ,23.049
+ ,25.076
+ ,24.037
+ ,24.430
+ ,24.667
+ ,26.451
+ ,25.618
+ ,25.014
+ ,25.110
+ ,22.964
+ ,23.981
+ ,23.798
+ ,22.270
+ ,24.775
+ ,22.646
+ ,23.988
+ ,24.737
+ ,26.276
+ ,25.816
+ ,25.210
+ ,25.199
+ ,23.162
+ ,24.707
+ ,24.364
+ ,22.644
+ ,25.565
+ ,24.062
+ ,25.431
+ ,24.635
+ ,27.009
+ ,26.606
+ ,26.268
+ ,26.462
+ ,25.246
+ ,25.180
+ ,24.657
+ ,23.304
+ ,26.982
+ ,26.199
+ ,27.210
+ ,26.122
+ ,26.706
+ ,26.878
+ ,26.152
+ ,26.379
+ ,24.712
+ ,25.688
+ ,24.990
+ ,24.239
+ ,26.721
+ ,23.475
+ ,24.767
+ ,26.219
+ ,28.361
+ ,28.599
+ ,27.914
+ ,27.784
+ ,25.693
+ ,26.881
+ ,26.217
+ ,24.218
+ ,27.914
+ ,26.975
+ ,28.527
+ ,27.139
+ ,28.982
+ ,28.169
+ ,28.056
+ ,29.136
+ ,26.291
+ ,26.987
+ ,26.589
+ ,24.848
+ ,27.543
+ ,26.896
+ ,28.878
+ ,27.390
+ ,28.065
+ ,28.141
+ ,29.048
+ ,28.484
+ ,26.634
+ ,27.735
+ ,27.132
+ ,24.924
+ ,28.963
+ ,26.589
+ ,27.931
+ ,28.009
+ ,29.229
+ ,28.759
+ ,28.405
+ ,27.945
+ ,25.912
+ ,26.619
+ ,26.076
+ ,25.286
+ ,27.660
+ ,25.951
+ ,26.398
+ ,25.565
+ ,28.865
+ ,30.000
+ ,29.261
+ ,29.012
+ ,26.992
+ ,27.897)
> 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] 168
> (np <- floor(n / par1))
[1] 14
> 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] 14 14 14 14 14 14 14 14 14 14 14 14
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 26.663 21.439 21.937 21.548 22.604 23.287 23.798 24.364 24.657 24.990
 [2,] 23.598 21.089 20.035 20.000 20.894 23.049 22.270 22.644 23.304 24.239
 [3,] 26.931 23.709 23.590 22.424 24.677 25.076 24.775 25.565 26.982 26.721
 [4,] 24.740 21.669 21.672 20.615 23.673 24.037 22.646 24.062 26.199 23.475
 [5,] 25.806 21.752 22.222 21.761 25.320 24.430 23.988 25.431 27.210 24.767
 [6,] 24.364 20.761 22.123 22.874 23.583 24.667 24.737 24.635 26.122 26.219
 [7,] 24.477 23.479 23.950 24.104 24.671 26.451 26.276 27.009 26.706 28.361
 [8,] 23.901 23.824 23.504 23.748 24.454 25.618 25.816 26.606 26.878 28.599
 [9,] 23.175 23.105 22.238 23.262 24.122 25.014 25.210 26.268 26.152 27.914
[10,] 23.227 23.110 23.142 22.907 24.252 25.110 25.199 26.462 26.379 27.784
[11,] 21.672 21.759 21.059 21.519 22.084 22.964 23.162 25.246 24.712 25.693
[12,] 21.870 22.073 21.573 22.025 22.991 23.981 24.707 25.180 25.688 26.881
       [,11]  [,12]  [,13]  [,14] [,15]
 [1,] 26.217 26.589 27.132 26.076    NA
 [2,] 24.218 24.848 24.924 25.286    NA
 [3,] 27.914 27.543 28.963 27.660    NA
 [4,] 26.975 26.896 26.589 25.951    NA
 [5,] 28.527 28.878 27.931 26.398    NA
 [6,] 27.139 27.390 28.009 25.565    NA
 [7,] 28.982 28.065 29.229 28.865    NA
 [8,] 28.169 28.141 28.759 30.000    NA
 [9,] 28.056 29.048 28.405 29.261    NA
[10,] 29.136 28.484 27.945 29.012    NA
[11,] 26.291 26.634 25.912 26.992    NA
[12,] 26.987 27.735 26.619 27.897    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/1ckdh1273836612.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/2ckdh1273836612.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/3ckdh1273836612.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/45uu21273836612.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]
[1,] 21.4390 20.0000 22.424 20.6150 21.7520 20.761 23.4790 23.504 22.238 22.907
[2,] 22.6040 21.0890 24.677 22.6460 23.9880 23.583 24.4770 23.901 23.262 23.227
[3,] 24.5105 23.1765 26.143 24.0495 25.3755 24.702 26.5785 26.211 25.681 25.789
[4,] 26.2170 24.2390 27.543 26.1990 27.2100 26.219 28.3610 28.169 28.056 27.945
[5,] 27.1320 25.2860 28.963 26.9750 28.8780 28.009 29.2290 30.000 29.261 29.136
      [,11]   [,12]
[1,] 21.059 21.5730
[2,] 21.759 22.0730
[3,] 23.937 24.9435
[4,] 25.912 26.8810
[5,] 26.992 27.8970

$n
 [1] 14 14 14 14 14 14 14 14 14 14 14 14

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 22.98483 21.84634 24.93277 22.54916 24.01494 23.58889 24.93839 24.40874
[2,] 26.03617 24.50666 27.35323 25.54984 26.73606 25.81511 28.21861 28.01326
         [,9]    [,10]   [,11]    [,12]
[1,] 23.65662 23.79672 22.1833 22.91321
[2,] 27.70538 27.78128 25.6907 26.97379

$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(21.439, 22.604, 24.5105, 26.217, 27.132, 20,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/55uu21273836612.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]
[1,] 21.6720 20.7610 20.0350 20.0000 20.8940 22.9640 22.2700 22.6440 23.3040
[2,] 23.2010 21.5540 21.6225 21.5335 22.7975 23.6340 23.4800 24.4995 25.2000
[3,] 24.1325 21.9160 22.1725 22.2245 23.8975 24.5485 24.7220 25.3385 26.1755
[4,] 25.2730 23.2945 23.3230 23.0845 24.5625 25.0930 25.2045 26.3650 26.7920
[5,] 26.9310 23.8240 23.9500 24.1040 25.3200 26.4510 26.2760 27.0090 27.2100
       [,10]   [,11]   [,12]  [,13]   [,14] [,15]
[1,] 23.4750 24.2180 24.8480 24.924 25.2860    NA
[2,] 24.8785 26.6330 26.7650 26.604 26.0135    NA
[3,] 26.4700 27.5265 27.6390 27.938 27.3260    NA
[4,] 27.8490 28.3480 28.3125 28.582 28.9385    NA
[5,] 28.5990 29.1360 29.0480 29.229 30.0000    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 23.18745 21.12215 21.39689 21.51708 23.09247 23.88304 23.93544 24.48763
[2,] 25.07755 22.70985 22.94811 22.93192 24.70253 25.21396 25.50856 26.18937
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14] [,15]
[1,] 25.44938 25.11514 26.74428 26.93317 27.03582 25.99189    NA
[2,] 26.90162 27.82486 28.30872 28.34483 28.84018 28.66011    NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(21.672, 23.201, 24.1325, 25.273, 26.931, 20.761,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6y3u51273836612.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.801035 5.2860 2.416500
[2,] 2.018570 5.8415 2.879000
[3,] 2.124151 6.3420 3.584250
[4,] 2.320454 6.7810 4.272125
[5,] 2.422331 7.2480 4.543500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1.986460 5.913488 2.948836
[2,] 2.261843 6.770512 4.219664

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1.80103486124744, 2.01857046813265, 2.12415134255889,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1ckdh1273836612.ps tmp/1ckdh1273836612.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ckdh1273836612.ps tmp/2ckdh1273836612.png",intern=TRUE))
character(0)
> try(system("convert tmp/3ckdh1273836612.ps tmp/3ckdh1273836612.png",intern=TRUE))
character(0)
> try(system("convert tmp/45uu21273836612.ps tmp/45uu21273836612.png",intern=TRUE))
character(0)
> try(system("convert tmp/55uu21273836612.ps tmp/55uu21273836612.png",intern=TRUE))
character(0)
> try(system("convert tmp/6y3u51273836612.ps tmp/6y3u51273836612.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.058   0.801   1.402