R version 2.8.1 (2008-12-22)
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(18.09
+ ,18.13
+ ,18
+ ,17.72
+ ,17.62
+ ,17.13
+ ,17.39
+ ,17.09
+ ,17.14
+ ,17.38
+ ,16.8
+ ,16.51
+ ,16.01
+ ,15.05
+ ,13.56
+ ,15.22
+ ,14.91
+ ,15.13
+ ,15.25
+ ,14.61
+ ,14.87
+ ,15.1
+ ,15.22
+ ,15.46
+ ,14.96
+ ,14
+ ,14.2
+ ,13.9
+ ,13.63
+ ,13.32
+ ,13.8
+ ,14.5
+ ,14.12
+ ,13.88
+ ,14.11
+ ,14.26
+ ,14.71
+ ,14.52
+ ,14.32
+ ,14.69
+ ,15.25
+ ,15.04
+ ,14.82
+ ,14.5
+ ,14.72
+ ,14.6
+ ,14.58
+ ,14
+ ,14.75
+ ,14.41
+ ,15.19
+ ,14.96
+ ,14.83
+ ,14.25
+ ,14.32
+ ,14.93
+ ,14.65
+ ,15.65
+ ,15.65
+ ,15.61
+ ,15.95
+ ,15.83
+ ,15.77
+ ,16.7
+ ,16.69
+ ,16.4
+ ,16.77
+ ,16.78
+ ,16.84
+ ,16.68
+ ,16.67
+ ,16.3
+ ,16.37
+ ,16.6
+ ,16.72
+ ,16.82
+ ,17.5
+ ,17.2
+ ,17.29
+ ,17.2
+ ,17.2
+ ,17.32
+ ,17.16
+ ,17.41
+ ,17.31
+ ,17.3
+ ,17.34
+ ,17.19
+ ,17.05
+ ,17.07
+ ,17.07
+ ,16.81
+ ,16.81
+ ,16.96
+ ,17.05
+ ,17
+ ,16.77
+ ,16.66
+ ,16.2
+ ,16.26
+ ,15.84
+ ,15.85
+ ,15.71
+ ,15.84
+ ,15.73
+ ,15.77
+ ,15.3
+ ,15.41
+ ,15.4
+ ,15.61
+ ,15
+ ,14.12
+ ,14.01
+ ,13.46
+ ,13.85
+ ,13.92
+ ,13.59
+ ,13.67
+ ,13.05
+ ,12.87
+ ,12.28
+ ,11.88
+ ,12.49
+ ,11.9
+ ,10.8
+ ,10.99
+ ,10.15
+ ,10.07
+ ,10.05
+ ,10.31
+ ,9.94
+ ,9.65
+ ,9.74
+ ,9.85
+ ,9.96
+ ,9.63
+ ,9.43
+ ,8.77
+ ,9.53
+ ,9.5
+ ,9.78
+ ,9.9
+ ,9.93
+ ,10.35
+ ,9.79
+ ,9.63
+ ,9.02
+ ,9.25
+ ,9.11
+ ,8.95
+ ,9.3
+ ,9.13
+ ,9.75
+ ,9.65
+ ,9.27
+ ,9.59
+ ,9.58
+ ,9.98
+ ,9.57
+ ,9.6
+ ,9.64
+ ,9.46
+ ,9.19
+ ,9.02
+ ,8.9
+ ,9.12
+ ,8.86
+ ,8.94
+ ,9
+ ,9.23
+ ,9.39
+ ,9.62
+ ,9.9
+ ,9.8
+ ,9.2
+ ,9.87
+ ,9.6
+ ,9.37
+ ,9.21
+ ,9.15
+ ,8.7
+ ,8.2
+ ,8.1
+ ,6.68
+ ,7.7
+ ,8.2
+ ,7.55
+ ,7.53
+ ,7.02
+ ,6.6
+ ,6
+ ,3.95
+ ,4.91
+ ,5.15
+ ,5.7
+ ,1.93
+ ,1.36
+ ,1.1
+ ,0.98
+ ,1
+ ,1.1
+ ,1.06
+ ,1.01
+ ,0.93
+ ,0.89
+ ,0.9
+ ,0.88
+ ,0.85
+ ,0.84
+ ,0.94
+ ,1
+ ,1.1
+ ,1.15
+ ,1.05
+ ,1.06
+ ,0.99
+ ,0.93
+ ,0.84
+ ,0.9
+ ,0.86
+ ,0.78
+ ,0.77
+ ,0.6
+ ,0.57
+ ,0.62
+ ,0.62
+ ,0.58
+ ,0.6
+ ,0.73
+ ,0.75
+ ,0.63
+ ,0.71
+ ,0.68
+ ,0.64
+ ,0.66
+ ,0.69
+ ,0.72
+ ,0.92
+ ,0.85
+ ,0.95
+ ,1
+ ,1.15
+ ,1.07
+ ,1.01
+ ,0.99
+ ,0.95
+ ,0.92
+ ,0.94
+ ,0.96
+ ,1.05
+ ,1.04
+ ,1.1
+ ,1.14
+ ,1.12
+ ,1.19
+ ,1.35
+ ,1.62
+ ,1.43
+ ,1.45
+ ,1.47
+ ,1.35
+ ,1.15
+ ,1.46
+ ,1.3
+ ,1.3
+ ,1.5
+ ,1.52
+ ,1.63
+ ,1.9
+ ,1.65
+ ,1.5
+ ,1.38
+ ,1.39)
> par3 = '0.01'
> par2 = '0.99'
> par1 = '0.01'
> ylab = 'Prijs aandelen'
> xlab = 'Percentile'
> 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/rcomp/tmp/1hgzl1236604730.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/2d9ou1236604730.tab") 
> 
> system("convert tmp/1hgzl1236604730.ps tmp/1hgzl1236604730.png")
> 
> 
> proc.time()
   user  system elapsed 
  0.918   0.195   1.001