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(90.8,96.4,90,92.1,97.2,95.1,88.5,91,90.5,75,66.3,66,68.4,70.6,83.9,90.1,90.6,87.1,90.8,94.1,99.8,96.8,87,96.3,107.1,115.2,106.1,89.5,91.3,97.6,100.7,104.6,94.7,101.8,102.5,105.3,110.3,109.8,117.3,118.8,131.3,125.9,133.1,147,145.8,164.4,149.8,137.7,151.7,156.8,180,180.4,170.4,191.6,199.5,218.2,217.5,205,194,199.3)
> 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] 60
> (np <- floor(n / par1))
[1] 5
> 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] 5 5 5 5 5 5 5 5 5 5 5 5
> arr
      [,1] [,2]  [,3]  [,4]  [,5] [,6]
 [1,] 90.8 68.4 107.1 110.3 151.7   NA
 [2,] 96.4 70.6 115.2 109.8 156.8   NA
 [3,] 90.0 83.9 106.1 117.3 180.0   NA
 [4,] 92.1 90.1  89.5 118.8 180.4   NA
 [5,] 97.2 90.6  91.3 131.3 170.4   NA
 [6,] 95.1 87.1  97.6 125.9 191.6   NA
 [7,] 88.5 90.8 100.7 133.1 199.5   NA
 [8,] 91.0 94.1 104.6 147.0 218.2   NA
 [9,] 90.5 99.8  94.7 145.8 217.5   NA
[10,] 75.0 96.8 101.8 164.4 205.0   NA
[11,] 66.3 87.0 102.5 149.8 194.0   NA
[12,] 66.0 96.3 105.3 137.7 199.3   NA
> darr
       [,1] [,2]  [,3]  [,4]  [,5] [,6]
 [1,]   5.6  2.2   8.1  -0.5   5.1   NA
 [2,]  -6.4 13.3  -9.1   7.5  23.2   NA
 [3,]   2.1  6.2 -16.6   1.5   0.4   NA
 [4,]   5.1  0.5   1.8  12.5 -10.0   NA
 [5,]  -2.1 -3.5   6.3  -5.4  21.2   NA
 [6,]  -6.6  3.7   3.1   7.2   7.9   NA
 [7,]   2.5  3.3   3.9  13.9  18.7   NA
 [8,]  -0.5  5.7  -9.9  -1.2  -0.7   NA
 [9,] -15.5 -3.0   7.1  18.6 -12.5   NA
[10,]  -8.7 -9.8   0.7 -14.6 -11.0   NA
[11,]  -0.3  9.3   2.8 -12.1   5.3   NA
[12,]   2.4 10.8   5.0  14.0    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/rcomp/tmp/1zyem1229595782.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/221291229595782.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/3ptb01229595782.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/40fsx1229595782.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,]  68.4  70.6  83.9  89.5  90.6  87.1  88.5  91.0  90.5  75.0  66.3  66.0
[2,]  90.8  96.4  90.0  90.1  91.3  95.1  90.8  94.1  94.7  96.8  87.0  96.3
[3,] 107.1 109.8 106.1  92.1  97.2  97.6 100.7 104.6  99.8 101.8 102.5 105.3
[4,] 110.3 115.2 117.3 118.8 131.3 125.9 133.1 147.0 145.8 164.4 149.8 137.7
[5,] 110.3 115.2 117.3 118.8 170.4 125.9 133.1 218.2 217.5 205.0 194.0 199.3

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

$conf
          [,1]      [,2]      [,3]      [,4]     [,5]     [,6]      [,7]
[1,]  93.32135  96.51597  86.80989  71.82065  68.9361  75.8368  70.81093
[2,] 120.87865 123.08403 125.39011 112.37935 125.4639 119.3632 130.58907
          [,8]      [,9]     [,10]     [,11]     [,12]
[1,]  67.22099  63.69287  54.03401  58.12568  76.04686
[2,] 141.97901 135.90713 149.56599 146.87432 134.55314

