R version 2.10.1 (2009-12-14)
Copyright (C) 2009 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(505.7,55.7,735.7,575.9,545.8,905.8,765.8,945.7,15.7,645.7,155.9,416,825.8,725.9,925.9,556,116.1,876.3,336.2,186.1,286.1,26,915.8,405.7,965.7,395.6,425.8,545.6,65.6,445.6,895.5,175.4,715.4,865.5,57.4,145.4,315.3,635.4,5.2,515.2,515.1,955,955,634.9,205,275,425,84.9,534.7,4.8,704.7,684.7,884.6,994.6,294.7,524.7,914.5,564.4,984.5,934.4,514.6,474.5,784.4,504.5,824.4,414.6,964.7,64.6,244.7,344.7,34.7,685,425,484.8,785.1,704.9,245.4,285.6,218.8,706.1,856.2,456.6,606.8,527.3,657.8,948.2,486.6,238.9,289.4,969.5,589.5,189.7,639.8,9710.1,969.9,939.9,859.7,679.9,879.9,329.8,349.6,39.5,849.5,449.6,749.6,249.7,649.8,619.4,939,778.9)
> 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] 110
> (np <- floor(n / par1))
[1] 9
> 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] 10 10  9  9  9  9  9  9  9  9  9  9
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]   [,8]  [,9] [,10]
 [1,] 505.7 825.8 965.7 315.3 534.7 514.6 425.0  657.8 859.7 939.0
 [2,]  55.7 725.9 395.6 635.4   4.8 474.5 484.8  948.2 679.9 778.9
 [3,] 735.7 925.9 425.8   5.2 704.7 784.4 785.1  486.6 879.9    NA
 [4,] 575.9 556.0 545.6 515.2 684.7 504.5 704.9  238.9 329.8    NA
 [5,] 545.8 116.1  65.6 515.1 884.6 824.4 245.4  289.4 349.6    NA
 [6,] 905.8 876.3 445.6 955.0 994.6 414.6 285.6  969.5  39.5    NA
 [7,] 765.8 336.2 895.5 955.0 294.7 964.7 218.8  589.5 849.5    NA
 [8,] 945.7 186.1 175.4 634.9 524.7  64.6 706.1  189.7 449.6    NA
 [9,]  15.7 286.1 715.4 205.0 914.5 244.7 856.2  639.8 749.6    NA
[10,] 645.7  26.0 865.5 275.0 564.4 344.7 456.6 9710.1 249.7    NA
[11,] 155.9 915.8  57.4 425.0 984.5  34.7 606.8  969.9 649.8    NA
[12,] 416.0 405.7 145.4  84.9 934.4 685.0 527.3  939.9 619.4    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/yougetitorg/rcomp/tmp/14vzf1306578279.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/yougetitorg/rcomp/tmp/2aqqv1306578279.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/yougetitorg/rcomp/tmp/3tseh1306578279.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/yougetitorg/rcomp/tmp/49eed1306578279.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] [,12]
[1,] 315.30   4.8 425.8 504.5  65.6  39.5 218.8  64.6  15.7  26.0  34.7  84.9
[2,] 505.70 395.6 486.6 504.5 245.4 414.6 336.2 186.1 244.7 275.0 155.9 405.7
[3,] 596.25 560.1 735.7 545.6 349.6 876.3 765.8 449.6 639.8 456.6 606.8 527.3
[4,] 859.70 725.9 785.1 575.9 545.8 955.0 895.5 634.9 749.6 645.7 915.8 685.0
[5,] 965.70 948.2 925.9 575.9 884.6 994.6 964.7 945.7 914.5 865.5 984.5 939.9

$n
 [1] 10 10  9  9  9  9  9  9  9  9  9  9

$conf
         [,1]     [,2]   [,3]    [,4]     [,5]      [,6]      [,7]    [,8]
[1,] 419.3775 395.0690 578.49 507.996 191.3893  591.6893  471.2353 213.232
[2,] 773.1225 725.1310 892.91 583.204 507.8107 1160.9107 1060.3647 685.968
        [,9]    [,10]    [,11]   [,12]
[1,] 373.886 261.3647  206.586 380.202
[2,] 905.714 651.8353 1007.014 674.398

$out
[1]    5.2  684.7  704.9  238.9  329.8 9710.1

$group
[1]  3  4  4  4  4 10

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

Warning message:
In bxp(list(stats = c(315.3, 505.7, 596.25, 859.7, 965.7, 4.8, 395.6,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/51bhy1306578279.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]
[1,]  15.70  26.00  57.40   5.20   4.80  34.7 218.80 189.70  39.50 778.90
[2,] 285.95 236.10 160.40 240.00 529.70 294.7 355.30 388.00 339.70 778.90
[3,] 560.85 480.85 435.70 470.05 694.70 489.5 506.05 648.80 634.60 858.95
[4,] 750.75 851.05 790.45 635.15 924.45 734.7 705.50 958.85 799.55 939.00
[5,] 945.70 925.90 965.70 955.00 994.60 964.7 856.20 969.90 879.90 939.00

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 348.8516 200.3672 148.3300 289.8194 514.6519 288.8130 346.3214 388.4315
[2,] 772.8484 761.3328 723.0700 650.2806 874.7481 690.1870 665.7786 909.1685
         [,9]     [,10]
[1,] 424.8593  680.0817
[2,] 844.3407 1037.8183

$out
[1] 9710.1

$group
[1] 8

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

Warning message:
In bxp(list(stats = c(15.7, 285.95, 560.85, 750.75, 945.7, 26, 236.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/yougetitorg/rcomp/tmp/6br2n1306578279.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,] 229.2595 650.40  71.4000
[2,] 289.5256 782.45 298.7875
[3,] 301.8471 889.95 357.0000
[4,] 343.3844 946.60 522.6500
[5,] 385.4744 955.10 759.9000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 277.2818 815.0801 254.8948
[2,] 326.4125 964.8199 459.1052

$out
[1]  150.7399 3103.6043  466.0000 9684.1000

$group
[1] 1 1 2 2

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

Warning message:
In bxp(list(stats = c(229.259532941259, 289.525607035137, 301.847119071763,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/14vzf1306578279.ps tmp/14vzf1306578279.png",intern=TRUE))
character(0)
> try(system("convert tmp/2aqqv1306578279.ps tmp/2aqqv1306578279.png",intern=TRUE))
character(0)
> try(system("convert tmp/3tseh1306578279.ps tmp/3tseh1306578279.png",intern=TRUE))
character(0)
> try(system("convert tmp/49eed1306578279.ps tmp/49eed1306578279.png",intern=TRUE))
character(0)
> try(system("convert tmp/51bhy1306578279.ps tmp/51bhy1306578279.png",intern=TRUE))
character(0)
> try(system("convert tmp/6br2n1306578279.ps tmp/6br2n1306578279.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.250   0.960   1.782