R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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.

  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(255
+ ,280.2
+ ,299.9
+ ,339.2
+ ,374.2
+ ,393.5
+ ,389.2
+ ,381.7
+ ,375.2
+ ,369
+ ,357.4
+ ,352.1
+ ,346.5
+ ,342.9
+ ,340.3
+ ,328.3
+ ,322.9
+ ,314.3
+ ,308.9
+ ,294
+ ,285.6
+ ,281.2
+ ,280.3
+ ,278.8
+ ,274.5
+ ,270.4
+ ,263.4
+ ,259.9
+ ,258
+ ,262.7
+ ,284.7
+ ,311.3
+ ,322.1
+ ,327
+ ,331.3
+ ,333.3
+ ,321.4
+ ,327
+ ,320
+ ,314.7
+ ,316.7
+ ,314.4
+ ,321.3
+ ,318.2
+ ,307.2
+ ,301.3
+ ,287.5
+ ,277.7
+ ,274.4
+ ,258.8
+ ,253.3
+ ,251
+ ,248.4
+ ,249.5
+ ,246.1
+ ,244.5
+ ,243.6
+ ,244
+ ,240.8
+ ,249.8
+ ,248
+ ,259.4
+ ,260.5
+ ,260.8
+ ,261.3
+ ,259.5
+ ,256.6
+ ,257.9
+ ,256.5
+ ,254.2
+ ,253.3
+ ,253.8
+ ,255.5
+ ,257.1
+ ,257.3
+ ,253.2
+ ,252.8
+ ,252
+ ,250.7
+ ,252.2
+ ,250
+ ,251
+ ,253.4
+ ,251.2
+ ,255.6
+ ,261.1
+ ,258.9
+ ,259.9
+ ,261.2
+ ,264.7
+ ,267.1
+ ,266.4
+ ,267.7
+ ,268.6
+ ,267.5
+ ,268.5
+ ,268.5
+ ,270.5
+ ,270.9
+ ,270.1
+ ,269.3
+ ,269.8
+ ,270.1
+ ,264.9
+ ,263.7
+ ,264.8
+ ,263.7
+ ,255.9
+ ,276.2
+ ,360.1
+ ,380.5
+ ,373.7
+ ,369.8
+ ,366.6
+ ,359.3
+ ,345.8
+ ,326.2
+ ,324.5
+ ,328.1
+ ,327.5
+ ,324.4
+ ,316.5
+ ,310.9
+ ,301.5
+ ,291.7
+ ,290.4
+ ,287.4
+ ,277.7
+ ,281.6
+ ,288
+ ,276
+ ,272.9
+ ,283
+ ,283.3
+ ,276.8
+ ,284.5
+ ,282.7
+ ,281.2
+ ,287.4
+ ,283.1
+ ,284
+ ,285.5
+ ,289.2
+ ,292.5
+ ,296.4
+ ,305.2
+ ,303.9
+ ,311.5
+ ,316.3
+ ,316.7
+ ,322.5
+ ,317.1
+ ,309.8
+ ,303.8
+ ,290.3
+ ,293.7
+ ,291.7
+ ,296.5
+ ,289.1
+ ,288.5
+ ,293.8
+ ,297.7
+ ,305.4
+ ,302.7
+ ,302.5
+ ,303
+ ,294.5
+ ,294.1
+ ,294.5
+ ,297.1
+ ,289.4
+ ,292.4
+ ,287.9
+ ,286.6
+ ,280.5
+ ,272.4
+ ,269.2
+ ,270.6
+ ,267.3
+ ,262.5
+ ,266.8
+ ,268.8
+ ,263.1
+ ,261.2
+ ,266
+ ,262.5
+ ,265.2
+ ,261.3
+ ,253.7
+ ,249.2
+ ,239.1
+ ,236.4
+ ,235.2
+ ,245.2
+ ,246.2
+ ,247.7
+ ,251.4
+ ,253.3
+ ,254.8
+ ,250
+ ,249.3
+ ,241.5
+ ,243.3
+ ,248
+ ,253
+ ,252.9
+ ,251.5
+ ,251.6
+ ,253.5
+ ,259.8
+ ,334.1
+ ,448
+ ,445.8
+ ,445
+ ,448.2
+ ,438.2
+ ,439.8
+ ,423.4
+ ,410.8
+ ,408.4
+ ,406.7
+ ,405.9
+ ,402.7
+ ,405.1
+ ,399.6
+ ,386.5
+ ,381.4
+ ,375.2
+ ,357.7
+ ,359
+ ,355
+ ,352.7
+ ,344.4
+ ,343.8
+ ,338
+ ,339
+ ,333.3
+ ,334.4
+ ,328.3
+ ,330.7
+ ,330
+ ,331.6
+ ,351.2
+ ,389.4
+ ,410.9
+ ,442.8
+ ,462.8
+ ,466.9
+ ,461.7
+ ,439.2
+ ,430.3
+ ,416.1
+ ,402.5
+ ,397.3
+ ,403.3
+ ,395.9
+ ,387.8
+ ,378.6
+ ,377.1
+ ,370.4
+ ,362
+ ,350.3
+ ,348.2
+ ,344.6
+ ,343.5
+ ,342.8
+ ,347.6
+ ,346.6
+ ,349.5
+ ,342.1
+ ,342
+ ,342.8
+ ,339.3
+ ,348.2
+ ,333.7
+ ,334.7
+ ,354
+ ,367.7
+ ,363.3
+ ,358.4
+ ,353.1
+ ,343.1
+ ,344.6
+ ,344.4
+ ,333.9
+ ,331.7
+ ,324.3
+ ,321.2
+ ,322.4
+ ,321.7
+ ,320.5
+ ,312.8
+ ,309.7
+ ,315.6
+ ,309.7
+ ,304.6
+ ,302.5
+ ,301.5
+ ,298.8
+ ,291.3
+ ,293.6
+ ,294.6
+ ,285.9
+ ,297.6
+ ,301.1
+ ,293.8
+ ,297.7
+ ,292.9
+ ,292.1
+ ,287.2
+ ,288.2
+ ,283.8
+ ,299.9
+ ,292.4
+ ,293.3
+ ,300.8
+ ,293.7
+ ,293.1
+ ,294.4
+ ,292.1
+ ,291.9
+ ,282.5
+ ,277.9
+ ,287.5
+ ,289.2
+ ,285.6
+ ,293.2
+ ,290.8
+ ,283.1
+ ,275
+ ,287.8
+ ,287.8
+ ,287.4
+ ,284
+ ,277.8
+ ,277.6
+ ,304.9
+ ,294
+ ,300.9
+ ,324
+ ,332.9
+ ,341.6
+ ,333.4
+ ,348.2
+ ,344.7
+ ,344.7
+ ,329.3
+ ,323.5
+ ,323.2
+ ,317.4
+ ,330.1
+ ,329.2
+ ,334.9
+ ,315.8
+ ,315.4
+ ,319.6
+ ,317.3
+ ,313.8
+ ,315.8
+ ,311.3)
> 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] 360
> (np <- floor(n / par1))
[1] 30
> arr <- array(NA,dim=c(par1,np+1))
> darr <- array(NA,dim=c(par1,np+1))
> ari <- array(0,dim=par1)
> dx <- diff(x)
> j <- 0
> for (i in 1:n)
+ {
+ j = j + 1
+ ari[j] = ari[j] + 1
+ arr[j,ari[j]] <- x[i]
+ darr[j,ari[j]] <- dx[i]
+ if (j == par1) j = 0
+ }
> ari
 [1] 30 30 30 30 30 30 30 30 30 30 30 30
