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(87.28
+ ,87.28
+ ,87.09
+ ,86.92
+ ,87.59
+ ,90.72
+ ,90.69
+ ,90.3
+ ,89.55
+ ,88.94
+ ,88.41
+ ,87.82
+ ,87.07
+ ,86.82
+ ,86.4
+ ,86.02
+ ,85.66
+ ,85.32
+ ,85
+ ,84.67
+ ,83.94
+ ,82.83
+ ,81.95
+ ,81.19
+ ,80.48
+ ,78.86
+ ,69.47
+ ,68.77
+ ,70.06
+ ,73.95
+ ,75.8
+ ,77.79
+ ,81.57
+ ,83.07
+ ,84.34
+ ,85.1
+ ,85.25
+ ,84.26
+ ,83.63
+ ,86.44
+ ,85.3
+ ,84.1
+ ,83.36
+ ,82.48
+ ,81.58
+ ,80.47
+ ,79.34
+ ,82.13
+ ,81.69
+ ,80.7
+ ,79.88
+ ,79.16
+ ,78.38
+ ,77.42
+ ,76.47
+ ,75.46
+ ,74.48
+ ,78.27
+ ,80.7
+ ,79.91
+ ,78.75
+ ,77.78
+ ,81.14
+ ,81.08
+ ,80.03
+ ,78.91
+ ,78.01
+ ,76.9
+ ,75.97
+ ,81.93
+ ,80.27
+ ,78.67
+ ,77.42
+ ,76.16
+ ,74.7
+ ,76.39
+ ,76.04
+ ,74.65
+ ,73.29
+ ,71.79
+ ,74.39
+ ,74.91
+ ,74.54
+ ,73.08
+ ,72.75
+ ,71.32
+ ,70.38
+ ,70.35
+ ,70.01
+ ,69.36
+ ,67.77
+ ,69.26
+ ,69.8
+ ,68.38
+ ,67.62
+ ,68.39
+ ,66.95
+ ,65.21
+ ,66.64
+ ,63.45
+ ,60.66
+ ,62.34
+ ,60.32
+ ,58.64
+ ,60.46
+ ,58.59
+ ,61.87
+ ,61.85
+ ,67.44
+ ,77.06
+ ,91.74
+ ,93.15
+ ,94.15
+ ,93.11
+ ,91.51
+ ,89.96
+ ,88.16
+ ,86.98
+ ,88.03
+ ,86.24
+ ,84.65
+ ,83.23
+ ,81.7
+ ,80.25
+ ,78.8
+ ,77.51
+ ,76.2
+ ,75.04
+ ,74
+ ,75.49
+ ,77.14
+ ,76.15
+ ,76.27
+ ,78.19
+ ,76.49
+ ,77.31
+ ,76.65
+ ,74.99
+ ,73.51
+ ,72.07
+ ,70.59
+ ,71.96
+ ,76.29
+ ,74.86
+ ,74.93
+ ,71.9
+ ,71.01
+ ,77.47
+ ,75.78
+ ,76.6
+ ,76.07
+ ,74.57
+ ,73.02
+ ,72.65
+ ,73.16
+ ,71.53
+ ,69.78
+ ,67.98
+ ,69.96
+ ,72.16
+ ,70.47
+ ,68.86
+ ,67.37
+ ,65.87
+ ,72.16
+ ,71.34
+ ,69.93
+ ,68.44
+ ,67.16
+ ,66.01
+ ,67.25
+ ,70.91
+ ,69.75
+ ,68.59
+ ,67.48
+ ,66.31
+ ,64.81
+ ,66.58
+ ,65.97
+ ,64.7
+ ,64.7
+ ,60.94
+ ,59.08
+ ,58.42
+ ,57.77
+ ,57.11
+ ,53.31
+ ,49.96
+ ,49.4
+ ,48.84
+ ,48.3
+ ,47.74
+ ,47.24
+ ,46.76
+ ,46.29
+ ,48.9
+ ,49.23
+ ,48.53
+ ,48.03
+ ,54.34
+ ,53.79
+ ,53.24
+ ,52.96
+ ,52.17
+ ,51.7
+ ,58.55
+ ,78.2
+ ,77.03
+ ,76.19
+ ,77.15
+ ,75.87
+ ,95.47
+ ,109.67
+ ,112.28
+ ,112.01
+ ,107.93
+ ,105.96
+ ,105.06
+ ,102.98
+ ,102.2
+ ,105.23
+ ,101.85
+ ,99.89
+ ,96.23
+ ,94.76
+ ,91.51
+ ,91.63
+ ,91.54
+ ,85.23
+ ,87.83
+ ,87.38
+ ,84.44
+ ,85.19
+ ,84.03
+ ,86.73
+ ,102.52
+ ,104.45
+ ,106.98
+ ,107.02
+ ,99.26
+ ,94.45
+ ,113.44
+ ,157.33
+ ,147.38
+ ,171.89
+ ,171.95
+ ,132.71
+ ,126.02
+ ,121.18
+ ,115.45
+ ,110.48
+ ,117.85
+ ,117.63
+ ,124.65
+ ,109.59
+ ,111.27
+ ,99.78
+ ,98.21
+ ,99.2
+ ,97.97
+ ,89.55
+ ,87.91
+ ,93.34
+ ,94.42
+ ,93.2
+ ,90.29
+ ,91.46
+ ,89.98
+ ,88.35
+ ,88.41
+ ,82.44
+ ,79.89
+ ,75.69
+ ,75.66
+ ,84.5
+ ,96.73
+ ,87.48
+ ,82.39
+ ,83.48
+ ,79.31
+ ,78.16
+ ,72.77
+ ,72.45
+ ,68.46
+ ,67.62
+ ,68.76
+ ,70.07
+ ,68.55
+ ,65.3
+ ,58.96
+ ,59.17
+ ,62.37
+ ,66.28
+ ,55.62
+ ,55.23
+ ,55.85
+ ,56.75
+ ,50.89
+ ,53.88
+ ,52.95
+ ,55.08
+ ,53.61
+ ,58.78
+ ,61.85
+ ,55.91
+ ,53.32
+ ,46.41
+ ,44.57
+ ,50
+ ,50
+ ,53.36
+ ,46.23
+ ,50.45
+ ,49.07
+ ,45.85
+ ,48.45
+ ,49.96
+ ,46.53
+ ,50.51
+ ,47.58
+ ,48.05
+ ,46.84
+ ,47.67
+ ,49.16
+ ,55.54
+ ,55.82
+ ,58.22
+ ,56.19
+ ,57.77
+ ,63.19
+ ,54.76
+ ,55.74
+ ,62.54
+ ,61.39
+ ,69.6
+ ,79.23
+ ,80
+ ,93.68
+ ,107.63
+ ,100.18
+ ,97.3
+ ,90.45
+ ,80.64
+ ,80.58
+ ,75.82
+ ,85.59
+ ,89.35
+ ,89.42
+ ,104.73
+ ,95.32
+ ,89.27
+ ,90.44
+ ,86.97
+ ,79.98
+ ,81.22
+ ,87.35
+ ,83.64
+ ,82.22
+ ,94.4
+ ,102.18)
> 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,] 87.28 87.07 80.48 85.25 81.69 78.75 77.42 72.75 66.95 67.44 84.65 76.27
 [2,] 87.28 86.82 78.86 84.26 80.70 77.78 76.16 71.32 65.21 77.06 83.23 78.19
 [3,] 87.09 86.40 69.47 83.63 79.88 81.14 74.70 70.38 66.64 91.74 81.70 76.49
 [4,] 86.92 86.02 68.77 86.44 79.16 81.08 76.39 70.35 63.45 93.15 80.25 77.31
 [5,] 87.59 85.66 70.06 85.30 78.38 80.03 76.04 70.01 60.66 94.15 78.80 76.65
 [6,] 90.72 85.32 73.95 84.10 77.42 78.91 74.65 69.36 62.34 93.11 77.51 74.99
 [7,] 90.69 85.00 75.80 83.36 76.47 78.01 73.29 67.77 60.32 91.51 76.20 73.51
 [8,] 90.30 84.67 77.79 82.48 75.46 76.90 71.79 69.26 58.64 89.96 75.04 72.07
 [9,] 89.55 83.94 81.57 81.58 74.48 75.97 74.39 69.80 60.46 88.16 74.00 70.59
