R version 2.9.0 (2009-04-17)
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(3030.29,2803.47,2767.63,2882.6,2863.36,2897.06,3012.61,3142.95,3032.93,3045.78,3110.52,3013.24,2987.1,2995.55,2833.18,2848.96,2794.83,2845.26,2915.02,2892.63,2604.42,2641.65,2659.81,2638.53,2720.25,2745.88,2735.7,2811.7,2799.43,2555.28,2304.98,2214.95,2065.81,1940.49,2042,1995.37,1946.81,1765.9,1635.25,1833.42,1910.43,1959.67,1969.6,2061.41,2091.48,2120.88,2174.56,2196.72,2350.44,2440.25,2408.64,2472.81,2407.6,2452.62,2448.05,2497.84,2645.64,2756.76,2849.27,2921.44,2981.85,3080.58,3106.22,3119.31,3061.26,3097.31,3161.69,3257.16,3277.01,3295.32,3363.99,3494.17,3667.03,3813.06,3917.96,3895.51,3801.06,3570.12,3701.61,3862.27,3970.1,4138.52,4199.75,4290.89,4443.91,4502.64,4356.98,4591.27,4696.96,4621.4,4562.84,4202.52,4296.49,4435.23,4105.18,4116.68,3844.49,3720.98,3674.4,3857.62,3801.06,3504.37,3032.6,3047.03,2962.34,2197.82,2014.45,1862.83,1905.41)
> 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] 109
> (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  9  9  9  9  9  9  9  9  9  9  9
> arr
         [,1]    [,2]    [,3]    [,4]    [,5]    [,6]    [,7]    [,8]    [,9]
 [1,] 3030.29 2987.10 2720.25 1946.81 2350.44 2981.85 3667.03 4443.91 3844.49
 [2,] 2803.47 2995.55 2745.88 1765.90 2440.25 3080.58 3813.06 4502.64 3720.98
 [3,] 2767.63 2833.18 2735.70 1635.25 2408.64 3106.22 3917.96 4356.98 3674.40
 [4,] 2882.60 2848.96 2811.70 1833.42 2472.81 3119.31 3895.51 4591.27 3857.62
 [5,] 2863.36 2794.83 2799.43 1910.43 2407.60 3061.26 3801.06 4696.96 3801.06
 [6,] 2897.06 2845.26 2555.28 1959.67 2452.62 3097.31 3570.12 4621.40 3504.37
 [7,] 3012.61 2915.02 2304.98 1969.60 2448.05 3161.69 3701.61 4562.84 3032.60
 [8,] 3142.95 2892.63 2214.95 2061.41 2497.84 3257.16 3862.27 4202.52 3047.03
 [9,] 3032.93 2604.42 2065.81 2091.48 2645.64 3277.01 3970.10 4296.49 2962.34
[10,] 3045.78 2641.65 1940.49 2120.88 2756.76 3295.32 4138.52 4435.23 2197.82
[11,] 3110.52 2659.81 2042.00 2174.56 2849.27 3363.99 4199.75 4105.18 2014.45
[12,] 3013.24 2638.53 1995.37 2196.72 2921.44 3494.17 4290.89 4116.68 1862.83
        [,10]
 [1,] 1905.41
 [2,]      NA
 [3,]      NA
 [4,]      NA
 [5,]      NA
 [6,]      NA
 [7,]      NA
 [8,]      NA
 [9,]      NA
[10,]      NA
[11,]      NA
[12,]      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/rcomp/tmp/197j31243508773.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/html/rcomp/tmp/2jd911243508773.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/html/rcomp/tmp/3w4dm1243508773.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/html/rcomp/tmp/4u17b1243508773.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]
[1,] 1905.410 1765.90 1635.25 1833.42 1910.43 1959.67 1969.60 2061.41 2065.81
[2,] 2350.440 2745.88 2735.70 2811.70 2794.83 2555.28 2448.05 2497.84 2604.42
[3,] 2984.475 2995.55 2833.18 2882.60 2863.36 2897.06 3012.61 3047.03 2962.34
[4,] 3667.030 3720.98 3674.40 3857.62 3801.06 3504.37 3161.69 3257.16 3277.01
[5,] 4443.910 4502.64 4356.98 4591.27 4696.96 4621.40 3701.61 4202.52 3970.10
       [,10]   [,11]   [,12]
[1,] 1940.49 2014.45 1862.83
[2,] 2197.82 2174.56 2196.72
[3,] 2756.76 2849.27 2921.44
[4,] 3295.32 3363.99 3494.17
[5,] 4435.23 4199.75 4290.89

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2326.654 2481.997 2338.798 2331.749 2333.412 2397.206 2636.760 2647.121
[2,] 3642.296 3509.103 3327.562 3433.451 3393.308 3396.914 3388.460 3446.939
         [,9]    [,10]    [,11]    [,12]
[1,] 2608.109 2178.743 2222.837 2238.116
[2,] 3316.571 3334.777 3475.703 3604.764

$out
[1] 4562.84 4296.49

$group
[1] 7 9

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

Warning message:
In bxp(list(stats = c(1905.41, 2350.44, 2984.475, 3667.03, 4443.91,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5vr7l1243508773.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,] 2767.630 2604.420 1940.490 1635.250 2350.440 2981.850 3570.120 4105.180
[2,] 2872.980 2650.730 2053.905 1871.925 2424.445 3088.945 3751.335 4249.505
[3,] 3012.925 2839.220 2430.130 1964.635 2462.715 3140.500 3878.890 4439.570
[4,] 3039.355 2903.825 2740.790 2106.180 2701.200 3286.165 4054.310 4577.055
[5,] 3142.950 2995.550 2811.700 2196.720 2921.440 3494.170 4290.890 4696.960
        [,9]   [,10]
[1,] 1862.83 1905.41
[2,] 2580.08 1905.41
[3,] 3275.70 1905.41
[4,] 3761.02 1905.41
[5,] 3857.62 1905.41

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

$conf
        [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2937.04 2723.782 2116.837 1857.790 2336.485 3050.547 3740.701 4290.172
[2,] 3088.81 2954.658 2743.423 2071.480 2588.945 3230.453 4017.079 4588.968
         [,9]   [,10]
[1,] 2737.065 1905.41
[2,] 3814.335 1905.41

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(2767.63, 2872.98, 3012.925, 3039.355, 3142.95,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6vxeg1243508773.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,] 710.2586 2141.110  672.590
[2,] 775.4576 2329.370  849.010
[3,] 824.1135 2565.870  990.665
[4,] 839.1050 2729.235 1081.226
[5,] 881.5980 2786.530 1297.450

$n
[1] 12 12 12

$conf
         [,1]     [,2]      [,3]
[1,] 795.0835 2383.489  884.7496
[2,] 853.1435 2748.251 1096.5804

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(710.258649591432, 775.457635675971, 824.113452168745,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/197j31243508773.ps tmp/197j31243508773.png")
> system("convert tmp/2jd911243508773.ps tmp/2jd911243508773.png")
> system("convert tmp/3w4dm1243508773.ps tmp/3w4dm1243508773.png")
> system("convert tmp/4u17b1243508773.ps tmp/4u17b1243508773.png")
> system("convert tmp/5vr7l1243508773.ps tmp/5vr7l1243508773.png")
> system("convert tmp/6vxeg1243508773.ps tmp/6vxeg1243508773.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.052   0.808   1.396