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(90,51,47,59,54,79,59,80,46,62,55,77,72,72,71,50,66,78,59,52,71,98,70,84,90,98,98,78,59,0,58,55,62,80,91,86,61,49,61,56,73,85,82,32,39,30,51,48,57,59,32,56,54,74,62,78,72,48,59,61,80,69,58,63,27,23,34,45,51,51,73,37,35,66,54,30,66,61,37,55,64,53,63,70,72,52,53,50,60,73,66,78)
> par1 = '7'
> #'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] 92
> (np <- floor(n / par1))
[1] 13
> 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] 14 13 13 13 13 13 13
> arr
     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14]
[1,]   90   80   71   98   59   86   82   59   72    63    73    61    72    78
[2,]   51   46   50   70    0   61   32   32   48    27    37    37    52    NA
[3,]   47   62   66   84   58   49   39   56   59    23    35    55    53    NA
[4,]   59   55   78   90   55   61   30   54   61    34    66    64    50    NA
[5,]   54   77   59   98   62   56   51   74   80    45    54    53    60    NA
[6,]   79   72   52   98   80   73   48   62   69    51    30    63    73    NA
[7,]   59   72   71   78   91   85   57   78   58    51    66    70    66    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/1dth21305852959.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/2kxa11305852959.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/3f4hz1305852959.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/439ri1305852959.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,] 59.0   27   35   50   45   30   51
[2,] 63.0   32   47   54   54   52   59
[3,] 72.5   46   55   59   59   69   70
[4,] 82.0   51   59   64   74   73   78
[5,] 98.0   70   66   78   98   98   91

$n
[1] 14 13 13 13 13 13 13

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]
[1,] 64.47682 37.67395 49.74144 54.61787 50.23574 59.79752 61.67395
[2,] 80.52318 54.32605 60.25856 63.38213 67.76426 78.20248 78.32605

$out
[1]  0 84 23 90 30 34

$group
[1] 2 3 3 4 4 4

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

Warning message:
In bxp(list(stats = c(59, 63, 72.5, 82, 98, 27, 32, 46, 51, 70,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/52zb71305852959.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] [,9] [,10] [,11] [,12] [,13] [,14]
[1,] 47.0 46.0 50.0   70 55.0 49.0 30.0   54 48.0  23.0    30  53.0  50.0    78
[2,] 52.5 58.5 55.5   81 56.5 58.5 35.5   55 58.5  30.5    36  54.0  52.5    78
[3,] 59.0 72.0 66.0   90 59.0 61.0 48.0   59 61.0  45.0    54  61.0  60.0    78
[4,] 69.0 74.5 71.0   98 71.0 79.0 54.0   68 70.5  51.0    66  63.5  69.0    78
[5,] 90.0 80.0 78.0   98 91.0 86.0 57.0   78 80.0  63.0    73  70.0  73.0    78

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

$conf
         [,1]     [,2]     [,3]      [,4]     [,5]     [,6]    [,7]     [,8]
[1,] 49.14647 62.44506 56.74365  79.84787 50.34083 48.75773 36.9521 51.23661
[2,] 68.85353 81.55494 75.25635 100.15213 67.65917 73.24227 59.0479 66.76339
         [,9]    [,10]    [,11]    [,12]    [,13] [,14]
[1,] 53.83379 32.75773 36.08448 55.32675 50.14647    78
[2,] 68.16621 57.24227 71.91552 66.67325 69.85353    78

$out
[1]  0 82 32 37

$group
[1]  5  7  8 12

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

Warning message:
In bxp(list(stats = c(47, 52.5, 59, 69, 90, 46, 58.5, 72, 74.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/68pzx1305852959.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,] 11.62281 39.0 10.00
[2,] 13.36611 46.5 14.25
[3,] 15.15560 60.0 19.00
[4,] 16.55802 64.5 19.50
[5,] 17.51263 70.0 21.00

$n
[1] 7 7 7

$conf
         [,1]     [,2]     [,3]
[1,] 13.24945 49.25069 15.86478
[2,] 17.06176 70.74931 22.13522

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(11.6228113462167, 13.3661131896434, 15.155603178109,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1dth21305852959.ps tmp/1dth21305852959.png",intern=TRUE))
character(0)
> try(system("convert tmp/2kxa11305852959.ps tmp/2kxa11305852959.png",intern=TRUE))
character(0)
> try(system("convert tmp/3f4hz1305852959.ps tmp/3f4hz1305852959.png",intern=TRUE))
character(0)
> try(system("convert tmp/439ri1305852959.ps tmp/439ri1305852959.png",intern=TRUE))
character(0)
> try(system("convert tmp/52zb71305852959.ps tmp/52zb71305852959.png",intern=TRUE))
character(0)
> try(system("convert tmp/68pzx1305852959.ps tmp/68pzx1305852959.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.613   0.178   1.801