> arr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,] 255.0 346.5 274.5 321.4 274.4 248.0 255.5 255.6 268.5 276.2 324.4 283.0
 [2,] 280.2 342.9 270.4 327.0 258.8 259.4 257.1 261.1 270.5 360.1 316.5 283.3
 [3,] 299.9 340.3 263.4 320.0 253.3 260.5 257.3 258.9 270.9 380.5 310.9 276.8
 [4,] 339.2 328.3 259.9 314.7 251.0 260.8 253.2 259.9 270.1 373.7 301.5 284.5
 [5,] 374.2 322.9 258.0 316.7 248.4 261.3 252.8 261.2 269.3 369.8 291.7 282.7
 [6,] 393.5 314.3 262.7 314.4 249.5 259.5 252.0 264.7 269.8 366.6 290.4 281.2
 [7,] 389.2 308.9 284.7 321.3 246.1 256.6 250.7 267.1 270.1 359.3 287.4 287.4
 [8,] 381.7 294.0 311.3 318.2 244.5 257.9 252.2 266.4 264.9 345.8 277.7 283.1
 [9,] 375.2 285.6 322.1 307.2 243.6 256.5 250.0 267.7 263.7 326.2 281.6 284.0
[10,] 369.0 281.2 327.0 301.3 244.0 254.2 251.0 268.6 264.8 324.5 288.0 285.5
[11,] 357.4 280.3 331.3 287.5 240.8 253.3 253.4 267.5 263.7 328.1 276.0 289.2
[12,] 352.1 278.8 333.3 277.7 249.8 253.8 251.2 268.5 255.9 327.5 272.9 292.5
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24]
 [1,] 296.4 291.7 294.5 266.8 235.2 253.0 439.8 357.7 330.0 402.5 343.5 354.0
 [2,] 305.2 296.5 297.1 268.8 245.2 252.9 423.4 359.0 331.6 397.3 342.8 367.7
 [3,] 303.9 289.1 289.4 263.1 246.2 251.5 410.8 355.0 351.2 403.3 347.6 363.3
 [4,] 311.5 288.5 292.4 261.2 247.7 251.6 408.4 352.7 389.4 395.9 346.6 358.4
 [5,] 316.3 293.8 287.9 266.0 251.4 253.5 406.7 344.4 410.9 387.8 349.5 353.1
 [6,] 316.7 297.7 286.6 262.5 253.3 259.8 405.9 343.8 442.8 378.6 342.1 343.1
 [7,] 322.5 305.4 280.5 265.2 254.8 334.1 402.7 338.0 462.8 377.1 342.0 344.6
 [8,] 317.1 302.7 272.4 261.3 250.0 448.0 405.1 339.0 466.9 370.4 342.8 344.4
 [9,] 309.8 302.5 269.2 253.7 249.3 445.8 399.6 333.3 461.7 362.0 339.3 333.9
