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(2430.47
+ ,2516.3
+ ,2633.63
+ ,2799.84
+ ,3001.93
+ ,3229.29
+ ,3173.02
+ ,3322.08
+ ,3417.88
+ ,3486.95
+ ,3016.22
+ ,2709.61
+ ,2914.87
+ ,3203.08
+ ,3320.25
+ ,3446.25
+ ,3456.85
+ ,3566.53
+ ,3763.67
+ ,3607.75
+ ,3747.38
+ ,3623.91
+ ,3699.76
+ ,3629.61
+ ,3911.52
+ ,4281.47
+ ,4742.42
+ ,4522.42
+ ,4879.79
+ ,5059.11
+ ,5093.19
+ ,4941.81
+ ,4832.67
+ ,4876.18
+ ,5018.07
+ ,4780.34
+ ,4953.59
+ ,4622.32
+ ,4557.13
+ ,4560.03
+ ,4105.66
+ ,4004.89
+ ,4277.26
+ ,4245.98
+ ,4057.64
+ ,3931.42
+ ,3637.15
+ ,3339.91
+ ,3465.74
+ ,3571.25
+ ,3706.93
+ ,3584.17
+ ,3552.11
+ ,3695.24
+ ,3510
+ ,3357.7
+ ,3060.91
+ ,2736.98
+ ,2709.45
+ ,2314.96
+ ,2561.29
+ ,2663.49
+ ,2407.87
+ ,2237.74
+ ,2165.44
+ ,2098.89
+ ,2318.54
+ ,2315.49
+ ,2395.47
+ ,2474.07
+ ,2479.57
+ ,2386.92
+ ,2537.84
+ ,2567.13
+ ,2660.37
+ ,2696.28
+ ,2748.5
+ ,2663.32
+ ,2707.69
+ ,2669.36
+ ,2687.68
+ ,2650.24
+ ,2620.03
+ ,2668.47
+ ,2692.06
+ ,2737.67
+ ,2774.77
+ ,2819.19
+ ,2892.56
+ ,2866.08
+ ,2817.41
+ ,2934.75
+ ,3036.54
+ ,3139.5
+ ,3114.31
+ ,3261.3
+ ,3201.79
+ ,3264.53
+ ,3349.1
+ ,3446.17
+ ,3469.48
+ ,3507.13
+ ,3536.2
+ ,3359.05
+ ,3378.85
+ ,3449.15
+ ,3522.89
+ ,3551.04
+ ,3669.15
+ ,3602
+ ,3697.22
+ ,3760.9
+ ,3665.08
+ ,3708.8
+ ,3858.21
+ ,3933.16
+ ,3946.98
+ ,3794.29
+ ,3765.56
+ ,3820.33
+ ,3885.12
+ ,3752.67
+ ,3683.79
+ ,3240.75
+ ,3188.82
+ ,3017.98
+ ,3237.2
+ ,3182.53
+ ,2906.42
+ ,2881.35
+ ,2915.64
+ ,2635.13
+ ,2331.43
+ ,2159.04
+ ,2065.46
+ ,1983.48
+ ,1770.41
+ ,1815.99
+ ,2026.97
+ ,2124.81
+ ,2098.28
+ ,2291.39
+ ,2401.57
+ ,2453.89
+ ,2409.53)
> 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]
 [1,] 2430.47 2914.87 3911.52 4953.59 3465.74 2561.29 2537.84 2692.06 3201.79
 [2,] 2516.30 3203.08 4281.47 4622.32 3571.25 2663.49 2567.13 2737.67 3264.53
 [3,] 2633.63 3320.25 4742.42 4557.13 3706.93 2407.87 2660.37 2774.77 3349.10
 [4,] 2799.84 3446.25 4522.42 4560.03 3584.17 2237.74 2696.28 2819.19 3446.17
 [5,] 3001.93 3456.85 4879.79 4105.66 3552.11 2165.44 2748.50 2892.56 3469.48
 [6,] 3229.29 3566.53 5059.11 4004.89 3695.24 2098.89 2663.32 2866.08 3507.13
 [7,] 3173.02 3763.67 5093.19 4277.26 3510.00 2318.54 2707.69 2817.41 3536.20
 [8,] 3322.08 3607.75 4941.81 4245.98 3357.70 2315.49 2669.36 2934.75 3359.05
 [9,] 3417.88 3747.38 4832.67 4057.64 3060.91 2395.47 2687.68 3036.54 3378.85
[10,] 3486.95 3623.91 4876.18 3931.42 2736.98 2474.07 2650.24 3139.50 3449.15
[11,] 3016.22 3699.76 5018.07 3637.15 2709.45 2479.57 2620.03 3114.31 3522.89
[12,] 2709.61 3629.61 4780.34 3339.91 2314.96 2386.92 2668.47 3261.30 3551.04
        [,10]   [,11]   [,12]   [,13]
 [1,] 3669.15 3885.12 2331.43 2409.53
 [2,] 3602.00 3752.67 2159.04      NA
 [3,] 3697.22 3683.79 2065.46      NA
 [4,] 3760.90 3240.75 1983.48      NA
 [5,] 3665.08 3188.82 1770.41      NA
 [6,] 3708.80 3017.98 1815.99      NA
 [7,] 3858.21 3237.20 2026.97      NA
 [8,] 3933.16 3182.53 2124.81      NA
 [9,] 3946.98 2906.42 2098.28      NA
