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(112.1294759
+ ,123.1474443
+ ,119.3806175
+ ,121.5465429
+ ,124.3716631
+ ,85.47917603
+ ,88.86932018
+ ,128.3268312
+ ,141.0398718
+ ,117.9680574
+ ,101.9590434
+ ,103.0890915
+ ,111.9411345
+ ,120.6990069
+ ,111.9411345
+ ,111.9411345
+ ,117.0263507
+ ,89.99936823
+ ,68.90513797
+ ,113.4478653
+ ,105.9142116
+ ,119.2864468
+ ,114.954596
+ ,106.7617476
+ ,110.2460624
+ ,143.0174559
+ ,125.5958818
+ ,106.7617476
+ ,147.2551361
+ ,79.16974108
+ ,92.63614701
+ ,137.6497276
+ ,141.0398718
+ ,123.4299564
+ ,112.5061585
+ ,127.5734659
+ ,130.5869273
+ ,152.1520109
+ ,131.3402927
+ ,106.1967236
+ ,144.6183573
+ ,88.02178414
+ ,99.60477665
+ ,142.8291145
+ ,116.8380094
+ ,125.0308578
+ ,111.2819398
+ ,115.0487667
+ ,134.2595835
+ ,152.3403523
+ ,107.8917957
+ ,139.3447997
+ ,147.7259894
+ ,89.24600286
+ ,97.43885122
+ ,145.6542347
+ ,144.9008693
+ ,144.4300159
+ ,106.7617476
+ ,130.7752687
+ ,125.0308578
+ ,139.250629
+ ,115.0487667
+ ,121.9232256
+ ,150.7394509
+ ,104.1249688
+ ,96.59131518
+ ,136.8021916
+ ,144.0533333
+ ,133.8829008
+ ,118.3447401
+ ,127.1967832
+ ,127.0084418
+ ,128.3268312
+ ,136.3313383
+ ,124.7483457
+ ,137.1788743
+ ,104.9725049
+ ,92.54197634
+ ,136.4255089
+ ,134.1654128
+ ,97.8155339
+ ,86.42088274
+ ,77.28632767
+ ,86.89173609)
> 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] 85
> (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] 8 7 7 7 7 7 7 7 7 7 7 7
> arr
           [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
 [1,] 112.12948 111.94113 110.24606 130.58693 134.25958 125.03086 127.00844
 [2,] 123.14744 120.69901 143.01746 152.15201 152.34035 139.25063 128.32683
 [3,] 119.38062 111.94113 125.59588 131.34029 107.89180 115.04877 136.33134
 [4,] 121.54654 111.94113 106.76175 106.19672 139.34480 121.92323 124.74835
 [5,] 124.37166 117.02635 147.25514 144.61836 147.72599 150.73945 137.17887
 [6,]  85.47918  89.99937  79.16974  88.02178  89.24600 104.12497 104.97250
 [7,]  88.86932  68.90514  92.63615  99.60478  97.43885  96.59132  92.54198
 [8,] 128.32683 113.44787 137.64973 142.82911 145.65423 136.80219 136.42551
 [9,] 141.03987 105.91421 141.03987 116.83801 144.90087 144.05333 134.16541
[10,] 117.96806 119.28645 123.42996 125.03086 144.43002 133.88290  97.81553
[11,] 101.95904 114.95460 112.50616 111.28194 106.76175 118.34474  86.42088
[12,] 103.08909 106.76175 127.57347 115.04877 130.77527 127.19678  77.28633
          [,8]
 [1,] 86.89174
 [2,]       NA
 [3,]       NA
 [4,]       NA
 [5,]       NA
 [6,]       NA
 [7,]       NA
 [8,]       NA
 [9,]       NA
[10,]       NA
[11,]       NA
[12,]       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/18xj61258128507.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/23jiy1258128507.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/3z5p11258128507.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/4pcem1258128507.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]
[1,]  86.89174 120.6990 107.8918 106.1967 117.0264  79.16974 88.86932 128.3268
[2,] 111.09360 125.7371 113.4950 109.3514 130.7753  86.75048 90.70565 132.3762
[3,] 118.58017 139.2506 119.3806 121.5465 144.6184  89.24600 92.63615 136.8022
[4,] 128.79768 147.5847 128.4681 123.3358 147.4906  97.06217 97.01508 140.2394
[5,] 134.25958 152.3404 136.3313 139.3448 150.7395 104.97250 99.60478 145.6542
         [,9]    [,10]    [,11]     [,12]
[1,] 105.9142 117.9681 101.9590  77.28633
[2,] 125.5017 118.6273 104.3604 104.92542
[3,] 141.0399 123.4300 111.2819 115.04877
[4,] 142.5466 129.4569 113.7304 127.38512
[5,] 144.9009 144.4300 118.3447 130.77527

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 108.6904 126.2036 110.4389 113.1953 134.6363 83.08803 88.86825 132.1064
[2,] 128.4699 152.2977 128.3223 129.8978 154.6005 95.40398 96.40404 141.4980
         [,9]    [,10]    [,11]    [,12]
[1,] 130.8609 116.9627 105.6863 101.6362
[2,] 151.2188 129.8972 116.8775 128.4613

$out
[1]  68.90514 113.44787  97.81553  86.42088

$group
[1]  7  8 10 11

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

Warning message:
In bxp(list(stats = c(86.89173609, 111.09359845, 118.58016685, 128.79768455,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5i37l1258128507.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,]  85.47918 105.9142  79.16974  88.02178  89.2460  96.59132  77.28633
[2,] 102.52407 106.3380 108.50391 108.73933 107.3268 116.69675  95.17876
[3,] 118.67434 111.9411 124.51292 120.93443 136.8022 126.11382 125.87839
[4,] 123.75955 115.9905 139.34480 137.08470 145.2776 138.02641 135.24838
[5,] 141.03987 120.6990 147.25514 152.15201 152.3404 150.73945 137.17887
         [,8]
[1,] 86.89174
[2,] 86.89174
[3,] 86.89174
[4,] 86.89174
[5,] 86.89174

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 108.9887 107.5386 110.4462 108.0059 119.4926 116.3852 107.6024 86.89174
[2,] 128.3600 116.3437 138.5797 133.8629 154.1118 135.8424 144.1544 86.89174

$out
[1] 89.99937 68.90514

$group
[1] 2 2

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

Warning message:
In bxp(list(stats = c(85.47917603, 102.52406745, 118.67433745, 123.7595537,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6ei6l1258128507.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,]  9.561873 25.80276  6.309435
[2,] 10.578622 31.17049  9.840835
[3,] 12.356458 32.67722 14.478741
[4,] 14.836775 42.80057 16.880093
[5,] 18.908438 53.48894 22.459705

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 10.41429 27.37267 11.26809
[2,] 14.29863 37.98178 17.68939

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(9.5618734263619, 10.5786220547616, 12.3564578208352,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/18xj61258128507.ps tmp/18xj61258128507.png")
> system("convert tmp/23jiy1258128507.ps tmp/23jiy1258128507.png")
> system("convert tmp/3z5p11258128507.ps tmp/3z5p11258128507.png")
> system("convert tmp/4pcem1258128507.ps tmp/4pcem1258128507.png")
> system("convert tmp/5i37l1258128507.ps tmp/5i37l1258128507.png")
> system("convert tmp/6ei6l1258128507.ps tmp/6ei6l1258128507.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.077   0.829   1.535