R version 2.11.1 (2010-05-31)
Copyright (C) 2010 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(158,157,156,154,152,151,152,154,155,155,156,158,155,151,152,156,147,159,152,148,146,148,142,141,139,143,146,150,141,153,142,139,134,152,142,139,137,132,137,134,127,140,127,120,112,128,119,117,115,104,108,104,98,111,104,98,91,94,82,81,70,64,71,73,71,80,72,68,60,71,62,69,55,58,67,69,71,85,76,79,69,95,94,99)
> 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,]  158  155  139  137  115   70   55   NA
 [2,]  157  151  143  132  104   64   58   NA
 [3,]  156  152  146  137  108   71   67   NA
 [4,]  154  156  150  134  104   73   69   NA
 [5,]  152  147  141  127   98   71   71   NA
 [6,]  151  159  153  140  111   80   85   NA
 [7,]  152  152  142  127  104   72   76   NA
 [8,]  154  148  139  120   98   68   79   NA
 [9,]  155  146  134  112   91   60   69   NA
[10,]  155  148  152  128   94   71   95   NA
[11,]  156  142  142  119   82   62   94   NA
[12,]  158  141  139  117   81   69   99   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/rcomp/tmp/1hftr1282093178.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/rcomp/tmp/2hftr1282093178.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/rcomp/tmp/3hftr1282093178.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/rcomp/tmp/4a6au1282093178.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] [,12]
[1,]  55.0   58  67.0  69.0  71.0   80   72  68.0   60  71.0    62    69
[2,]  92.5   84  89.5  88.5  84.5   98   90  88.5   80  94.5    88    90
[3,] 137.0  132 137.0 134.0 127.0  140  127 120.0  112 128.0   119   117
[4,] 147.0  147 149.0 152.0 144.0  152  147 143.5  140 150.0   142   140
[5,] 158.0  157 156.0 156.0 152.0  159  152 154.0  155 155.0   156   158

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

$conf
         [,1]      [,2]     [,3]      [,4]      [,5]     [,6]      [,7]
[1,] 104.4535  94.37742 101.4676  96.07882  91.46756 107.7521  92.96052
[2,] 169.5465 169.62258 172.5324 171.92118 162.53244 172.2479 161.03948
          [,8]      [,9]    [,10]     [,11]    [,12]
[1,]  87.15489  76.16897  94.8563  86.75207  87.1408
[2,] 152.84511 147.83103 161.1437 151.24793 146.8592

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(55, 92.5, 137, 147, 158, 58, 84, 132, 147,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/5a6au1282093178.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]
[1,] 151.0 141.0  134 112.0  81.0 60.0 55.0   NA
[2,] 153.0 146.5  139 119.5  92.5 66.0 68.0   NA
[3,] 155.0 149.5  142 127.5 101.0 70.5 73.5   NA
[4,] 156.5 153.5  148 135.5 106.0 71.5 89.5   NA
[5,] 158.0 159.0  153 140.0 115.0 73.0 99.0   NA

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

$conf
         [,1]     [,2]     [,3]     [,4]      [,5]     [,6]    [,7] [,8]
[1,] 153.4036 146.3073 137.8950 120.2023  94.84256 67.99141 63.6937   NA
[2,] 156.5964 152.6927 146.1050 134.7977 107.15744 73.00859 83.3063   NA

$out
[1] 80

$group
[1] 6

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

Warning message:
In bxp(list(stats = c(151, 153, 155, 156.5, 158, 141, 146.5, 149.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/rcomp/tmp/6lg9f1282093178.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,] 33.24870  79.0 50.00
[2,] 33.91385  82.5 54.25
[3,] 35.18768  88.0 56.25
[4,] 37.91625  94.5 59.75
[5,] 41.03193 103.0 63.50

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 33.36216 82.52672 53.74141
[2,] 37.01320 93.47328 58.75859

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(33.2487020870919, 33.9138543565698, 35.1876834499031,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1hftr1282093178.ps tmp/1hftr1282093178.png",intern=TRUE))
character(0)
> try(system("convert tmp/2hftr1282093178.ps tmp/2hftr1282093178.png",intern=TRUE))
character(0)
> try(system("convert tmp/3hftr1282093178.ps tmp/3hftr1282093178.png",intern=TRUE))
character(0)
> try(system("convert tmp/4a6au1282093178.ps tmp/4a6au1282093178.png",intern=TRUE))
character(0)
> try(system("convert tmp/5a6au1282093178.ps tmp/5a6au1282093178.png",intern=TRUE))
character(0)
> try(system("convert tmp/6lg9f1282093178.ps tmp/6lg9f1282093178.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.780   1.080   2.826