$out
[1] 151.7 156.8 180.0 180.4 191.6 199.5

$group
[1] 1 2 3 4 6 7

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

Warning message:
In bxp(list(stats = c(68.4, 90.8, 107.1, 110.3, 110.3, 70.6, 96.4,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5rbt01229595782.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,] -0.5 -9.1  0.4  0.5 -5.4  3.1  2.5 -1.2 -15.5 -11.0  -0.3   2.4
[2,]  2.2 -6.4  0.4  0.5 -3.5  3.1  3.3 -1.2 -12.5 -11.0  -0.3   3.7
[3,]  5.1  7.5  1.5  1.8 -2.1  3.7  3.9 -0.7  -3.0  -9.8   2.8   7.9
[4,]  5.6 13.3  2.1  5.1  6.3  7.2 13.9 -0.5   7.1  -8.7   5.3  12.4
[5,]  8.1 23.2  2.1  5.1  6.3  7.9 18.7 -0.5  18.6  -8.7   9.3  14.0

$n
 [1] 5 5 5 5 5 5 5 5 5 5 5 4

$conf
         [,1]     [,2]      [,3]      [,4]      [,5]      [,6]      [,7]
[1,] 2.697569 -6.41997 0.2987843 -1.450348 -9.024655 0.8029503 -3.589933
[2,] 7.502431 21.41997 2.7012157  5.050348  4.824655 6.5970497 11.389933
           [,8]      [,9]      [,10]     [,11]  [,12]
[1,] -1.1946182 -16.84931 -11.425174 -1.156946  1.027
[2,] -0.2053818  10.84931  -8.174826  6.756946 14.773

$out
 [1]   6.2 -16.6  12.5 -10.0  21.2  -6.6   5.7  -9.9   0.7 -14.6 -12.1

$group
 [1]  3  3  4  4  5  6  8  8 10 10 11

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

Warning message:
In bxp(list(stats = c(-0.5, 2.19999999999999, 5.10000000000002,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6b73v1229595782.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]
[1,] 66.00 83.90  89.50 109.80 151.70   NA
[2,] 81.75 85.45  96.15 118.05 175.20   NA
[3,] 90.65 90.35 102.15 132.20 192.80   NA
[4,] 93.60 95.20 105.70 146.40 202.25   NA
[5,] 97.20 99.80 115.20 164.40 218.20   NA

$n
[1] 12 12 12 12 12  0

$conf
         [,1]     [,2]      [,3]     [,4]     [,5] [,6]
[1,] 85.24514 85.90296  97.79418 119.2694 180.4623   NA
[2,] 96.05486 94.79704 106.50582 145.1306 205.1377   NA

$out
[1] 68.4 70.6

$group
[1] 2 2

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

Warning message:
In bxp(list(stats = c(66, 81.75, 90.65, 93.6, 97.2, 83.9, 85.45,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/7u8tn1229595782.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,] 105.66  92.10 100.550
[2,] 114.82  98.70 105.125
[3,] 119.69 102.15 111.625
[4,] 125.56 105.70 119.325
[5,] 130.98 109.80 130.600

$n
[1] 12 12 12

$conf
         [,1]      [,2]     [,3]
[1,] 114.7914  98.95725 105.1483
[2,] 124.5886 105.34275 118.1017

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(105.66, 114.82, 119.69, 125.56, 130.98, 92.1,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1zyem1229595782.ps tmp/1zyem1229595782.png")
> system("convert tmp/221291229595782.ps tmp/221291229595782.png")
> system("convert tmp/3ptb01229595782.ps tmp/3ptb01229595782.png")
> system("convert tmp/40fsx1229595782.ps tmp/40fsx1229595782.png")
> system("convert tmp/5rbt01229595782.ps tmp/5rbt01229595782.png")
> system("convert tmp/6b73v1229595782.ps tmp/6b73v1229595782.png")
> system("convert tmp/7u8tn1229595782.ps tmp/7u8tn1229595782.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.249   0.975   1.985