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.350799688370536
+ ,-0.140769858461987
+ ,-0.176994002962568
+ ,0.562290817956577
+ ,0.144201690134944
+ ,-0.111030533067628
+ ,-0.243045038133344
+ ,-0.151508019282699
+ ,-0.396513842133135
+ ,0.0329753101728757
+ ,-0.0833088349918484
+ ,0.156956311022783
+ ,0.523621600794454
+ ,-0.216145977751341
+ ,-0.208014642793551
+ ,-0.341520041804527
+ ,-0.114593780635957
+ ,-0.558823753124076
+ ,0.155427432046773
+ ,0.138052779163656
+ ,-0.4090756185683
+ ,0.110316654295769
+ ,-0.250998024914370
+ ,0.329481615878391
+ ,0.346878506478169
+ ,0.323299830417798
+ ,-0.0310312875854553
+ ,-0.84926656008146
+ ,-0.77680329705359
+ ,0.250909338924582
+ ,0.0831309127438593
+ ,0.435693250316565
+ ,-0.130263895758751
+ ,-0.232129238451711
+ ,0.165342569504054
+ ,0.142248180166828
+ ,-0.280455237088636
+ ,0.155402490351535
+ ,-0.142784384867873
+ ,0.28524670566859
+ ,0.157418555269204
+ ,-0.172706911536110
+ ,-0.0149928813013760
+ ,0.306296873369522
+ ,-0.345207720981392
+ ,0.722409122191829
+ ,-0.158636568880998
+ ,-0.431031704131918
+ ,-0.308095207623409
+ ,0.447320732827413
+ ,0.447799821943335
+ ,0.0214235496574372
+ ,0.152471875346584
+ ,0.163942423990300
+ ,-0.00498382593447193
+ ,0.308356992481182
+ ,-0.230794586105762
+ ,-0.0542129596699259
+ ,-0.137751473869932
+ ,0.137775171895457
+ ,-0.208611181469940
+ ,0.530877388565287
+ ,-0.307691326359942
+ ,0.287643632777539
+ ,0.239698181689028
+ ,-0.0324235808931527
+ ,-0.110760321621385
+ ,0.187857496084540
+ ,-0.629717422499043
+ ,-0.464897710438974
+ ,0.381446921054368
+ ,0.388355710424064
+ ,-0.59380461268784
+ ,0.456215309938507
+ ,-0.35419239371042
+ ,0.278101258603134
+ ,-0.32386448817015
+ ,-0.115433965509661
+ ,-0.0483653154188056
+ ,-0.0424966476262035
+ ,-0.179962353470623
+ ,0.634520247547924
+ ,1.07894182435992
+ ,0.395987866990447
+ ,-0.219758520860607
+ ,0.361998505460662
+ ,0.660014172078838
+ ,-0.331005583090054
+ ,0.695781527796754
+ ,0.646124350414131
+ ,-0.0835481139905882
+ ,-0.763493350940545
+ ,0.0627068862169823
+ ,-0.0788278540379063
+ ,-1.01076620067800
+ ,-0.208872220642831
+ ,-0.515267639977992
+ ,0.319253920126987
+ ,-0.739240050461461
+ ,-0.0632470686715278
+ ,-0.182648154761123
+ ,-0.242076757818195
+ ,-0.704413262583499
+ ,0.851791451893186
+ ,-0.183365626294880)
> 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/1576g1260910805.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/2ty8q1260910805.tab") 
> try(system("convert tmp/1576g1260910805.ps tmp/1576g1260910805.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.191   0.225   5.511