[10,] 88.94 82.83 83.07 80.47 78.27 81.93 74.91 68.38 58.59 86.98 75.49 71.96
[11,] 88.41 81.95 84.34 79.34 80.70 80.27 74.54 67.62 61.87 88.03 77.14 76.29
[12,] 87.82 81.19 85.10 82.13 79.91 78.67 73.08 68.39 61.85 86.24 76.15 74.86
      [,13] [,14] [,15] [,16] [,17]  [,18]  [,19]  [,20]  [,21]  [,22] [,23]
 [1,] 74.93 69.78 67.16 64.70 47.24  51.70 105.96  85.23  94.45 117.63 93.20
 [2,] 71.90 67.98 66.01 60.94 46.76  58.55 105.06  87.83 113.44 124.65 90.29
 [3,] 71.01 69.96 67.25 59.08 46.29  78.20 102.98  87.38 157.33 109.59 91.46
 [4,] 77.47 72.16 70.91 58.42 48.90  77.03 102.20  84.44 147.38 111.27 89.98
 [5,] 75.78 70.47 69.75 57.77 49.23  76.19 105.23  85.19 171.89  99.78 88.35
 [6,] 76.60 68.86 68.59 57.11 48.53  77.15 101.85  84.03 171.95  98.21 88.41
 [7,] 76.07 67.37 67.48 53.31 48.03  75.87  99.89  86.73 132.71  99.20 82.44
 [8,] 74.57 65.87 66.31 49.96 54.34  95.47  96.23 102.52 126.02  97.97 79.89
 [9,] 73.02 72.16 64.81 49.40 53.79 109.67  94.76 104.45 121.18  89.55 75.69
