R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(213008
+ ,120982
+ ,176537
+ ,184969
+ ,128552
+ ,52746
+ ,385534
+ ,33170
+ ,149061
+ ,165446
+ ,255558
+ ,180915
+ ,138816
+ ,279314
+ ,185362
+ ,330754
+ ,242201
+ ,248074
+ ,137731
+ ,217015
+ ,238060
+ ,352819
+ ,153935
+ ,132943
+ ,176297
+ ,174415
+ ,227308
+ ,223632
+ ,128603
+ ,238315
+ ,220067
+ ,171645
+ ,269655
+ ,84853
+ ,319589
+ ,215641
+ ,331650
+ ,167548
+ ,114673
+ ,275650
+ ,275861
+ ,155363
+ ,155774
+ ,128649
+ ,129301
+ ,368166
+ ,78800
+ ,183167
+ ,289575
+ ,154368
+ ,168809
+ ,24188
+ ,329267
+ ,65029
+ ,101097
+ ,223355
+ ,256087
+ ,248342
+ ,269536
+ ,214872
+ ,324427
+ ,254411
+ ,184360
+ ,214011
+ ,196553
+ ,178886
+ ,143250
+ ,189237
+ ,187681
+ ,120478
+ ,182192
+ ,73566
+ ,194979
+ ,167488
+ ,147958
+ ,289924
+ ,251351
+ ,185082
+ ,138570
+ ,152450
+ ,123905
+ ,364565
+ ,145790
+ ,201714
+ ,80953
+ ,122774
+ ,130585
+ ,298851
+ ,255899
+ ,148553
+ ,211322
+ ,187632
+ ,149018
+ ,220959
+ ,247964
+ ,162765
+ ,200476
+ ,232138
+ ,290838
+ ,315389
+ ,233269
+ ,232317
+ ,151811
+ ,43287
+ ,157799
+ ,168180
+ ,220195
+ ,252225
+ ,99466
+ ,100750
+ ,230772
+ ,262070
+ ,22938
+ ,152474
+ ,61857
+ ,132487
+ ,317394
+ ,21054
+ ,209641
+ ,31414
+ ,249237
+ ,184528
+ ,140189
+ ,112178
+ ,38214
+ ,158671
+ ,287703
+ ,172558
+ ,340803
+ ,253210
+ ,352257
+ ,158015
+ ,85439
+ ,236202
+ ,366492
+ ,84925
+ ,331866
+ ,131069
+ ,211775
+ ,170863
+ ,157949
+ ,314332
+ ,202656
+ ,173260
+ ,268676
+ ,112739
+ ,200511
+ ,228606
+ ,1
+ ,14688
+ ,98
+ ,455
+ ,0
+ ,0
+ ,184207
+ ,293412
+ ,0
+ ,203
+ ,7199
+ ,46660
+ ,17547
+ ,95227
+ ,969
+ ,156828)
> 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] 164
> (np <- floor(n / par1))
[1] 13
> 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 13 13 13 13
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 213008 138816 176297 331650 289575 324427 194979  80953 200476  99466
 [2,] 120982 279314 174415 167548 154368 254411 167488 122774 232138 100750
 [3,] 176537 185362 227308 114673 168809 184360 147958 130585 290838 230772
 [4,] 184969 330754 223632 275650  24188 214011 289924 298851 315389 262070
 [5,] 128552 242201 128603 275861 329267 196553 251351 255899 233269  22938
 [6,]  52746 248074 238315 155363  65029 178886 185082 148553 232317 152474
 [7,] 385534 137731 220067 155774 101097 143250 138570 211322 151811  61857
 [8,]  33170 217015 171645 128649 223355 189237 152450 187632  43287 132487
 [9,] 149061 238060 269655 129301 256087 187681 123905 149018 157799 317394
[10,] 165446 352819  84853 368166 248342 120478 364565 220959 168180  21054
[11,] 255558 153935 319589  78800 269536 182192 145790 247964 220195 209641
[12,] 180915 132943 215641 183167 214872  73566 201714 162765 252225  31414
       [,11]  [,12]  [,13]  [,14]
 [1,] 249237  85439 268676      0
 [2,] 184528 236202 112739    203
 [3,] 140189 366492 200511   7199
 [4,] 112178  84925 228606  46660
 [5,]  38214 331866      1  17547
 [6,] 158671 131069  14688  95227
 [7,] 287703 211775     98    969
 [8,] 172558 170863    455 156828
 [9,] 340803 157949      0     NA
