R version 2.15.2 (2012-10-26) -- "Trick or Treat"
Copyright (C) 2012 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
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(235.1
+ ,280.7
+ ,264.6
+ ,240.7
+ ,201.4
+ ,240.8
+ ,241.1
+ ,223.8
+ ,206.1
+ ,174.7
+ ,203.3
+ ,220.5
+ ,299.5
+ ,347.4
+ ,338.3
+ ,327.7
+ ,351.6
+ ,396.6
+ ,438.8
+ ,395.6
+ ,363.5
+ ,378.8
+ ,357
+ ,369
+ ,464.8
+ ,479.1
+ ,431.3
+ ,366.5
+ ,326.3
+ ,355.1
+ ,331.6
+ ,261.3
+ ,249
+ ,205.5
+ ,235.6
+ ,240.9
+ ,264.9
+ ,253.8
+ ,232.3
+ ,193.8
+ ,177
+ ,213.2
+ ,207.2
+ ,180.6
+ ,188.6
+ ,175.4
+ ,199
+ ,179.6
+ ,225.8
+ ,234
+ ,200.2
+ ,183.6
+ ,178.2
+ ,203.2
+ ,208.5
+ ,191.8
+ ,172.8
+ ,148
+ ,159.4
+ ,154.5
+ ,213.2
+ ,196.4
+ ,182.8
+ ,176.4
+ ,153.6
+ ,173.2
+ ,171
+ ,151.2
+ ,161.9
+ ,157.2
+ ,201.7
+ ,236.4
+ ,356.1
+ ,398.3
+ ,403.7
+ ,384.6
+ ,365.8
+ ,368.1
+ ,367.9
+ ,347
+ ,343.3
+ ,292.9
+ ,311.5
+ ,300.9
+ ,366.9
+ ,356.9
+ ,329.7
+ ,316.2
+ ,269
+ ,289.3
+ ,266.2
+ ,253.6
+ ,233.8
+ ,228.4
+ ,253.6
+ ,260.1
+ ,306.6
+ ,309.2
+ ,309.5
+ ,271
+ ,279.9
+ ,317.9
+ ,298.4
+ ,246.7
+ ,227.3
+ ,209.1
+ ,259.9
+ ,266
+ ,320.6
+ ,308.5
+ ,282.2
+ ,262.7
+ ,263.5
+ ,313.1
+ ,284.3
+ ,252.6
+ ,250.3
+ ,246.5
+ ,312.7
+ ,333.2
+ ,446.4
+ ,511.6
+ ,515.5
+ ,506.4
+ ,483.2
+ ,522.3
+ ,509.8
+ ,460.7
+ ,405.8
+ ,375
+ ,378.5
+ ,406.8
+ ,467.8
+ ,469.8
+ ,429.8
+ ,355.8
+ ,332.7
+ ,378
+ ,360.5
+ ,334.7
+ ,319.5
+ ,323.1
+ ,363.6
+ ,352.1
+ ,411.9
+ ,388.6
+ ,416.4
+ ,360.7
+ ,338
+ ,417.2
+ ,388.4
+ ,371.1
+ ,331.5
+ ,353.7
+ ,396.7
+ ,447
+ ,533.5
+ ,565.4
+ ,542.3
+ ,488.7
+ ,467.1
+ ,531.3
+ ,496.1
+ ,444
+ ,403.4
+ ,386.3
+ ,394.1
+ ,404.1
+ ,462.1
+ ,448.1
+ ,432.3
+ ,386.3
+ ,395.2
+ ,421.9
+ ,382.9
+ ,384.2
+ ,345.5
+ ,323.4
+ ,372.6
+ ,376
+ ,462.7
+ ,487
+ ,444.2
+ ,399.3
+ ,394.9
+ ,455.4
+ ,414
+ ,375.5
+ ,347
+ ,339.4
+ ,385.8
+ ,378.8
+ ,451.8
+ ,446.1
+ ,422.5
+ ,383.1
+ ,352.8
+ ,445.3
+ ,367.5
+ ,355.1
+ ,326.2
+ ,319.8
+ ,331.8
+ ,340.9
+ ,394.1
+ ,417.2
+ ,369.9
+ ,349.2
+ ,321.4
+ ,405.7
+ ,342.9
+ ,316.5
+ ,284.2
+ ,270.9
+ ,288.8
+ ,278.8
+ ,324.4
+ ,310.9
+ ,299
+ ,273
+ ,279.3
+ ,359.2
+ ,305
+ ,282.1
+ ,250.3
+ ,246.5
+ ,257.9
+ ,266.5
+ ,315.9
+ ,318.4
+ ,295.4
+ ,266.4
+ ,245.8
+ ,362.8
+ ,324.9
+ ,294.2
+ ,289.5
+ ,295.2
+ ,290.3
+ ,272
+ ,307.4
+ ,328.7
+ ,292.9
+ ,249.1
+ ,230.4
+ ,361.5
+ ,321.7
+ ,277.2
+ ,260.7
+ ,251
+ ,257.6
+ ,241.8
+ ,287.5
+ ,292.3
+ ,274.7
+ ,254.2
+ ,230
+ ,339
+ ,318.2
+ ,287
+ ,295.8
+ ,284
+ ,271
+ ,262.7
+ ,340.6
+ ,379.4
+ ,373.3
+ ,355.2
+ ,338.4
+ ,466.9
+ ,451
+ ,422
+ ,429.2
+ ,425.9
+ ,460.7
+ ,463.6
+ ,541.4
+ ,544.2
+ ,517.5
+ ,469.4
+ ,439.4
+ ,549
+ ,533
+ ,506.1
+ ,484
+ ,457
+ ,481.5
+ ,469.5
+ ,544.7
+ ,541.2
+ ,521.5
+ ,469.7
+ ,434.4
+ ,542.6
+ ,517.3
+ ,485.7
+ ,465.8
+ ,447
+ ,426.6
+ ,411.6
+ ,467.5
+ ,484.5
+ ,451.2
+ ,417.4
+ ,379.9
+ ,484.7
+ ,455
+ ,420.8
+ ,416.5
+ ,376.3
+ ,405.6
+ ,405.8
+ ,500.8
+ ,514
+ ,475.5
+ ,430.1
+ ,414.4
+ ,538
+ ,526
+ ,488.5
+ ,520.2
+ ,504.4
+ ,568.5
+ ,610.6
+ ,818
+ ,830.9
+ ,835.9
+ ,782
+ ,762.3
+ ,856.9
+ ,820.9
+ ,769.6
+ ,752.2
+ ,724.4
+ ,723.1
+ ,719.5
+ ,817.4
+ ,803.3
+ ,752.5
+ ,689
+ ,630.4
+ ,765.5
+ ,757.7
+ ,732.2
+ ,702.6
+ ,683.3
+ ,709.5
+ ,702.2
+ ,784.8
+ ,810.9
+ ,755.6
+ ,656.8
+ ,615.1
+ ,745.3
+ ,694.1
+ ,675.7
+ ,643.7
+ ,622.1
+ ,634.6
+ ,588
+ ,689.7
+ ,673.9
+ ,647.9
+ ,568.8
+ ,545.7
+ ,632.6
+ ,643.8
+ ,593.1
+ ,579.7
+ ,546
+ ,562.9
+ ,572.5)
> 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] 372
> (np <- floor(n / par1))
[1] 31
> 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] 31 31 31 31 31 31 31 31 31 31 31 31
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 235.1 299.5 464.8 264.9 225.8 213.2 356.1 366.9 306.6 320.6 446.4 467.8
 [2,] 280.7 347.4 479.1 253.8 234.0 196.4 398.3 356.9 309.2 308.5 511.6 469.8
 [3,] 264.6 338.3 431.3 232.3 200.2 182.8 403.7 329.7 309.5 282.2 515.5 429.8
 [4,] 240.7 327.7 366.5 193.8 183.6 176.4 384.6 316.2 271.0 262.7 506.4 355.8
 [5,] 201.4 351.6 326.3 177.0 178.2 153.6 365.8 269.0 279.9 263.5 483.2 332.7
 [6,] 240.8 396.6 355.1 213.2 203.2 173.2 368.1 289.3 317.9 313.1 522.3 378.0
 [7,] 241.1 438.8 331.6 207.2 208.5 171.0 367.9 266.2 298.4 284.3 509.8 360.5
 [8,] 223.8 395.6 261.3 180.6 191.8 151.2 347.0 253.6 246.7 252.6 460.7 334.7
 [9,] 206.1 363.5 249.0 188.6 172.8 161.9 343.3 233.8 227.3 250.3 405.8 319.5
