R version 2.7.2 (2008-08-25)
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(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 = '15'
> x <- x[x>quantile(x,0.05) & x<quantile(x,0.95)] 
> #'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] 53
> (np <- floor(n / par1))
[1] 3
> 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] 4 4 4 4 4 4 4 4 3 3 3 3 3 3 3
> arr
      [,1]  [,2] [,3]  [,4]
 [1,] 88.6 103.0 84.2  74.5
 [2,] 94.3  82.1 75.1  71.2
 [3,] 98.3  96.4 92.0  92.4
 [4,] 86.4  94.3 80.8  81.5
 [5,] 80.6  90.0 73.1  85.3
 [6,] 93.4  88.0 99.8  84.2
 [7,] 71.9  76.1 90.0  90.7
 [8,] 94.1  82.5 83.1 100.3
 [9,] 94.9  81.4 72.4    NA
[10,] 96.4  97.2 78.8    NA
[11,] 91.1  94.1 87.3    NA
[12,] 84.4  80.7 91.0    NA
[13,] 86.4  87.8 80.1    NA
[14,] 88.0  89.5 73.6    NA
[15,] 75.1  99.6 86.4    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/1mxa01225413621.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/2dlxt1225413621.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/3dv1u1225413621.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/4x5x91225413621.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,]  74.50 71.20 92.00 80.80 73.10 84.2 71.90  82.5 72.40  78.8  87.3 80.70
[2,]  79.35 73.15 92.20 81.15 76.85 86.1 74.00  82.8 76.90  87.6  89.2 82.55
[3,]  86.40 78.60 94.40 83.95 82.95 90.7 83.05  88.6 81.40  96.4  91.1 84.40
[4,]  95.80 88.20 97.35 90.35 87.65 96.6 90.35  97.2 88.15  96.8  92.6 87.70
[5,] 103.00 94.30 98.30 94.30 90.00 99.8 90.70 100.3 94.90  97.2  94.1 91.00
     [,13] [,14] [,15]
[1,] 80.10 73.60 75.10
[2,] 83.25 80.80 80.75
[3,] 86.40 88.00 86.40
[4,] 87.10 88.75 93.00
[5,] 87.80 89.50 99.60

$n
 [1] 4 4 4 4 4 4 4 4 3 3 3 3 3 3 3

$conf
        [,1]    [,2]    [,3]   [,4]   [,5]   [,6]    [,7]   [,8]    [,9]
[1,] 73.4045 66.7105 90.3315 76.682 74.418 82.405 70.1335 77.224 71.1376
[2,] 99.3955 90.4895 98.4685 91.218 91.482 98.995 95.9665 99.976 91.6624
         [,10]    [,11]   [,12]    [,13]   [,14]    [,15]
[1,]  88.00764 87.99847 79.7021 82.88798 80.7479 75.22539
[2,] 104.79236 94.20153 89.0979 89.91202 95.2521 97.57461

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(74.5, 79.35, 86.4, 95.8, 103, 71.2, 73.15,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/54hfy1225413621.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]
[1,] 75.1  76.10 72.40  71.20
[2,] 85.4  82.30 76.95  78.00
[3,] 88.6  89.50 83.10  84.75
[4,] 94.2  95.35 88.65  91.55
[5,] 98.3 103.00 99.80 100.30

$n
[1] 15 15 15  8

$conf
      [,1]    [,2]     [,3]     [,4]
[1,] 85.01 84.1762 78.32694 77.18078
[2,] 92.19 94.8238 87.87306 92.31922

$out
[1] 71.9

$group
[1] 1

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

Warning message:
In bxp(list(stats = c(75.1, 85.4, 88.6, 94.2, 98.3, 76.1, 82.3,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6ll9j1225413621.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,] 80.675 78.60 79.63750
[2,] 83.300 83.50 83.69375
[3,] 85.750 86.40 85.17500
[4,] 90.400 89.65 90.10000
[5,] 94.775 96.40 94.58750

$n
[1] 15 15 15

$conf
         [,1]     [,2]     [,3]
[1,] 82.85352 83.89108 82.56154
[2,] 88.64648 88.90892 87.78846

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(80.675, 83.3, 85.75, 90.4, 94.775, 78.6, 83.5,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> 
> system("convert tmp/1mxa01225413621.ps tmp/1mxa01225413621.png")
> system("convert tmp/2dlxt1225413621.ps tmp/2dlxt1225413621.png")
> system("convert tmp/3dv1u1225413621.ps tmp/3dv1u1225413621.png")
> system("convert tmp/4x5x91225413621.ps tmp/4x5x91225413621.png")
> system("convert tmp/54hfy1225413621.ps tmp/54hfy1225413621.png")
> system("convert tmp/6ll9j1225413621.ps tmp/6ll9j1225413621.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.288   0.820   2.806