[10,] 253210 314332      0     NA
[11,] 352257 202656 184207     NA
[12,] 158015 173260 293412     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/wessaorg/rcomp/tmp/18g7y1323126683.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/250pm1323126683.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/30b311323126683.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/40ytr1323126683.ps",horizontal=F,onefile=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,]      0.0    203 114673.0  24188      1  14688.0     98.0  43287.0      0
[2,]  99466.0 120982 140189.0 112178  38214  95227.0 101097.0 128649.0 149018
[3,] 197727.5 167518 180448.5 226119 214911 153918.5 147530.5 163845.5 157949
[4,] 268676.0 232138 227308.0 289924 255899 185082.0 211775.0 187632.0 256087
[5,] 331650.0 279314 290838.0 330754 331866 248074.0 287703.0 223355.0 340803
      [,10]  [,11]  [,12]
[1,]      0  78800  73566
[2,] 120478 182192 158015
[3,] 220959 209641 180915
[4,] 314332 255558 214872
[5,] 368166 352257 293412

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 126274.7 120579.8 143660.5 151061.7 122988.6 115975.2 100794.2 138938.6
[2,] 269180.3 214456.2 217236.5 301176.3 306833.4 191861.8 194266.8 188752.4
       [,9]    [,10]    [,11]    [,12]
[1,] 111030 136009.6 177491.1 155999.5
[2,] 204868 305908.4 241790.9 205830.5

$out
[1] 366492   7199 385534  33170    455  31414

$group
[1]  3  3  7  8  8 12

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

Warning message:
In bxp(list(stats = c(0, 99466, 197727.5, 268676, 331650, 203, 120982,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5bowc1323126683.ps",horizontal=F,onefile=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,]  33170.0 132943.0  84853.0  78800.0  24188.0 120478.0 123905.0  80953.0
[2,] 124767.0 146375.5 173030.0 128975.0 127732.5 161068.0 146874.0 139569.0
[3,] 170991.5 227537.5 217854.0 161661.0 219113.5 186020.5 176285.0 175198.5
[4,] 198988.5 263694.0 232811.5 275755.5 262811.5 205282.0 226532.5 234461.5
[5,] 255558.0 352819.0 319589.0 368166.0 329267.0 254411.0 289924.0 298851.0
         [,9]    [,10]    [,11]  [,12]    [,13]    [,14]
[1,] 151811.0  21054.0  38214.0  84925      0.0      0.0
[2,] 162989.5  46635.5 149102.0 144509     49.5    586.0
[3,] 226166.5 116618.5 178543.0 187958  63713.5  12373.0
[4,] 242747.0 220206.5 270456.5 275267 214558.5  70943.5
[5,] 315389.0 317394.0 352257.0 366492 293412.0 156828.0

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

$conf
         [,1]     [,2]     [,3]      [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 137138.6 174027.7 190587.3  94713.43 157503.1 165854.2 139952.2 131917.4
[2,] 204844.4 281047.3 245120.7 228608.57 280723.9 206186.8 212617.8 218479.6
         [,9]    [,10]    [,11]    [,12]     [,13]     [,14]
[1,] 189788.6  37451.6 123192.4 128318.4 -34125.49 -26929.71
[2,] 262544.4 195785.4 233893.6 247597.6 161552.49  51675.71

$out
[1] 385534 324427  73566 364565  43287

$group
[1] 1 6 6 7 9

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

Warning message:
In bxp(list(stats = c(33170, 124767, 170991.5, 198988.5, 255558,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6fxxi1323126683.ps",horizontal=F,onefile=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,]  68736.52 222900.0  54255.00
[2,]  71661.32 267727.5  75921.75
[3,]  87942.37 319108.0 100105.88
[4,] 102971.16 350048.0 155246.25
[5,] 126056.48 385436.0 193963.50

$n
[1] 12 12 12

$conf
          [,1]     [,2]      [,3]
[1,]  73661.75 281561.1  63925.44
[2,] 102223.00 356654.9 136286.31

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(68736.5191708147, 71661.3171847433, 87942.3746935508,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/18g7y1323126683.ps tmp/18g7y1323126683.png",intern=TRUE))
character(0)
> try(system("convert tmp/250pm1323126683.ps tmp/250pm1323126683.png",intern=TRUE))
character(0)
> try(system("convert tmp/30b311323126683.ps tmp/30b311323126683.png",intern=TRUE))
character(0)
> try(system("convert tmp/40ytr1323126683.ps tmp/40ytr1323126683.png",intern=TRUE))
character(0)
> try(system("convert tmp/5bowc1323126683.ps tmp/5bowc1323126683.png",intern=TRUE))
character(0)
> try(system("convert tmp/6fxxi1323126683.ps tmp/6fxxi1323126683.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.428   0.282   1.723