R version 3.3.1 (2016-06-21) -- "Bug in Your Hair"
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(106474.50
+ ,106559.75
+ ,106636.75
+ ,106722.00
+ ,106804.50
+ ,106889.75
+ ,106972.25
+ ,107057.50
+ ,107142.75
+ ,107225.25
+ ,107310.50
+ ,107393.00
+ ,107478.25
+ ,107563.50
+ ,107640.50
+ ,107725.75
+ ,107808.25
+ ,107893.50
+ ,107976.00
+ ,108061.25
+ ,108146.50
+ ,108229.00
+ ,108314.25
+ ,108396.75
+ ,108482.00
+ ,108567.25
+ ,108647.00
+ ,108732.25
+ ,108814.75
+ ,108900.00
+ ,108982.50
+ ,109067.75
+ ,109153.00
+ ,109235.50
+ ,109320.75
+ ,109403.25
+ ,109488.50
+ ,109573.75
+ ,109650.75
+ ,109736.00
+ ,109818.50
+ ,109903.75
+ ,109986.25
+ ,110071.50
+ ,110156.75
+ ,110239.25
+ ,110324.50
+ ,110407.00
+ ,110492.25
+ ,110577.50
+ ,110654.50
+ ,110739.75
+ ,110822.25
+ ,110907.50
+ ,110990.00
+ ,111075.25
+ ,111160.50
+ ,111243.00
+ ,111328.25
+ ,111410.75
+ ,111496.00
+ ,111581.25
+ ,111658.25
+ ,111743.50
+ ,111826.00
+ ,111911.25
+ ,111993.75
+ ,112079.00
+ ,112164.25
+ ,112246.75
+ ,112332.00
+ ,112414.50
+ ,112499.75
+ ,112585.00
+ ,112664.75
+ ,112750.00
+ ,112832.50
+ ,112917.75
+ ,113000.25
+ ,113085.50
+ ,113170.75
+ ,113253.25
+ ,113338.50
+ ,113421.00
+ ,113506.25
+ ,113591.50
+ ,113668.50
+ ,113753.75
+ ,113836.25
+ ,113921.50
+ ,114004.00
+ ,114089.25
+ ,114174.50
+ ,114257.00
+ ,114342.25
+ ,114424.75
+ ,114510.00
+ ,114595.25
+ ,114672.25
+ ,114757.50
+ ,114840.00
+ ,114925.25
+ ,115007.75
+ ,115093.00
+ ,115178.25
+ ,115260.75
+ ,115346.00
+ ,115428.50
+ ,115513.75
+ ,115599.00
+ ,115676.00
+ ,115761.25
+ ,115843.75
+ ,115929.00
+ ,116011.50
+ ,116096.75
+ ,116182.00
+ ,116264.50
+ ,116349.75
+ ,116432.25)
> par1 = '10'
> par1 <- '10'
> #'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] 120
> (np <- floor(n / par1))
[1] 12
> 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] 12 12 12 12 12 12 12 12 12 12
> arr
          [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
 [1,] 106474.5 107310.5 108146.5 108982.5 109818.5 110654.5 111496.0 112332.0
 [2,] 106559.8 107393.0 108229.0 109067.8 109903.8 110739.8 111581.2 112414.5
 [3,] 106636.8 107478.2 108314.2 109153.0 109986.2 110822.2 111658.2 112499.8
 [4,] 106722.0 107563.5 108396.8 109235.5 110071.5 110907.5 111743.5 112585.0
 [5,] 106804.5 107640.5 108482.0 109320.8 110156.8 110990.0 111826.0 112664.8
 [6,] 106889.8 107725.8 108567.2 109403.2 110239.2 111075.2 111911.2 112750.0
 [7,] 106972.2 107808.2 108647.0 109488.5 110324.5 111160.5 111993.8 112832.5
 [8,] 107057.5 107893.5 108732.2 109573.8 110407.0 111243.0 112079.0 112917.8
 [9,] 107142.8 107976.0 108814.8 109650.8 110492.2 111328.2 112164.2 113000.2
[10,] 107225.2 108061.2 108900.0 109736.0 110577.5 111410.8 112246.8 113085.5
          [,9]    [,10]    [,11]    [,12] [,13]
 [1,] 113170.8 114004.0 114840.0 115676.0    NA
 [2,] 113253.2 114089.2 114925.2 115761.2    NA
 [3,] 113338.5 114174.5 115007.8 115843.8    NA
 [4,] 113421.0 114257.0 115093.0 115929.0    NA
 [5,] 113506.2 114342.2 115178.2 116011.5    NA
 [6,] 113591.5 114424.8 115260.8 116096.8    NA
 [7,] 113668.5 114510.0 115346.0 116182.0    NA
 [8,] 113753.8 114595.2 115428.5 116264.5    NA
 [9,] 113836.2 114672.2 115513.8 116349.8    NA
[10,] 113921.5 114757.5 115599.0 116432.2    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/1ssxa1471205903.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/2ddc61471205903.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/3uc6x1471205903.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/4m87y1471205903.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]
[1,] 106474.5 106559.8 106636.8 106722.0 106804.5 106889.8 106972.2 107057.5
[2,] 108564.5 108648.4 108733.6 108816.1 108901.4 108985.2 109067.8 109153.0
[3,] 111075.2 111160.5 111240.2 111325.5 111408.0 111493.2 111577.1 111661.0
[4,] 113587.4 113671.2 113756.5 113839.0 113924.2 114008.1 114089.2 114174.5
[5,] 115676.0 115761.2 115843.8 115929.0 116011.5 116096.8 116182.0 116264.5
         [,9]    [,10]
