x <- c(37729,48191,52498,57319,44377,48081,52597,53331,39587,46278,50365,57176,39251,47946,50427,54317,41210,50592,55728,59099,47519,53203,53882,55163,45255,50423,52161,54562,40971,48014,48440,44967,27218,30269,33234,36811,27745,31891,32398,34093,28358,29532,30769,32080,23951,34628,22978,35704,23090,22111,28925,35968,28963,34074,39160,51314,34527,40722,50609,52435)
#'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/1s68w1476392575.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/2h3wn1476392575.tab") 

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