[10,] 72.65 71.34 66.58 48.84 53.24 112.28  91.51 106.98 115.45  87.91 75.66
[11,] 73.16 69.93 65.97 48.30 52.96 112.01  91.63 107.02 110.48  93.34 84.50
[12,] 71.53 68.44 64.70 47.74 52.17 107.93  91.54  99.26 117.85  94.42 96.73
      [,24] [,25] [,26] [,27] [,28]  [,29]  [,30] [,31]
 [1,] 87.48 65.30 55.08 50.45 55.54  80.00 104.73    NA
 [2,] 82.39 58.96 53.61 49.07 55.82  93.68  95.32    NA
 [3,] 83.48 59.17 58.78 45.85 58.22 107.63  89.27    NA
 [4,] 79.31 62.37 61.85 48.45 56.19 100.18  90.44    NA
 [5,] 78.16 66.28 55.91 49.96 57.77  97.30  86.97    NA
 [6,] 72.77 55.62 53.32 46.53 63.19  90.45  79.98    NA
 [7,] 72.45 55.23 46.41 50.51 54.76  80.64  81.22    NA
 [8,] 68.46 55.85 44.57 47.58 55.74  80.58  87.35    NA
 [9,] 67.62 56.75 50.00 48.05 62.54  75.82  83.64    NA
[10,] 68.76 50.89 50.00 46.84 61.39  85.59  82.22    NA
[11,] 70.07 53.88 53.36 47.67 69.60  89.35  94.40    NA
[12,] 68.55 52.95 46.23 49.16 79.23  89.42 102.18    NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]  [,8]  [,9] [,10] [,11] [,12]
 [1,]  0.00 -0.25 -1.62 -0.99 -0.99 -0.97 -1.26 -1.43 -1.74  9.62 -1.42  1.92
 [2,] -0.19 -0.42 -9.39 -0.63 -0.82  3.36 -1.46 -0.94  1.43 14.68 -1.53 -1.70
 [3,] -0.17 -0.38 -0.70  2.81 -0.72 -0.06  1.69 -0.03 -3.19  1.41 -1.45  0.82
 [4,]  0.67 -0.36  1.29 -1.14 -0.78 -1.05 -0.35 -0.34 -2.79  1.00 -1.45 -0.66
 [5,]  3.13 -0.34  3.89 -1.20 -0.96 -1.12 -1.39 -0.65  1.68 -1.04 -1.29 -1.66
 [6,] -0.03 -0.32  1.85 -0.74 -0.95 -0.90 -1.36 -1.59 -2.02 -1.60 -1.31 -1.48
 [7,] -0.39 -0.33  1.99 -0.88 -1.01 -1.11 -1.50  1.49 -1.68 -1.55 -1.16 -1.44
 [8,] -0.75 -0.73  3.78 -0.90 -0.98 -0.93  2.60  0.54  1.82 -1.80 -1.04 -1.48
 [9,] -0.61 -1.11  1.50 -1.11  3.79  5.96  0.52 -1.42 -1.87 -1.18  1.49  1.37
