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(93.26,100.61,109.57,107.08,110.33,110.36,106.5,104.3,107.21,109.34,108.2,109.86,108.68,113.38,117.12,116.23,114.75,115.81,115.86,117.8,117.11,116.31,118.38,121.57,121.65,124.2,126.12,128.6,128.16,130.12,135.83,138.05,134.99,132.38,128.94,128.12,127.84,132.43,134.13,134.78,133.13,129.08,134.48,132.86,134.08,134.54,134.51,135.97,136.09,139.14,135.63,136.55,138.83,138.84,135.37,132.22,134.75,135.98,136.06,138.05,139.59,140.58,139.81,140.77,140.96,143.59,142.7,145.11,146.7,148.53,148.99,149.65,151.11,154.82,156.56,157.6,155.24,160.68,163.22,164.55,166.76,159.05,159.82,164.95,162.89,163.55,158.68,157.97,156.59,161.56,162.31,166.26,168.45,163.63,153.2,133.52,123.28)
> 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] 97
> (np <- floor(n / par1))
[1] 8
> 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] 9 8 8 8 8 8 8 8 8 8 8 8
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6]   [,7]   [,8]   [,9]
 [1,]  93.26 108.68 121.65 127.84 136.09 139.59 151.11 162.89 123.28
 [2,] 100.61 113.38 124.20 132.43 139.14 140.58 154.82 163.55     NA
 [3,] 109.57 117.12 126.12 134.13 135.63 139.81 156.56 158.68     NA
 [4,] 107.08 116.23 128.60 134.78 136.55 140.77 157.60 157.97     NA
 [5,] 110.33 114.75 128.16 133.13 138.83 140.96 155.24 156.59     NA
 [6,] 110.36 115.81 130.12 129.08 138.84 143.59 160.68 161.56     NA
 [7,] 106.50 115.86 135.83 134.48 135.37 142.70 163.22 162.31     NA
 [8,] 104.30 117.80 138.05 132.86 132.22 145.11 164.55 166.26     NA
 [9,] 107.21 117.11 134.99 134.08 134.75 146.70 166.76 168.45     NA
[10,] 109.34 116.31 132.38 134.54 135.98 148.53 159.05 163.63     NA
[11,] 108.20 118.38 128.94 134.51 136.06 148.99 159.82 153.20     NA
[12,] 109.86 121.57 128.12 135.97 138.05 149.65 164.95 133.52     NA
> darr
       [,1]  [,2]  [,3]  [,4]  [,5]  [,6]  [,7]   [,8] [,9]
 [1,]  7.35  4.70  2.55  4.59  3.05  0.99  3.71   0.66   NA
 [2,]  8.96  3.74  1.92  1.70 -3.51 -0.77  1.74  -4.87   NA
 [3,] -2.49 -0.89  2.48  0.65  0.92  0.96  1.04  -0.71   NA
 [4,]  3.25 -1.48 -0.44 -1.65  2.28  0.19 -2.36  -1.38   NA
 [5,]  0.03  1.06  1.96 -4.05  0.01  2.63  5.44   4.97   NA
 [6,] -3.86  0.05  5.71  5.40 -3.47 -0.89  2.54   0.75   NA
 [7,] -2.20  1.94  2.22 -1.62 -3.15  2.41  1.33   3.95   NA
 [8,]  2.91 -0.69 -3.06  1.22  2.53  1.59  2.21   2.19   NA
 [9,]  2.13 -0.80 -2.61  0.46  1.23  1.83 -7.71  -4.82   NA
[10,] -1.14  2.07 -3.44 -0.03  0.08  0.46  0.77 -10.43   NA
[11,]  1.66  3.19 -0.82  1.46  1.99  0.66  5.13 -19.68   NA
[12,] -1.18  0.08 -0.28  0.12  1.54  1.46 -2.06 -10.24   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/1041l1229541820.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/210yy1229541820.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/3588j1229541820.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/4eao31229541820.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,] 108.68 100.610 109.570 107.080 110.330 110.360 106.500 104.300 107.210
[2,] 121.65 118.790 121.620 122.415 121.455 122.445 125.170 125.010 125.595
[3,] 127.84 135.785 134.880 135.665 135.980 134.480 135.600 135.455 134.870
[4,] 139.59 147.700 148.185 149.185 148.100 152.135 152.505 154.830 156.730
[5,] 162.89 163.550 158.680 157.970 156.590 161.560 163.220 166.260 168.450
       [,10]   [,11]   [,12]
[1,] 109.340 108.200 109.860
[2,] 124.345 123.660 124.845
[3,] 135.260 135.285 134.745
[4,] 153.790 151.095 143.850
[5,] 163.630 159.820 164.950