[10,] 174.7 378.8 205.5 175.4 148.0 157.2 292.9 228.4 209.1 246.5 375.0 323.1
[11,] 203.3 357.0 235.6 199.0 159.4 201.7 311.5 253.6 259.9 312.7 378.5 363.6
[12,] 220.5 369.0 240.9 179.6 154.5 236.4 300.9 260.1 266.0 333.2 406.8 352.1
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
 [1,] 411.9 533.5 462.1 462.7 451.8 394.1 324.4 315.9 307.4 287.5 340.6 541.4
 [2,] 388.6 565.4 448.1 487.0 446.1 417.2 310.9 318.4 328.7 292.3 379.4 544.2
 [3,] 416.4 542.3 432.3 444.2 422.5 369.9 299.0 295.4 292.9 274.7 373.3 517.5
 [4,] 360.7 488.7 386.3 399.3 383.1 349.2 273.0 266.4 249.1 254.2 355.2 469.4
 [5,] 338.0 467.1 395.2 394.9 352.8 321.4 279.3 245.8 230.4 230.0 338.4 439.4
 [6,] 417.2 531.3 421.9 455.4 445.3 405.7 359.2 362.8 361.5 339.0 466.9 549.0
 [7,] 388.4 496.1 382.9 414.0 367.5 342.9 305.0 324.9 321.7 318.2 451.0 533.0
 [8,] 371.1 444.0 384.2 375.5 355.1 316.5 282.1 294.2 277.2 287.0 422.0 506.1
 [9,] 331.5 403.4 345.5 347.0 326.2 284.2 250.3 289.5 260.7 295.8 429.2 484.0
