R version 3.0.3 (2014-03-06) -- "Warm Puppy" Copyright (C) 2014 The R Foundation for Statistical Computing Platform: i686-pc-linux-gnu (32-bit) 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(-2,0,1,-3,-3,-5,-7,-7,-5,-13,-16,-20,-18,-21,-20,-16,-14,-12,-10,-3,-4,-4,-1,-8,-10,-11,-7,-2,-6,-4,0,2,2,5,8,8,5,10,6,6,9,5,5,-4,-5,-1,-8,-8,-13,-18,-8,-8,-6,-5,-11,-14,-12,-13,-19,-21,-22,-13,-21,-17,-15,-14,-11,-8,-3,-2,-1,1) > par1 = '12' > par1 <- '12' > #'GNU S' R Code compiled by R2WASP v. 1.2.291 () > #Author: root > #To cite this work: Wessa P., (2012), Mean Plot (v1.0.4) 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 > # > par1 <- as.numeric(par1) > (n <- length(x)) [1] 72 > (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] 6 6 6 6 6 6 6 6 6 6 6 6 > arr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -2 -18 -10 5 -13 -22 NA [2,] 0 -21 -11 10 -18 -13 NA [3,] 1 -20 -7 6 -8 -21 NA [4,] -3 -16 -2 6 -8 -17 NA [5,] -3 -14 -6 9 -6 -15 NA [6,] -5 -12 -4 5 -5 -14 NA [7,] -7 -10 0 5 -11 -11 NA [8,] -7 -3 2 -4 -14 -8 NA [9,] -5 -4 2 -5 -12 -3 NA [10,] -13 -4 5 -1 -13 -2 NA [11,] -16 -1 8 -8 -19 -1 NA [12,] -20 -8 8 -8 -21 1 NA > darr [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] 2 -3 -1 5 -5 9 NA [2,] 1 1 4 -4 10 -8 NA [3,] -4 4 5 0 0 4 NA [4,] 0 2 -4 3 2 2 NA [5,] -2 2 2 -4 1 1 NA [6,] -2 2 4 0 -6 3 NA [7,] 0 7 2 -9 -3 3 NA [8,] 2 -1 0 -1 2 5 NA [9,] -8 0 3 4 -1 1 NA [10,] -3 3 3 -7 -6 1 NA [11,] -4 -7 0 0 -2 2 NA [12,] 2 -2 -3 -5 -1 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/wessaorg/rcomp/tmp/1ntsd1394796752.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/276hs1394796752.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/3ruyk1394796752.ps",horizontal=F,onefile=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/wessaorg/rcomp/tmp/4n7qg1394796752.ps",horizontal=F,onefile=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,] -22.0 -21 -21.0 -17.0 -15 -14 -11.0 -14.0 -5.0 -13 -19.0 -21 [2,] -18.0 -18 -20.0 -16.0 -14 -12 -11.0 -8.0 -5.0 -13 -16.0 -20 [3,] -11.5 -12 -7.5 -5.5 -6 -5 -8.5 -5.5 -4.5 -3 -4.5 -8 [4,] -2.0 0 1.0 -2.0 -3 -4 0.0 -3.0 -3.0 -1 -1.0 1 [5,] 5.0 10 6.0 6.0 9 5 5.0 2.0 -3.0 5 8.0 8 $n [1] 6 6 6 6 6 6 6 6 6 6 6 6 $conf [,1] [,2] [,3] [,4] [,5] [,6] [1,] -21.820517 -23.6105814 -21.045678 -14.530452 -13.095355 -10.1602584 [2,] -1.179483 -0.3894186 6.045678 3.530452 1.095355 0.1602584 [,7] [,8] [,9] [,10] [,11] [,12] [1,] -15.595355 -8.725161 -5.790065 -10.740388 -14.175484 -21.545678 [2,] -1.404645 -2.274839 -3.209935 4.740388 5.175484 5.545678 $out [1] 2 -12 $group [1] 9 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-22, -18, -11.5, -2, 5, -21, -18, -12, 0, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/5r7hu1394796752.ps",horizontal=F,onefile=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,] -5.0 -8 -4 0 -4 -6 -9 -1 -1.0 -7 -7 -5 [2,] -3.0 -4 0 0 -2 -2 -3 -1 -1.0 -6 -4 -3 [3,] 0.5 1 2 2 1 1 1 1 0.5 -1 -1 -2 [4,] 5.0 4 4 2 2 3 3 2 3.0 3 0 -1 [5,] 9.0 10 5 3 2 4 7 5 4.0 3 2 2 $n [1] 6 6 6 6 6 6 6 6 6 6 6 5 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -4.660258 -4.160258 -0.5801292 0.7099354 -1.580129 -2.225161 -2.870194 [2,] 5.660258 6.160258 4.5801292 3.2900646 3.580129 4.225161 4.870194 [,8] [,9] [,10] [,11] [,12] [1,] -0.9350969 -2.080129 -6.805291 -3.580129 -3.413195 [2,] 2.9350969 3.080129 4.805291 1.580129 -0.586805 $out [1] -4 -8 $group [1] 4 9 $names [1] "1" "2" "3" "4" "5" "6" "7" "8" "9" "10" "11" "12" Warning message: In bxp(list(stats = c(-5, -3, 0.5, 5, 9, -8, -4, 1, 4, 10, -4, 0, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/6opfa1394796752.ps",horizontal=F,onefile=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,] -20.0 -21 -11.0 -8.0 -21.0 -22.0 NA [2,] -10.0 -17 -6.5 -4.5 -16.0 -16.0 NA [3,] -5.0 -11 -1.0 5.0 -12.5 -12.0 NA [4,] -2.5 -4 3.5 6.0 -8.0 -2.5 NA [5,] 1.0 -1 8.0 10.0 -5.0 1.0 NA $n [1] 12 12 12 12 12 12 0 $conf [,1] [,2] [,3] [,4] [,5] [,6] [,7] [1,] -8.4208 -16.929387 -5.561067 0.2108795 -16.148854 -18.157441 NA [2,] -1.5792 -5.070613 3.561067 9.7891205 -8.851146 -5.842559 NA $out numeric(0) $group numeric(0) $names [1] "1" "2" "3" "4" "5" "6" NA Warning message: In bxp(list(stats = c(-20, -10, -5, -2.5, 1, -21, -17, -11, -4, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > postscript(file="/var/wessaorg/rcomp/tmp/7leap1394796752.ps",horizontal=F,onefile=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,] -10.000000 -12.00 -10.3750 [2,] -8.083333 -8.25 -9.0625 [3,] -6.000000 -5.75 -7.6875 [4,] -5.666667 -4.75 -6.1250 [5,] -4.500000 -3.00 -4.1250 $n [1] 12 12 12 $conf [,1] [,2] [,3] [1,] -7.102258 -7.346373 -9.027313 [2,] -4.897742 -4.153627 -6.347687 $out numeric(0) $group numeric(0) $names [1] "mean" "median" "midrange" Warning message: In bxp(list(stats = c(-10, -8.08333333333333, -6, -5.66666666666667, : some notches went outside hinges ('box'): maybe set notch=FALSE > dev.off() null device 1 > > try(system("convert tmp/1ntsd1394796752.ps tmp/1ntsd1394796752.png",intern=TRUE)) character(0) > try(system("convert tmp/276hs1394796752.ps tmp/276hs1394796752.png",intern=TRUE)) character(0) > try(system("convert tmp/3ruyk1394796752.ps tmp/3ruyk1394796752.png",intern=TRUE)) character(0) > try(system("convert tmp/4n7qg1394796752.ps tmp/4n7qg1394796752.png",intern=TRUE)) character(0) > try(system("convert tmp/5r7hu1394796752.ps tmp/5r7hu1394796752.png",intern=TRUE)) character(0) > try(system("convert tmp/6opfa1394796752.ps tmp/6opfa1394796752.png",intern=TRUE)) character(0) > try(system("convert tmp/7leap1394796752.ps tmp/7leap1394796752.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 4.892 0.958 5.828