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(227.86
+ ,198.24
+ ,194.97
+ ,184.88
+ ,196.79
+ ,205.36
+ ,226.72
+ ,226.05
+ ,202.50
+ ,194.79
+ ,192.43
+ ,219.25
+ ,217.47
+ ,192.34
+ ,196.83
+ ,186.07
+ ,197.31
+ ,215.02
+ ,242.67
+ ,225.17
+ ,206.69
+ ,197.75
+ ,196.43
+ ,213.55
+ ,222.75
+ ,194.03
+ ,201.85
+ ,189.50
+ ,206.07
+ ,225.59
+ ,247.91
+ ,247.64
+ ,213.01
+ ,203.01
+ ,200.26
+ ,220.50
+ ,237.90
+ ,216.94
+ ,214.01
+ ,196.00
+ ,208.37
+ ,232.75
+ ,257.46
+ ,267.69
+ ,220.18
+ ,210.61
+ ,209.59
+ ,232.75
+ ,232.75
+ ,219.82
+ ,226.74
+ ,208.04
+ ,220.12
+ ,235.69
+ ,257.05
+ ,258.69
+ ,227.15
+ ,219.91
+ ,219.30
+ ,259.04
+ ,237.29
+ ,212.88
+ ,226.03
+ ,211.07
+ ,222.91
+ ,249.18
+ ,266.38
+ ,268.53
+ ,238.02
+ ,224.69
+ ,213.75
+ ,237.43
+ ,248.46
+ ,210.82
+ ,221.40
+ ,209.00
+ ,234.37
+ ,248.43
+ ,271.98
+ ,268.11
+ ,233.88
+ ,223.43
+ ,221.38
+ ,233.76
+ ,243.97
+ ,217.76
+ ,224.66
+ ,210.84
+ ,220.35
+ ,236.84
+ ,266.15
+ ,255.20
+ ,234.76
+ ,221.29
+ ,221.26
+ ,244.13
+ ,245.78
+ ,224.62
+ ,234.80
+ ,211.37
+ ,222.39
+ ,249.63
+ ,282.29
+ ,279.13
+ ,236.60
+ ,223.62
+ ,225.86
+ ,246.41
+ ,261.70
+ ,225.01
+ ,231.54
+ ,214.82
+ ,227.70
+ ,263.86
+ ,278.15
+ ,274.64
+ ,237.66
+ ,227.97
+ ,224.75
+ ,242.91
+ ,253.08
+ ,228.13
+ ,233.68
+ ,217.38
+ ,236.38
+ ,256.08
+ ,292.83
+ ,304.71
+ ,245.57
+ ,234.41
+ ,234.12
+ ,258.17
+ ,268.66
+ ,245.31
+ ,247.47
+ ,226.25
+ ,251.67
+ ,268.79
+ ,288.94
+ ,290.16
+ ,250.69
+ ,240.80)
> 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] 142
> (np <- floor(n / par1))
[1] 11
> 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] 12 12 12 12 12 12 12 12 12 12 11 11
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 227.86 217.47 222.75 237.90 232.75 237.29 248.46 243.97 245.78 261.70
 [2,] 198.24 192.34 194.03 216.94 219.82 212.88 210.82 217.76 224.62 225.01
 [3,] 194.97 196.83 201.85 214.01 226.74 226.03 221.40 224.66 234.80 231.54
 [4,] 184.88 186.07 189.50 196.00 208.04 211.07 209.00 210.84 211.37 214.82
 [5,] 196.79 197.31 206.07 208.37 220.12 222.91 234.37 220.35 222.39 227.70
 [6,] 205.36 215.02 225.59 232.75 235.69 249.18 248.43 236.84 249.63 263.86
 [7,] 226.72 242.67 247.91 257.46 257.05 266.38 271.98 266.15 282.29 278.15
 [8,] 226.05 225.17 247.64 267.69 258.69 268.53 268.11 255.20 279.13 274.64
 [9,] 202.50 206.69 213.01 220.18 227.15 238.02 233.88 234.76 236.60 237.66
