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(26281,23899,25727,30733,28599,16723,43738,45272,46532,41032,37967,35366,33892,21560,26588,33527,24859,17952,45504,40129,40357,41913,33730,37842,33025,24050,30429,34507,25189,20253,48527,44446,46380,48950,38883,42928,37107,30186,32602,39892,32194,21629,59968,45694,55756,48554,41052,49822,39191,31994,35735,38930,33658,23849,58972,59249,63955,53785,52760,44795,37348,32370,32717,40974,33591,21124,58608,46865,51378,46235,47206,45382,41227,33795,31295,42625,33625,21538,56421,53152,53536,52408,41454,38271,35306,26414,31917,38030,27534,18387)
> 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] 90
> (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] 8 8 8 8 8 8 7 7 7 7 7 7
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]
 [1,] 26281 33892 33025 37107 39191 37348 41227 35306
 [2,] 23899 21560 24050 30186 31994 32370 33795 26414
 [3,] 25727 26588 30429 32602 35735 32717 31295 31917
 [4,] 30733 33527 34507 39892 38930 40974 42625 38030
 [5,] 28599 24859 25189 32194 33658 33591 33625 27534
 [6,] 16723 17952 20253 21629 23849 21124 21538 18387
 [7,] 43738 45504 48527 59968 58972 58608 56421    NA
 [8,] 45272 40129 44446 45694 59249 46865 53152    NA
 [9,] 46532 40357 46380 55756 63955 51378 53536    NA
[10,] 41032 41913 48950 48554 53785 46235 52408    NA
[11,] 37967 33730 38883 41052 52760 47206 41454    NA
[12,] 35366 37842 42928 49822 44795 45382 38271    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/1fczk1304964024.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/2al0y1304964024.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/3st531304964024.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/45wdf1304964024.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,] 26281.0 21560.0 25727.0 30733 24859.0 16723.0 43738.0 40129.0 40357 41032
[2,] 33458.5 23974.5 28508.5 34017 26361.5 18169.5 47015.5 44859.0 46456 44074
[3,] 36206.5 28300.0 31606.0 38480 30396.5 20688.5 56421.0 45694.0 51378 48554
[4,] 38269.5 32182.0 32659.5 40433 33608.0 21583.5 58790.0 50008.5 54646 50679
[5,] 41227.0 33795.0 35735.0 42625 33658.0 23849.0 59968.0 53152.0 63955 53785
     [,11]   [,12]
[1,] 33730 35366.0
[2,] 38425 38056.5
[3,] 41052 42928.0
[4,] 44330 45088.5
[5,] 52760 49822.0

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6]     [,7]    [,8]
[1,] 33519.01 23715.17 29287.19 34895.93 26348.5 18781.39 49389.46 42618.8
[2,] 38893.99 32884.83 33924.81 42064.07 34444.5 22595.61 63452.54 48769.2
         [,9]   [,10]    [,11]   [,12]
[1,] 46487.06 44609.6 37525.63 38728.6
[2,] 56268.94 52498.4 44578.37 47127.4

$out
[1] 59249

$group
[1] 8

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

Warning message:
In bxp(list(stats = c(26281, 33458.5, 36206.5, 38269.5, 41227, 21560,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/54mzy1304964024.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,] 16723.0 17952.0 20253 21629 23849.0 21124.0 21538.0 18387.0
[2,] 26004.0 25723.5 27809 32398 34696.5 33154.0 33710.0 26414.0
[3,] 33049.5 33811.0 36695 40472 41993.0 43178.0 41340.5 29725.5
[4,] 42385.0 40243.0 45413 49188 56378.5 47035.5 52780.0 35306.0
[5,] 46532.0 45504.0 48950 59968 63955.0 58608.0 56421.0 38030.0

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

$conf
         [,1]     [,2]    [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 25578.02 27188.56 28665.7 32813.97 32103.69 36846.55 32642.54 23989.87
[2,] 40520.98 40433.44 44724.3 48130.03 51882.31 49509.45 50038.46 35461.13

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(16723, 26004, 33049.5, 42385, 46532, 17952,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/67hna1304964024.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,] 2343.119  7126.0  3162.000
[2,] 3949.324 10950.0  4641.500
[3,] 4740.642 13604.5  6255.000
[4,] 6323.422 17630.0  7553.875
[5,] 7653.893 23598.0 11774.500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 3657.799 10557.71 4926.646
[2,] 5823.484 16651.29 7583.354

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2343.11894933106, 3949.32418373775, 4740.64155035683,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1fczk1304964024.ps tmp/1fczk1304964024.png",intern=TRUE))
character(0)
> try(system("convert tmp/2al0y1304964024.ps tmp/2al0y1304964024.png",intern=TRUE))
character(0)
> try(system("convert tmp/3st531304964024.ps tmp/3st531304964024.png",intern=TRUE))
character(0)
> try(system("convert tmp/45wdf1304964024.ps tmp/45wdf1304964024.png",intern=TRUE))
character(0)
> try(system("convert tmp/54mzy1304964024.ps tmp/54mzy1304964024.png",intern=TRUE))
character(0)
> try(system("convert tmp/67hna1304964024.ps tmp/67hna1304964024.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.627   0.191   1.888