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(1213.8,1245.6,1306.3,1255.8,1257.6,1287.8,1300.4,1320.9,1370.8,1327.3,1320,1345.3,1346.7,1395.4,1462,1491.6,1461.8,1477.9,1490.3,1521.1,1561.9,1552.6,1523.6,1548.3,1552.4,1587,1621.3,1648.7,1641.8,1650.6,1688.6,1670.7,1682.2,1678.9,1650.6,1662.4,1664.5,1683.2,1736.2,1747.6,1749,1759.7,1793.6,1817.4,1858.4,1839.9,1809.1,1877.7,1880.3,1930.9,2039.3,1992.7,1987.8,1984.4,2016.5,2016.7,2064.1,2031.5,2000.3,2057.8,2041.2,2093.2,2158.3,2128.8,2131.9,2170.3,2190.8,2217.7,2254.4,2223.3,2210.5,2250.8,2249.1,2288.6,2329.2,2313.8,2309.8,2345.9,2361.3,2372,2410.4,2398.5,2362.3,2419.1,2421.6,2465,2480.5,2506.1,2506.6,2525.8,2550,2578.3,2807.8,2815.3,2767.7,2815.4,2838.8,2864,2948.6,2922.8,2917.2,2936.8,2993.4,3007.8,3046.3,3011.5,2958.6,3019.8,2998.5,3040.4,3166,3110,3099.2,3150.3,3163.6,3182.6,3244.4,3223.2,3143.6,3217,3182.3,3217.2,3262.5,3227.9,3171.6,3219,3195.4,3221.6,3262.1,3179.5,3133.6,3219.2,3245,3265.3,3312.5,3383.6,3386.3,3411.1,3467.2,3487.7,3575.5,3571.5,3582.3,3637.1,3685)
> 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] 145
> (np <- floor(n / par1))
[1] 12
> 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] 13 12 12 12 12 12 12 12 12 12 12 12
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 1213.8 1346.7 1552.4 1664.5 1880.3 2041.2 2249.1 2421.6 2838.8 2998.5
 [2,] 1245.6 1395.4 1587.0 1683.2 1930.9 2093.2 2288.6 2465.0 2864.0 3040.4
 [3,] 1306.3 1462.0 1621.3 1736.2 2039.3 2158.3 2329.2 2480.5 2948.6 3166.0
 [4,] 1255.8 1491.6 1648.7 1747.6 1992.7 2128.8 2313.8 2506.1 2922.8 3110.0
 [5,] 1257.6 1461.8 1641.8 1749.0 1987.8 2131.9 2309.8 2506.6 2917.2 3099.2
 [6,] 1287.8 1477.9 1650.6 1759.7 1984.4 2170.3 2345.9 2525.8 2936.8 3150.3
 [7,] 1300.4 1490.3 1688.6 1793.6 2016.5 2190.8 2361.3 2550.0 2993.4 3163.6
 [8,] 1320.9 1521.1 1670.7 1817.4 2016.7 2217.7 2372.0 2578.3 3007.8 3182.6
 [9,] 1370.8 1561.9 1682.2 1858.4 2064.1 2254.4 2410.4 2807.8 3046.3 3244.4
[10,] 1327.3 1552.6 1678.9 1839.9 2031.5 2223.3 2398.5 2815.3 3011.5 3223.2
[11,] 1320.0 1523.6 1650.6 1809.1 2000.3 2210.5 2362.3 2767.7 2958.6 3143.6
[12,] 1345.3 1548.3 1662.4 1877.7 2057.8 2250.8 2419.1 2815.4 3019.8 3217.0
       [,11]  [,12] [,13]
 [1,] 3182.3 3245.0  3685
 [2,] 3217.2 3265.3    NA
 [3,] 3262.5 3312.5    NA
 [4,] 3227.9 3383.6    NA
 [5,] 3171.6 3386.3    NA
 [6,] 3219.0 3411.1    NA
 [7,] 3195.4 3467.2    NA
 [8,] 3221.6 3487.7    NA
 [9,] 3262.1 3575.5    NA