[10,] 303.8 303.0 270.6 249.2 241.5 445.0 386.5 334.4 439.2 350.3 348.2 331.7
[11,] 290.3 294.5 267.3 239.1 243.3 448.2 381.4 328.3 430.3 348.2 333.7 324.3
[12,] 293.7 294.1 262.5 236.4 248.0 438.2 375.2 330.7 416.1 344.6 334.7 321.2
      [,25] [,26] [,27] [,28] [,29] [,30] [,31]
 [1,] 322.4 293.6 299.9 289.2 304.9 323.2    NA
 [2,] 321.7 294.6 292.4 285.6 294.0 317.4    NA
 [3,] 320.5 285.9 293.3 293.2 300.9 330.1    NA
 [4,] 312.8 297.6 300.8 290.8 324.0 329.2    NA
 [5,] 309.7 301.1 293.7 283.1 332.9 334.9    NA
 [6,] 315.6 293.8 293.1 275.0 341.6 315.8    NA
 [7,] 309.7 297.7 294.4 287.8 333.4 315.4    NA
 [8,] 304.6 292.9 292.1 287.8 348.2 319.6    NA
 [9,] 302.5 292.1 291.9 287.4 344.7 317.3    NA
[10,] 301.5 287.2 282.5 284.0 344.7 313.8    NA
[11,] 298.8 288.2 277.9 277.8 329.3 315.8    NA
[12,] 291.3 283.8 287.5 277.6 323.5 311.3    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13]
 [1,]  25.2  -3.6  -4.1   5.6 -15.6 11.4  1.6  5.5  2.0  83.9  -7.9   0.3   8.8
 [2,]  19.7  -2.6  -7.0  -7.0  -5.5  1.1  0.2 -2.2  0.4  20.4  -5.6  -6.5  -1.3
 [3,]  39.3 -12.0  -3.5  -5.3  -2.3  0.3 -4.1  1.0 -0.8  -6.8  -9.4   7.7   7.6
 [4,]  35.0  -5.4  -1.9   2.0  -2.6  0.5 -0.4  1.3 -0.8  -3.9  -9.8  -1.8   4.8
 [5,]  19.3  -8.6   4.7  -2.3   1.1 -1.8 -0.8  3.5  0.5  -3.2  -1.3  -1.5   0.4
 [6,]  -4.3  -5.4  22.0   6.9  -3.4 -2.9 -1.3  2.4  0.3  -7.3  -3.0   6.2   5.8
 [7,]  -7.5 -14.9  26.6  -3.1  -1.6  1.3  1.5 -0.7 -5.2 -13.5  -9.7  -4.3  -5.4
 [8,]  -6.5  -8.4  10.8 -11.0  -0.9 -1.4 -2.2  1.3 -1.2 -19.6   3.9   0.9  -7.3
 [9,]  -6.2  -4.4   4.9  -5.9   0.4 -2.3  1.0  0.9  1.1  -1.7   6.4   1.5  -6.0
