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(155,176,180,181,208,183,185,205,180,196,183,147,128,158,186,165,191,168,171,169,157,175,156,129,89,138,146,151,156,129,146,141,137,155,147,128,92,136,159,131,134,148,146,144,161,140,141,139,94,136,164,141,159,162,154,166,156,147,161,135,98,150,173,144,167,161,156,175,163,159,167,148,119,150,161,136,166,155,140,141)
> 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] 80
> (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] 7 7 7 7 7 7 7 7 6 6 6 6
> arr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,]  155  128   89   92   94   98  119
 [2,]  176  158  138  136  136  150  150
 [3,]  180  186  146  159  164  173  161
 [4,]  181  165  151  131  141  144  136
 [5,]  208  191  156  134  159  167  166
 [6,]  183  168  129  148  162  161  155
 [7,]  185  171  146  146  154  156  140
 [8,]  205  169  141  144  166  175  141
 [9,]  180  157  137  161  156  163   NA
[10,]  196  175  155  140  147  159   NA
[11,]  183  156  147  141  161  167   NA
[12,]  147  129  128  139  135  148   NA
> darr
      [,1] [,2] [,3] [,4] [,5] [,6] [,7]
 [1,]   21   30   49   44   42   52   31
 [2,]    4   28    8   23   28   23   11
 [3,]    1  -21    5  -28  -23  -29  -25
 [4,]   27   26    5    3   18   23   30
 [5,]  -25  -23  -27   14    3   -6  -11
 [6,]    2    3   17   -2   -8   -5  -15
 [7,]   20   -2   -5   -2   12   19    1
 [8,]  -25  -12   -4   17  -10  -12   NA
 [9,]   16   18   18  -21   -9   -4   NA
[10,]  -13  -19   -8    1   14    8   NA
[11,]  -36  -27  -19   -2  -26  -19   NA
[12,]  -19  -40  -36  -45  -37  -29   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/19dmo1244199120.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/2hzay1244199120.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/3k0hx1244199120.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/4qo3u1244199120.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] [,11] [,12]
[1,]  89.0  136 146.0 131.0 134.0 148.0 140.0 141.0  156   140 141.0   128
[2,]  93.0  137 160.0 138.5 157.5 151.5 146.0 142.5  156   147 147.0   129
[3,]  98.0  150 164.0 144.0 166.0 161.0 154.0 166.0  159   157 158.5   137
[4,] 123.5  154 176.5 158.0 179.0 165.0 163.5 172.0  163   175 167.0   147
[5,] 155.0  176 186.0 181.0 208.0 183.0 185.0 205.0  163   196 183.0   148

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,]  79.7859 139.8479 154.1465 132.3549 153.1605 152.9380 143.5493 148.3831
[2,] 116.2141 160.1521 173.8535 155.6451 178.8395 169.0620 164.4507 183.6169
         [,9]    [,10]    [,11]    [,12]
[1,] 154.4848 138.9391 145.5994 125.3894
[2,] 163.5152 175.0609 171.4006 148.6106

$out
[1] 129 180 137

$group
[1] 6 9 9

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

Warning message:
In bxp(list(stats = c(89, 93, 98, 123.5, 155, 136, 137, 150, 154,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5gxvr1244199120.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] [,12]
[1,] 21.0  4.0 -29.0  3.0 -27.0 -15.0 -5.0  -12  -21 -19.0 -36.0 -45.0
[2,] 30.5  9.5 -26.5 11.5 -24.0  -6.5 -2.0  -12   -9 -13.0 -27.0 -40.0
[3,] 42.0 23.0 -23.0 23.0 -11.0  -2.0  1.0  -11    6  -3.5 -22.5 -36.5
[4,] 46.5 25.5 -10.0 26.5  -1.5   2.5 15.5   -4   18   8.0 -19.0 -29.0
[5,] 52.0 28.0   5.0 30.0  14.0   3.0 20.0   -4   18  14.0 -19.0 -19.0

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

$conf
         [,1]     [,2]      [,3]     [,4]       [,5]      [,6]      [,7]
[1,] 32.44506 13.44506 -32.85353 14.04224 -24.436637 -7.374655 -9.450718
[2,] 51.55494 32.55494 -13.14647 31.95776   2.436637  3.374655 11.450718
           [,8]      [,9]     [,10]     [,11]     [,12]
[1,] -16.160258 -11.41587 -17.04568 -27.66026 -43.59536
[2,]  -5.839742  23.41587  10.04568 -17.33974 -29.40464

$out
[1]  17 -25  17  -2

$group
[1]  6  8  8 11

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

Warning message:
In bxp(list(stats = c(21, 30.5, 42, 46.5, 52, 4, 9.5, 23, 25.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6gm4h1244199120.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]
[1,] 176.0 156.0 128.0 131.0 135.0  144 119.0
[2,] 178.0 156.5 133.0 135.0 138.5  149 138.0
[3,] 182.0 166.5 143.5 140.5 155.0  160 145.5
[4,] 190.5 173.0 149.0 147.0 161.5  167 158.0
[5,] 208.0 191.0 156.0 161.0 166.0  175 166.0

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]     [,7]
[1,] 176.2987 158.9742 136.2023 135.0267 144.5095 151.7901 134.3277
[2,] 187.7013 174.0258 150.7977 145.9733 165.4905 168.2099 156.6723

$out
[1] 155 147 128 129  89  92  94  98

$group
[1] 1 1 2 2 3 4 5 6

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

Warning message:
In bxp(list(stats = c(176, 178, 182, 190.5, 208, 156, 156.5, 166.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/76fm21244199120.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,] 137.6667 137.00 137.7500
[2,] 149.5000 147.00 146.8750
[3,] 158.5000 157.75 157.3125
[4,] 162.5000 162.50 159.6875
[5,] 168.7143 166.00 168.2500

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 152.5706 150.6803 151.4686
[2,] 164.4294 164.8197 163.1564

$out
[1] 110.7143  98.0000 108.2500

$group
[1] 1 2 3

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

Warning message:
In bxp(list(stats = c(137.666666666667, 149.5, 158.5, 162.5, 168.714285714286,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/19dmo1244199120.ps tmp/19dmo1244199120.png")
> system("convert tmp/2hzay1244199120.ps tmp/2hzay1244199120.png")
> system("convert tmp/3k0hx1244199120.ps tmp/3k0hx1244199120.png")
> system("convert tmp/4qo3u1244199120.ps tmp/4qo3u1244199120.png")
> system("convert tmp/5gxvr1244199120.ps tmp/5gxvr1244199120.png")
> system("convert tmp/6gm4h1244199120.ps tmp/6gm4h1244199120.png")
> system("convert tmp/76fm21244199120.ps tmp/76fm21244199120.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.183   0.934   1.486