R version 3.0.2 (2013-09-25) -- "Frisbee Sailing"
Copyright (C) 2013 The R Foundation for Statistical Computing
Platform: x86_64-pc-linux-gnu (64-bit)

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(-350840
+ ,354950
+ ,342750
+ ,318610
+ ,303890
+ ,293480
+ ,286820
+ ,283650
+ ,276190
+ ,268670
+ ,275960
+ ,287040
+ ,291500
+ ,287540
+ ,272340
+ ,256570
+ ,246190
+ ,237340
+ ,235890
+ ,232790
+ ,226100
+ ,221110
+ ,221950
+ ,222990
+ ,232840
+ ,231380
+ ,221690
+ ,210190
+ ,202210
+ ,195210
+ ,194050
+ ,192360
+ ,187170
+ ,185930
+ ,188840
+ ,189350
+ ,194320
+ ,195850
+ ,188840
+ ,182380
+ ,177200
+ ,169140
+ ,171900
+ ,167320
+ ,164130
+ ,162550
+ ,166260
+ ,166010
+ ,177070
+ ,180800
+ ,179570
+ ,178130
+ ,178180
+ ,179430
+ ,179130
+ ,182220
+ ,183650
+ ,187110
+ ,197780
+ ,204000
+ ,221330
+ ,232500
+ ,237570
+ ,241100
+ ,244420
+ ,247240
+ ,254300
+ ,258970
+ ,262230
+ ,268880
+ ,277950
+ ,286630
+ ,299830
+ ,308090
+ ,311400
+ ,313520
+ ,310980
+ ,314910
+ ,310950
+ ,310720
+ ,310380
+ ,310570
+ ,314790
+ ,323400
+ ,335010
+ ,341630
+ ,340820
+ ,336280
+ ,325490
+ ,323750
+ ,317510
+ ,313890
+ ,308610
+ ,303720
+ ,303090
+ ,305140
+ ,304040
+ ,307100
+ ,304330
+ ,294710
+ ,286890
+ ,279050
+ ,271860
+ ,266710
+ ,259590
+ ,253830
+ ,250640
+ ,249140
+ ,250840
+ ,247590
+ ,237830
+ ,226380
+ ,217230
+ ,211420
+ ,207620
+ ,204310
+ ,197490
+ ,193580
+ ,192330
+ ,191970
+ ,196070
+ ,191940
+ ,185620
+ ,179410
+ ,173920
+ ,169190
+ ,166840
+ ,165170
+ ,161450
+ ,160830
+ ,163670
+ ,170830
+ ,182690
+ ,190940
+ ,197770
+ ,205090
+ ,210720
+ ,220210
+ ,229730
+ ,237070
+ ,241620
+ ,250370
+ ,258570
+ ,269860
+ ,283220
+ ,289610
+ ,281770
+ ,274700
+ ,267650
+ ,261380
+ ,260500
+ ,260730
+ ,254200
+ ,250450
+ ,253380
+ ,263740
+ ,276240
+ ,273820
+ ,265890
+ ,258400
+ ,253520
+ ,250710
+ ,252850
+ ,255260
+ ,251170
+ ,252500
+ ,257780
+ ,269900
+ ,291590
+ ,298870
+ ,295570
+ ,292100
+ ,290870
+ ,290580
+ ,297970
+ ,304010
+ ,304340
+ ,309850
+ ,322320
+ ,340170
+ ,369280
+ ,376690
+ ,379700
+ ,379520
+ ,377770
+ ,381560
+ ,394580
+ ,399320
+ ,400370
+ ,408200
+ ,419070
+ ,437730)
> par3 = '0.01'
> par2 = '0.99'
> par1 = '0.01'
> ylab = 'value'
> xlab = 'quantile'
> main = 'Harrell-Davis Quantiles'
> par3 <- '0.01'
> par2 <- '0.99'
> par1 <- '0.01'
> #'GNU S' R Code compiled by R2WASP v. 1.2.291 ()
> #Author: root
> #To cite this work: Wessa, (2012), Harrell-Davis Quantile Estimator (v1.0.13) in Free Statistics Software (v$_version), Office for Research Development and Education, URL http://www.wessa.net/rwasp_harrell_davis.wasp/
> #Source of accompanying publication: Office for Research, Development, and Education
> #
> par1 <- as(par1,'numeric')
> par2 <- as(par2,'numeric')
> par3 <- as(par3,'numeric')
> library(Hmisc)
Loading required package: grid
Loading required package: lattice
Loading required package: survival
Loading required package: splines
Loading required package: Formula

Attaching package: 'Hmisc'

The following objects 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/wessaorg/rcomp/tmp/1oyb91413030376.ps",horizontal=F,onefile=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/wessaorg/rcomp/createtable file can be downloaded at http://www.wessa.net/cretab
> load(file="/var/wessaorg/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/wessaorg/rcomp/tmp/24s091413030376.tab") 
> 
> try(system("convert tmp/1oyb91413030376.ps tmp/1oyb91413030376.png",intern=TRUE))
character(0)
> 
> 
> proc.time()
   user  system elapsed 
  0.989   0.102   1.090