[10,] 353.7 386.3 323.4 339.4 319.8 270.9 246.5 295.2 251.0 284.0 425.9 457.0
[11,] 396.7 394.1 372.6 385.8 331.8 288.8 257.9 290.3 257.6 271.0 460.7 481.5
[12,] 447.0 404.1 376.0 378.8 340.9 278.8 266.5 272.0 241.8 262.7 463.6 469.5
      [,25] [,26] [,27] [,28] [,29] [,30] [,31] [,32]
 [1,] 544.7 467.5 500.8 818.0 817.4 784.8 689.7    NA
 [2,] 541.2 484.5 514.0 830.9 803.3 810.9 673.9    NA
 [3,] 521.5 451.2 475.5 835.9 752.5 755.6 647.9    NA
 [4,] 469.7 417.4 430.1 782.0 689.0 656.8 568.8    NA
 [5,] 434.4 379.9 414.4 762.3 630.4 615.1 545.7    NA
 [6,] 542.6 484.7 538.0 856.9 765.5 745.3 632.6    NA
 [7,] 517.3 455.0 526.0 820.9 757.7 694.1 643.8    NA
 [8,] 485.7 420.8 488.5 769.6 732.2 675.7 593.1    NA
 [9,] 465.8 416.5 520.2 752.2 702.6 643.7 579.7    NA
[10,] 447.0 376.3 504.4 724.4 683.3 622.1 546.0    NA
[11,] 426.6 405.6 568.5 723.1 709.5 634.6 562.9    NA
[12,] 411.6 405.8 610.6 719.5 702.2 588.0 572.5    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/1pw9m1355248261.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/2d5pl1355248261.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/3f6ls1355248261.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/47adr1355248261.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,] 213.20 196.40 182.8 176.40 153.60 173.20 171.00 151.20 161.90 148.0 159.4
[2,] 311.65 314.65 297.2 268.70 266.25 347.05 311.60 269.25 250.30 246.5 258.9
[3,] 411.90 417.20 416.4 360.70 338.40 405.70 367.90 355.10 331.50 323.1 357.0
[4,] 484.30 512.80 495.5 449.75 424.40 526.80 502.95 452.35 422.85 406.1 416.1
[5,] 689.70 803.30 755.6 689.00 630.40 765.50 757.70 675.70 643.70 622.1 634.6
      [,12]
[1,] 154.50
[2,] 264.35
[3,] 352.10
[4,] 429.30
[5,] 610.60

$n
 [1] 31 31 31 31 31 31 31 31 31 31 31 31

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 362.906 360.9697 360.1271 309.3223 293.5208 354.6912 313.5994 303.1405
[2,] 460.894 473.4303 472.6729 412.0777 383.2792 456.7088 422.2006 407.0595
         [,9]    [,10]    [,11]    [,12]
[1,] 282.5344 277.8093 312.3904 305.2911
[2,] 380.4656 368.3907 401.6096 398.9089

$out
 [1] 818.0 817.4 784.8 830.9 810.9 835.9 782.0 762.3 856.9 820.9 769.6 732.2
