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(3.98079081276633e-05
+ ,-0.000973976133257226
+ ,9.54105981046696e-05
+ ,0.0025088460252085
+ ,-0.00317190862576495
+ ,-0.000294031330328119
+ ,0.000968982944777646
+ ,-0.000851197318024926
+ ,-0.00138793660546551
+ ,0.000768157042933171
+ ,-0.00043345520475556
+ ,0.00366870846041696
+ ,-0.00115605912816505
+ ,-0.000726387151364303
+ ,0.00158304848369840
+ ,-0.000705551333861491
+ ,-0.000686845735372371
+ ,0.000856524956680535
+ ,0.00112951149233031
+ ,-0.000633986568076334
+ ,0.000485377794173421
+ ,0.00350968792492686
+ ,0.00110875144650834
+ ,-0.000454222753280748
+ ,-0.00070519397282043
+ ,-0.000564138898205549
+ ,-0.00323625658579197
+ ,-0.000207775895461774
+ ,0.000121573704094360
+ ,0.00043576463403868
+ ,-0.00105292199456297
+ ,-0.0002403967976497
+ ,-0.00118907538442198
+ ,0.000625841684237297
+ ,0.00143341373589394
+ ,-0.00192988553325272
+ ,-0.00157492400816634
+ ,-0.000826348774317081
+ ,0.00200273506082489
+ ,0.00234881104580967
+ ,-0.00111885056021927
+ ,-0.00081068847426387
+ ,-0.000153003860193700
+ ,-0.000584636837683116
+ ,0.00160325078754813
+ ,-0.00183253836358638
+ ,-2.44819684689276e-06
+ ,-0.000771460372823094
+ ,-0.00104001978218138
+ ,4.48050729071727e-05
+ ,-0.00169203184569203
+ ,0.000174462273812442
+ ,-0.00114098465483701
+ ,0.000769566716964055
+ ,-0.000818911240646756
+ ,-0.000828750325242718
+ ,0.00018261275748915
+ ,-0.00159799022419187
+ ,0.00152164245558621
+ ,0.000884169332101304
+ ,-0.00180485032467462
+ ,-0.000145630936549182
+ ,-0.00147324017244386
+ ,0.00060792316144348
+ ,0.000905629038218515
+ ,-0.00169040652912523
+ ,-3.12105291171649e-05
+ ,-0.000906676164921031
+ ,0.000514094616829656
+ ,0.000837486441965654
+ ,0.000374199742847721
+ ,-0.000878094545087892
+ ,-0.001182893501353
+ ,0.000646111491672511
+ ,-0.000599857591795931
+ ,-0.00144802173609301
+ ,-0.000104582917017275
+ ,0.000492759150213751
+ ,-0.000931164919957729
+ ,0.000371424202406294
+ ,0.00140975891549500
+ ,6.6750425077598e-05
+ ,-0.000277340783067839
+ ,-0.000455194022166956
+ ,0.00149545439655812
+ ,-0.00131297831713598
+ ,-0.000172568590443257
+ ,-0.000606874399482126
+ ,0.000336885535519071
+ ,-0.000642830360379811
+ ,-0.000512599342982913
+ ,0.000704580414167671
+ ,-0.000945254989964433
+ ,-0.000475584476934928
+ ,-0.000848558711677428
+ ,0.000314088890732562
+ ,-0.000617706260337561
+ ,-0.000400848410183710
+ ,-0.000838652504832539
+ ,-0.000333202351853059
+ ,-0.00104162067143271
+ ,-0.000146313497017913
+ ,-8.59352300135553e-05
+ ,0.00137121135042611
+ ,0.000958241811236665
+ ,0.00339310785749926
+ ,0.00259450104170567
+ ,0.00237818486786294
+ ,-0.00116307131503507
+ ,0.000681404236029412
+ ,-0.000892729488794187
+ ,-0.000758226783340658
+ ,-0.00148062583439785
+ ,-0.00170642385731143
+ ,0.00129224214825562
+ ,-0.00154330598635992
+ ,0.000855094334823046)
> 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/17dbw1260910259.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/2owtj1260910259.tab") 
> try(system("convert tmp/17dbw1260910259.ps tmp/17dbw1260910259.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  1.201   0.236   1.861