[10,] -0.53 -0.88  1.27 -1.13  2.43 -1.66 -0.37 -0.76  3.28  1.05  1.65  4.33
[11,] -0.59 -0.76  0.76  2.79 -0.79 -1.60 -1.46  0.77 -0.02 -1.79 -0.99 -1.43
[12,] -0.75 -0.71  0.15 -0.44 -1.16 -1.25 -0.33 -1.44  5.59 -1.59  0.12  0.07
      [,13] [,14] [,15] [,16] [,17] [,18] [,19] [,20]  [,21]  [,22] [,23] [,24]
 [1,] -3.03 -1.80 -1.15 -3.76 -0.48  6.85 -0.90  2.60  18.99   7.02 -2.91 -5.09
 [2,] -0.89  1.98  1.24 -1.86 -0.47 19.65 -2.08 -0.45  43.89 -15.06  1.17  1.09
 [3,]  6.46  2.20  3.66 -0.66  2.61 -1.17 -0.78 -2.94  -9.95   1.68 -1.48 -4.17
 [4,] -1.69 -1.69 -1.16 -0.65  0.33 -0.84  3.03  0.75  24.51 -11.49 -1.63 -1.15
 [5,]  0.82 -1.61 -1.16 -0.66 -0.70  0.96 -3.38 -1.16   0.06  -1.57  0.06 -5.39
 [6,] -0.53 -1.49 -1.11 -3.80 -0.50 -1.28 -1.96  2.70 -39.24   0.99 -5.97 -0.32
 [7,] -1.50 -1.50 -1.17 -3.35  6.31 19.60 -3.66 15.79  -6.69  -1.23 -2.55 -3.99
 [8,] -1.55  6.29 -1.50 -0.56 -0.55 14.20 -1.47  1.93  -4.84  -8.42 -4.20 -0.84
 [9,] -0.37 -0.82  1.77 -0.56 -0.55  2.61 -3.25  2.53  -5.73  -1.64 -0.03  1.14
[10,]  0.51 -1.41 -0.61 -0.54 -0.28 -0.27  0.12  0.04  -4.97   5.43  8.84  1.31
[11,] -1.63 -1.49 -1.27 -0.56 -0.79 -4.08 -0.09 -7.76   7.37   1.08 12.23 -1.52
[12,] -1.75 -1.28  0.00 -0.50 -0.47 -1.97 -6.31 -4.81  -0.22  -1.22 -9.25 -3.25
       [,25] [,26] [,27] [,28] [,29] [,30] [,31]
 [1,]  -6.34 -1.47 -1.38  0.28 13.68 -9.41    NA
 [2,]   0.21  5.17 -3.22  2.40 13.95 -6.05    NA
 [3,]   3.20  3.07  2.60 -2.03 -7.45  1.17    NA
 [4,]   3.91 -5.94  1.51  1.58 -2.88 -3.47    NA
 [5,] -10.66 -2.59 -3.43  5.42 -6.85 -6.99    NA
 [6,]  -0.39 -6.91  3.98 -8.43 -9.81  1.24    NA
 [7,]   0.62 -1.84 -2.93  0.98 -0.06  6.13    NA
 [8,]   0.90  5.43  0.47  6.80 -4.76 -3.71    NA
 [9,]  -5.86  0.00 -1.21 -1.15  9.77 -1.42    NA
[10,]   2.99  3.36  0.83  8.21  3.76 12.18    NA
[11,]  -0.93 -7.13  1.49  9.63  0.07  7.78    NA
[12,]   2.13  4.22  6.38  0.77 15.31    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/1mkv01290693121.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/2mkv01290693121.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/3xbcl1290693121.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/4xbcl1290693121.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,]  47.240  46.760  45.85  48.450  49.230  46.530 46.41  44.57  48.050
[2,]  66.950  65.210  67.25  68.770  69.750  68.590 67.37  65.87  64.810
[3,]  78.085  77.985  79.04  78.315  77.405  76.875 75.97  75.25  74.435
[4,]  87.070  87.280  87.38  86.920  86.970  85.320 83.36  87.35  83.940
[5,] 105.960 113.440 109.59 111.270 105.230 101.850 99.89 102.52 109.670
       [,10]  [,11]  [,12]
