R version 2.8.0 (2008-10-20)
Copyright (C) 2008 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(-2.08459015
+ ,-2.084139482
+ ,-2.088701546
+ ,-2.091544359
+ ,-2.089390446
+ ,-2.086901271
+ ,-2.085599746
+ ,-2.08742381
+ ,-2.093887342
+ ,-2.095893691
+ ,-2.091975472
+ ,-2.09244865
+ ,-2.09565
+ ,-2.095845274
+ ,-2.093161272
+ ,-2.092853846
+ ,-2.0928748
+ ,-2.093683784
+ ,-2.094184837
+ ,-2.097200613
+ ,-2.095753659
+ ,-2.096785246
+ ,-2.096368447
+ ,-2.095345725
+ ,-2.092687879
+ ,-2.093193124
+ ,-2.092260024
+ ,-2.092495694
+ ,-2.093859451
+ ,-2.095758586
+ ,-2.094597509
+ ,-2.091013085
+ ,-2.093223477
+ ,-2.093852802
+ ,-2.094982547
+ ,-2.09398957
+ ,-2.095647939
+ ,-2.095228638
+ ,-2.096383324
+ ,-2.095880221
+ ,-2.095007621
+ ,-2.093853488
+ ,-2.094144186
+ ,-2.095746868
+ ,-2.095680347
+ ,-2.097695049
+ ,-2.098450058
+ ,-2.099198044
+ ,-2.099089189
+ ,-2.099476174
+ ,-2.099957143
+ ,-2.100178567
+ ,-2.102059091
+ ,-2.102647002
+ ,-2.103123235
+ ,-2.102623138
+ ,-2.103165422
+ ,-2.102972118
+ ,-2.10155247
+ ,-2.101613961
+ ,-2.098782587
+ ,-2.099326322
+ ,-2.098527451
+ ,-2.098528763
+ ,-2.099315042
+ ,-2.099877728
+ ,-2.099215526
+ ,-2.099215526
+ ,-2.098943382
+ ,-2.096225737
+ ,-2.096803441
+ ,-2.09587561
+ ,-2.097095122
+ ,-2.098978875
+ ,-2.101873612
+ ,-2.10299589
+ ,-2.101251943
+ ,-2.100396996
+ ,-2.099605882
+ ,-2.099829232
+ ,-2.101086072
+ ,-2.102830941
+ ,-2.102186801
+ ,-2.103663557
+ ,-2.107623577
+ ,-2.105919365
+ ,-2.104433762
+ ,-2.101850959
+ ,-2.099995847
+ ,-2.101438052
+ ,-2.104421277
+ ,-2.107563303
+ ,-2.105547558
+ ,-2.105337666
+ ,-2.106823077
+ ,-2.106771686
+ ,-2.104862884
+ ,-2.104725766
+ ,-2.103426316
+ ,-2.103803075
+ ,-2.105536364
+ ,-2.106065803
+ ,-2.106540827
+ ,-2.105206859
+ ,-2.103772297
+ ,-2.103229913
+ ,-2.10204982
+ ,-2.099946178
+ ,-2.100715803
+ ,-2.099925575
+ ,-2.098369388
+ ,-2.099502446
+ ,-2.098890826
+ ,-2.09588369
+ ,-2.095679817
+ ,-2.095280711
+ ,-2.092904552
+ ,-2.0949
+ ,-2.094259233
+ ,-2.097355919
+ ,-2.098822457
+ ,-2.098621805
+ ,-2.09945
+ ,-2.096040089
+ ,-2.094370821
+ ,-2.095459389
+ ,-2.096877887
+ ,-2.096779607
+ ,-2.098821569
+ ,-2.098386036
+ ,-2.099474156
+ ,-2.099239853
+ ,-2.095739844
+ ,-2.097076987
+ ,-2.096919932
+ ,-2.101134867
+ ,-2.103471774
+ ,-2.105650845
+ ,-2.103740231
+ ,-2.103842308
+ ,-2.103856835
+ ,-2.103715916
+ ,-2.104636842
+ ,-2.103348111
+ ,-2.104129852
+ ,-2.101606387
+ ,-2.101515055
+ ,-2.102852381
+ ,-2.104966414
+ ,-2.102236802
+ ,-2.102861738
+ ,-2.101456604
+ ,-2.101417436
+ ,-2.102461205
+ ,-2.101163504
+ ,-2.101746298
+ ,-2.102459663
+ ,-2.10401215
+ ,-2.104627018
+ ,-2.100961453
+ ,-2.097760465
+ ,-2.099432191
+ ,-2.098763216
+ ,-2.099480452
+ ,-2.100553775
+ ,-2.096750575
+ ,-2.098197384
+ ,-2.09949193
+ ,-2.099561482
+ ,-2.098159479
+ ,-2.095900912
+ ,-2.097377314
+ ,-2.098537392
+ ,-2.096913575
+ ,-2.095743102
+ ,-2.095890991
+ ,-2.095678578
+ ,-2.098798876
+ ,-2.098370166
+ ,-2.096946168
+ ,-2.098127395
+ ,-2.099275
+ ,-2.102392226
+ ,-2.101208204
+ ,-2.100205445
+ ,-2.101636633
+ ,-2.102835982
+ ,-2.102719383
+ ,-2.104532644
+ ,-2.104588323
+ ,-2.103367777
+ ,-2.102719615
+ ,-2.10263523
+ ,-2.102163986
+ ,-2.101952402
+ ,-2.102130467
+ ,-2.104317793
+ ,-2.104281533
+ ,-2.103940887
+ ,-2.105030435
+ ,-2.103340487
+ ,-2.100194627
+ ,-2.097172727
+ ,-2.096563789
+ ,-2.094979776
+ ,-2.094277682
+ ,-2.091227095
+ ,-2.092315254
+ ,-2.087907598
+ ,-2.085437748
+ ,-2.085015789
+ ,-2.087271956
+ ,-2.087045249
+ ,-2.090159306
+ ,-2.096816667
+ ,-2.098950791
+ ,-2.10025343
+ ,-2.10188172
+ ,-2.107328453
+ ,-2.107414499
+ ,-2.111427663
+ ,-2.114428613
+ ,-2.116807895
+ ,-2.113098032)
> par3 = '0.05'
> par2 = '0.95'
> par1 = '0.05'
> 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)

Attaching package: 'Hmisc'


	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/freestat/rcomp/tmp/1ziqh1255689656.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/freestat/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/www/html/freestat/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] 19
> 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/freestat/rcomp/tmp/2ma2j1255689656.tab") 
> 
> system("convert tmp/1ziqh1255689656.ps tmp/1ziqh1255689656.png")
> 
> 
> proc.time()
   user  system elapsed 
  1.135   0.259   1.239