[10,] 194.79 197.75 203.01 210.61 219.91 224.69 223.43 221.29 223.62 227.97
[11,] 192.43 196.43 200.26 209.59 219.30 213.75 221.38 221.26 225.86 224.75
[12,] 219.25 213.55 220.50 232.75 259.04 237.43 233.76 244.13 246.41 242.91
       [,11]  [,12]
 [1,] 253.08 268.66
 [2,] 228.13 245.31
 [3,] 233.68 247.47
 [4,] 217.38 226.25
 [5,] 236.38 251.67
 [6,] 256.08 268.79
 [7,] 292.83 288.94
 [8,] 304.71 290.16
 [9,] 245.57 250.69
[10,] 234.41 240.80
[11,] 234.12     NA
[12,] 258.17     NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] -29.62 -25.13 -28.72 -20.96 -12.93 -24.41 -37.64 -26.21 -21.16 -36.69
 [2,]  -3.27   4.49   7.82  -2.93   6.92  13.15  10.58   6.90  10.18   6.53
 [3,] -10.09 -10.76 -12.35 -18.01 -18.70 -14.96 -12.40 -13.82 -23.43 -16.72
 [4,]  11.91  11.24  16.57  12.37  12.08  11.84  25.37   9.51  11.02  12.88
 [5,]   8.57  17.71  19.52  24.38  15.57  26.27  14.06  16.49  27.24  36.16
 [6,]  21.36  27.65  22.32  24.71  21.36  17.20  23.55  29.31  32.66  14.29
 [7,]  -0.67 -17.50  -0.27  10.23   1.64   2.15  -3.87 -10.95  -3.16  -3.51
 [8,] -23.55 -18.48 -34.63 -47.51 -31.54 -30.51 -34.23 -20.44 -42.53 -36.98
 [9,]  -7.71  -8.94 -10.00  -9.57  -7.24 -13.33 -10.45 -13.47 -12.98  -9.69
[10,]  -2.36  -1.32  -2.75  -1.02  -0.61 -10.94  -2.05  -0.03   2.24  -3.22
[11,]  26.82  17.12  20.24  23.16  39.74  23.68  12.38  22.87  20.55  18.16
[12,]  -1.78   9.20  17.40   0.00 -21.75  11.03  10.21   1.65  15.29  10.17
       [,11]  [,12]
 [1,] -24.95 -23.35
 [2,]   5.55   2.16
 [3,] -16.30 -21.22
 [4,]  19.00  25.42
 [5,]  19.70  17.12
 [6,]  36.75  20.15
 [7,]  11.88   1.22
 [8,] -59.14 -39.47
 [9,] -11.16  -9.89
[10,]  -0.29     NA
[11,]  24.05     NA
[12,]  10.49     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/rcomp/tmp/17gwa1272299002.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/rcomp/tmp/2i7ev1272299002.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/rcomp/tmp/3i7ev1272299002.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/rcomp/tmp/4byvy1272299002.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,] 217.470 192.340 194.970 184.880 196.790 205.360 226.720 225.170 202.500
[2,] 230.305 204.530 207.930 192.750 207.220 229.170 252.480 251.420 216.595
[3,] 240.935 217.350 225.345 209.920 221.370 242.635 266.265 267.900 234.320
[4,] 250.770 224.815 232.610 213.095 231.035 252.855 280.220 276.885 237.840
[5,] 268.660 245.310 247.470 226.250 251.670 268.790 292.830 304.710 250.690
      [,10]   [,11]   [,12]
[1,] 194.79 192.430 213.550
[2,] 206.81 204.925 226.625
[3,] 222.36 219.300 237.430
[4,] 226.33 223.065 245.270
[5,] 240.80 234.120 259.040

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 231.6008 208.0979 214.0883 200.6405 210.5078 231.8321 253.6126 256.2852
[2,] 250.2692 226.6021 236.6017 219.1995 232.2322 253.4379 278.9174 279.5148
       [,9]    [,10]    [,11]    [,12]
[1,] 224.63 213.4568 210.6583 228.5477
[2,] 244.01 231.2632 227.9417 246.3123