[10,] 3794.29 2881.35 2291.39      NA
[11,] 3765.56 2915.64 2401.57      NA
[12,] 3820.33 2635.13 2453.89      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/1605f1260967149.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/2reqt1260967149.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/3mdop1260967149.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/4e9zk1260967149.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]
[1,] 2331.43 2159.040 2065.460 1983.480 1770.410 1815.99 2026.97 2124.810
[2,] 2537.84 2615.310 2647.000 2748.060 2820.530 2764.70 2762.55 2802.055
[3,] 2914.87 3233.805 3334.675 3343.460 3322.835 3368.21 3373.60 3339.890
[4,] 3669.15 3677.335 3702.075 3672.535 3608.595 3702.02 3810.94 3770.455
[5,] 4953.59 4622.320 4742.420 4560.030 4105.660 5059.11 5093.19 4941.810
        [,9]    [,10]    [,11]    [,12]
[1,] 2098.28 2291.390 2401.570 2314.960
[2,] 2797.05 2693.610 2664.740 2544.510
[3,] 3219.88 3294.325 3065.265 2985.455
[4,] 3847.18 3709.100 3668.455 3590.325
[5,] 4832.67 4876.180 5018.070 4780.340

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 2419.115 2749.408 2853.448 2921.801 2963.393 2940.692 2895.422 2898.196
[2,] 3410.625 3718.202 3815.902 3765.119 3682.277 3795.728 3851.778 3781.584
         [,9]    [,10]    [,11]    [,12]
[1,] 2740.909 2831.153 2607.464 2508.452
[2,] 3698.851 3757.497 3523.066 3462.458

$out
[1] 4879.79

$group
[1] 5

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

Warning message:
In bxp(list(stats = c(2331.43, 2537.84, 2914.87, 3669.15, 4953.59,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/502af1260967149.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,] 2430.470 3203.080 4281.470 3339.910 2314.960 2098.890 2567.130 2692.060
[2,] 2671.620 3383.250 4632.420 3968.155 2898.945 2276.615 2635.135 2796.090
[3,] 3009.075 3587.140 4854.425 4175.820 3487.870 2391.195 2665.895 2879.320
[4,] 3275.685 3664.685 4979.940 4558.580 3577.710 2476.820 2691.980 3075.425
[5,] 3486.950 3763.670 5093.190 4953.590 3706.930 2663.490 2748.500 3261.300
         [,9]    [,10]    [,11]    [,12]   [,13]
[1,] 3201.790 3602.000 2635.130 1770.410 2409.53
[2,] 3354.075 3683.185 2911.030 2005.225 2409.53
[3,] 3447.660 3763.230 3185.675 2111.545 2409.53
[4,] 3515.010 3839.270 3462.270 2311.410 2409.53
[5,] 3551.040 3946.980 3885.120 2453.890 2409.53

$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,] 2733.557 3458.776 4695.919 3906.523 3178.281 2299.88 2639.968 2751.913
[2,] 3284.593 3715.504 5012.931 4445.117 3797.459 2482.51 2691.822 3006.727
         [,9]    [,10]    [,11]    [,12]   [,13]
[1,] 3374.256 3692.039 2934.251 1971.892 2409.53
[2,] 3521.064 3834.421 3437.099 2251.198 2409.53

$out
[1] 2914.87 3911.52 2537.84

$group
[1] 2 3 7

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

Warning message:
In bxp(list(stats = c(2430.47, 2671.62, 3009.075, 3275.685, 3486.95,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6wmud1260967149.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,] 733.8132 2463.28  723.8075
[2,] 745.7258 2580.67  868.8212
[3,] 792.5941 2649.56  958.4588
[4,] 829.4428 2941.61  998.7963
[5,] 870.7644 3243.12 1131.3100

$n
[1] 12 12 12

$conf
         [,1]     [,2]      [,3]
[1,] 754.4102 2484.933  899.1763
[2,] 830.7779 2814.187 1017.7412

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(733.813201795624, 745.725756990264, 792.59405059232,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1605f1260967149.ps tmp/1605f1260967149.png",intern=TRUE))
character(0)
> try(system("convert tmp/2reqt1260967149.ps tmp/2reqt1260967149.png",intern=TRUE))
character(0)
> try(system("convert tmp/3mdop1260967149.ps tmp/3mdop1260967149.png",intern=TRUE))
character(0)
> try(system("convert tmp/4e9zk1260967149.ps tmp/4e9zk1260967149.png",intern=TRUE))
character(0)
> try(system("convert tmp/502af1260967149.ps tmp/502af1260967149.png",intern=TRUE))
character(0)
> try(system("convert tmp/6wmud1260967149.ps tmp/6wmud1260967149.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.044   0.784   1.773