R version 2.9.0 (2009-04-17)
Copyright (C) 2009 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(0.312330033
+ ,0.009995004
+ ,0.295646257
+ ,-0.303089192
+ ,0.317705894
+ ,-0.774667791
+ ,-0.792455285
+ ,-0.737256151
+ ,-0.792615022
+ ,-0.103254692
+ ,-0.483337454
+ ,-0.956913885
+ ,-0.846353349
+ ,-0.765080851
+ ,-0.181880976
+ ,-0.567384004
+ ,0.458949967
+ ,-1.084592482
+ ,-0.253979902
+ ,0.701791272
+ ,0.818236832
+ ,0.456052157
+ ,0.074104935
+ ,0.275791233
+ ,-0.477487044
+ ,-0.332111908
+ ,-0.706076846
+ ,-0.342106113
+ ,-0.538193936
+ ,-0.826917102
+ ,-0.039307061
+ ,-0.343163219
+ ,-0.761713269
+ ,0.243897876
+ ,-0.544587845
+ ,-0.368157046
+ ,-0.456027302
+ ,0.098597925
+ ,-0.663962256
+ ,-1.004123798
+ ,-0.257525654
+ ,0.031550618
+ ,0.882878401
+ ,0.63655095
+ ,0.571850909
+ ,0.254923185
+ ,1.128351059
+ ,0.595043508
+ ,0.18626974
+ ,0.583583804
+ ,0.161376723
+ ,0.067941595
+ ,0.853194934
+ ,0.772618858
+ ,0.016638009
+ ,-0.187641356
+ ,-0.133382955
+ ,1.142275157
+ ,0.777932181
+ ,0.416700154
+ ,0.467389716
+ ,0.488730153
+ ,0.569296476
+ ,0.663207928
+ ,0.388595347
+ ,0.693972659
+ ,0.142277688
+ ,0.681523087
+ ,0.255301554
+ ,0.511672353
+ ,0.654016047
+ ,0.251032312
+ ,0.530130259
+ ,0.152422278
+ ,0.348828865
+ ,0.353194248
+ ,0.129707075
+ ,-0.15822426
+ ,-0.319798168
+ ,0.949650332
+ ,-0.205002364
+ ,0.531089433
+ ,0.325019689
+ ,0.388826704
+ ,-0.160516627
+ ,0.062700579
+ ,0.780427355
+ ,0.22376635
+ ,-0.01722534
+ ,-0.124478324
+ ,-0.167209445
+ ,-0.127058382
+ ,-0.43590623
+ ,-0.316612043
+ ,0.260307372
+ ,-0.454875497
+ ,-0.842489883
+ ,-0.392208712
+ ,-0.201080773
+ ,-0.798172208
+ ,-0.366462777
+ ,-0.74323943
+ ,-0.127907657
+ ,-0.48577647
+ ,0.172894288
+ ,0.033915608
+ ,-0.672168605
+ ,-0.410915983)
> 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/html/rcomp/tmp/1prhy1260565443.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/html/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/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/html/rcomp/tmp/2enat1260565443.tab") 
> 
> system("convert tmp/1prhy1260565443.ps tmp/1prhy1260565443.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.171   0.212   1.272