R version 2.13.0 (2011-04-13)
Copyright (C) 2011 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-pc-linux-gnu (32-bit)

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(105.32,106.93,107.99,109.41,111.79,112.7,113.41,115.14,115.6,116.21,116.9,117.31,117.62,118.66,120,120.41,120.8,120.93,121.54,122.04,122.5,123.01,123.79,123.99,124.58,125.77,127.6,130.97,131.79,132.29,132.29,132.65,133.88,134.79,135.9,136.33,136.33,137.91,138.92,140.12,141.92,142.57,143.69,144.61,145.02,146.07,146.77,147.55,148.14,148.62,151.16,152.56,154.55,156.17,158.8,159.39,162.44,166.48,168.3,170.32,170.64,172.85,174.38,177.2,178.96,179.62,180.27,183.38,190.81,193.72,196.86,197.73)
> 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] 72
> (np <- floor(n / par1))
[1] 6
> 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] 6 6 6 6 6 6 6 6 6 6 6 6
> arr
        [,1]   [,2]   [,3]   [,4]   [,5]   [,6] [,7]
 [1,] 105.32 117.62 124.58 136.33 148.14 170.64   NA
 [2,] 106.93 118.66 125.77 137.91 148.62 172.85   NA
 [3,] 107.99 120.00 127.60 138.92 151.16 174.38   NA
 [4,] 109.41 120.41 130.97 140.12 152.56 177.20   NA
 [5,] 111.79 120.80 131.79 141.92 154.55 178.96   NA
 [6,] 112.70 120.93 132.29 142.57 156.17 179.62   NA
 [7,] 113.41 121.54 132.29 143.69 158.80 180.27   NA
 [8,] 115.14 122.04 132.65 144.61 159.39 183.38   NA
 [9,] 115.60 122.50 133.88 145.02 162.44 190.81   NA
[10,] 116.21 123.01 134.79 146.07 166.48 193.72   NA
[11,] 116.90 123.79 135.90 146.77 168.30 196.86   NA
[12,] 117.31 123.99 136.33 147.55 170.32 197.73   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,] 1.61 1.04 1.19 1.58 0.48 2.21   NA
 [2,] 1.06 1.34 1.83 1.01 2.54 1.53   NA
 [3,] 1.42 0.41 3.37 1.20 1.40 2.82   NA
 [4,] 2.38 0.39 0.82 1.80 1.99 1.76   NA
 [5,] 0.91 0.13 0.50 0.65 1.62 0.66   NA
 [6,] 0.71 0.61 0.00 1.12 2.63 0.65   NA
 [7,] 1.73 0.50 0.36 0.92 0.59 3.11   NA
 [8,] 0.46 0.46 1.23 0.41 3.05 7.43   NA
 [9,] 0.61 0.51 0.91 1.05 4.04 2.91   NA
[10,] 0.69 0.78 1.11 0.70 1.82 3.14   NA
[11,] 0.41 0.20 0.43 0.78 2.02 0.87   NA
[12,] 0.31 0.59 0.00 0.59 0.32   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/wessaorg/rcomp/tmp/1j7wa1331730574.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/2n5i31331730574.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3cary1331730574.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4a3eb1331730574.ps",horizontal=F,onefile=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,] 105.320 106.93 107.99 109.410 111.790 112.70 113.41 115.14 115.60 116.21
[2,] 117.620 118.66 120.00 120.410 120.800 120.93 121.54 122.04 122.50 123.01
[3,] 130.455 131.84 133.26 135.545 136.855 137.43 137.99 138.63 139.45 140.43
[4,] 148.140 148.62 151.16 152.560 154.550 156.17 158.80 159.39 162.44 166.48
[5,] 170.640 172.85 174.38 177.200 178.960 179.62 180.27 183.38 190.81 193.72
       [,11]  [,12]
[1,] 116.900 117.31
[2,] 123.790 123.99
[3,] 141.335 141.94
[4,] 168.300 170.32
[5,] 196.860 197.73

$n
 [1] 6 6 6 6 6 6 6 6 6 6 6 6

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]    [,8]
[1,] 110.7686 112.5148 113.1608 114.8072 115.0852 114.6991 113.9561 114.538
[2,] 150.1414 151.1652 153.3592 156.2828 158.6248 160.1609 162.0239 162.722
         [,9]    [,10]    [,11]    [,12]
