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. Natural language support but running in an English locale 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(-1.8207657603853E-14 + ,-0.011286801534618 + ,0.011286801534612 + ,-1.8207657603853E-14 + ,-1.8207657603853E-14 + ,-0.011286801534618 + ,1.1990408665952E-14 + ,1.1990408665952E-14 + ,0.011286801534612 + ,-1.8207657603853E-14 + ,-1.8207657603853E-14 + ,-0.011286801534618 + ,1.1990408665952E-14 + ,1.1990408665952E-14 + ,-0.013714500670788 + ,-0.011574203281101 + ,0.011574203281136 + ,-8.8817841970013E-16 + ,-0.011574203281101 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,0.011574203281136 + ,-8.8817841970013E-16 + ,-8.8817841970013E-16 + ,-8.8817841970013E-16 + ,-8.8817841970013E-16 + ,-0.011574203281101 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,0.011574203281136 + ,-8.8817841970013E-16 + ,-8.8817841970013E-16 + ,-8.8817841970013E-16 + ,-8.8817841970013E-16 + ,-0.011574203281101 + ,0.011574203281136 + ,-8.8817841970013E-16 + ,-8.8817841970013E-16 + ,-8.8817841970013E-16 + ,-0.011574203281101 + ,3.5971225997855E-14 + ,3.5971225997855E-14 + ,0.011574203281136 + ,-8.8817841970013E-16 + ,-0.005770356461301 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,-0.011641575015464 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,0.0070011954590217 + ,3.8191672047105E-14 + ,3.8191672047105E-14 + ,0.004640379556538 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,-0.0046403795564642 + ,0.004640379556538 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,-0.0046403795564642 + ,3.8191672047105E-14 + ,0.004640379556538 + ,3.5527136788005E-14 + ,3.5527136788005E-14 + ,-0.54696467038186 + ,-2.2204460492503E-16 + ,0.3364722366212 + ,0.16430305129129 + ,1.0658141036402E-14 + ,0.029852963149711 + ,2.9309887850104E-14 + ,2.9309887850104E-14 + ,-0.0082693914824707 + ,3.5971225997855E-14 + ,-0.007142887512364 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,1.5987211554602E-14 + ,0.007142887512416 + ,-0.007142887512364 + ,0.0011940299926159 + ,-0.0071856596608919 + ,0.0071856596608826 + ,7.9936057773011E-15 + ,-0.0071856596608919 + ,0.0071856596608826 + ,-0.0071856596608919 + ,-1.7319479184152E-14 + ,-1.7319479184152E-14 + ,-1.7319479184152E-14 + ,-1.7319479184152E-14 + ,0.0071856596608826 + ,-0.0059844583151918 + ,4.8849813083507E-14 + ,4.8849813083507E-14 + ,4.8849813083507E-14 + ,-0.012077441489751 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,2.1760371282653E-14 + ,-0.0085418457129784 + ,-1.5543122344752E-14 + ,-1.5543122344752E-14 + ,-0.0073801072976156 + ,0.0073801072976067 + ,-1.5543122344752E-14 + ,-1.5543122344752E-14 + ,-1.5543122344752E-14 + ,-1.5543122344752E-14 + ,-0.0073801072976156 + ,6.6613381477509E-15 + ,6.6613381477509E-15 + ,6.6613381477509E-15 + ,6.6613381477509E-15 + ,0.0073801072976067 + ,-0.0073801072976156 + ,6.6613381477509E-15 + ,6.6613381477509E-15 + ,6.6613381477509E-15 + ,0.0073801072976067 + ,-1.5543122344752E-14 + ,-0.0061463008487155 + ,-0.012407107047022 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,0.012407107047032 + ,-2.1316282072803E-14 + ,-2.1316282072803E-14 + ,-2.1316282072803E-14 + ,-0.012407107047022 + ,0.012407107047032 + ,-2.1316282072803E-14 + ,-2.1316282072803E-14 + ,-0.012407107047022 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,3.1974423109205E-14 + ,0.012407107047032 + ,-2.1316282072803E-14 + ,-2.1316282072803E-14 + ,-2.1316282072803E-14 + ,-2.1316282072803E-14 + ,0.0036923118870789 + ,-4.8849813083507E-14 + ,-4.8849813083507E-14 + ,-4.8849813083507E-14 + ,-4.8849813083507E-14 + ,-0.012361096823549 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,0.012361096823525 + ,-4.8849813083507E-14 + ,-4.8849813083507E-14 + ,-0.012361096823549 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,0.012361096823525 + ,-4.8849813083507E-14 + ,-4.8849813083507E-14 + ,-0.0061614489440491 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-0.0124379712922 + ,0.012437971292217 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-0.0124379712922 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,0.012437971292217 + ,-0.0124379712922 + ,1.6875389974302E-14 + ,0.012437971292217 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-0.0124379712922 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,0.012437971292217 + ,-4.4408920985006E-16 + ,-0.0124379712922 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,1.6875389974302E-14 + ,0.012437971292217 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-0.0124379712922 + ,1.6875389974302E-14 + ,0.012437971292217 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-0.0124379712922 + ,-0.0050188311119834 + ,-4.0856207306206E-14 + ,0.0050188311119594 + ,0.012437971292217 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-0.0111871893906 + ,0.011187189390564 + ,-4.4408920985006E-16 + ,-0.0111871893906 + ,-3.5527136788005E-14 + ,-3.5527136788005E-14 + ,-3.5527136788005E-14 + ,-3.5527136788005E-14 + ,0.011187189390564 + ,-0.0111871893906 + ,-3.5527136788005E-14 + ,-3.5527136788005E-14 + ,0.011187189390564 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-4.4408920985006E-16 + ,-0.0111871893906 + ,0.011187189390564 + ,-0.0061996478795003 + ,2.5313084961454E-14 + ,-0.47499117075677 + ,0.4749911707568 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,-0.012515807931875 + ,0.012515807931856 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,-0.012515807931875 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,0.012515807931856 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,-0.012515807931875 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,0.012515807931856 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,-0.012515807931875 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,0.012515807931856 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,2.5313084961454E-14 + ,-0.012515807931875 + ,0.012515807931856 + ,2.5313084961454E-14 + ,-0.012515807931875 + ,-4.4408920985006E-14 + ,-0.084038927104744 + ,-4.5519144009631E-14 + ,-4.5519144009631E-14 + ,0.084038927104654 + ,-4.4408920985006E-14 + ,-4.4408920985006E-14 + ,-0.0012602396122445 + ,4.3076653355456E-14 + ,-0.010139503852457 + ,-1.6875389974302E-14 + ,0.010139503852483 + ,4.3076653355456E-14 + ,4.3076653355456E-14 + ,4.3076653355456E-14 + ,4.3076653355456E-14 + ,4.3076653355456E-14) > 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) 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/19c5f1290173011.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] 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/freestat/rcomp/tmp/2cvml1290173011.tab") > > try(system("convert tmp/19c5f1290173011.ps tmp/19c5f1290173011.png",intern=TRUE)) character(0) > > > proc.time() user system elapsed 2.283 0.295 15.613