R version 3.1.0 (2014-04-10) -- "Spring Dance"
Copyright (C) 2014 The R Foundation for Statistical Computing
Platform: i686-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(71.97
+ ,72.32
+ ,74.07
+ ,77.95
+ ,81.75
+ ,80.81
+ ,74.1
+ ,71.37
+ ,75.21
+ ,76.9
+ ,74.44
+ ,74.76
+ ,76.23
+ ,76.97
+ ,78.4
+ ,78.6
+ ,80.08
+ ,81.12
+ ,80.31
+ ,84.59
+ ,81.34
+ ,80.95
+ ,80.48
+ ,75.26
+ ,76.32
+ ,78.92
+ ,80.47
+ ,83.14
+ ,85.42
+ ,81.53
+ ,87.31
+ ,86.01
+ ,85.1
+ ,79.91
+ ,78.6
+ ,78.6
+ ,79.37
+ ,82.89
+ ,84.43
+ ,85.32
+ ,87.71
+ ,84.68
+ ,80.62
+ ,84.79
+ ,85.49
+ ,81.68
+ ,77.69
+ ,78.31
+ ,79.18
+ ,80.91
+ ,83.91
+ ,86.3
+ ,89.76
+ ,85.11
+ ,83.81
+ ,85.36
+ ,85.89
+ ,82.59
+ ,80.87
+ ,80.27
+ ,81.36
+ ,84.81
+ ,90.3
+ ,95.43
+ ,97.59
+ ,97.8
+ ,99.48
+ ,97.52
+ ,104.39
+ ,97.74
+ ,91.37
+ ,92.42
+ ,96.9
+ ,101.58
+ ,105.46
+ ,110.06
+ ,107.9
+ ,102.87
+ ,96.28
+ ,98.59
+ ,103.22
+ ,98.6
+ ,91.79
+ ,93.83
+ ,95.17
+ ,95.19
+ ,99.44
+ ,109.18
+ ,109.15
+ ,109.72
+ ,108.41
+ ,102.96
+ ,107.64
+ ,97.28
+ ,97.25
+ ,91.84
+ ,94.12
+ ,97.86
+ ,98.83
+ ,102.29
+ ,104.49
+ ,102.11
+ ,102.14
+ ,101.28
+ ,101.21
+ ,94.2
+ ,88.47
+ ,88.08
+ ,88.02
+ ,92.95
+ ,97.05
+ ,101.44
+ ,100.34
+ ,99.98
+ ,94.17
+ ,94.54
+ ,95.12
+ ,98.04
+ ,93.72
+ ,93.83
+ ,93.03
+ ,95.81
+ ,99.1
+ ,100.12
+ ,100.67
+ ,103.87
+ ,102.39
+ ,107.21
+ ,105.71
+ ,99.79
+ ,96.12
+ ,96.17
+ ,97.23
+ ,98.08
+ ,99.84
+ ,99.72
+ ,99.92
+ ,102.7
+ ,102.06
+ ,102.36
+ ,102.43
+ ,100.6
+ ,98.4
+ ,98.61
+ ,103.03
+ ,104.7
+ ,107.45
+ ,109.67
+ ,110.54
+ ,112.05
+ ,113.19
+ ,114.2
+ ,112.56
+ ,107.36
+ ,103.93
+ ,103.83
+ ,104.74
+ ,107.5
+ ,109.53
+ ,109.42
+ ,108.6
+ ,110.72
+ ,105.1
+ ,105.19
+ ,102.55
+ ,101.25
+ ,101.56
+ ,101.62
+ ,101.7
+ ,102.94
+ ,104.37
+ ,106.93
+ ,107.82
+ ,110.83
+ ,106.86
+ ,109.46
+ ,108.8
+ ,108.69
+ ,107.77
+ ,108.64
+ ,108.5
+ ,113.84
+ ,114.59
+ ,116.27
+ ,113.63
+ ,112.29
+ ,110.31
+ ,108.47
+ ,110.67
+ ,109.1
+ ,107.02
+ ,108.12
+ ,106.69
+ ,109.87
+ ,110.82
+ ,114.14
+ ,113.31
+ ,115.16
+ ,111.06
+ ,111.13
+ ,115.96
+ ,117.57
+ ,114.69
+ ,119.42
+ ,118.4
+ ,123.32
+ ,123.39
+ ,127.04
+ ,129.35
+ ,127.12
+ ,122.1
+ ,120.22
+ ,121.53
+ ,119.01
+ ,114.27
+ ,114.46)
> 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] 216
> (np <- floor(n / par1))
[1] 18
> 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] 18 18 18 18 18 18 18 18 18 18 18 18
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]   [,6]   [,7]   [,8]   [,9]  [,10]  [,11]
 [1,] 71.97 76.23 76.32 79.37 79.18  81.36  96.90  95.17  94.12  88.02  93.03
 [2,] 72.32 76.97 78.92 82.89 80.91  84.81 101.58  95.19  97.86  92.95  95.81
 [3,] 74.07 78.40 80.47 84.43 83.91  90.30 105.46  99.44  98.83  97.05  99.10
 [4,] 77.95 78.60 83.14 85.32 86.30  95.43 110.06 109.18 102.29 101.44 100.12
 [5,] 81.75 80.08 85.42 87.71 89.76  97.59 107.90 109.15 104.49 100.34 100.67
 [6,] 80.81 81.12 81.53 84.68 85.11  97.80 102.87 109.72 102.11  99.98 103.87
 [7,] 74.10 80.31 87.31 80.62 83.81  99.48  96.28 108.41 102.14  94.17 102.39
 [8,] 71.37 84.59 86.01 84.79 85.36  97.52  98.59 102.96 101.28  94.54 107.21
 [9,] 75.21 81.34 85.10 85.49 85.89 104.39 103.22 107.64 101.21  95.12 105.71
