R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-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(2341
+ ,2115
+ ,2402
+ ,2180
+ ,2453
+ ,2507
+ ,2679
+ ,2622
+ ,2618
+ ,2648
+ ,2523
+ ,2473
+ ,2513
+ ,2466
+ ,2544
+ ,2537
+ ,2564
+ ,2582
+ ,2716
+ ,2904
+ ,2851
+ ,2932
+ ,2772
+ ,2811
+ ,2935
+ ,2783
+ ,3003
+ ,2995
+ ,3127
+ ,2985
+ ,3287
+ ,3236
+ ,3252
+ ,3228
+ ,2856
+ ,3176
+ ,3362
+ ,3036
+ ,3330
+ ,3251
+ ,3318
+ ,3238
+ ,3597
+ ,3708
+ ,3902
+ ,3745
+ ,3426
+ ,3526
+ ,3483
+ ,3458
+ ,3824
+ ,3696
+ ,3518
+ ,3814
+ ,3996
+ ,4136
+ ,4037
+ ,3915
+ ,3760
+ ,3955
+ ,4160
+ ,4115
+ ,4202
+ ,4018
+ ,4233
+ ,4029
+ ,4401
+ ,4645
+ ,4491
+ ,4379
+ ,4394
+ ,4472
+ ,4614
+ ,4160
+ ,4328
+ ,4202
+ ,4635
+ ,4542
+ ,4920
+ ,4774
+ ,4698
+ ,4916
+ ,4703
+ ,4616
+ ,4873
+ ,4375
+ ,4801
+ ,4427
+ ,4684
+ ,4648
+ ,5225
+ ,5174
+ ,5181
+ ,5266
+ ,4839
+ ,5032
+ ,5221
+ ,4658
+ ,5014
+ ,4980
+ ,4952
+ ,4946
+ ,5365
+ ,5456
+ ,5397
+ ,5436
+ ,4995
+ ,5019
+ ,5249
+ ,4799
+ ,5137
+ ,4979
+ ,4951
+ ,5265
+ ,5612
+ ,5572
+ ,5403
+ ,5373
+ ,5252
+ ,5437
+ ,5296
+ ,5011
+ ,5294
+ ,5335
+ ,5398
+ ,5396
+ ,5724
+ ,5898
+ ,5718
+ ,5625
+ ,5380
+ ,5488
+ ,5678
+ ,5224
+ ,5596
+ ,5184
+ ,5620
+ ,5531
+ ,5816
+ ,6086
+ ,6175
+ ,6112
+ ,5813
+ ,5740
+ ,5821
+ ,5294
+ ,5881
+ ,5589
+ ,5845
+ ,5706
+ ,6355
+ ,6404
+ ,6426
+ ,6375
+ ,5869
+ ,5994
+ ,6105
+ ,5792
+ ,6011
+ ,5968
+ ,6255
+ ,6208
+ ,6897
+ ,6814
+ ,6897
+ ,6596
+ ,6188
+ ,6406
+ ,6548
+ ,5842
+ ,6555
+ ,6424
+ ,6596
+ ,6645
+ ,7203
+ ,7128
+ ,7133
+ ,6778
+ ,6593
+ ,6591
+ ,6120
+ ,5612
+ ,6070
+ ,5983
+ ,6145
+ ,6303
+ ,6588
+ ,6640
+ ,6719
+ ,6575
+ ,6487
+ ,6510
+ ,6365
+ ,5844
+ ,5974
+ ,5880
+ ,6279
+ ,6342
+ ,6598
+ ,6801
+ ,6529
+ ,6369
+ ,6028
+ ,6187
+ ,6164
+ ,5866
+ ,6198
+ ,5898
+ ,6462
+ ,6063
+ ,6496
+ ,6678
+ ,6554
+ ,6513
+ ,6210
+ ,5928
+ ,6268
+ ,5582
+ ,5869
+ ,5764
+ ,6082
+ ,6062
+ ,6810
+ ,6727
+ ,6537
+ ,6175
+ ,6014
+ ,6109)
> par1 = '12'
> par1 <- '12'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa P., (2012), Standard Deviation Plot (v1.0.2) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_sdplot.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 228
> (np <- floor(n / par1))
[1] 19
> 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] 19 19 19 19 19 19 19 19 19 19 19 19
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
 [1,] 2341 2513 2935 3362 3483 4160 4614 4873 5221  5249  5296  5678  5821
 [2,] 2115 2466 2783 3036 3458 4115 4160 4375 4658  4799  5011  5224  5294
 [3,] 2402 2544 3003 3330 3824 4202 4328 4801 5014  5137  5294  5596  5881
 [4,] 2180 2537 2995 3251 3696 4018 4202 4427 4980  4979  5335  5184  5589
 [5,] 2453 2564 3127 3318 3518 4233 4635 4684 4952  4951  5398  5620  5845
 [6,] 2507 2582 2985 3238 3814 4029 4542 4648 4946  5265  5396  5531  5706
 [7,] 2679 2716 3287 3597 3996 4401 4920 5225 5365  5612  5724  5816  6355
 [8,] 2622 2904 3236 3708 4136 4645 4774 5174 5456  5572  5898  6086  6404
 [9,] 2618 2851 3252 3902 4037 4491 4698 5181 5397  5403  5718  6175  6426
