R version 2.6.0 (2007-10-03)
Copyright (C) 2007 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.

  Natural language support but running in an English locale

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(109.20,88.60,94.30,98.30,86.40,80.60,104.10,108.20,93.40,71.90,94.10,94.90,96.40,91.10,84.40,86.40,88.00,75.10,109.70,103.00,82.10,68.00,96.40,94.30,90.00,88.00,76.10,82.50,81.40,66.50,97.20,94.10,80.70,70.50,87.80,89.50,99.60,84.20,75.10,92.00,80.80,73.10,99.80,90.00,83.10,72.40,78.80,87.30,91.00,80.10,73.60,86.40,74.50,71.20,92.40,81.50,85.30,69.90,84.20,90.70,100.30)
> par1 = '12'
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: Wessa P., (2007), Mean Plot (v1.0.1) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_meanplot.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description
> par1 <- as.numeric(par1)
> (n <- length(x))
[1] 61
> (np <- floor(n / par1))
[1] 5
> arr <- array(NA,dim=c(par1,np+1))
> ari <- array(0,dim=par1)
> j <- 0
> for (i in 1:n)
+ {
+ j = j + 1
+ ari[j] = ari[j] + 1
+ arr[j,ari[j]] <- x[i]
+ if (j == par1) j = 0
+ }
> ari
 [1] 6 5 5 5 5 5 5 5 5 5 5 5
> arr
       [,1]  [,2] [,3] [,4] [,5]  [,6]
 [1,] 109.2  96.4 90.0 99.6 91.0 100.3
 [2,]  88.6  91.1 88.0 84.2 80.1    NA
 [3,]  94.3  84.4 76.1 75.1 73.6    NA
 [4,]  98.3  86.4 82.5 92.0 86.4    NA
 [5,]  86.4  88.0 81.4 80.8 74.5    NA
 [6,]  80.6  75.1 66.5 73.1 71.2    NA
 [7,] 104.1 109.7 97.2 99.8 92.4    NA
 [8,] 108.2 103.0 94.1 90.0 81.5    NA
 [9,]  93.4  82.1 80.7 83.1 85.3    NA
[10,]  71.9  68.0 70.5 72.4 69.9    NA
[11,]  94.1  96.4 87.8 78.8 84.2    NA
[12,]  94.9  94.3 89.5 87.3 90.7    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/1iooq1193397414.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/2hkml1193397414.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/36sp51193397414.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/4c97k1193397414.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,]  90.0 80.1 73.6 82.5 74.5 66.5  92.4  81.5 80.7  68.0  78.8  87.3
[2,]  91.0 84.2 75.1 86.4 80.8 71.2  97.2  90.0 82.1  69.9  84.2  89.5
[3,]  98.0 88.0 76.1 86.4 81.4 73.1  99.8  94.1 83.1  70.5  87.8  90.7
[4,] 100.3 88.6 84.4 92.0 86.4 75.1 104.1 103.0 85.3  71.9  94.1  94.3
[5,] 109.2 91.1 94.3 98.3 88.0 80.6 109.7 108.2 85.3  72.4  96.4  94.9

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

$conf
         [,1]     [,2]     [,3]     [,4]     [,5]     [,6]      [,7]      [,8]
[1,]  92.0012 84.89097 69.52864 82.44305 77.44305 70.34427  94.92448  84.91423
[2,] 103.9988 91.10903 82.67136 90.35695 85.35695 75.85573 104.67552 103.28577
         [,9]   [,10]    [,11]    [,12]
[1,] 80.83889 69.0868 80.80468 87.30833
[2,] 85.36111 71.9132 94.79532 94.09167

$out
[1] 93.4

$group
[1] 9

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

Warning message:
In bxp(list(stats = c(90, 91, 98, 100.3, 109.2, 80.1, 84.2, 88,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/5xuye1193397414.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,]  71.9  68.00 66.50 72.40 69.90 100.3
[2,]  87.5  83.25 78.40 76.95 74.05 100.3
[3,]  94.2  89.55 85.15 83.65 82.85 100.3
[4,] 101.2  96.40 89.75 91.00 88.55 100.3
[5,] 109.2 109.70 97.20 99.80 92.40 100.3

$n
[1] 12 12 12 12 12  1

$conf
          [,1]    [,2]     [,3]    [,4]     [,5]  [,6]
[1,]  87.95134 83.5522 79.97319 77.2417 76.23645 100.3
[2,] 100.44866 95.5478 90.32681 90.0583 89.46355 100.3

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(71.9, 87.5, 94.2, 101.2, 109.2, 68, 83.25,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6836e1193397414.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,]  70.54 70.50  70.90000
[2,]  81.46 78.75  81.67500
[3,]  87.33 87.10  87.77500
[4,]  93.35 92.40  94.06875
[5,] 100.64 99.80 100.65000

$n
[1] 12 12 12

$conf
         [,1]     [,2]     [,3]
[1,] 81.90689 80.87414 82.12213
[2,] 92.75311 93.32586 93.42787

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(70.54, 81.46, 87.33, 93.35, 100.64, 70.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1iooq1193397414.ps tmp/1iooq1193397414.png")
> system("convert tmp/2hkml1193397414.ps tmp/2hkml1193397414.png")
> system("convert tmp/36sp51193397414.ps tmp/36sp51193397414.png")
> system("convert tmp/4c97k1193397414.ps tmp/4c97k1193397414.png")
> system("convert tmp/5xuye1193397414.ps tmp/5xuye1193397414.png")
> system("convert tmp/6836e1193397414.ps tmp/6836e1193397414.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.394   1.356   2.613