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 = '20'
> #'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] 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 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3
> arr
       [,1] [,2] [,3]  [,4]
 [1,] 109.2 82.1 80.8 100.3
 [2,]  88.6 68.0 73.1    NA
 [3,]  94.3 96.4 99.8    NA
 [4,]  98.3 94.3 90.0    NA
 [5,]  86.4 90.0 83.1    NA
 [6,]  80.6 88.0 72.4    NA
 [7,] 104.1 76.1 78.8    NA
 [8,] 108.2 82.5 87.3    NA
 [9,]  93.4 81.4 91.0    NA
[10,]  71.9 66.5 80.1    NA
[11,]  94.1 97.2 73.6    NA
[12,]  94.9 94.1 86.4    NA
[13,]  96.4 80.7 74.5    NA
[14,]  91.1 70.5 71.2    NA
[15,]  84.4 87.8 92.4    NA
[16,]  86.4 89.5 81.5    NA
[17,]  88.0 99.6 85.3    NA
[18,]  75.1 84.2 69.9    NA
[19,] 109.7 75.1 84.2    NA
[20,] 103.0 92.0 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/1jkzq1193403405.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/25a651193403405.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/311l11193403405.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/4lb441193403405.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,]  80.80 68.00 94.30 90.00 83.10 72.4  76.10  82.50 81.4  66.5 73.60 86.40
[2,]  81.45 70.55 95.35 92.15 84.75 76.5  77.45  84.90 86.2  69.2 83.85 90.25
[3,]  91.20 73.10 96.40 94.30 86.40 80.6  78.80  87.30 91.0  71.9 94.10 94.10
[4,] 104.75 80.85 98.10 96.30 88.20 84.3  91.45  97.75 92.2  76.0 95.65 94.50
[5,] 109.20 88.60 99.80 98.30 90.00 88.0 104.10 108.20 93.4  80.1 97.20 94.90
     [,13] [,14] [,15] [,16] [,17] [,18]  [,19]  [,20]
[1,] 74.50 70.50  84.4 81.50 85.30 69.90  75.10  90.70
[2,] 77.60 70.85  86.1 83.95 86.65 72.50  79.65  91.35
[3,] 80.70 71.20  87.8 86.40 88.00 75.10  84.20  92.00
[4,] 88.55 81.15  90.1 87.95 93.80 79.65  96.95  97.50
[5,] 96.40 91.10  92.4 89.50 99.60 84.20 109.70 103.00

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

$conf
        [,1]    [,2]     [,3]     [,4]     [,5]     [,6]     [,7]     [,8]
[1,]  72.793 63.7042 93.89141 90.51431 83.25286 73.48474 66.02901 75.57806
[2,] 109.607 82.4958 98.90859 98.08569 89.54714 87.71526 91.57099 99.02194
         [,9]    [,10]     [,11]   [,12]    [,13]   [,14]    [,15]    [,16]
[1,] 85.52672 65.69695  83.33588 90.2231 70.71126 61.8042 84.15115 82.75115
[2,] 96.47328 78.10305 104.86412 97.9769 90.68874 80.5958 91.44885 90.04885
        [,17]    [,18]    [,19]    [,20]
[1,] 81.47767 68.57767 68.41871 86.38989
[2,] 94.52233 81.62233 99.98129 97.61011

$out
numeric(0)

$group
numeric(0)

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

Warning message:
In bxp(list(stats = c(80.8, 81.45, 91.2, 104.75, 109.2, 68, 70.55,  :
  some notches went outside hinges ('box'): maybe set notch=FALSE
> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/56ec81193403406.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,]  71.90 66.50 69.90 100.3
[2,]  86.40 78.40 74.05 100.3
[3,]  93.75 86.00 82.30 100.3
[4,] 100.65 93.05 88.65 100.3
[5,] 109.70 99.60 99.80 100.3

$n
[1] 20 20 20  1

$conf
        [,1]     [,2]     [,3]  [,4]
[1,] 88.7155 80.82417 77.14184 100.3
[2,] 98.7845 91.17583 87.45816 100.3

$out
numeric(0)

$group
numeric(0)

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

> dev.off()
null device 
          1 
> postscript(file="/var/www/html/rcomp/tmp/6gm241193403406.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,] 72.83333 71.20 72.6000
[2,] 82.10000 79.70 81.7375
[3,] 88.25000 86.85 88.2000
[4,] 92.23333 91.60 91.7375
[5,] 96.83333 96.40 96.7250

$n
[1] 20 20 20

$conf
        [,1]     [,2]     [,3]
[1,] 84.6699 82.64574 84.66701
[2,] 91.8301 91.05426 91.73299

$out
numeric(0)

$group
numeric(0)

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

> dev.off()
null device 
          1 
> 
> system("convert tmp/1jkzq1193403405.ps tmp/1jkzq1193403405.png")
> system("convert tmp/25a651193403405.ps tmp/25a651193403405.png")
> system("convert tmp/311l11193403405.ps tmp/311l11193403405.png")
> system("convert tmp/4lb441193403405.ps tmp/4lb441193403405.png")
> system("convert tmp/56ec81193403406.ps tmp/56ec81193403406.png")
> system("convert tmp/6gm241193403406.ps tmp/6gm241193403406.png")
> 
> 
> proc.time()
   user  system elapsed 
  2.462   1.379   2.660