[10,] -11.6  -0.9   4.3 -13.8  -3.2 -0.9  2.4 -1.1 -1.1   3.6 -12.0   3.7 -13.5
[11,]  -5.3  -1.5   2.0  -9.8   9.0  0.5 -2.2  1.0 -7.8  -0.6  -3.1   3.3   3.4
[12,]  -5.6  -4.3 -11.9  -3.3  -1.8  1.7  4.4  0.0 20.3  -3.1  10.1   3.9  -2.0
      [,14] [,15] [,16] [,17] [,18] [,19] [,20] [,21] [,22] [,23] [,24] [,25]
 [1,]   4.8   2.6   2.0  10.0  -0.1 -16.4   1.3   1.6  -5.2  -0.7  13.7  -0.7
 [2,]  -7.4  -7.7  -5.7   1.0  -1.4 -12.6  -4.0  19.6   6.0   4.8  -4.4  -1.2
 [3,]  -0.6   3.0  -1.9   1.5   0.1  -2.4  -2.3  38.2  -7.4  -1.0  -4.9  -7.7
 [4,]   5.3  -4.5   4.8   3.7   1.9  -1.7  -8.3  21.5  -8.1   2.9  -5.3  -3.1
 [5,]   3.9  -1.3  -3.5   1.9   6.3  -0.8  -0.6  31.9  -9.2  -7.4 -10.0   5.9
 [6,]   7.7  -6.1   2.7   1.5  74.3  -3.2  -5.8  20.0  -1.5  -0.1   1.5  -5.9
 [7,]  -2.7  -8.1  -3.9  -4.8 113.9   2.4   1.0   4.1  -6.7   0.8  -0.2  -5.1
 [8,]  -0.2  -3.2  -7.6  -0.7  -2.2  -5.5  -5.7  -5.2  -8.4  -3.5 -10.5  -2.1
 [9,]   0.5   1.4  -4.5  -7.8  -0.8 -13.1   1.1 -22.5 -11.7   8.9  -2.2  -1.0