$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(217.47, 230.305, 240.935, 250.77, 268.66,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5byvy1272299002.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]
[1,] -37.640 -3.270 -23.430  9.510  8.570 14.290 -10.950 -59.14 -13.470 -3.220
[2,] -29.170  3.325 -18.355 11.540 16.030 20.755  -3.690 -41.00 -12.070 -2.555
[3,] -25.040  6.715 -15.630 12.225 18.615 22.935  -0.470 -34.43  -9.945 -1.320
[4,] -22.255  9.000 -12.375 17.785 25.325 28.480   1.895 -27.03  -9.255 -0.450
[5,] -12.930 13.150 -10.090 25.420 36.160 36.750  10.230 -18.48  -7.240  2.240
      [,11]  [,12]
[1,] 12.380 -1.780
[2,] 19.200  0.825
[3,] 22.870 10.170
[4,] 23.865 10.760
[5,] 26.820 17.400

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

$conf
          [,1]     [,2]      [,3]      [,4]     [,5]     [,6]      [,7]
[1,] -28.19398 4.126594 -18.35752  9.376614 14.37549 19.41158 -3.017356
[2,] -21.88602 9.303406 -12.90248 15.073386 22.85451 26.45842  2.077356
          [,8]      [,9]      [,10]    [,11]     [,12]
[1,] -40.80181 -11.22894 -2.3227966 20.64765  5.437086
[2,] -28.05819  -8.66106 -0.3172034 25.09235 14.902914

$out
[1] -17.50  11.88 -10.94  39.74 -21.75

$group
[1]  7  7 10 11 12

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

Warning message:
In bxp(list(stats = c(-37.64, -29.17, -25.04, -22.255, -12.93, -3.27000000000001,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6byvy1272299002.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,] 184.88 186.070 189.500 196.000 208.040 211.07 209.000 210.840 211.37
[2,] 194.88 196.630 201.055 210.100 219.865 218.33 221.390 220.805 224.12
[3,] 200.37 202.220 209.540 218.560 226.945 231.66 233.820 229.710 235.70
[4,] 222.65 216.245 224.170 235.325 246.370 243.60 248.445 244.050 248.02
[5,] 227.86 242.670 247.910 267.690 259.040 268.53 271.980 266.150 282.29
       [,10]   [,11]  [,12]
[1,] 214.820 217.380 226.25
[2,] 226.355 233.900 245.31
[3,] 234.600 240.975 251.18
[4,] 262.780 257.125 268.79
[5,] 278.150 258.170 290.16

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 187.7039 193.2735 198.9971 207.0547 214.8559 220.1342 221.4800 219.1078
[2,] 213.0361 211.1665 220.0829 230.0653 239.0341 243.1858 246.1600 240.3122
         [,9]    [,10]    [,11]    [,12]
[1,] 224.7990 217.9863 230.3819 239.4485
[2,] 246.6010 251.2137 251.5681 262.9115

$out
[1] 292.83 304.71

$group
[1] 11 11

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

Warning message:
In bxp(list(stats = c(184.88, 194.88, 200.37, 222.65, 227.86, 186.07,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7l7c11272299002.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,] 205.4350 209.9200 203.3038
[2,] 217.0075 220.3350 216.6531
[3,] 225.0288 229.8325 224.7956
[4,] 241.0371 241.7850 240.8363
[5,] 264.8775 267.9000 266.9750

$n
[1] 12 12 12

$conf
         [,1]    [,2]     [,3]
[1,] 214.0687 220.049 213.7655
[2,] 235.9888 239.616 235.8257

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(205.435, 217.0075, 225.02875, 241.037083333333,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/17gwa1272299002.ps tmp/17gwa1272299002.png",intern=TRUE))
character(0)
> try(system("convert tmp/2i7ev1272299002.ps tmp/2i7ev1272299002.png",intern=TRUE))
character(0)
> try(system("convert tmp/3i7ev1272299002.ps tmp/3i7ev1272299002.png",intern=TRUE))
character(0)
> try(system("convert tmp/4byvy1272299002.ps tmp/4byvy1272299002.png",intern=TRUE))
character(0)
> try(system("convert tmp/5byvy1272299002.ps tmp/5byvy1272299002.png",intern=TRUE))
character(0)
> try(system("convert tmp/6byvy1272299002.ps tmp/6byvy1272299002.png",intern=TRUE))
character(0)
> try(system("convert tmp/7l7c11272299002.ps tmp/7l7c11272299002.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.203   1.004   1.596