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(6195800,6172725,6149325,6100900,6579950,6554600,6195800,5957250,5980325,5980325,6006000,6052150,6123975,6123975,6077825,5957250,6579950,6674850,6531525,6195800,6339450,6123975,6221150,6267625,6316050,6195800,6221150,6052150,6579950,6746675,6603350,6339450,6626425,6316050,6603350,6579950,6651775,6387875,6674850,6651775,7082400,6985225,6603350,6410950,6674850,6316050,6579950,6626425,6723600,6508450,6626425,6698250,6962150,6746675,6459700,6149325,6436625,5646875,6029075,6244225,6459700,6149325,6149325,6149325,6316050,6077825,5765175,5503550,5693350,4952350,5406375,5670275,5718700,5454800,5477875,5406375,5646875,5477875,5144750,4903925,5311150,4426825,5001100,5262725,5262725,4952350,4665375,4642300,4903925,4665375,4211675,3899025,4234750,3445325,4162925,4544800,4665375,4401475,4068025,4306575,4401475,4329650,3611725,3278600,3516825,2799225,3540225,3804125,4019275,3660475,3324750,3516825,3611725,3421925,2704325,2391675,2678650,1889225,2750475,3278600) > par3 = '0.1' > par2 = '0.9' > par1 = '0.1' > ylab = 'Omzet (EUR)' > xlab = 'Decielen' > main = 'Omzet product banden Michelin' > par3 <- '0.1' > par2 <- '0.9' > par1 <- '0.1' > #'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/1ff1l1437996243.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] 9 > 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/2gl401437996243.tab") > > try(system("convert tmp/1ff1l1437996243.ps tmp/1ff1l1437996243.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 0.756 0.168 0.923