R version 2.8.0 (2008-10-20)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

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(23100
+ ,22650
+ ,22440
+ ,22910
+ ,22980
+ ,22535
+ ,22300
+ ,22780
+ ,22780
+ ,23300
+ ,23800
+ ,24510
+ ,24660
+ ,24730
+ ,25070
+ ,24690
+ ,24880
+ ,23920
+ ,23880
+ ,23990
+ ,24590
+ ,23610
+ ,23580
+ ,23360
+ ,23910
+ ,23940
+ ,23060
+ ,22800
+ ,23020
+ ,22890
+ ,22780
+ ,22530
+ ,22290
+ ,22820
+ ,22480
+ ,22110
+ ,22000
+ ,22230
+ ,22260
+ ,22590
+ ,22820
+ ,22420
+ ,22230
+ ,21600
+ ,21000
+ ,21360
+ ,21640
+ ,21450
+ ,21710
+ ,21620
+ ,21800
+ ,21490
+ ,21670
+ ,22130
+ ,22050
+ ,22050
+ ,22140
+ ,22390
+ ,22220
+ ,21790
+ ,21510
+ ,21670
+ ,21745
+ ,21850
+ ,22105
+ ,22050
+ ,21670
+ ,21680
+ ,21800
+ ,21920
+ ,21980
+ ,22270
+ ,21740
+ ,21950
+ ,22010
+ ,21890
+ ,21920
+ ,22110
+ ,22340
+ ,22210
+ ,22240
+ ,21960
+ ,22220
+ ,22060
+ ,22090
+ ,21960
+ ,21940
+ ,21790
+ ,21710
+ ,21690
+ ,21710
+ ,21670
+ ,21640
+ ,21500
+ ,21290
+ ,21250
+ ,21580
+ ,21670
+ ,21620
+ ,21510
+ ,21360
+ ,21420
+ ,21470
+ ,21370
+ ,21370
+ ,21340
+ ,21130
+ ,21130
+ ,20990
+ ,21240
+ ,21320
+ ,21430
+ ,21390
+ ,21530
+ ,21510
+ ,21630
+ ,21560
+ ,21610
+ ,21560
+ ,21310
+ ,21340
+ ,21410
+ ,21550
+ ,21380
+ ,21600
+ ,21530
+ ,21560
+ ,21670
+ ,21540
+ ,21540
+ ,21550
+ ,21590
+ ,21420
+ ,21420
+ ,21370
+ ,21380
+ ,21210
+ ,21505
+ ,21365
+ ,21385
+ ,21350
+ ,21360
+ ,21530
+ ,21380
+ ,21630
+ ,22145
+ ,22315
+ ,22340
+ ,22440
+ ,22135
+ ,21955
+ ,22060
+ ,22050
+ ,22035
+ ,22280
+ ,22315
+ ,22205
+ ,21970
+ ,22075
+ ,22115
+ ,22105
+ ,21885
+ ,21805
+ ,21910
+ ,21995
+ ,22245
+ ,22100
+ ,22130
+ ,22300
+ ,22915
+ ,23040
+ ,22880
+ ,23000
+ ,23160
+ ,23020
+ ,22770
+ ,22660
+ ,22740
+ ,22905
+ ,22720
+ ,22705
+ ,22735
+ ,22600
+ ,22510
+ ,22560
+ ,22575
+ ,22685
+ ,22980
+ ,23275
+ ,23845
+ ,23640
+ ,23640
+ ,23835
+ ,23625
+ ,24055
+ ,24005
+ ,24325
+ ,24445
+ ,24670
+ ,24615
+ ,24700
+ ,25065
+ ,25185
+ ,25220
+ ,25235
+ ,24975
+ ,25055
+ ,25520
+ ,25880
+ ,25960
+ ,25740
+ ,24965
+ ,25235
+ ,24895
+ ,24635
+ ,24835
+ ,24635
+ ,24695
+ ,25090
+ ,25220
+ ,24740
+ ,25005
+ ,24650
+ ,24460
+ ,24680
+ ,24840
+ ,24630
+ ,24490
+ ,24695)
> par1 = '12'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 229
> (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] 20 19 19 19 19 19 19 19 19 19 19 19
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 23100 24660 23910 22000 21710 21510 21740 22090 21580 20990 21340 21420
 [2,] 22650 24730 23940 22230 21620 21670 21950 21960 21670 21240 21410 21420
 [3,] 22440 25070 23060 22260 21800 21745 22010 21940 21620 21320 21550 21370
 [4,] 22910 24690 22800 22590 21490 21850 21890 21790 21510 21430 21380 21380
 [5,] 22980 24880 23020 22820 21670 22105 21920 21710 21360 21390 21600 21210
 [6,] 22535 23920 22890 22420 22130 22050 22110 21690 21420 21530 21530 21505
 [7,] 22300 23880 22780 22230 22050 21670 22340 21710 21470 21510 21560 21365
 [8,] 22780 23990 22530 21600 22050 21680 22210 21670 21370 21630 21670 21385
 [9,] 22780 24590 22290 21000 22140 21800 22240 21640 21370 21560 21540 21350
[10,] 23300 23610 22820 21360 22390 21920 21960 21500 21340 21610 21540 21360
[11,] 23800 23580 22480 21640 22220 21980 22220 21290 21130 21560 21550 21530
[12,] 24510 23360 22110 21450 21790 22270 22060 21250 21130 21310 21590 21380
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20]
 [1,] 21630 22205 22300 22705 23835 25235 24635 24695
 [2,] 22145 21970 22915 22735 23625 24975 24695    NA
 [3,] 22315 22075 23040 22600 24055 25055 25090    NA
 [4,] 22340 22115 22880 22510 24005 25520 25220    NA
 [5,] 22440 22105 23000 22560 24325 25880 24740    NA
 [6,] 22135 21885 23160 22575 24445 25960 25005    NA
 [7,] 21955 21805 23020 22685 24670 25740 24650    NA
 [8,] 22060 21910 22770 22980 24615 24965 24460    NA
 [9,] 22050 21995 22660 23275 24700 25235 24680    NA