[10,] 3179.5 3571.5    NA
[11,] 3133.6 3582.3    NA
[12,] 3219.2 3637.1    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/165ff1260696155.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/2xcqs1260696155.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/38uw51260696155.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/4mxl91260696155.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,] 1213.8 1245.6 1306.30 1255.80 1257.60 1287.80 1300.40 1320.90 1370.80
[2,] 1664.5 1635.1 1678.75 1698.15 1695.40 1705.15 1741.10 1744.05 1770.30
[3,] 2249.1 2190.9 2243.75 2221.30 2220.85 2258.10 2276.05 2294.85 2332.40
[4,] 2998.5 2952.2 3057.30 3016.40 3008.20 3043.55 3078.50 3095.20 3145.35
[5,] 3685.0 3265.3 3312.50 3383.60 3386.30 3411.10 3467.20 3487.70 3575.50
      [,10]   [,11]   [,12]
[1,] 1327.3 1320.00 1345.30
[2,] 1759.4 1729.85 1770.05
[3,] 2310.9 2286.40 2334.95
[4,] 3095.5 3046.10 3118.40
[5,] 3571.5 3582.30 3637.10

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 1664.524 1590.162 1614.984 1620.037 1622.073 1647.647 1666.053 1678.581
[2,] 2833.676 2791.638 2872.516 2822.563 2819.627 2868.553 2886.047 2911.119
         [,9]    [,10]    [,11]    [,12]
[1,] 1705.230 1701.496 1686.050 1719.959
[2,] 2959.570 2920.304 2886.750 2949.941

$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(1213.8, 1664.5, 2249.1, 2998.5, 3685, 1245.6,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5zrku1260696155.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]    [,9]
[1,] 1213.80 1395.40 1587.00 1664.50 1930.9 2041.20 2249.10 2421.60 2838.80
[2,] 1256.70 1461.90 1631.55 1741.90 1986.1 2130.35 2311.80 2493.30 2920.00
[3,] 1303.35 1490.95 1650.60 1776.65 2008.4 2180.55 2353.60 2537.90 2953.60
[4,] 1324.10 1535.95 1674.80 1828.65 2035.4 2220.50 2385.25 2787.75 3009.65
[5,] 1370.80 1561.90 1688.60 1877.70 2064.1 2254.40 2419.10 2815.40 3046.30
       [,10]   [,11]   [,12] [,13]
[1,] 2998.50 3133.60 3245.00  3685
[2,] 3104.60 3180.90 3348.05  3685
[3,] 3156.95 3218.10 3439.15  3685
[4,] 3199.80 3224.75 3573.50  3685
[5,] 3244.40 3262.50 3637.10  3685

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 1272.608 1457.175 1630.873 1737.083 1985.914 2139.432 2320.099 2403.599
[2,] 1334.092 1524.725 1670.327 1816.217 2030.886 2221.668 2387.101 2672.201
        [,9]    [,10]  [,11]    [,12] [,13]
[1,] 2912.71 3113.529 3198.1 3336.321  3685
[2,] 2994.49 3200.371 3238.1 3541.979  3685

$out
[1] 1346.7 1552.4 1880.3

$group
[1] 2 3 5

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

Warning message:
In bxp(list(stats = c(1213.8, 1256.7, 1303.35, 1324.1, 1370.8, 1395.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/67fxn1260696155.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,] 716.1562 2006.20 1232.875
[2,] 720.6409 2125.55 1245.975
[3,] 727.7527 2166.80 1255.800
[4,] 744.3781 2253.25 1276.125
[5,] 751.6040 2291.80 1295.475

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 716.9260 2108.555 1242.048
[2,] 738.5794 2225.045 1269.552

$out
[1]  799.2848 2471.2000 1334.0000

$group
[1] 1 2 3

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

Warning message:
In bxp(list(stats = c(716.156220683897, 720.640922033498, 727.752719413209,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/165ff1260696155.ps tmp/165ff1260696155.png",intern=TRUE))
character(0)
> try(system("convert tmp/2xcqs1260696155.ps tmp/2xcqs1260696155.png",intern=TRUE))
character(0)
> try(system("convert tmp/38uw51260696155.ps tmp/38uw51260696155.png",intern=TRUE))
character(0)
> try(system("convert tmp/4mxl91260696155.ps tmp/4mxl91260696155.png",intern=TRUE))
character(0)
> try(system("convert tmp/5zrku1260696155.ps tmp/5zrku1260696155.png",intern=TRUE))
character(0)
> try(system("convert tmp/67fxn1260696155.ps tmp/67fxn1260696155.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.083   0.833   1.448