[1,]  46.840  47.67  46.23
[2,]  66.580  67.62  68.39
[3,]  75.575  78.24  78.95
[4,]  85.590  88.41  89.42
[5,] 112.280 112.01 117.85

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

$conf
         [,1]     [,2]     [,3]     [,4]    [,5]     [,6]     [,7]     [,8]
[1,] 72.28104 71.61853 73.23315 73.07932 72.4376 72.04894 71.35741 69.05372
[2,] 83.88896 84.35147 84.84685 83.55068 82.3724 81.70106 80.58259 81.44628
         [,9]    [,10]    [,11]    [,12]
[1,] 68.91662 70.09124 72.24277 72.88353
[2,] 79.95338 81.05876 84.23723 85.01647

$out
 [1] 117.63 124.65 157.33 147.38 171.89 171.95 132.71 126.02 121.18 115.45

$group
 [1]  1  2  3  4  5  6  7  8  9 10

$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/5xbcl1290693121.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,] -3.76 -6.050 -4.170 -3.47 -3.43 -3.800 -3.990 -4.840 -3.250 -4.97 -4.080
[2,] -1.74 -1.530 -1.450 -1.63 -1.66 -1.960 -1.680 -1.500 -1.210 -0.54 -1.490
[3,] -1.07 -0.435 -0.115 -0.72 -1.14 -1.195 -1.165 -0.795 -0.555  0.67 -0.775
[4,]  0.28  1.980  2.200  0.75  0.06 -0.320  0.620  1.820  1.490  3.28  0.770
[5,]  2.60  5.170  6.460  3.91  1.68  1.850  1.990  6.800  3.790  8.84  2.790
     [,12]
[1,] -3.25
[2,] -1.44
[3,] -0.50
[4,]  0.12
[5,]  2.13

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

$conf
           [,1]       [,2]       [,3]        [,4]       [,5]       [,6]
[1,] -1.6527038 -1.4475199 -1.1679053 -1.40655197 -1.6361636 -1.6680862
[2,] -0.4872962  0.5775199  0.9379053 -0.03344803 -0.6438364 -0.7219138
           [,7]       [,8]       [,9]      [,10]      [,11]       [,12]
[1,] -1.8284746 -1.7527111 -1.3338615 -0.4319448 -1.4269359 -0.95770187
[2,] -0.5015254  0.1627111  0.2238615  1.7719448 -0.1230641 -0.04229813

$out
 [1]   9.62   6.85  18.99   7.02  -5.09  -6.34  13.68  -9.41  -9.39  14.68
[11]  19.65  43.89 -15.06  13.95  -9.95  -7.45  24.51 -11.49  -5.94   3.13
[21]   3.89  -5.39 -10.66   5.42  -6.85  -6.99   2.70 -39.24  -5.97  -6.91
[31]   3.98  -8.43  -9.81   6.31  19.60  15.79  -6.69   6.13  14.20  -8.42
[41]   5.96  -5.73  -5.86   9.77  12.18  -7.76   7.37  12.23  -7.13   9.63
[51]   7.78   5.59  -6.31  -4.81  -9.25   4.22   6.38  15.31

$group
 [1]  1  1  1  1  1  1  1  1  2  2  2  2  2  2  3  3  4  4  4  5  5  5  5  5  5
[26]  5  6  6  6  6  6  6  6  7  7  7  7  7  8  8  9  9  9  9 10 11 11 11 11 11
[51] 11 12 12 12 12 12 12 12

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

