R version 2.12.0 (2010-10-15)
Copyright (C) 2010 The R Foundation for Statistical Computing
ISBN 3-900051-07-0
Platform: i486-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(35.323
+ ,35.478
+ ,4.39
+ ,41.667
+ ,22.173
+ ,28.021
+ ,18.109
+ ,13.962
+ ,40.174
+ ,16.065
+ ,18.145
+ ,18.439
+ ,10.603
+ ,34.811
+ ,69.064
+ ,51.202
+ ,14.786
+ ,33.01
+ ,81.101
+ ,89.232
+ ,21.223
+ ,15.173
+ ,241.66
+ ,26.848
+ ,8.752
+ ,60.535
+ ,60.535
+ ,26.052
+ ,49.218
+ ,30.669
+ ,18.673
+ ,86
+ ,10.632
+ ,35.802
+ ,33.974
+ ,36.972
+ ,4.928
+ ,53.976
+ ,15.467
+ ,35.723
+ ,40.424
+ ,9.706
+ ,26.532
+ ,23.843
+ ,18.062
+ ,35.681
+ ,68.125
+ ,23.937
+ ,31.479
+ ,66.659
+ ,250.234
+ ,49.469
+ ,42.951
+ ,43.402
+ ,24.112
+ ,56.95
+ ,17.313
+ ,25.658
+ ,48.172
+ ,13.891
+ ,32.048
+ ,19.797
+ ,31.317
+ ,20.966
+ ,22.708
+ ,26.81
+ ,52.004
+ ,32.354
+ ,27.128
+ ,26.529
+ ,28.392
+ ,57.393
+ ,194.731
+ ,9.415
+ ,91.076
+ ,57.751
+ ,8.236
+ ,20.407
+ ,13.681
+ ,79.659
+ ,53.48
+ ,6.906
+ ,50.202
+ ,37.877
+ ,85.903
+ ,35.351
+ ,283.801
+ ,5.974
+ ,3.441
+ ,51.987
+ ,13.22
+ ,1.455
+ ,18.187
+ ,21.29
+ ,5.686
+ ,4.944
+ ,32.789
+ ,50.494
+ ,35.162
+ ,38.095
+ ,19.172
+ ,24.5
+ ,20.573
+ ,42.042
+ ,302.912
+ ,25.027
+ ,16.488
+ ,32.36
+ ,6.193
+ ,37.7
+ ,6.343
+ ,23.025
+ ,48.578
+ ,21.564
+ ,33.697
+ ,10.831
+ ,19.172
+ ,21.075
+ ,33.189
+ ,60.5
+ ,33.686
+ ,40.838
+ ,13.491
+ ,106.637
+ ,35.897
+ ,7.314
+ ,49.094
+ ,14.667
+ ,54.179
+ ,145.846
+ ,18.56
+ ,23.525
+ ,21.804
+ ,26.301
+ ,41.33
+ ,10.5
+ ,13.338
+ ,60.31
+ ,34.256
+ ,48.267
+ ,41.559
+ ,32.45
+ ,10.951
+ ,22.561
+ ,57.095
+ ,19.105
+ ,13.151
+ ,27.426
+ ,15.355
+ ,13.82
+ ,47.21
+ ,110.349
+ ,34.985
+ ,27.257
+ ,23.556
+ ,50.108
+ ,18.158
+ ,87.357
+ ,18.187
+ ,28.33
+ ,13.474
+ ,26.244)
> #'GNU S' R Code compiled by R2WASP v. 1.0.44 ()
> #Author: Prof. Dr. P. Wessa
> #To cite this work: AUTHOR(S), (YEAR), YOUR SOFTWARE TITLE (vNUMBER) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_YOURPAGE.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #Technical description: Write here your technical program description (don't use hard returns!)
> library(Hmisc)
Loading required package: survival
Loading required package: splines

Attaching package: 'Hmisc'

The following object(s) are masked from 'package:survival':

    untangle.specials

The following object(s) are masked from 'package:base':

    format.pval, round.POSIXt, trunc.POSIXt, units

> m <- mean(x)
> e <- median(x)
> postscript(file="/var/www/rcomp/tmp/170h11323192391.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()
null device 
          1 
> 
> #Note: the /var/www/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/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/www/rcomp/tmp/2mim11323192391.tab") 
> 
> try(system("convert tmp/170h11323192391.ps tmp/170h11323192391.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.720   0.020   0.726