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.0180678527088341
+ ,0.0562901166564763
+ ,-0.0383920228971653
+ ,0.0197362615228291
+ ,-0.0349783616211472
+ ,0.0068518670200992
+ ,0.0280842746777466
+ ,0.0446752709325745
+ ,0.0346027844881839
+ ,-0.04235896700803
+ ,-0.0704530877146328
+ ,-0.141169759207391
+ ,0.190464885503065
+ ,0.0160178050891860
+ ,0.136736106377406
+ ,0.00968031609540383
+ ,-0.0209143775407505
+ ,-0.0925769658093795
+ ,0.0384619588840446
+ ,0.0682195700863305
+ ,-0.134498105554547
+ ,0.0218255727200960
+ ,0.0776457687368569
+ ,0.0689225629397706
+ ,0.0740965027033676
+ ,-0.103338356927779
+ ,-0.0429672132927176
+ ,0.0365409556880891
+ ,-0.0414354395090189
+ ,-0.012330339132303
+ ,0.0286531468201482
+ ,-0.0263319233764757
+ ,0.0785732307259898
+ ,0.0237994722126526
+ ,0.0377483003142014
+ ,0.0507363385159117
+ ,-0.0223289969571086
+ ,-0.0483741687167987
+ ,-0.138702439764697
+ ,0.0658943416950095
+ ,-0.015257998756574
+ ,-0.0966648642708767
+ ,-0.021423000325382
+ ,-0.00651105347738536
+ ,0.00512035105957504
+ ,0.073029300952841
+ ,0.0573348071614868
+ ,0.172161649472528
+ ,-0.0969619108653857
+ ,-0.0984926589389618
+ ,0.0508824323815179
+ ,-0.0623193761665625
+ ,-0.137439210181595
+ ,0.13754570197851
+ ,-0.0237979817564337
+ ,0.00561691044717635
+ ,0.0328043530195327
+ ,0.00654130655901269
+ ,0.0601321062076455
+ ,-0.0639295914075303
+ ,0.0224173380303157
+ ,0.317941659607052
+ ,0.0225647602888595
+ ,-0.0206822092677761
+ ,0.00654057567239308
+ ,-0.0190310812241401
+ ,0.0506293024636681
+ ,0.0309034012148992
+ ,0.161159786003573
+ ,-0.0350860450190002
+ ,0.0766629141540984
+ ,0.104018416656661
+ ,-0.0253219589883051
+ ,0.0607550354735824
+ ,-0.115326979930159
+ ,0.0242361757084084
+ ,-0.0348254913157667)
> 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/1r8yp1259934947.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/2e3dq1259934948.tab") 
> 
> system("convert tmp/1r8yp1259934947.ps tmp/1r8yp1259934947.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.131   0.195   1.197