$n
 [1] 9 8 8 8 8 8 8 8 8 8 8 8

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,] 118.3916 119.6355 120.0404 120.7109 121.0957 117.8947 120.3303 118.7971
[2,] 137.2884 151.9345 149.7196 150.6191 150.8643 151.0653 150.8697 152.1129
         [,9]    [,10]    [,11]    [,12]
[1,] 117.4775 118.8116 119.9594 124.1285
[2,] 152.2625 151.7084 150.6106 145.3615

$out
[1] 93.26

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(108.68, 121.65, 127.84, 139.59, 162.89, 100.61,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/5q2o81229541820.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,] 0.660 -4.87 -2.490 -2.360 -4.05 -3.86 -3.150 -0.690 -7.710 -3.440 -0.82
[2,] 1.770 -2.14 -0.800 -1.565  0.02 -2.18 -1.910  0.265 -3.715 -2.290 -0.08
[3,] 3.380  1.72  0.785 -0.910  1.51  0.40  1.635  1.890 -0.170  0.025  1.56
[4,] 4.645  2.83  1.000  1.235  3.80  3.97  2.315  2.370  1.530  0.615  2.59
[5,] 7.350  8.96  2.480  3.250  5.44  5.71  3.950  2.910  2.130  2.070  5.13
     [,12]
[1,] -2.06
[2,] -1.62
[3,] -0.10
[4,]  0.79
[5,]  1.54

$n
 [1] 8 8 8 8 8 8 8 8 8 8 8 8

$conf
         [,1]      [,2]       [,3]       [,4]       [,5]      [,6]       [,7]
[1,] 1.773984 -1.056313 -0.2205058 -2.4741202 -0.6015623 -3.035478 -0.7251457
[2,] 4.986016  4.496313  1.7905058  0.6541202  3.6215623  3.835478  3.9951457
          [,8]      [,9]     [,10]      [,11]     [,12]
[1,] 0.7141168 -3.099932 -1.597775 0.06849967 -1.446261
[2,] 3.0658832  2.759932  1.647775 3.05150033  1.246261

$out
[1]  -3.06 -10.43 -19.68 -10.24

$group
[1]  8 10 11 12

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

Warning message:
In bxp(list(stats = c(0.660000000000025, 1.77000000000000, 3.37999999999998,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/6b8pc1229541820.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,] 100.610 113.38 121.650 132.430 132.220 139.590 151.110 153.200 123.28
[2,] 105.400 115.28 127.120 132.645 135.500 140.675 155.900 157.280 123.28
[3,] 107.705 116.27 128.770 134.105 136.075 143.145 159.435 161.935 123.28
[4,] 109.715 117.46 133.685 134.525 138.440 147.615 163.885 163.590 123.28
[5,] 110.360 118.38 138.050 135.970 139.140 149.650 166.760 168.450 123.28

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]    [,7]     [,8]
[1,] 105.7369 115.2757 125.7757 133.2475 134.7340 139.9796 155.793 159.0570
[2,] 109.6731 117.2643 131.7643 134.9625 137.4160 146.3104 163.077 164.8130
       [,9]
[1,] 123.28
[2,] 123.28

$out
[1]  93.26 108.68 121.57 127.84 129.08 133.52

$group
[1] 1 2 2 4 4 8

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

Warning message:
In bxp(list(stats = c(100.61, 105.4, 107.705, 109.715, 110.36, 113.38,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/freestat/rcomp/tmp/7w0wz1229541820.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,] 133.5888 134.4800 130.6200
[2,] 134.7256 134.8075 133.8250
[3,] 135.6119 135.2725 136.0275
[4,] 137.2519 135.6325 139.0031
[5,] 138.7562 135.9800 140.7763

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 134.4596 134.8962 133.6657
[2,] 136.7641 135.6488 138.3893

$out
[1] 129.3767 127.8400

$group
[1] 1 2

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

Warning message:
In bxp(list(stats = c(133.58875, 134.725625, 135.611875, 137.251875,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1041l1229541820.ps tmp/1041l1229541820.png")
> system("convert tmp/210yy1229541820.ps tmp/210yy1229541820.png")
> system("convert tmp/3588j1229541820.ps tmp/3588j1229541820.png")
> system("convert tmp/4eao31229541820.ps tmp/4eao31229541820.png")
> system("convert tmp/5q2o81229541820.ps tmp/5q2o81229541820.png")
> system("convert tmp/6b8pc1229541820.ps tmp/6b8pc1229541820.png")
> system("convert tmp/7w0wz1229541820.ps tmp/7w0wz1229541820.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.850   1.454   2.084