[10,] 2648 2932 3228 3745 3915 4379 4916 5266 5436  5373  5625  6112  6375
[11,] 2523 2772 2856 3426 3760 4394 4703 4839 4995  5252  5380  5813  5869
[12,] 2473 2811 3176 3526 3955 4472 4616 5032 5019  5437  5488  5740  5994
      [,14] [,15] [,16] [,17] [,18] [,19] [,20]
 [1,]  6105  6548  6120  6365  6164  6268    NA
 [2,]  5792  5842  5612  5844  5866  5582    NA
 [3,]  6011  6555  6070  5974  6198  5869    NA
 [4,]  5968  6424  5983  5880  5898  5764    NA
 [5,]  6255  6596  6145  6279  6462  6082    NA
 [6,]  6208  6645  6303  6342  6063  6062    NA
 [7,]  6897  7203  6588  6598  6496  6810    NA
 [8,]  6814  7128  6640  6801  6678  6727    NA
 [9,]  6897  7133  6719  6529  6554  6537    NA
[10,]  6596  6778  6575  6369  6513  6175    NA
[11,]  6188  6593  6487  6028  6210  6014    NA
[12,]  6406  6591  6510  6187  5928  6109    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/1fvaz1416344457.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/2ttmo1416344457.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/3u2ej1416344457.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/4lnym1416344457.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] [,10] [,11]
[1,] 2341.0 2115.0 2402.0 2180 2453.0 2507.0 2679.0 2622.0 2618  2648  2523
[2,] 3821.5 3786.5 4013.0 3857 3875.5 3921.5 4198.5 4390.5 4264  4147  4077
[3,] 5249.0 4799.0 5137.0 4980 4952.0 5265.0 5612.0 5572.0 5403  5436  5252
[4,] 6112.5 5597.0 5927.5 5822 6113.5 6062.5 6542.0 6659.0 6533  6372  6021
[5,] 6548.0 5866.0 6555.0 6424 6596.0 6645.0 7203.0 7128.0 7133  6778  6593
      [,12]
[1,] 2473.0
[2,] 4213.5
[3,] 5437.0
[4,] 6051.5
[5,] 6591.0

$n
 [1] 19 19 19 19 19 19 19 19 19 19 19 19

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 4418.566 4142.736 4443.038 4267.733 4140.777 4488.937 4762.535 4749.721
[2,] 6079.434 5455.264 5830.962 5692.267 5763.223 6041.063 6461.465 6394.279
        [,9]    [,10]    [,11]    [,12]
[1,] 4580.54 4629.489 4547.345 4770.768
[2,] 6225.46 6242.511 5956.655 6103.232

$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(2341, 3821.5, 5249, 6112.5, 6548, 2115, 3786.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5vryx1416344457.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]   [,9]  [,10]  [,11]
[1,] 2115.0 2466.0 2783 3036.0 3458.0 4018.0 4160.0 4375.0 4658.0 4799.0 5011.0
[2,] 2371.5 2540.5 2960 3284.5 3607.0 4137.5 4435.0 4666.0 4966.0 5058.0 5315.5
[3,] 2490.0 2649.0 3065 3394.0 3819.0 4306.0 4625.5 4856.0 5016.5 5258.5 5397.0
[4,] 2620.0 2831.0 3232 3652.5 3975.5 4436.5 4738.5 5177.5 5381.0 5420.0 5671.5
[5,] 2679.0 2932.0 3287 3902.0 4136.0 4645.0 4920.0 5266.0 5456.0 5612.0 5898.0
      [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19] [,20]
[1,] 5184.0 5294.0 5792.0 6424.0 5612.0 5844.0 5866.0 5582.0    NA
[2,] 5563.5 5763.5 6058.0 6551.5 6095.0 6001.0 5995.5 5941.5    NA
[3,] 5709.0 5875.0 6231.5 6594.5 6395.0 6310.5 6204.0 6095.5    NA
[4,] 5951.0 6365.0 6705.0 6953.0 6581.5 6449.0 6504.5 6402.5    NA
[5,] 6175.0 6426.0 6897.0 7203.0 6719.0 6801.0 6678.0 6810.0    NA

$n
 [1] 12 12 12 12 12 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,] 2376.657 2516.501 2940.939 3226.153 3650.925 4169.624 4487.072 4622.701
[2,] 2603.343 2781.499 3189.061 3561.847 3987.075 4442.376 4763.928 5089.299
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 4827.216 5093.389 5234.626 5532.259 5600.652 5936.399 6411.373 6173.104
[2,] 5205.784 5423.611 5559.374 5885.741 6149.348 6526.601 6777.627 6616.896
        [,17]    [,18]    [,19] [,20]
[1,] 6106.164 5971.842 5885.235    NA
[2,] 6514.836 6436.158 6305.765    NA

$out
[1] 5842

$group
[1] 15

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

Warning message:
In bxp(list(stats = c(2115, 2371.5, 2490, 2620, 2679, 2466, 2540.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6fuyr1416344457.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,] 1232.589 3751 1810.50
[2,] 1301.762 4124 1929.25
[3,] 1348.973 4148 2183.00
[4,] 1403.613 4375 2268.75
[5,] 1453.642 4524 2343.50

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1302.518 4033.517 2028.152
[2,] 1395.428 4262.483 2337.848

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1232.58918992141, 1301.76179968122, 1348.97313343793,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1fvaz1416344457.ps tmp/1fvaz1416344457.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ttmo1416344457.ps tmp/2ttmo1416344457.png",intern=TRUE))
character(0)
> try(system("convert tmp/3u2ej1416344457.ps tmp/3u2ej1416344457.png",intern=TRUE))
character(0)
> try(system("convert tmp/4lnym1416344457.ps tmp/4lnym1416344457.png",intern=TRUE))
character(0)
> try(system("convert tmp/5vryx1416344457.ps tmp/5vryx1416344457.png",intern=TRUE))
character(0)
> try(system("convert tmp/6fuyr1416344457.ps tmp/6fuyr1416344457.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.166   0.435   2.629