R version 3.2.4 Revised (2016-03-16 r70336) -- "Very Secure Dishes"
Copyright (C) 2016 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(340.4,343.2,345,346.6,348.7,351.1,352.7,354.8,359.8,364.4,366.2,368.8,369.6,370.6,374.2,378.1,381,383.2,387.3,391.4,395.1,399.1,403,406.3,410.2,413.3,418.4,421.4,422.5,425.5,427.3,430.7,433.2,437.5,439.9,443,445.6,446.2,449.3,453.9,458,461.2,463.7,466,468.3,471.7,474.7,477.3,479.8,482.6,485.6,488.5,492,494.8,498.3,502.1,505.8,511.7,516.6,521.3,526.1,530.4,534.7,538.4,544.6,547.7,551.4,554.3,557.5,560.7,563.8,566.2,567.2,569.3,570.9,573,575.1,578.1,581,584.4,340.4,343.2,345,346.6,348.7,351.1,352.7,354.8,359.8,364.4,366.2,368.8,369.6,370.6,374.2,378.1,381,383.2,387.3,391.4,395.1,399.1,403,406.3,410.2,413.3,418.4,421.4,422.5,425.5,427.3,430.7,433.2,437.5,439.9,443,445.6,446.2,449.3,453.9,458,461.2,463.7,466,468.3,471.7,474.7,477.3,479.8,482.6,485.6,488.5,492,494.8,498.3,502.1,505.8,511.7,516.6,521.3,526.1,530.4,534.7,538.4,544.6,547.7,551.4,554.3,557.5,560.7,563.8,566.2,567.2,569.3,570.9,573,575.1,578.1,581,584.4)
> 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] 160
> (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 13 13 13 13 13 13 13 13
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 340.4 369.6 410.2 445.6 479.8 526.1 567.2 348.7 381.0 422.5 458.0 492.0
 [2,] 343.2 370.6 413.3 446.2 482.6 530.4 569.3 351.1 383.2 425.5 461.2 494.8
 [3,] 345.0 374.2 418.4 449.3 485.6 534.7 570.9 352.7 387.3 427.3 463.7 498.3
 [4,] 346.6 378.1 421.4 453.9 488.5 538.4 573.0 354.8 391.4 430.7 466.0 502.1
 [5,] 348.7 381.0 422.5 458.0 492.0 544.6 575.1 359.8 395.1 433.2 468.3 505.8
 [6,] 351.1 383.2 425.5 461.2 494.8 547.7 578.1 364.4 399.1 437.5 471.7 511.7
 [7,] 352.7 387.3 427.3 463.7 498.3 551.4 581.0 366.2 403.0 439.9 474.7 516.6
 [8,] 354.8 391.4 430.7 466.0 502.1 554.3 584.4 368.8 406.3 443.0 477.3 521.3
 [9,] 359.8 395.1 433.2 468.3 505.8 557.5 340.4 369.6 410.2 445.6 479.8 526.1
[10,] 364.4 399.1 437.5 471.7 511.7 560.7 343.2 370.6 413.3 446.2 482.6 530.4
[11,] 366.2 403.0 439.9 474.7 516.6 563.8 345.0 374.2 418.4 449.3 485.6 534.7
[12,] 368.8 406.3 443.0 477.3 521.3 566.2 346.6 378.1 421.4 453.9 488.5 538.4
      [,13] [,14]
 [1,] 544.6 575.1
 [2,] 547.7 578.1
 [3,] 551.4 581.0
 [4,] 554.3 584.4
 [5,] 557.5    NA
 [6,] 560.7    NA
 [7,] 563.8    NA
 [8,] 566.2    NA
 [9,] 567.2    NA
[10,] 569.3    NA
[11,] 570.9    NA
[12,] 573.0    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/1buxk1458677739.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/220uj1458677739.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/34nxg1458677739.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/4dte11458677739.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] [,12]
[1,] 340.4 343.2 345.0 346.60 348.7 351.1 352.7 354.8 340.4 343.2 345.0 346.6
[2,] 381.0 383.2 387.3 391.40 395.1 399.1 403.0 406.3 395.1 399.1 403.0 406.3
[3,] 451.8 453.7 456.5 459.95 458.0 461.2 463.7 466.0 445.6 446.2 449.3 453.9
[4,] 526.1 530.4 534.7 538.40 505.8 511.7 516.6 521.3 505.8 511.7 516.6 521.3
[5,] 575.1 578.1 581.0 584.40 575.1 578.1 581.0 584.4 567.2 569.3 570.9 573.0

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

$conf
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6]    [,7]     [,8]
[1,] 390.5282 391.5414 394.257 397.8759 409.4898 411.8572 413.919 415.6055
[2,] 513.0718 515.8586 518.743 522.0241 506.5102 510.5428 513.481 516.3945
         [,9]    [,10]   [,11]    [,12]
[1,] 397.0898 396.8572 399.519 403.5055
[2,] 494.1102 495.5428 499.081 504.2945

$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(340.4, 381, 451.8, 526.1, 575.1, 343.2, 383.2,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5h8b71458677739.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,] 340.4 369.60 410.20 445.60 479.80 526.10 340.4 348.70 381.00 422.5 458.00
[2,] 345.8 376.15 419.90 451.60 487.05 536.55 345.8 353.75 389.35 429.0 464.85
[3,] 351.9 385.25 426.40 462.45 496.55 549.55 570.1 365.30 401.05 438.7 473.20
[4,] 362.1 397.10 435.35 470.00 508.75 559.10 576.6 370.10 411.75 445.9 481.20
[5,] 368.8 406.30 443.00 477.30 521.30 566.20 584.4 378.10 421.40 453.9 488.50
      [,12]  [,13]  [,14]
[1,] 492.00 544.60 575.10
[2,] 500.20 552.85 576.60
[3,] 514.15 562.25 579.55
[4,] 528.25 568.25 582.70
[5,] 538.40 573.00 584.40

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 344.4655 375.6946 419.3532 454.0576 486.6525 539.2648 464.8306 357.8427
[2,] 359.3345 394.8054 433.4468 470.8424 506.4475 559.8352 675.3694 372.7573
         [,9]    [,10]    [,11]    [,12]   [,13]   [,14]
[1,] 390.8332 430.9918 465.7427 501.3562 555.226 574.731
[2,] 411.2668 446.4082 480.6573 526.9438 569.274 584.369

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(340.4, 345.8, 351.9, 362.1, 368.8, 369.6,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6hqjy1458677739.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,] 74.68054 225.90 110.700
[2,] 75.04144 226.40 112.600
[3,] 75.53317 227.65 114.300
[4,] 79.78910 234.80 129.850
[5,] 80.14128 237.80 130.775

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 73.36773 223.8187 106.4322
[2,] 77.69861 231.4813 122.1678

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(74.6805367154158, 75.0414380150255, 75.5331717321426,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1buxk1458677739.ps tmp/1buxk1458677739.png",intern=TRUE))
character(0)
> try(system("convert tmp/220uj1458677739.ps tmp/220uj1458677739.png",intern=TRUE))
character(0)
> try(system("convert tmp/34nxg1458677739.ps tmp/34nxg1458677739.png",intern=TRUE))
character(0)
> try(system("convert tmp/4dte11458677739.ps tmp/4dte11458677739.png",intern=TRUE))
character(0)
> try(system("convert tmp/5h8b71458677739.ps tmp/5h8b71458677739.png",intern=TRUE))
character(0)
> try(system("convert tmp/6hqjy1458677739.ps tmp/6hqjy1458677739.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.315   0.412   2.750