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.

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(26.663
+ ,23.598
+ ,26.931
+ ,24.740
+ ,25.806
+ ,24.364
+ ,24.477
+ ,23.901
+ ,23.175
+ ,23.227
+ ,21.672
+ ,21.870
+ ,21.439
+ ,21.089
+ ,23.709
+ ,21.669
+ ,21.752
+ ,20.761
+ ,23.479
+ ,23.824
+ ,23.105
+ ,23.110
+ ,21.759
+ ,22.073
+ ,21.937
+ ,20.035
+ ,23.590
+ ,21.672
+ ,22.222
+ ,22.123
+ ,23.950
+ ,23.504
+ ,22.238
+ ,23.142
+ ,21.059
+ ,21.573
+ ,21.548
+ ,20.000
+ ,22.424
+ ,20.615
+ ,21.761
+ ,22.874
+ ,24.104
+ ,23.748
+ ,23.262
+ ,22.907
+ ,21.519
+ ,22.025
+ ,22.604
+ ,20.894
+ ,24.677
+ ,23.673
+ ,25.320
+ ,23.583
+ ,24.671
+ ,24.454
+ ,24.122
+ ,24.252
+ ,22.084
+ ,22.991
+ ,23.287
+ ,23.049
+ ,25.076
+ ,24.037
+ ,24.430
+ ,24.667
+ ,26.451
+ ,25.618
+ ,25.014
+ ,25.110
+ ,22.964
+ ,23.981
+ ,23.798
+ ,22.270
+ ,24.775
+ ,22.646
+ ,23.988
+ ,24.737
+ ,26.276
+ ,25.816
+ ,25.210
+ ,25.199
+ ,23.162
+ ,24.707
+ ,24.364
+ ,22.644
+ ,25.565
+ ,24.062
+ ,25.431
+ ,24.635
+ ,27.009
+ ,26.606
+ ,26.268
+ ,26.462
+ ,25.246
+ ,25.180
+ ,24.657
+ ,23.304
+ ,26.982
+ ,26.199
+ ,27.210
+ ,26.122
+ ,26.706
+ ,26.878
+ ,26.152
+ ,26.379
+ ,24.712
+ ,25.688
+ ,24.990
+ ,24.239
+ ,26.721
+ ,23.475
+ ,24.767
+ ,26.219
+ ,28.361
+ ,28.599
+ ,27.914
+ ,27.784
+ ,25.693
+ ,26.881
+ ,26.217
+ ,24.218
+ ,27.914
+ ,26.975
+ ,28.527
+ ,27.139
+ ,28.982
+ ,28.169
+ ,28.056
+ ,29.136
+ ,26.291
+ ,26.987
+ ,26.589
+ ,24.848
+ ,27.543
+ ,26.896
+ ,28.878
+ ,27.390
+ ,28.065
+ ,28.141
+ ,29.048
+ ,28.484
+ ,26.634
+ ,27.735
+ ,27.132
+ ,24.924
+ ,28.963
+ ,26.589
+ ,27.931
+ ,28.009
+ ,29.229
+ ,28.759
+ ,28.405
+ ,27.945
+ ,25.912
+ ,26.619
+ ,26.076
+ ,25.286
+ ,27.660
+ ,25.951
+ ,26.398
+ ,25.565
+ ,28.865
+ ,30.000
+ ,29.261
+ ,29.012
+ ,26.992
+ ,27.897)
> 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] 168
> (np <- floor(n / par1))
[1] 14
> 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] 14 14 14 14 14 14 14 14 14 14 14 14
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] 26.663 21.439 21.937 21.548 22.604 23.287 23.798 24.364 24.657 24.990
 [2,] 23.598 21.089 20.035 20.000 20.894 23.049 22.270 22.644 23.304 24.239
 [3,] 26.931 23.709 23.590 22.424 24.677 25.076 24.775 25.565 26.982 26.721
 [4,] 24.740 21.669 21.672 20.615 23.673 24.037 22.646 24.062 26.199 23.475
 [5,] 25.806 21.752 22.222 21.761 25.320 24.430 23.988 25.431 27.210 24.767
 [6,] 24.364 20.761 22.123 22.874 23.583 24.667 24.737 24.635 26.122 26.219
 [7,] 24.477 23.479 23.950 24.104 24.671 26.451 26.276 27.009 26.706 28.361
 [8,] 23.901 23.824 23.504 23.748 24.454 25.618 25.816 26.606 26.878 28.599
 [9,] 23.175 23.105 22.238 23.262 24.122 25.014 25.210 26.268 26.152 27.914
