x <- c(97.91,98.51,98.54,98.52,98.66,98.53,98.71,98.92,98.96,99.25,99.32,99.41,99.36,99.58,99.77,99.77,100.03,100.2,100.24,100.1,100.03,100.18,100.29,100.41,100.6,100.75,100.79,100.44,100.29,100.34,100.46,100.12,100.06,100.28,100.28,100.4,100.61,100.89,100.73,101.12,101.16,101.33,101.37,101.61,101.85,102.27,102.28,102.23,102.42,102.53,103.47,103.53,103.77,103.74,103.93,103.97,103.68,103.86,103.97,104.05)
#'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
#Author: root
#To cite this work: Wessa P., (2012), Mean versus Median (v1.0.2) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_meanversusmedian.wasp/
#Source of accompanying publication: Office for Research, Development, and Education
#
library(Hmisc)
m <- mean(x)
e <- median(x)
postscript(file="/var/wessaorg/rcomp/tmp/1agjv1489050218.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
op <- par(mfrow=c(2,1))
mydensity1 <- density(x,kernel='gaussian',na.rm=TRUE)
plot(mydensity1,main='Density Plot - Gaussian Kernel',xlab='Median (0 -> full line) | Mean (0 -> dashed line)',ylab='density')
abline(v=e,lty=1)
abline(v=m,lty=5)
grid()
myseq <- seq(0.01, 0.99, 0.01)
hd <- hdquantile(x, probs = myseq, se = TRUE, na.rm = FALSE, names = TRUE, weights=FALSE)
plot(myseq,hd,col=2,main='Harrell-Davis Quantiles',xlab='quantiles',ylab='Median (0 -> full) | Mean (0 -> dashed)')
abline(h=m,lty=5)
abline(h=e,lty=1)
grid()
par(op)
dev.off()

#Note: the /var/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
load(file="/var/wessaorg/rcomp/createtable")

a<-table.start()
a<-table.row.start(a)
a<-table.element(a,'Median versus Mean',2,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'mean',header=TRUE)
a<-table.element(a,mean(x))
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'median',header=TRUE)
a<-table.element(a,median(x))
a<-table.row.end(a)
a<-table.end(a)
table.save(a,file="/var/wessaorg/rcomp/tmp/2k0y91489050218.tab") 

try(system("convert tmp/1agjv1489050218.ps tmp/1agjv1489050218.png",intern=TRUE))