R version 2.12.1 (2010-12-16)
Copyright (C) 2010 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.

  Natural language support but running in an English locale

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(70938,34077,45409,40809,37013,44953,19848,32745,43412,34931,33008,8620,68906,39556,50669,36432,40891,48428,36222,33425,39401,37967,34801,12657,69116,41519,51321,38529,41547,52073,38401,40898,40439,41888,37898,8771,68184,50530,47221,41756,45633,48138,39486,39341,41117,41629,29722,7054,56676,34870,35117,30169,30936,35699,33228,27733,33666,35429,27438,8170,63410,38040,45389,37353,37024,50957,37994,36454,46080,43373,37395,10963,75001)
> 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] 73
> (np <- floor(n / par1))
[1] 6
> 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] 7 6 6 6 6 6 6 6 6 6 6 6
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]
 [1,] 70938 68906 69116 68184 56676 63410 75001
 [2,] 34077 39556 41519 50530 34870 38040    NA
 [3,] 45409 50669 51321 47221 35117 45389    NA
 [4,] 40809 36432 38529 41756 30169 37353    NA
 [5,] 37013 40891 41547 45633 30936 37024    NA
 [6,] 44953 48428 52073 48138 35699 50957    NA
 [7,] 19848 36222 38401 39486 33228 37994    NA
 [8,] 32745 33425 40898 39341 27733 36454    NA
 [9,] 43412 39401 40439 41117 33666 46080    NA
[10,] 34931 37967 41888 41629 35429 43373    NA
[11,] 33008 34801 37898 29722 27438 37395    NA
[12,]  8620 12657  8771  7054  8170 10963    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/1lddl1304948374.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/www/html/freestat/rcomp/tmp/2anxi1304948374.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/www/html/freestat/rcomp/tmp/3lqdb1304948374.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/www/html/freestat/rcomp/tmp/4gtnf1304948374.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,] 63410 34077 45389 30169 30936.0 44953 33228 27733.0 33666 34931 27438.0
[2,] 65797 34870 45389 36432 37013.0 44953 33228 32745.0 39401 35429 29722.0
[3,] 68906 38798 46315 37941 38957.5 48283 37108 34939.5 40778 39798 33904.5
[4,] 70027 41519 50669 40809 41547.0 50957 38401 39341.0 43412 41888 37395.0
[5,] 75001 50530 51321 41756 45633.0 52073 39486 40898.0 46080 43373 37898.0
       [,12]
[1,]  7054.0
[2,]  8170.0
[3,]  8695.5
[4,] 10963.0
[5,] 12657.0

$n
 [1] 7 6 6 6 6 6 6 6 6 6 6 6

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 66379.91 34509.18 42909.23 35117.69 36032.92 44410.23 33771.25 30684.87
[2,] 71432.09 43086.82 49720.77 40764.31 41882.08 52155.77 40444.75 39194.13
         [,9]    [,10]    [,11]     [,12]
[1,] 38190.78 35631.74 28955.17  6893.925
[2,] 43365.22 43964.26 38853.83 10497.075

$out
[1] 56676 35117 35699 19848

$group
[1] 1 3 6 7

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

Warning message:
In bxp(list(stats = c(63410, 65797, 68906, 70027, 75001, 34077,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/52wjk1304948374.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]
[1,] 19848.0 33425.0 37898.0 29722.0 27438 36454.0 75001
[2,] 32876.5 35511.5 38465.0 39413.5 28951 37188.5 75001
[3,] 35972.0 38684.0 41208.5 41692.5 33447 38017.0 75001
[4,] 44182.5 44659.5 46604.5 47679.5 35273 45734.5 75001
[5,] 45409.0 50669.0 52073.0 50530.0 35699 50957.0 75001

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]  [,7]
[1,] 30815.26 34511.54 37496.02 37922.32 30563.49 34119.11 75001
[2,] 41128.74 42856.46 44920.98 45462.68 36330.51 41914.89 75001

$out
 [1] 70938  8620 68906 12657 69116  8771 68184  7054 56676  8170 63410 10963

$group
 [1] 1 1 2 2 3 3 4 4 5 5 6 6

$names
[1] "1" "2" "3" "4" "5" "6" NA 

Warning message:
In bxp(list(stats = c(19848, 32876.5, 35972, 44182.5, 45409, 33425,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6naai1304948374.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,] 2052.664  5603.0 2132.500
[2,] 4152.080 11023.5 3903.375
[3,] 4926.767 13931.0 4390.125
[4,] 5904.440 16413.5 5535.000
[5,] 7362.034 19638.0 6203.000

$n
[1] 12 12 12

$conf
         [,1]     [,2]    [,3]
[1,] 4127.504 11472.58 3645.93
[2,] 5726.030 16389.42 5134.32

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2052.66351358424, 4152.08019651607, 4926.76671943296,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1lddl1304948374.ps tmp/1lddl1304948374.png",intern=TRUE))
character(0)
> try(system("convert tmp/2anxi1304948374.ps tmp/2anxi1304948374.png",intern=TRUE))
character(0)
> try(system("convert tmp/3lqdb1304948374.ps tmp/3lqdb1304948374.png",intern=TRUE))
character(0)
> try(system("convert tmp/4gtnf1304948374.ps tmp/4gtnf1304948374.png",intern=TRUE))
character(0)
> try(system("convert tmp/52wjk1304948374.ps tmp/52wjk1304948374.png",intern=TRUE))
character(0)
> try(system("convert tmp/6naai1304948374.ps tmp/6naai1304948374.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.619   0.187   1.967