[10,] 23.227 23.110 23.142 22.907 24.252 25.110 25.199 26.462 26.379 27.784
[11,] 21.672 21.759 21.059 21.519 22.084 22.964 23.162 25.246 24.712 25.693
[12,] 21.870 22.073 21.573 22.025 22.991 23.981 24.707 25.180 25.688 26.881
       [,11]  [,12]  [,13]  [,14] [,15]
 [1,] 26.217 26.589 27.132 26.076    NA
 [2,] 24.218 24.848 24.924 25.286    NA
 [3,] 27.914 27.543 28.963 27.660    NA
 [4,] 26.975 26.896 26.589 25.951    NA
 [5,] 28.527 28.878 27.931 26.398    NA
 [6,] 27.139 27.390 28.009 25.565    NA
 [7,] 28.982 28.065 29.229 28.865    NA
 [8,] 28.169 28.141 28.759 30.000    NA
 [9,] 28.056 29.048 28.405 29.261    NA
[10,] 29.136 28.484 27.945 29.012    NA
[11,] 26.291 26.634 25.912 26.992    NA
[12,] 26.987 27.735 26.619 27.897    NA
> darr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]  [,10]
 [1,] -3.065 -0.350 -1.902 -1.548 -1.710 -0.238 -1.528 -1.720 -1.353 -0.751
 [2,]  3.333  2.620  3.555  2.424  3.783  2.027  2.505  2.921  3.678  2.482
 [3,] -2.191 -2.040 -1.918 -1.809 -1.004 -1.039 -2.129 -1.503 -0.783 -3.246
 [4,]  1.066  0.083  0.550  1.146  1.647  0.393  1.342  1.369  1.011  1.292
 [5,] -1.442 -0.991 -0.099  1.113 -1.737  0.237  0.749 -0.796 -1.088  1.452
 [6,]  0.113  2.718  1.827  1.230  1.088  1.784  1.539  2.374  0.584  2.142
 [7,] -0.576  0.345 -0.446 -0.356 -0.217 -0.833 -0.460 -0.403  0.172  0.238
 [8,] -0.726 -0.719 -1.266 -0.486 -0.332 -0.604 -0.606 -0.338 -0.726 -0.685
 [9,]  0.052  0.005  0.904 -0.355  0.130  0.096 -0.011  0.194  0.227 -0.130
[10,] -1.555 -1.351 -2.083 -1.388 -2.168 -2.146 -2.037 -1.216 -1.667 -2.091
[11,]  0.198  0.314  0.514  0.506  0.907  1.017  1.545 -0.066  0.976  1.188
[12,] -0.431 -0.136 -0.025  0.579  0.296 -0.183 -0.343 -0.523 -0.698 -0.664
       [,11]  [,12]  [,13]  [,14] [,15]
 [1,] -1.999 -1.741 -2.208 -0.790    NA
 [2,]  3.696  2.695  4.039  2.374    NA
 [3,] -0.939 -0.647 -2.374 -1.709    NA
 [4,]  1.552  1.982  1.342  0.447    NA
 [5,] -1.388 -1.488  0.078 -0.833    NA
 [6,]  1.843  0.675  1.220  3.300    NA
 [7,] -0.813  0.076 -0.470  1.135    NA
 [8,] -0.113  0.907 -0.354 -0.739    NA
 [9,]  1.080 -0.564 -0.460 -0.249    NA
[10,] -2.845 -1.850 -2.033 -2.020    NA
[11,]  0.696  1.101  0.707  0.905    NA
[12,] -0.398 -0.603 -0.543     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/10pnv1272300931.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/20pnv1272300931.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/30pnv1272300931.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/4bgng1272300931.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]
[1,] 21.4390 20.0000 22.424 20.6150 21.7520 20.761 23.4790 23.504 22.238 22.907
[2,] 22.6040 21.0890 24.677 22.6460 23.9880 23.583 24.4770 23.901 23.262 23.227
[3,] 24.5105 23.1765 26.143 24.0495 25.3755 24.702 26.5785 26.211 25.681 25.789
[4,] 26.2170 24.2390 27.543 26.1990 27.2100 26.219 28.3610 28.169 28.056 27.945
[5,] 27.1320 25.2860 28.963 26.9750 28.8780 28.009 29.2290 30.000 29.261 29.136
      [,11]   [,12]
[1,] 21.059 21.5730
[2,] 21.759 22.0730
[3,] 23.937 24.9435
[4,] 25.912 26.8810
[5,] 26.992 27.8970

$n
 [1] 14 14 14 14 14 14 14 14 14 14 14 14

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 22.98483 21.84634 24.93277 22.54916 24.01494 23.58889 24.93839 24.40874
[2,] 26.03617 24.50666 27.35323 25.54984 26.73606 25.81511 28.21861 28.01326
         [,9]    [,10]   [,11]    [,12]
[1,] 23.65662 23.79672 22.1833 22.91321
[2,] 27.70538 27.78128 25.6907 26.97379

