R version 2.8.1 (2008-12-22)
Copyright (C) 2008 The R Foundation for Statistical Computing
ISBN 3-900051-07-0

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(10.505
+ ,10.505
+ ,10.315
+ ,10.145
+ ,10.815
+ ,13.945
+ ,13.915
+ ,13.525
+ ,12.775
+ ,12.165
+ ,11.635
+ ,11.045
+ ,10.295
+ ,10.045
+ ,9.625
+ ,9.245
+ ,8.885
+ ,8.545
+ ,8.225
+ ,7.895
+ ,7.165
+ ,6.055
+ ,5.175
+ ,4.415
+ ,3.705
+ ,2.085
+ ,-7.305
+ ,-8.005
+ ,-6.715
+ ,-2.825
+ ,-0.975
+ ,1.015
+ ,4.795
+ ,6.295
+ ,7.565
+ ,8.325
+ ,8.475
+ ,7.485
+ ,6.855
+ ,9.665
+ ,8.525
+ ,7.325
+ ,6.585
+ ,5.705
+ ,4.805
+ ,3.695
+ ,2.565
+ ,5.355
+ ,4.915
+ ,3.925
+ ,3.105
+ ,2.385
+ ,1.605
+ ,0.645
+ ,-0.305
+ ,-1.315
+ ,-2.295
+ ,1.495
+ ,3.925
+ ,3.135
+ ,1.975
+ ,1.005
+ ,4.365
+ ,4.305
+ ,3.255
+ ,2.135
+ ,1.235
+ ,0.125
+ ,-0.805
+ ,5.155
+ ,3.495
+ ,1.895
+ ,0.645
+ ,-0.615
+ ,-2.075
+ ,-0.385
+ ,-0.735
+ ,-2.125
+ ,-3.485
+ ,-4.985
+ ,-2.385
+ ,-1.865
+ ,-2.235
+ ,-3.695
+ ,-4.025
+ ,-5.455
+ ,-6.395
+ ,-6.425
+ ,-6.765
+ ,-7.415
+ ,-9.005
+ ,-7.515
+ ,-6.975
+ ,-8.395
+ ,-9.155
+ ,-8.385
+ ,-9.825
+ ,-11.565
+ ,-10.135
+ ,-13.325
+ ,-16.115
+ ,-14.435
+ ,-16.455
+ ,-18.135
+ ,-16.315
+ ,-18.185
+ ,-14.905
+ ,-14.925
+ ,-9.335
+ ,0.285
+ ,14.965
+ ,16.375
+ ,17.375
+ ,16.335
+ ,14.735
+ ,13.185
+ ,11.385
+ ,10.205
+ ,11.255
+ ,9.465
+ ,7.875
+ ,6.455
+ ,4.925
+ ,3.475
+ ,2.025
+ ,0.735
+ ,-0.575
+ ,-1.735
+ ,-2.775
+ ,-1.285
+ ,0.365
+ ,-0.625
+ ,-0.505
+ ,1.415
+ ,-0.285
+ ,0.535
+ ,-0.125
+ ,-1.785
+ ,-3.265
+ ,-4.705
+ ,-6.185
+ ,-4.815
+ ,-0.485
+ ,-1.915
+ ,-1.845
+ ,-4.875
+ ,-5.765
+ ,0.695
+ ,-0.995
+ ,-0.175
+ ,-0.705
+ ,-2.205
+ ,-3.755
+ ,-4.125
+ ,-3.615
+ ,-5.245
+ ,-6.995
+ ,-8.795
+ ,-6.815
+ ,-4.615
+ ,-6.305
+ ,-7.915
+ ,-9.405
+ ,-10.905
+ ,-4.615
+ ,-5.435
+ ,-6.845
+ ,-8.335
+ ,-9.615
+ ,-10.765
+ ,-9.525
+ ,-5.865
+ ,-7.025
+ ,-8.185
+ ,-9.295
+ ,-10.465
+ ,-11.965
+ ,-10.195
+ ,-10.805
+ ,-12.075
+ ,-12.075
+ ,-15.835
+ ,-17.695
+ ,-18.355
+ ,-19.005
+ ,-19.665
+ ,-23.465
+ ,-26.815
+ ,-27.375
+ ,-27.935
+ ,-28.475
+ ,-29.035
+ ,-29.535
+ ,-30.015
+ ,-30.485
+ ,-27.875
+ ,-27.545
+ ,-28.245
+ ,-28.745
+ ,-22.435
+ ,-22.985
+ ,-23.535
+ ,-23.815
+ ,-24.605
+ ,-25.075
+ ,-18.225
+ ,1.425
+ ,0.255
+ ,-0.585
+ ,0.375
+ ,-0.905
+ ,18.695
+ ,32.895
+ ,35.505
+ ,35.235
+ ,31.155
+ ,29.185
+ ,28.285
+ ,26.205
+ ,25.425
+ ,28.455
+ ,25.075
+ ,23.115
+ ,19.455
+ ,17.985
+ ,14.735
+ ,14.855
+ ,14.765
+ ,8.455
+ ,11.055
+ ,10.605
+ ,7.665
+ ,8.415
+ ,7.255
+ ,9.955
+ ,25.745
+ ,27.675
+ ,30.205
+ ,30.245
+ ,22.485
+ ,17.675
+ ,36.665
+ ,80.555
+ ,70.605
+ ,95.115
+ ,95.175
+ ,55.935
+ ,49.245
+ ,44.405
+ ,38.675
+ ,33.705
+ ,41.075
+ ,40.855
+ ,47.875
+ ,32.815
+ ,34.495
+ ,23.005
+ ,21.435
+ ,22.425
+ ,21.195
+ ,12.775
+ ,11.135
+ ,16.565
+ ,17.645
+ ,16.425
+ ,13.515
+ ,14.685
+ ,13.205
+ ,11.575
+ ,11.635
+ ,5.665
+ ,3.115
+ ,-1.085
+ ,-1.115
+ ,7.725
+ ,19.955
+ ,10.705
+ ,5.615
+ ,6.705
+ ,2.535
+ ,1.385
+ ,-4.005
+ ,-4.325
+ ,-8.315
+ ,-9.155
+ ,-8.015
+ ,-6.705
+ ,-8.225
+ ,-11.475
+ ,-17.815
+ ,-17.605
+ ,-14.405
+ ,-10.495
+ ,-21.155
+ ,-21.545
+ ,-20.925
+ ,-20.025
+ ,-25.885
+ ,-22.895
+ ,-23.825
+ ,-21.695
+ ,-23.165
+ ,-17.995
+ ,-14.925
+ ,-20.865
+ ,-23.455
+ ,-30.365
+ ,-32.205
+ ,-26.775
+ ,-26.775
+ ,-23.415
+ ,-30.545
+ ,-26.325
+ ,-27.705
+ ,-30.925
+ ,-28.325
+ ,-26.815
+ ,-30.245
+ ,-26.265
+ ,-29.195
+ ,-28.725
+ ,-29.935
+ ,-29.105
+ ,-27.615
+ ,-21.235
+ ,-20.955
+ ,-18.555
+ ,-20.585
+ ,-19.005
+ ,-13.585
+ ,-22.015
+ ,-21.035
+ ,-14.235
+ ,-15.385
+ ,-7.175
+ ,2.455
+ ,3.225
+ ,16.905
+ ,30.855
+ ,23.405
+ ,20.525
+ ,13.675
+ ,3.865
+ ,3.805
+ ,-0.955
+ ,8.815
+ ,12.575
+ ,12.645
+ ,27.955
+ ,18.545
+ ,12.495
+ ,13.665
+ ,10.195
+ ,3.205
+ ,4.445
+ ,10.575
+ ,6.865
+ ,5.445
+ ,17.625
+ ,25.405)
> par3 = '0.005'
> 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/www/rcomp/tmp/1y8zi1256159667.ps",horizontal=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/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,'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] 197
> 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/www/rcomp/tmp/2dc9c1256159667.tab") 
> 
> system("convert tmp/1y8zi1256159667.ps tmp/1y8zi1256159667.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.740   0.640   2.382