Warning message:
In bxp(list(stats = c(-3.76000000000001, -1.74000000000001, -1.06999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6qkb61290693121.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,] 86.920 81.190 68.770 79.340 74.480 75.970 71.790 67.620 58.59 86.240
[2,] 87.280 83.385 72.005 81.855 76.945 77.895 73.840 68.385 60.39 86.610
[3,] 88.115 85.160 78.325 83.495 78.770 78.830 74.675 69.580 61.86 89.060
[4,] 89.925 86.210 82.320 84.755 80.305 80.675 76.100 70.365 64.33 92.425
[5,] 90.720 87.070 85.100 86.440 81.690 81.930 77.420 72.750 66.95 94.150
      [,11] [,12]  [,13]  [,14]  [,15] [,16]  [,17]   [,18]   [,19]   [,20]
[1,] 74.000 70.59 71.010 65.870 64.700 47.74 46.290  51.700  91.510  84.030
[2,] 75.820 72.79 72.275 68.210 65.990 49.12 47.635  76.030  93.195  85.210
[3,] 77.325 75.63 73.865 69.855 66.870 55.21 49.065  77.675 100.870  87.605
[4,] 80.975 76.57 75.925 70.905 68.035 58.75 53.100 108.800 104.020 103.485
[5,] 84.650 78.19 77.470 72.160 70.910 64.70 54.340 112.280 105.960 107.020
       [,21]   [,22]  [,23]  [,24]  [,25]  [,26] [,27]  [,28]   [,29]  [,30]
[1,]  94.450  87.910 75.660 67.620 50.890 44.570 45.85 54.760  75.820  79.98
[2,] 114.445  93.880 81.165 68.655 54.555 48.205 47.21 55.780  80.610  82.93
[3,] 123.600  98.705 88.380 72.610 56.300 53.340 48.25 57.995  89.385  88.31
[4,] 152.355 110.430 90.875 80.850 60.770 55.495 49.56 62.865  95.490  94.86
[5,] 171.950 124.650 96.730 87.480 66.280 61.850 50.51 69.600 107.630 104.73
     [,31]
[1,]    NA
[2,]    NA
[3,]    NA
[4,]    NA
[5,]    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,] 86.9086 83.8715 73.62026 82.17229 77.23748 77.56202 73.6442 68.67691
[2,] 89.3214 86.4485 83.02974 84.81771 80.30252 80.09798 75.7058 70.48309
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14]    [,15]    [,16]
[1,] 60.06294 86.40774 74.97377 73.90592 72.20021 68.62579 65.93726 50.81769
[2,] 63.65706 91.71226 79.67623 77.35408 75.52979 71.08421 67.80274 59.60231
        [,17]    [,18]     [,19]    [,20]   [,21]     [,22]   [,23]    [,24]
[1,] 46.57238 62.72838  95.93264 79.26965 106.309  91.15643 83.9512 67.04778
[2,] 51.55762 92.62162 105.80736 95.94035 140.891 106.25357 92.8088 78.17222
       [,25]    [,26]    [,27]    [,28]    [,29]    [,30] [,31]
[1,] 53.4653 50.01498 47.17815 54.76348 82.59813 82.86865    NA
[2,] 59.1347 56.66502 49.32185 61.22652 96.17187 93.75135    NA

$out
[1] 67.44 77.06 79.23

$group
[1] 10 10 28

$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(86.92, 87.28, 88.115, 89.925, 90.72, 81.19,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7qkb61290693121.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,] 75.40833 74.4350 74.68875
[2,] 75.95567 75.7725 76.14125
[3,] 77.62100 77.6950 76.82250
[4,] 78.84733 78.2775 78.09875
[5,] 79.73967 79.0400 78.71125

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 76.30209 76.55245 75.92967
[2,] 78.93991 78.83755 77.71533

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(75.4083333333333, 75.9556666666667, 77.621,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1mkv01290693121.ps tmp/1mkv01290693121.png",intern=TRUE))
character(0)
> try(system("convert tmp/2mkv01290693121.ps tmp/2mkv01290693121.png",intern=TRUE))
character(0)
> try(system("convert tmp/3xbcl1290693121.ps tmp/3xbcl1290693121.png",intern=TRUE))
character(0)
> try(system("convert tmp/4xbcl1290693121.ps tmp/4xbcl1290693121.png",intern=TRUE))
character(0)
> try(system("convert tmp/5xbcl1290693121.ps tmp/5xbcl1290693121.png",intern=TRUE))
character(0)
> try(system("convert tmp/6qkb61290693121.ps tmp/6qkb61290693121.png",intern=TRUE))
character(0)
> try(system("convert tmp/7qkb61290693121.ps tmp/7qkb61290693121.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.028   1.481   2.247