[10,] 22035 22245 22740 23845 25065 24895 24840    NA
[11,] 22280 22100 22905 23640 25185 24635 24630    NA
[12,] 22315 22130 22720 23640 25220 24835 24490    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/www/html/freestat/rcomp/tmp/1k3y91274912408.ps",horizontal=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/www/html/freestat/rcomp/tmp/2k3y91274912408.ps",horizontal=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/www/html/freestat/rcomp/tmp/3vcyu1274912408.ps",horizontal=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/www/html/freestat/rcomp/tmp/4vcyu1274912408.ps",horizontal=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,] 20990.0 21240 21320.0 21380 21210 21420.0 21365 21370 21000.0 21340 21130
[2,] 21605.0 21670 21772.5 21650 21690 21787.5 21690 21670 21600.0 21575 21600
[3,] 22147.5 22145 22260.0 22340 22440 22135.0 22230 22060 22140.0 22245 22220
[4,] 23872.5 23270 23050.0 22895 23010 23025.0 22900 22880 23027.5 23455 23610
[5,] 25235.0 24975 24055.0 24690 24880 24445.0 24670 24615 24700.0 25065 25185
     [,12]
[1,] 21130
[2,] 21520
[3,] 22130
[4,] 23500
[5,] 25220

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 21346.40 21565.04 21796.94 21888.72 21961.53 21686.43 21791.40 21621.40
[2,] 22948.60 22724.96 22723.06 22791.28 22918.47 22583.57 22668.60 22498.60
         [,9]    [,10]    [,11]    [,12]
[1,] 21622.56 21563.54 21491.42 21412.30
[2,] 22657.44 22926.46 22948.58 22847.70

$out
 [1] 25070 25055 25090 25520 25220 25880 25960 25005 25740 24965 25235

$group
 [1] 3 3 3 4 4 5 6 6 7 8 9

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

Warning message:
In bxp(list(stats = c(20990, 21605, 22147.5, 23872.5, 25235, 21240,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5vcyu1274912408.ps",horizontal=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,] 22300.0 23360 22110 21000 21490 21510 21740 21250 21130 20990 21340
[2,] 22592.5 23745 22505 21525 21690 21675 21935 21570 21350 21315 21470
[3,] 22845.0 24290 22810 22115 21925 21825 22035 21700 21395 21470 21545
[4,] 23200.0 24710 23040 22340 22135 22015 22215 21865 21545 21560 21575
[5,] 23800.0 25070 23060 22820 22390 22270 22340 22090 21670 21630 21670
       [,12]   [,13]   [,14]   [,15]   [,16]   [,17] [,18]   [,19] [,20]
[1,] 21350.0 21955.0 21805.0 22660.0 22510.0 23625.0 24635 24460.0 24695
[2,] 21362.5 22042.5 21940.0 22730.0 22587.5 24030.0 24930 24632.5 24695
[3,] 21380.0 22140.0 22087.5 22892.5 22720.0 24530.0 25145 24687.5 24695
[4,] 21420.0 22315.0 22122.5 23010.0 23457.5 24882.5 25630 24922.5 24695
[5,] 21505.0 22440.0 22245.0 23160.0 23845.0 25220.0 25960 25220.0 24695

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 22567.92 23849.86 22565.98 21743.27 21722.03 21669.92 21907.29 21565.45
[2,] 23122.08 24730.14 23054.02 22486.73 22127.97 21980.08 22162.71 21834.55
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 21306.06 21358.25 21497.11 21353.77 22015.71 22004.26 22764.79 22323.19
[2,] 21483.94 21581.75 21592.89 21406.23 22264.29 22170.74 23020.21 23116.81
        [,17]    [,18]    [,19] [,20]
[1,] 24141.17 24825.73 24555.23 24695
[2,] 24918.83 25464.27 24819.77 24695

$out
[1] 24510 23910 23940 21210 21530 21630 22300

$group
[1]  1  3  3 12 12 13 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(22300, 22592.5, 22845, 23200, 23800, 23360,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/60gnh1274912408.ps",horizontal=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,] 1152.858 3595.0 1210.00
[2,] 1249.560 3752.5 1241.25
[3,] 1274.869 4115.0 1373.75
[4,] 1317.919 4310.0 1930.00
[5,] 1343.122 4670.0 2236.25

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 1243.690 3860.721 1059.607
[2,] 1306.048 4369.279 1687.893

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(1152.85762423397, 1249.55958065991, 1274.86915972299,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1k3y91274912408.ps tmp/1k3y91274912408.png",intern=TRUE))
character(0)
> try(system("convert tmp/2k3y91274912408.ps tmp/2k3y91274912408.png",intern=TRUE))
character(0)
> try(system("convert tmp/3vcyu1274912408.ps tmp/3vcyu1274912408.png",intern=TRUE))
character(0)
> try(system("convert tmp/4vcyu1274912408.ps tmp/4vcyu1274912408.png",intern=TRUE))
character(0)
> try(system("convert tmp/5vcyu1274912408.ps tmp/5vcyu1274912408.png",intern=TRUE))
character(0)
> try(system("convert tmp/60gnh1274912408.ps tmp/60gnh1274912408.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.710   1.295   3.532