[1,] 107142.8 107225.2
[2,] 109232.8 109318.0
[3,] 111746.2 111828.8
[4,] 114254.2 114339.5
[5,] 116349.8 116432.2

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

$conf
         [,1]     [,2]     [,3]     [,4]   [,5]     [,6]     [,7]     [,8]
[1,] 108784.3 108869.5 108949.3 109034.5 109117 109202.3 109286.8 109370.7
[2,] 113366.2 113451.5 113531.2 113616.5 113699 113784.2 113867.5 113951.3
         [,9]    [,10]
[1,] 109455.9 109538.4
[2,] 114036.6 114119.1

$out
numeric(0)

$group
numeric(0)

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

> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5cp0k1471205903.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,] 106474.5 107310.5 108146.5 108982.5 109818.5 110654.5 111496.0 112332.0
[2,] 106636.8 107478.2 108314.2 109153.0 109986.2 110822.2 111658.2 112499.8
[3,] 106847.1 107683.1 108524.6 109362.0 110198.0 111032.6 111868.6 112707.4
[4,] 107057.5 107893.5 108732.2 109573.8 110407.0 111243.0 112079.0 112917.8
[5,] 107225.2 108061.2 108900.0 109736.0 110577.5 111410.8 112246.8 113085.5
         [,9]    [,10]    [,11]    [,12] [,13]
[1,] 113170.8 114004.0 114840.0 115676.0    NA
[2,] 113338.5 114174.5 115007.8 115843.8    NA
[3,] 113548.9 114383.5 115219.5 116054.1    NA
[4,] 113753.8 114595.2 115428.5 116264.5    NA
[5,] 113921.5 114757.5 115599.0 116432.2    NA

$n
 [1] 10 10 10 10 10 10 10 10 10 10 10 10  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 106636.9 107475.6 108315.8 109151.8 109987.8 110822.4 111658.4 112498.5
[2,] 107057.3 107890.6 108733.5 109572.2 110408.2 111242.8 112078.8 112916.2
         [,9]    [,10]    [,11]    [,12] [,13]
[1,] 113341.4 114173.3 115009.3 115843.9    NA
[2,] 113756.4 114593.7 115429.7 116264.3    NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(106474.5, 106636.75, 106847.125, 107057.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6pkey1471205903.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,] 3017.188 9207 4603.500
[2,] 3017.674 9207 4603.500
[3,] 3018.107 9207 4604.188
[4,] 3018.159 9207 4604.188
[5,] 3018.644 9207 4604.188

$n
[1] 10 10 10

$conf
         [,1] [,2]     [,3]
[1,] 3017.864 9207 4603.844
[2,] 3018.350 9207 4604.531

$out
[1] 3016.911 9201.500 9201.500 9209.750 4605.562 4605.562 4600.750 4600.750

$group
[1] 1 2 2 2 3 3 3 3

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

Warning message:
In bxp(list(stats = c(3017.1881233816, 3017.67350457407, 3018.10702818079,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1ssxa1471205903.ps tmp/1ssxa1471205903.png",intern=TRUE))
character(0)
> try(system("convert tmp/2ddc61471205903.ps tmp/2ddc61471205903.png",intern=TRUE))
character(0)
> try(system("convert tmp/3uc6x1471205903.ps tmp/3uc6x1471205903.png",intern=TRUE))
character(0)
> try(system("convert tmp/4m87y1471205903.ps tmp/4m87y1471205903.png",intern=TRUE))
character(0)
> try(system("convert tmp/5cp0k1471205903.ps tmp/5cp0k1471205903.png",intern=TRUE))
character(0)
> try(system("convert tmp/6pkey1471205903.ps tmp/6pkey1471205903.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.175   0.149   2.365