R version 2.13.0 (2011-04-13)
Copyright (C) 2011 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)
> 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)
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

> 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/1a6wk1322790169.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/21fqr1322790169.tab") 
> 
> try(system("convert tmp/1a6wk1322790169.ps tmp/1a6wk1322790169.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.138   0.085   1.545