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(1664.81,2397.53,2840.71,3547.29,3752.96,3714.74,4349.61,3566.34,5021.82,6423.48,7600.60,19756.21,2499.81,5198.24,7225.14,4806.03,5900.88,4951.34,6179.12,4752.15,5496.43,5835.10,12600.08,28541.72,4717.02,5702.63,9957.58,5304.78,6492.43,6630.80,7349.62,8176.62,8573.17,9690.50,15151.84,34061.01,5921.10,5814.58,12421.25,6369.77,7609.12,7224.75,8121.22,7979.25,8093.06,8476.70,17914.66,30114.41,4826.64,6470.23,9638.77,8821.17,8722.37,10209.48,11276.55,12552.22,11637.39,13606.89,21822.11,45060.69,7615.03,9849.69,14558.40,11587.33,9332.56,13082.09,16732.78,19888.61,23933.38,25391.35,36024.80,80721.71,10243.24,11266.88,21826.84,17357.33,15997.79,18601.53,26155.15,28586.52,30505.41,30821.33,46634.38,104660.67)
> 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] 84
> (np <- floor(n / par1))
[1] 7
> 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 7 7 7 7 7 7 7 7 7 7 7
> arr
          [,1]     [,2]     [,3]     [,4]     [,5]     [,6]      [,7] [,8]
 [1,]  1664.81  2499.81  4717.02  5921.10  4826.64  7615.03  10243.24   NA
 [2,]  2397.53  5198.24  5702.63  5814.58  6470.23  9849.69  11266.88   NA
 [3,]  2840.71  7225.14  9957.58 12421.25  9638.77 14558.40  21826.84   NA
 [4,]  3547.29  4806.03  5304.78  6369.77  8821.17 11587.33  17357.33   NA
 [5,]  3752.96  5900.88  6492.43  7609.12  8722.37  9332.56  15997.79   NA
 [6,]  3714.74  4951.34  6630.80  7224.75 10209.48 13082.09  18601.53   NA
 [7,]  4349.61  6179.12  7349.62  8121.22 11276.55 16732.78  26155.15   NA
 [8,]  3566.34  4752.15  8176.62  7979.25 12552.22 19888.61  28586.52   NA
 [9,]  5021.82  5496.43  8573.17  8093.06 11637.39 23933.38  30505.41   NA
[10,]  6423.48  5835.10  9690.50  8476.70 13606.89 25391.35  30821.33   NA
[11,]  7600.60 12600.08 15151.84 17914.66 21822.11 36024.80  46634.38   NA
[12,] 19756.21 28541.72 34061.01 30114.41 45060.69 80721.71 104660.67   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/1vp1o1304935661.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/2m2we1304935661.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/3cyyz1304935661.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/488o51304935661.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]
[1,]  1664.810  2397.530  2840.710  3547.290 3752.960  3714.74  4349.61
[2,]  3608.415  5450.435  8431.955  5055.405 6196.655  5791.07  6764.37
[3,]  4826.640  5814.580  9957.580  6369.770 7609.120  7224.75  8121.22
[4,]  6768.065  8159.960 13489.825 10204.250 9027.465 11645.78 14004.66
[5,] 10243.240 11266.880 14558.400 17357.330 9332.560 18601.53 16732.78
         [,8]      [,9]    [,10]    [,11]     [,12]
[1,]  3566.34  5021.820  5835.10  7600.60  19756.21
[2,]  6365.70  6794.745  7450.09 13875.96  29328.07
[3,]  8176.62  8573.170  9690.50 17914.66  34061.01
[4,] 16220.42 17785.385 19499.12 28923.46  62891.20
[5,] 28586.52 30505.410 30821.33 46634.38 104660.67

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

$conf
         [,1]     [,2]      [,3]     [,4]     [,5]      [,6]      [,7]
[1,] 2939.748 4196.495  6937.102 3294.963 5918.606  3728.409  3797.433
[2,] 6713.532 7432.665 12978.058 9444.577 9299.634 10721.091 12445.007
          [,8]      [,9]     [,10]     [,11]    [,12]
[1,]  2291.543  2009.737  2495.014  8928.539 14017.65
[2,] 14061.697 15136.603 16885.986 26900.781 54104.37

$out
[1] 21826.84 15997.79 26155.15

$group
[1] 3 5 7

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

Warning message:
In bxp(list(stats = c(1664.81, 3608.415, 4826.64, 6768.065, 10243.24,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5hs591304935661.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,] 1664.81 2499.810  4717.02  5814.580  4826.64  7615.03 10243.24   NA
[2,] 3194.00 4878.685  6097.53  6797.260  8771.77 10718.51 16677.56   NA
[3,] 3733.85 5665.765  7763.12  8036.155 10743.01 15645.59 23991.00   NA
[4,] 5722.65 6702.130  9824.04 10448.975 13079.56 24662.36 30663.37   NA
[5,] 7600.60 7225.140 15151.84 12421.250 13606.89 36024.80 46634.38   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]      [,5]      [,6]     [,7] [,8]
[1,] 2580.516 4834.079 6063.434 6370.583  8778.205  9285.704 17611.97   NA
[2,] 4887.184 6497.451 9462.806 9701.727 12707.825 22005.476 30370.02   NA

$out
 [1]  19756.21  12600.08  28541.72  34061.01  17914.66  30114.41  21822.11
 [8]  45060.69  80721.71 104660.67

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

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

Warning message:
In bxp(list(stats = c(1664.81, 3194, 3733.85, 5722.65, 7600.6, 2499.81,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6jr6p1304935661.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,]  2937.655  8578.43  2709.525
[2,]  4363.967 13027.44  4108.760
[3,]  6815.826 20395.83  6547.505
[4,]  9899.623 25251.89 11519.835
[5,] 13909.592 39033.78 15047.495

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 4290.976 14820.18 3167.264
[2,] 9340.677 25971.49 9927.746

$out
[1] 31552.80 84904.46 33563.14

$group
[1] 1 2 3

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

Warning message:
In bxp(list(stats = c(2937.65464975848, 4363.96659986192, 6815.8263548304,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1vp1o1304935661.ps tmp/1vp1o1304935661.png",intern=TRUE))
character(0)
> try(system("convert tmp/2m2we1304935661.ps tmp/2m2we1304935661.png",intern=TRUE))
character(0)
> try(system("convert tmp/3cyyz1304935661.ps tmp/3cyyz1304935661.png",intern=TRUE))
character(0)
> try(system("convert tmp/488o51304935661.ps tmp/488o51304935661.png",intern=TRUE))
character(0)
> try(system("convert tmp/5hs591304935661.ps tmp/5hs591304935661.png",intern=TRUE))
character(0)
> try(system("convert tmp/6jr6p1304935661.ps tmp/6jr6p1304935661.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.606   0.198   1.885