[13] 752.2 702.6 724.4 683.3 723.1 709.5 719.5 702.2

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

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

> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/55wkl1355248261.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,] 174.70 299.50 205.50 175.40 148.00 151.20 292.90 228.40 209.10 246.50
[2,] 204.70 342.85 244.95 180.10 166.10 159.55 327.40 253.60 253.30 257.65
[3,] 229.45 360.25 328.95 196.40 187.70 174.80 360.95 267.60 275.45 283.25
[4,] 240.95 387.20 398.90 222.75 205.85 199.05 376.35 322.95 307.90 312.90
[5,] 280.70 438.80 479.10 264.90 234.00 236.40 403.70 366.90 317.90 333.20
      [,11]  [,12]  [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
[1,] 375.00 319.50 331.50 386.30 323.40 339.40 319.80 270.90 246.50 245.80
[2,] 406.30 333.70 357.20 403.75 374.30 377.15 336.35 286.50 262.20 280.75
[3,] 471.95 358.15 388.50 477.90 385.25 397.10 361.30 332.15 280.70 294.70
[4,] 510.70 403.90 414.15 532.40 427.10 449.80 433.90 382.00 307.95 317.15
[5,] 522.30 469.80 447.00 565.40 462.10 487.00 451.80 417.20 359.20 362.80
      [,21]  [,22]  [,23]  [,24]  [,25]  [,26] [,27] [,28]  [,29]  [,30]  [,31]
[1,] 230.40 230.00 338.40 439.40 411.60 376.30 414.4 719.5 630.40 588.00 545.70
[2,] 250.05 266.85 364.25 469.45 440.70 405.70 482.0 738.3 695.60 628.35 565.85
[3,] 268.95 285.50 423.95 495.05 477.70 419.10 509.2 775.8 720.85 666.25 586.40
[4,] 314.55 294.05 455.85 537.20 531.35 461.25 532.0 825.9 761.60 750.45 645.85
[5,] 361.50 318.20 466.90 549.00 544.70 484.70 568.5 856.9 817.40 810.90 689.70
     [,32]
[1,]    NA
[2,]    NA
[3,]    NA
[4,]    NA
[5,]    NA

$n
 [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12
[26] 12 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 212.9161 340.0217 258.7324 176.947 169.5698 156.7838 338.6236 235.969
[2,] 245.9839 380.4783 399.1676 215.853 205.8302 192.8162 383.2764 299.231
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 250.5466 258.0501 424.3325 326.1313 362.5247 419.2219 361.1676 363.9638
[2,] 300.3534 308.4499 519.5675 390.1687 414.4753 536.5781 409.3324 430.2362
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 316.8068 288.5918 259.8331 278.0977 239.5311 273.0939 382.1706 464.1488
[2,] 405.7932 375.7082 301.5669 311.3023 298.3689 297.9061 465.7294 525.9512
        [,25]    [,26]    [,27]    [,28]   [,29]    [,30]    [,31] [,32]
[1,] 436.3539 393.7633 486.3947 735.8451 690.747 610.5594 549.9115    NA
[2,] 519.0461 444.4367 532.0053 815.7549 750.953 721.9406 622.8885    NA

$out
[1] 339.0 610.6

$group
[1] 22 27

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30"
[31] "31" NA  

Warning message:
In bxp(list(stats = c(174.7, 204.7, 229.45, 240.95, 280.7, 299.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6sdmq1355248261.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,] 140.8085 563.70 157.200
[2,] 146.4626 583.35 162.275
[3,] 153.3961 607.15 176.200
[4,] 161.2180 642.20 187.225
[5,] 165.1843 683.70 198.300

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 146.6660 580.3081 164.8201
[2,] 160.1261 633.9919 187.5799

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(140.808450229075, 146.462567853424, 153.396079766919,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1pw9m1355248261.ps tmp/1pw9m1355248261.png",intern=TRUE))
character(0)
> try(system("convert tmp/2d5pl1355248261.ps tmp/2d5pl1355248261.png",intern=TRUE))
character(0)
> try(system("convert tmp/3f6ls1355248261.ps tmp/3f6ls1355248261.png",intern=TRUE))
character(0)
> try(system("convert tmp/47adr1355248261.ps tmp/47adr1355248261.png",intern=TRUE))
character(0)
> try(system("convert tmp/55wkl1355248261.ps tmp/55wkl1355248261.png",intern=TRUE))
character(0)
> try(system("convert tmp/6sdmq1355248261.ps tmp/6sdmq1355248261.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  4.367   0.754   5.091