x <- c(1280	1024
,1024	768
,1120	700
,1024	768
,1280	800
,1280	1024
,1280	800
,1024	768
,1280	800
,1280	1024
,1280	800
,1280	800
,1280	1024
,1688	949
,1440	900
,1600	1200
,1280	800
,1280	800
,1280	768
,1176	735
,1280	800
,1503	845
,1440	900
,1366	768
,1280	768
,1024	768
,1280	800
,2560	1440
,1280	768
,1024	768
,1280	1024
,1280	800
,1440	900
,1280	800
,1440	900
,1024	768
,1440	900
,1143	857
,1280	800
,1440	900
,1280	800
,1366	768
,1024	768
,1408	880
,1366	768
,1176	735
,1920	1200
,1257	785
,1280	800
,1280	800
,1440	900
,1680	1050
,1440	900
,1024	768
,1140	641
,1280	1024
,1280	800
,1280	800
,1280	800
,1280	800
,1440	900
,1280	800
,1152	864
,1280	1024
,1280	800
,1440	900
,1280	800
,1280	1024
,1440	900
,1280	800
,1280	800
,1440	900
,1280	800
,1280	1024
,1600	900
,1024	768
,1366	768
,1280	800
,1280	800
,1440	900
,1366	768
,1280	800
,1024	768
,1280	800
,1440	900
,1280	800
,1280	800
,1408	880
,1280	800
,1600	900
,1600	900
,1680	1050
,1440	900
,1440	900
,917	550
,1280	800
,1760	990
,1280	800
,1280	800
,1280	800
,1024	768
,1366	768
,1440	900
,1280	800
,1280	1024
,1920	1080
,1024	768
,1024	768
,1600	900
,1117	698
,1440	900
,983	737
,1024	768
,1024	640
,1280	800
,1440	900
,1280	800
,1280	800
,1280	800
,1440	900
,1280	800
,1024	768
,1024	768
,1152	864
,1280	768
,1024	768
,1366	768
,1680	1050
,1680	1050
,1280	800
,1366	768
,1024	768
,1440	900
,1024	768
,1280	800
,1280	800
,1280	800
,1024	768
,1280	800)
par3 = '0.01'
par2 = '0.99'
par1 = '0.01'
ylab = 'value'
xlab = 'quantile'
main = 'Harrell-Davis Quantiles'
#'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!)
par1 <- as(par1,'numeric')
par2 <- as(par2,'numeric')
par3 <- as(par3,'numeric')
library(Hmisc)
myseq <- seq(par1, par2, par3)
hd <- hdquantile(x, probs = myseq, se = TRUE, na.rm = FALSE, names = TRUE, weights=FALSE)
postscript(file="/var/wessaorg/rcomp/tmp/1khhk1317463575.ps",horizontal=F,onefile=F,pagecentre=F,paper="special",width=8.3333333333333,height=5.5555555555556) 
plot(myseq,hd,col=2,main=main,xlab=xlab,ylab=ylab)
grid()
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,'Harrell-Davis Quantiles',3,TRUE)
a<-table.row.end(a)
a<-table.row.start(a)
a<-table.element(a,'quantiles',header=TRUE)
a<-table.element(a,'value',header=TRUE)
a<-table.element(a,'standard error',header=TRUE)
a<-table.row.end(a)
length(hd)
for (i in 1:length(hd))
{
a<-table.row.start(a)
a<-table.element(a,as(labels(hd)[i],'numeric'),header=TRUE)
a<-table.element(a,as.matrix(hd[i])[1,1])
a<-table.element(a,as.matrix(attr(hd,'se')[i])[1,1])
a<-table.row.end(a)
}
a<-table.end(a)
table.save(a,file="/var/wessaorg/rcomp/tmp/2vnnx1317463575.tab") 

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