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.651124187
+ ,0.677462043
+ ,0.487584459
+ ,0.35573232
+ ,0.467396208
+ ,0.5847797
+ ,0.638410776
+ ,0.562504993
+ ,0.226959834
+ ,0.096608848
+ ,0.346780629
+ ,0.32212801
+ ,0.124353661
+ ,0.116113081
+ ,0.298476107
+ ,0.320797584
+ ,0.320938168
+ ,0.272939992
+ ,0.242238517
+ ,0.029850007
+ ,0.144567637
+ ,0.073595478
+ ,0.105558217
+ ,0.179094701
+ ,0.347499306
+ ,0.323095391
+ ,0.379457113
+ ,0.373362194
+ ,0.295091587
+ ,0.173059076
+ ,0.25395154
+ ,0.466386248
+ ,0.343237476
+ ,0.307402512
+ ,0.235057345
+ ,0.300157187
+ ,0.190740285
+ ,0.223666633
+ ,0.144958306
+ ,0.188004361
+ ,0.248384872
+ ,0.323674763
+ ,0.305392718
+ ,0.200707615
+ ,0.209268012
+ ,0.060984282
+ ,-0.002231387
+ ,-0.069789538
+ ,-0.059028506
+ ,-0.091004634
+ ,-0.135694478
+ ,-0.152499926
+ ,-0.366844839
+ ,-0.450497399
+ ,-0.515235087
+ ,-0.442133219
+ ,-0.514521352
+ ,-0.490012558
+ ,-0.303601016
+ ,-0.310362096
+ ,-0.01578921
+ ,-0.056493576
+ ,0.020719536
+ ,0.026208378
+ ,-0.044823852
+ ,-0.098827417
+ ,-0.034225275
+ ,-0.034225275
+ ,-0.008453333
+ ,0.213738024
+ ,0.172428731
+ ,0.242113198
+ ,0.15114142
+ ,-0.002440472
+ ,-0.306723255
+ ,-0.453671587
+ ,-0.227233275
+ ,-0.130523048
+ ,-0.047519648
+ ,-0.068848603
+ ,-0.201545891
+ ,-0.421545609
+ ,-0.333128768
+ ,-0.5438376
+ ,-1.551038145
+ ,-0.983778668
+ ,-0.662441384
+ ,-0.269523484
+ ,-0.060934863
+ ,-0.217716491
+ ,-0.651830893
+ ,-1.502188008
+ ,-0.879597713
+ ,-0.831909014
+ ,-1.225813299
+ ,-1.208658248
+ ,-0.732318801
+ ,-0.700744426
+ ,-0.473166048
+ ,-0.522039199
+ ,-0.853039535
+ ,-0.979494757
+ ,-1.108876863
+ ,-0.772168602
+ ,-0.503727917
+ ,-0.414905721
+ ,-0.25200225
+ ,-0.014582013
+ ,-0.095054366
+ ,-0.009089297
+ ,0.138738084
+ ,0.036439646
+ ,0.09359806
+ ,0.33509075
+ ,0.349531434
+ ,0.381745218
+ ,0.536301301
+ ,0.413095585
+ ,0.456051135
+ ,0.237104602
+ ,0.113336816
+ ,0.136326454
+ ,0.062846398
+ ,0.352464439
+ ,0.471155508
+ ,0.402330127
+ ,0.302224828
+ ,0.309743661
+ ,0.143310479
+ ,0.182417007
+ ,0.084235603
+ ,0.107550334
+ ,0.393103741
+ ,0.295544151
+ ,0.310972647
+ ,-0.077335332
+ ,-0.385262055
+ ,-0.798163939
+ ,-0.423084824
+ ,-0.43872323
+ ,-0.438448813
+ ,-0.414863511
+ ,-0.573976594
+ ,-0.354421361
+ ,-0.478284249
+ ,-0.114965004
+ ,-0.103481844
+ ,-0.276699126
+ ,-0.62908592
+ ,-0.185555629
+ ,-0.269043576
+ ,-0.085442224
+ ,-0.079763405
+ ,-0.21001755
+ ,-0.04659717
+ ,-0.113020436
+ ,-0.203010969
+ ,-0.435472926
+ ,-0.5442069
+ ,-0.013832988
+ ,0.293298418
+ ,0.149790454
+ ,0.214316775
+ ,0.150318618
+ ,0.043147573
+ ,0.386091161
+ ,0.2715729
+ ,0.155364229
+ ,0.151210935
+ ,0.279908767
+ ,0.457889092
+ ,0.352014075
+ ,0.257325144
+ ,0.391406149
+ ,0.481729908
+ ,0.472232849
+ ,0.488285072
+ ,0.242085404
+ ,0.281978442
+ ,0.399431551
+ ,0.305312483
+ ,0.205735682
+ ,-0.136076382
+ ,0.012342482
+ ,0.123502113
+ ,-0.034258798
+ ,-0.190092633
+ ,-0.171264219
+ ,-0.460842425
+ ,-0.468821782
+ ,-0.260515601
+ ,-0.164445256
+ ,-0.152273865
+ ,-0.087961961
+ ,-0.060005142
+ ,-0.082178379
+ ,-0.411831932
+ ,-0.403900354
+ ,-0.343582591
+ ,-0.545014941
+ ,-0.246134632
+ ,0.148529096
+ ,0.420561811
+ ,0.470669144
+ ,0.583865257
+ ,0.635121216
+ ,0.823492271
+ ,0.76770808
+ ,0.995614642
+ ,1.102989155
+ ,1.121969728
+ ,1.030205654
+ ,1.043444445
+ ,0.897900099
+ ,0.49074351
+ ,0.313599762
+ ,0.188325762
+ ,0.004373587
+ ,-1.130992438
+ ,-1.163640132)
> par3 = '0.01'
> 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/1w79y1256476135.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] 99
> 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/2s7rw1256476136.tab") 
> 
> system("convert tmp/1w79y1256476135.ps tmp/1w79y1256476135.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.917   0.191   1.287