[10,]  -8.5  -3.3 -10.1   1.8   3.2  -5.1  -6.1  -8.9  -2.1 -14.5  -7.4  -2.7
[11,]  -0.4  -4.8  -2.7   4.7 -10.0  -6.2   2.4 -14.2  -3.6   1.0  -3.1  -7.5
[12,]   0.4   4.3  -1.2   5.0   1.6 -17.5  -0.7 -13.6  -1.1  19.3   1.2   2.3
      [,26] [,27] [,28] [,29] [,30] [,31]
 [1,]   1.0  -7.5  -3.6 -10.9  -5.8    NA
 [2,]  -8.7   0.9   7.6   6.9  12.7    NA
 [3,]  11.7   7.5  -2.4  23.1  -0.9    NA
 [4,]   3.5  -7.1  -7.7   8.9   5.7    NA
 [5,]  -7.3  -0.6  -8.1   8.7 -19.1    NA
 [6,]   3.9   1.3  12.8  -8.2  -0.4    NA
 [7,]  -4.8  -2.3   0.0  14.8   4.2    NA
 [8,]  -0.8  -0.2  -0.4  -3.5  -2.3    NA
 [9,]  -4.9  -9.4  -3.4   0.0  -3.5    NA
[10,]   1.0  -4.6  -6.2 -15.4   2.0    NA
[11,]  -4.4   9.6  -0.2  -5.8  -4.5    NA
[12,]  16.1   1.7  27.3  -0.3    NA    NA
> arr.mean <- array(NA,dim=par1)
> arr.median <- array(NA,dim=par1)
> arr.midrange <- array(NA,dim=par1)
> for (j in 1:par1)
+ {
+ arr.mean[j] <- mean(arr[j,],na.rm=TRUE)
+ arr.median[j] <- median(arr[j,],na.rm=TRUE)
+ arr.midrange[j] <- (quantile(arr[j,],0.75,na.rm=TRUE) + quantile(arr[j,],0.25,na.rm=TRUE)) / 2
+ }
> overall.mean <- mean(x)
> overall.median <- median(x)
> overall.midrange <- (quantile(x,0.75) + quantile(x,0.25)) / 2
> postscript(file="/var/www/html/freestat/rcomp/tmp/1d1aw1290714274.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.mean,type='b',ylab='mean',main='Mean Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.mean,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/2naaz1290714274.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.median,type='b',ylab='median',main='Median Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.median,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/3naaz1290714274.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> plot(arr.midrange,type='b',ylab='midrange',main='Midrange Plot',xlab='Periodic Index')
> mtext(paste('#blocks = ',np))
> abline(overall.midrange,0)
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/4g19k1290714274.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] [,10]  [,11]
[1,] 235.20 245.20 246.2 247.70 248.40 249.50 246.10 244.50 243.6 241.5 239.10
[2,] 268.50 270.40 263.4 261.20 266.00 264.70 280.50 272.40 269.2 270.6 267.50
[3,] 294.05 295.55 296.6 301.15 297.45 295.75 307.15 303.65 302.5 301.4 289.75
[4,] 324.40 331.60 340.3 339.20 344.40 342.10 338.00 344.40 333.9 334.4 329.30
[5,] 402.50 423.40 410.8 408.40 410.90 442.80 402.70 448.00 399.6 386.5 381.40
     [,12]
[1,] 236.4
[2,] 268.5
[3,] 291.9
[4,] 330.7
[5,] 416.1

$n
 [1] 30 30 30 30 30 30 30 30 30 30 30 30

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 277.9247 277.8958 274.4169 278.6496 274.8342 273.4226 290.5631 282.8804
[2,] 310.1753 313.2042 318.7831 323.6504 320.0658 318.0774 323.7369 324.4196
         [,9]    [,10]    [,11]    [,12]
[1,] 283.8362 282.9958 271.9227 273.9573
[2,] 321.1638 319.8042 307.5773 309.8427

$out
 [1] 439.8 462.8 466.9 445.8 461.7 445.0 439.2 448.2 430.3 438.2

$group
 [1]  1  7  8  9  9 10 10 11 11 12

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

> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5g19k1290714274.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> z <- data.frame(t(darr))
> names(z) <- c(1:par1)
> (boxplot(z,notch=TRUE,col='grey',xlab='Periodic Index',ylab='Value',main='Notched Box Plots - Differenced Periodic Subseries'))
$stats
       [,1]   [,2]   [,3] [,4]  [,5] [,6]  [,7]   [,8]   [,9]  [,10]  [,11]
[1,] -16.40 -12.60 -12.00 -9.8 -10.0 -8.2 -14.9 -11.00 -13.10 -15.40 -14.20
[2,]  -4.10  -5.70  -4.10 -4.5  -3.5 -3.4  -5.2  -6.50  -5.90  -8.90  -5.30
[3,]   1.15  -1.35  -0.95 -0.6  -0.8  0.1  -2.5  -2.25  -1.95  -3.25  -2.45
[4,]   5.50   4.80   3.00  3.7   3.5  5.8   1.3  -0.70   1.00   1.00   1.00
[5,]  13.70  20.40  11.70  8.9   8.7 12.8   4.2   3.90   8.90   4.30   9.60
     [,12]
[1,]  -5.6
[2,]  -2.0
[3,]   0.4
[4,]   4.3
[5,]  10.1

$n
 [1] 30 30 30 30 30 30 30 30 30 30 30 29

$conf
          [,1]      [,2]      [,3]      [,4]      [,5]      [,6]       [,7]
[1,] -1.619285 -4.378906 -2.998117 -2.965431 -2.819270 -2.553898 -4.3750369
[2,]  3.919285  1.678906  1.098117  1.765431  1.219270  2.753898 -0.6249631
           [,8]        [,9]      [,10]      [,11]     [,12]
[1,] -3.9231098 -3.94042377 -6.1058254 -4.2673434 -1.448411
[2,] -0.5768902  0.04042377 -0.3941746 -0.6326566  2.248411

$out
 [1]  25.2  83.9  39.3  38.2  23.1  35.0  21.5  19.3  31.9 -19.1  22.0  74.3
[13]  20.0  26.6 113.9  14.8  10.8 -19.6 -22.5 -11.9  20.3 -17.5 -13.6  19.3
[25]  16.1  27.3

$group
 [1]  1  1  3  3  3  4  4  5  5  5  6  6  6  7  7  7  8  8  9 12 12 12 12 12 12
[26] 12

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