$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(21.439, 22.604, 24.5105, 26.217, 27.132, 20,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5bgng1272300931.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,] -3.065 2.027 -3.246 0.083 -1.7370 0.1130 -0.8330 -1.266 -0.5640 -2.8450
[2,] -1.902 2.482 -2.129 0.550 -1.3880 1.0880 -0.4700 -0.726 -0.2490 -2.0910
[3,] -1.629 2.808 -1.759 1.219 -0.8145 1.6615 -0.3795 -0.605  0.0285 -2.0265
[4,] -0.790 3.678 -1.004 1.369  0.2370 2.1420  0.1720 -0.338  0.1940 -1.5550
[5,] -0.238 4.039 -0.647 1.982  1.4520 3.3000  0.3450 -0.113  0.2270 -1.2160
      [,11]  [,12]
[1,] -0.066 -0.698
[2,]  0.506 -0.543
[3,]  0.806 -0.398
[4,]  1.017 -0.136
[5,]  1.545  0.296

$n
 [1] 14 14 14 14 14 14 14 14 14 14 14 13

$conf
          [,1]     [,2]      [,3]      [,4]       [,5]     [,6]       [,7]
[1,] -2.098567 2.302962 -2.234057 0.8731586 -1.5006932 1.216425 -0.6505991
[2,] -1.159433 3.313038 -1.283943 1.5648414 -0.1283068 2.106575 -0.1084009
           [,8]       [,9]     [,10]     [,11]      [,12]
[1,] -0.7688418 -0.1585668 -2.252838 0.5902186 -0.5763528
[2,] -0.4411582  0.2155668 -1.800162 1.0217814 -0.2196472

$out
[1] 1.135 0.907 0.904 1.080 0.579

$group
[1]  7  8  9  9 12

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

Warning message:
In bxp(list(stats = c(-3.065, -1.902, -1.629, -0.789999999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/64p4j1272300931.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,] 21.6720 20.7610 20.0350 20.0000 20.8940 22.9640 22.2700 22.6440 23.3040
[2,] 23.2010 21.5540 21.6225 21.5335 22.7975 23.6340 23.4800 24.4995 25.2000
[3,] 24.1325 21.9160 22.1725 22.2245 23.8975 24.5485 24.7220 25.3385 26.1755
[4,] 25.2730 23.2945 23.3230 23.0845 24.5625 25.0930 25.2045 26.3650 26.7920
[5,] 26.9310 23.8240 23.9500 24.1040 25.3200 26.4510 26.2760 27.0090 27.2100
       [,10]   [,11]   [,12]  [,13]   [,14] [,15]
[1,] 23.4750 24.2180 24.8480 24.924 25.2860    NA
[2,] 24.8785 26.6330 26.7650 26.604 26.0135    NA
[3,] 26.4700 27.5265 27.6390 27.938 27.3260    NA
[4,] 27.8490 28.3480 28.3125 28.582 28.9385    NA
[5,] 28.5990 29.1360 29.0480 29.229 30.0000    NA

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 23.18745 21.12215 21.39689 21.51708 23.09247 23.88304 23.93544 24.48763
[2,] 25.07755 22.70985 22.94811 22.93192 24.70253 25.21396 25.50856 26.18937
         [,9]    [,10]    [,11]    [,12]    [,13]    [,14] [,15]
[1,] 25.44938 25.11514 26.74428 26.93317 27.03582 25.99189    NA
[2,] 26.90162 27.82486 28.30872 28.34483 28.84018 28.66011    NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(21.672, 23.201, 24.1325, 25.273, 26.931, 20.761,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/74p4j1272300931.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,] 22.88557 23.1765 22.80900
[2,] 24.30357 24.2800 24.48669
[3,] 25.09318 25.1595 25.26962
[4,] 25.88139 25.9660 25.90044
[5,] 26.47321 26.5785 26.40625

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 24.37352 24.39050 24.62480
[2,] 25.81283 25.92850 25.91445

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(22.8855714285714, 24.3035714285714, 25.0931785714286,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/10pnv1272300931.ps tmp/10pnv1272300931.png",intern=TRUE))
character(0)
> try(system("convert tmp/20pnv1272300931.ps tmp/20pnv1272300931.png",intern=TRUE))
character(0)
> try(system("convert tmp/30pnv1272300931.ps tmp/30pnv1272300931.png",intern=TRUE))
character(0)
> try(system("convert tmp/4bgng1272300931.ps tmp/4bgng1272300931.png",intern=TRUE))
character(0)
> try(system("convert tmp/5bgng1272300931.ps tmp/5bgng1272300931.png",intern=TRUE))
character(0)
> try(system("convert tmp/64p4j1272300931.ps tmp/64p4j1272300931.png",intern=TRUE))
character(0)
> try(system("convert tmp/74p4j1272300931.ps tmp/74p4j1272300931.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  2.013   1.543   2.498