[10,] 76.90 80.95 79.91 81.68 82.59  97.74  98.60  97.28  94.20  98.04  99.79
[11,] 74.44 80.48 78.60 77.69 80.87  91.37  91.79  97.25  88.47  93.72  96.12
[12,] 74.76 75.26 78.60 78.31 80.27  92.42  93.83  91.84  88.08  93.83  96.17
       [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18] [,19]
 [1,]  97.23 103.03 104.74 101.70 108.50 106.69 118.40    NA
 [2,]  98.08 104.70 107.50 102.94 113.84 109.87 123.32    NA
 [3,]  99.84 107.45 109.53 104.37 114.59 110.82 123.39    NA
 [4,]  99.72 109.67 109.42 106.93 116.27 114.14 127.04    NA
 [5,]  99.92 110.54 108.60 107.82 113.63 113.31 129.35    NA
 [6,] 102.70 112.05 110.72 110.83 112.29 115.16 127.12    NA
 [7,] 102.06 113.19 105.10 106.86 110.31 111.06 122.10    NA
 [8,] 102.36 114.20 105.19 109.46 108.47 111.13 120.22    NA
 [9,] 102.43 112.56 102.55 108.80 110.67 115.96 121.53    NA
[10,] 100.60 107.36 101.25 108.69 109.10 117.57 119.01    NA
[11,]  98.40 103.93 101.56 107.77 107.02 114.69 114.27    NA
[12,]  98.61 103.83 101.62 108.64 108.12 119.42 114.46    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/1n0vd1398688050.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/280ec1398688051.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/3zdp81398688051.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/44z8t1398688051.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]
[1,]  71.970  72.320  74.07  77.950  80.08  80.810  74.10  71.37  75.210  76.90
[2,]  79.370  82.890  84.43  86.300  89.76  85.110  87.31  86.01  85.890  82.59
[3,]  94.645  96.835  99.27 101.865 102.58 102.785 102.10 101.82 102.885  98.32
[4,] 103.030 104.700 107.45 109.670 109.15 110.830 108.41 108.47 108.800 107.36
[5,] 118.400 123.320 123.39 127.040 129.35 127.120 122.10 120.22 121.530 119.01
      [,11]  [,12]