Warning message:
In bxp(list(stats = c(-16.4000000000000, -4.10000000000002, 1.15000000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/69bq51290714274.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]  [,10]
[1,] 255.00 278.8 258.00 287.50 240.80 248.00 250.00 255.60 255.90 276.20
[2,] 319.55 283.4 263.05 304.25 244.25 254.00 251.10 260.50 264.25 326.85
[3,] 363.20 311.6 279.60 315.70 248.95 257.25 252.50 265.55 268.90 352.55
[4,] 378.45 334.3 324.55 320.65 252.15 260.00 254.45 267.60 270.10 368.20
[5,] 393.50 346.5 333.30 327.00 258.80 261.30 257.30 268.60 270.90 380.50
      [,11]  [,12] [,13]  [,14]  [,15]  [,16]  [,17]  [,18]  [,19]  [,20]
[1,] 272.90 281.20 290.3 288.50 262.50 236.40 235.20 251.50 375.20 328.30
[2,] 279.65 282.85 300.1 292.75 269.90 251.45 244.25 252.95 393.05 333.85
[3,] 289.20 283.65 307.5 295.50 283.55 261.90 247.85 296.95 405.50 341.40
[4,] 306.20 286.45 316.5 302.60 290.90 265.60 250.70 445.40 409.60 353.85
[5,] 324.40 289.20 322.5 305.40 297.10 268.80 254.80 448.20 423.40 359.00
      [,21]  [,22]  [,23] [,24]  [,25]  [,26]  [,27]  [,28]  [,29]  [,30] [,31]
[1,] 330.00 344.60 333.70 321.2 291.30 283.80 287.50 275.00 294.00 311.30    NA
[2,] 370.30 356.15 340.65 332.8 302.00 287.70 289.70 280.45 314.20 315.60    NA
[3,] 423.20 377.85 342.80 344.5 309.70 293.25 292.75 286.50 331.10 317.35    NA
[4,] 452.25 396.60 347.10 356.2 318.05 296.10 294.05 288.50 343.15 326.20    NA
[5,] 466.90 403.30 349.50 367.7 322.40 301.10 299.90 293.20 348.20 334.90    NA

$n
 [1] 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12 12
[26] 12 12 12 12 12  0

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 336.3353 288.3842 251.5494 308.2198 245.3468 254.5134 250.9720 262.3116
[2,] 390.0647 334.8158 307.6506 323.1802 252.5532 259.9866 254.0280 268.7884
         [,9]  [,10]    [,11]   [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 266.2318 333.69 277.0904 282.008 300.0198 291.0073 273.9718 255.4461
[2,] 271.5682 371.41 301.3096 285.292 314.9802 299.9927 293.1282 268.3539
        [,17]    [,18]    [,19]    [,20]    [,21]    [,22]    [,23]    [,24]
[1,] 244.9081 209.1723 397.9514 332.2779 385.8221 359.4005 339.8581 333.8271
[2,] 250.7919 384.7277 413.0486 350.5221 460.5779 396.2995 345.7419 355.1729
        [,25]    [,26]    [,27]    [,28]    [,29]    [,30] [,31]
[1,] 302.3795 289.4187 290.7659 282.8283 317.8957 312.5153    NA
[2,] 317.0205 297.0813 294.7341 290.1717 344.3043 322.1847    NA

$out
[1] 277.7 274.4 276.8 292.5 439.8 300.8 282.5 277.9

$group
[1]  4  5 12 12 19 27 27 27

$names
 [1] "1"  "2"  "3"  "4"  "5"  "6"  "7"  "8"  "9"  "10" "11" "12" "13" "14" "15"
[16] "16" "17" "18" "19" "20" "21" "22" "23" "24" "25" "26" "27" "28" "29" "30"
[31] NA  

Warning message:
In bxp(list(stats = c(255, 319.55, 363.2, 378.45, 393.5, 278.8,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/79bq51290714274.ps",horizontal=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
> z <- data.frame(cbind(arr.mean,arr.median,arr.midrange))
> names(z) <- list('mean','median','midrange')
> (boxplot(z,notch=TRUE,col='grey',ylab='Overall Central Tendency',main='Notched Box Plots'))
$stats
         [,1]    [,2]     [,3]
[1,] 302.7100 289.750 297.0375
[2,] 305.4283 294.800 299.9062
[3,] 308.8750 297.025 302.2687
[4,] 310.7633 301.950 304.2000
[5,] 315.4333 307.150 309.2875

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 306.4417 293.7638 300.3103
[2,] 311.3083 300.2862 304.2272

$out
numeric(0)

$group
numeric(0)

$names
[1] "mean"     "median"   "midrange"

Warning message:
In bxp(list(stats = c(302.71, 305.428333333333, 308.875, 310.763333333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1d1aw1290714274.ps tmp/1d1aw1290714274.png",intern=TRUE))
character(0)
> try(system("convert tmp/2naaz1290714274.ps tmp/2naaz1290714274.png",intern=TRUE))
character(0)
> try(system("convert tmp/3naaz1290714274.ps tmp/3naaz1290714274.png",intern=TRUE))
character(0)
> try(system("convert tmp/4g19k1290714274.ps tmp/4g19k1290714274.png",intern=TRUE))
character(0)
> try(system("convert tmp/5g19k1290714274.ps tmp/5g19k1290714274.png",intern=TRUE))
character(0)
> try(system("convert tmp/69bq51290714274.ps tmp/69bq51290714274.png",intern=TRUE))
character(0)
> try(system("convert tmp/79bq51290714274.ps tmp/79bq51290714274.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.047   1.488   2.290