source('/home/pw/wessanet/cretab')



myrfcuid = 'r0637275'

x <- c(14,19,17,17,15,20,15,19,15,15,19,16,20,18,15,14,20,16,16,16,10,19,19,16,15,18,17,19,17,14,19,20,19,16,15,16,18,16,15,17,14,20,19,13,16,16,18,18,16,17,19,16,19,13,16,13,12,17,17,17,16,16,14,16,13,16,14,20,12,18,14,19,18,14,18,15,14,17,19,13,19,18,20,15,15,15,20,15,19,18,18,15,17,12,18,19,20,13,17,16,18,18,14,15,12,17,14,18,17,17,20,16,14,15,18,20,17,17,17,15,17,18,17,20,15,16,18,15,18,20,19,14,16,15,17,18,20,17,18,15,16,11,15,18,17,16,19,18,15,17,19,18,19,16,16,16,14)
#'GNU S' R Code compiled by R2WASP v. 1.2.327 (Tue, 16 Aug 2016 15:27:00 +0200)
#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="/home/pw/wessanet/rcomp/tmp/1b04w1482304438.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()

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="/home/pw/wessanet/rcomp/tmp/2498g1482304438.tab") 

try(system("convert /home/pw/wessanet/rcomp/tmp/1b04w1482304438.ps /home/pw/wessanet/rcomp/tmp/1b04w1482304438.png",intern=TRUE))