[1,]  74.44  74.76
[2,]  80.87  80.27
[3,]  94.92  93.83
[4,] 103.93 103.83
[5,] 114.69 119.42

$n
 [1] 18 18 18 18 18 18 18 18 18 18 18 18

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,]  85.83379  88.71275  90.69713  93.16179  95.35898  93.20663  94.24216
[2,] 103.45621 104.95725 107.84287 110.56821 109.80102 112.36337 109.95784
          [,8]     [,9]     [,10]     [,11]     [,12]
[1,]  93.45568  94.3531  89.09541  86.33224  85.05603
[2,] 110.18432 111.4169 107.54459 103.50776 102.60397

$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(71.97, 79.37, 94.645, 103.03, 118.4, 72.32,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5h50w1398688051.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]
[1,] 71.370 75.260 76.32 77.690 79.180  81.360  91.79  91.84  88.080  88.020
[2,] 73.195 77.685 78.76 79.995 80.890  90.835  96.59  96.22  94.160  93.775
[3,] 74.600 80.195 81.00 83.660 83.860  96.475 100.09 101.20 100.020  94.830
[4,] 77.425 81.035 85.26 85.055 85.625  97.770 104.34 108.78 102.125  99.010
[5,] 81.750 84.590 87.31 87.710 89.760 104.390 110.06 109.72 104.490 101.440
       [,11]   [,12]   [,13]   [,14]   [,15]   [,16]   [,17]   [,18] [,19]
[1,]  93.030  97.230 103.030 101.250 101.700 107.020 106.690 114.270    NA
[2,]  96.145  98.505 104.315 102.085 105.615 108.485 110.940 118.705    NA
[3,]  99.955  99.880 108.560 105.145 107.795 110.490 113.725 121.815    NA
[4,] 103.130 102.210 112.305 109.010 108.745 113.735 115.560 125.215    NA
[5,] 107.210 102.700 114.200 110.720 110.830 116.270 119.420 129.350    NA

$n
 [1] 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,] 72.67067 78.66704 78.03531 81.3521 81.70033 93.3119  96.55517  95.4713
[2,] 76.52933 81.72296 83.96469 85.9679 86.01967 99.6381 103.62483 106.9287
          [,9]    [,10]     [,11]     [,12]    [,13]    [,14]    [,15]    [,16]
[1,]  96.38711 92.44228  96.76909  98.19012 104.9157 101.9865 106.3674 108.0954
[2,] 103.65289 97.21772 103.14091 101.56988 112.2043 108.3035 109.2226 112.8846
        [,17]    [,18] [,19]
[1,] 111.6178 118.8457    NA
[2,] 115.8322 124.7843    NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(71.37, 73.195, 74.6, 77.425, 81.75, 75.26,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6ucug1398688051.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,] 12.53083 40.250 18.99750
[2,] 12.76618 45.485 20.16250
[3,] 13.16892 47.215 20.72875
[4,] 13.55132 49.180 21.05500
[5,] 13.92901 51.000 21.05500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 12.81081 45.52969 20.32167
[2,] 13.52703 48.90031 21.13583

$out
[1] 22.830 17.295 22.520

$group
[1] 3 3 3

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

Warning message:
In bxp(list(stats = c(12.530834700875, 12.7661820808901, 13.1689203826566,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1n0vd1398688050.ps tmp/1n0vd1398688050.png",intern=TRUE))
character(0)
> try(system("convert tmp/280ec1398688051.ps tmp/280ec1398688051.png",intern=TRUE))
character(0)
> try(system("convert tmp/3zdp81398688051.ps tmp/3zdp81398688051.png",intern=TRUE))
character(0)
> try(system("convert tmp/44z8t1398688051.ps tmp/44z8t1398688051.png",intern=TRUE))
character(0)
> try(system("convert tmp/5h50w1398688051.ps tmp/5h50w1398688051.png",intern=TRUE))
character(0)
> try(system("convert tmp/6ucug1398688051.ps tmp/6ucug1398688051.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  4.036   0.685   4.751