R version 3.2.3 (2015-12-10) -- "Wooden Christmas-Tree"
Copyright (C) 2015 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-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(283411
+ ,234585
+ ,207695
+ ,204000
+ ,220791
+ ,222335
+ ,303510
+ ,201602
+ ,218279
+ ,207253
+ ,217639
+ ,231043
+ ,302988
+ ,221415
+ ,213989
+ ,190507
+ ,199397
+ ,201077
+ ,261048
+ ,241508
+ ,194637
+ ,240044
+ ,253534
+ ,269220
+ ,270820
+ ,354121
+ ,315233
+ ,297098
+ ,258289
+ ,276246
+ ,246652
+ ,242945
+ ,210838
+ ,246343
+ ,234636
+ ,229011
+ ,238756
+ ,283696
+ ,226656
+ ,229790
+ ,219831
+ ,232331
+ ,315447
+ ,251031
+ ,297019
+ ,208424
+ ,311626
+ ,352320
+ ,375214
+ ,419501
+ ,664867
+ ,483142
+ ,312717
+ ,228228
+ ,230978
+ ,218033
+ ,225566
+ ,207708
+ ,241861
+ ,208381
+ ,209071
+ ,214514
+ ,195868
+ ,190208
+ ,187651
+ ,215934
+ ,213012
+ ,236845
+ ,154595
+ ,193485
+ ,231875
+ ,192259
+ ,191609
+ ,212911
+ ,238596
+ ,203033
+ ,258272
+ ,314681
+ ,289789
+ ,297541
+ ,239578
+ ,207798
+ ,254091
+ ,242544
+ ,218067
+ ,217463
+ ,229438
+ ,216485
+ ,238410
+ ,223108
+ ,221267
+ ,224802
+ ,233630
+ ,235134
+ ,309660
+ ,253229
+ ,293530
+ ,273323
+ ,297578
+ ,341589
+ ,389911
+ ,398685
+ ,629637
+ ,522502
+ ,358941
+ ,252783
+ ,236585
+ ,224995
+ ,225721
+ ,223338
+ ,228952
+ ,212759
+ ,200270
+ ,211925
+ ,191206
+ ,200511
+ ,205944
+ ,289288
+ ,186209
+ ,197515
+ ,182038
+ ,181606
+ ,223337
+ ,201213
+ ,220700
+ ,194043
+ ,224593
+ ,249907
+ ,238613
+ ,266200
+ ,274197
+ ,288601
+ ,242662
+ ,232105
+ ,267842
+ ,219037
+ ,198404
+ ,203910
+ ,209685
+ ,214274
+ ,215376
+ ,217711
+ ,224819
+ ,217142
+ ,221445
+ ,213483
+ ,282063
+ ,236194
+ ,297296
+ ,255885
+ ,264502
+ ,311580
+ ,416595
+ ,356762
+ ,537286
+ ,532855)
> par3 = '0.01'
> par2 = '0.99'
> par1 = '0.01'
> ylab = 'value'
> xlab = 'quantile'
> main = 'Harrell-Davis Quantiles'
> par3 <- '0.01'
> par2 <- '0.99'
> par1 <- '0.01'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa, (2012), Harrell-Davis Quantile Estimator (v1.0.13) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_harrell_davis.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> par1 <- as(par1,'numeric')
> par2 <- as(par2,'numeric')
> par3 <- as(par3,'numeric')
> library(Hmisc)
Loading required package: grid
Loading required package: lattice
Loading required package: survival
Loading required package: Formula
Loading required package: ggplot2

Attaching package: 'Hmisc'

The following objects are masked from 'package:base':

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

> 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/1vxvr1456782681.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()
null device 
          1 
> 
> #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)
[1] 99
> 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/2xj931456782681.tab") 
> 
> try(system("convert tmp/1vxvr1456782681.ps tmp/1vxvr1456782681.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.178   0.169   1.345