[1,] 113.6874 112.3904 112.6246 112.0557
[2,] 165.2126 168.4696 170.0454 171.8243

$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(105.32, 117.62, 130.455, 148.14, 170.64, 106.93,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/5lqdu1331730574.ps",horizontal=F,onefile=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] [,12]
[1,] 0.480 1.010 0.41 0.39 0.130 0.00 0.360 0.410 0.51 0.690 0.200  0.00
[2,] 1.040 1.060 1.20 0.82 0.500 0.61 0.500 0.460 0.61 0.700 0.410  0.31
[3,] 1.385 1.435 1.41 1.78 0.655 0.68 0.755 0.845 0.98 0.945 0.605  0.32
[4,] 1.610 1.830 2.82 1.99 0.910 1.12 1.730 3.050 2.91 1.820 0.870  0.59
[5,] 2.210 2.540 3.37 2.38 0.910 1.12 3.110 3.050 4.04 3.140 0.870  0.59

$n
 [1] 6 6 6 6 6 6 6 6 6 6 6 5

$conf
         [,1]      [,2]      [,3]     [,4]      [,5]      [,6]        [,7]
[1,] 1.017332 0.9383251 0.3650477 1.025312 0.3905368 0.3510335 -0.03838973
[2,] 1.752668 1.9316749 2.4549523 2.534688 0.9194632 1.0089665  1.54838973
           [,8]       [,9]     [,10]     [,11]     [,12]
[1,] -0.8256337 -0.5035743 0.2225638 0.3082851 0.1221527
[2,]  2.5156337  2.4635743 1.6674362 0.9017149 0.5178473

$out
[1] 1.62 2.63 7.43 2.02

$group
[1]  5  6  8 11

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

Warning message:
In bxp(list(stats = c(0.480000000000018, 1.03999999999999, 1.38499999999999,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/6x1by1331730574.ps",horizontal=F,onefile=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]
[1,] 105.320 117.620 124.580 136.330 148.140 170.640   NA
[2,] 108.700 120.205 129.285 139.520 151.860 175.790   NA
[3,] 113.055 121.235 132.290 143.130 157.485 179.945   NA
[4,] 115.905 122.755 134.335 145.545 164.460 192.265   NA
[5,] 117.310 123.990 136.330 147.550 170.320 197.730   NA

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

$conf
         [,1]     [,2]     [,3]    [,4]     [,5]     [,6] [,7]
[1,] 109.7688 120.0719 129.9867 140.382 151.7381 172.4306   NA
[2,] 116.3412 122.3981 134.5933 145.878 163.2319 187.4594   NA

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(105.32, 108.7, 113.055, 115.905, 117.31, 117.62,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/wessaorg/rcomp/tmp/71nkm1331730574.ps",horizontal=F,onefile=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.7717 130.4550 132.2738
[2,] 137.5600 134.4025 135.6250
[3,] 141.1900 137.7100 138.9475
[4,] 145.8775 139.9400 142.6906
[5,] 148.8717 141.9400 145.8512

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 137.3963 135.1843 135.7248
[2,] 144.9837 140.2357 142.1702

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(133.771666666667, 137.56, 141.19, 145.8775,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> try(system("convert tmp/1j7wa1331730574.ps tmp/1j7wa1331730574.png",intern=TRUE))
character(0)
> try(system("convert tmp/2n5i31331730574.ps tmp/2n5i31331730574.png",intern=TRUE))
character(0)
> try(system("convert tmp/3cary1331730574.ps tmp/3cary1331730574.png",intern=TRUE))
character(0)
> try(system("convert tmp/4a3eb1331730574.ps tmp/4a3eb1331730574.png",intern=TRUE))
character(0)
> try(system("convert tmp/5lqdu1331730574.ps tmp/5lqdu1331730574.png",intern=TRUE))
character(0)
> try(system("convert tmp/6x1by1331730574.ps tmp/6x1by1331730574.png",intern=TRUE))
character(0)
> try(system("convert tmp/71nkm1331730574.ps tmp/71nkm1331730574.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.724   0.400   2.388