R version 2.9.0 (2009-04-17) Copyright (C) 2009 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(-1.05389946067187 + ,8.43509140231974 + ,0.188905996583695 + ,3.60508040509410 + ,2.33890045064559 + ,3.5342847449073 + ,-2.35616784621537 + ,-2.58434087716323 + ,-3.33170015520434 + ,-2.88374983949467 + ,2.10904821238508 + ,2.34225423383044 + ,4.42007659454271 + ,3.25957134888508 + ,-4.29990795342411 + ,-7.50552236135862 + ,-3.85555216991755 + ,-1.25670836774721 + ,1.20468595761797 + ,-1.19516777025716 + ,-0.092352328092581 + ,2.85271681269754 + ,3.13317544841988 + ,-0.555250367187646 + ,-2.27118992591041 + ,-1.69562706989690 + ,-3.71012427223098 + ,-1.15008382338577 + ,-4.49982473326374 + ,0.788475952905313 + ,-1.65946218597853 + ,6.6203226344579 + ,-0.00452307730535890 + ,-3.22908667923240 + ,-11.9904635818427 + ,-9.53914407361147 + ,-0.497805850609024 + ,-2.6370901267526 + ,-2.36309505601764 + ,-6.99795791906111 + ,5.68874133680357 + ,-2.21400877345519 + ,2.40687767646164 + ,10.0202018713738 + ,-7.31369072515343 + ,-3.19489319711621 + ,3.90851094431611 + ,-8.08886314768259 + ,16.8436006273332 + ,-12.4941601309334 + ,-12.3872717022138 + ,-10.7834024725768 + ,-6.61524611429378 + ,-12.7318764430555 + ,-6.77900911212372 + ,-3.46954646116616 + ,4.08429373071817 + ,-0.666330544069827 + ,-1.5947107700646 + ,11.7970761300939 + ,-7.68674401462327 + ,1.58759796889509 + ,4.34675919202847 + ,3.88192565451631 + ,4.90515629658590 + ,7.87166553362415 + ,-10.7039267277302 + ,-7.70068020194925 + ,3.16737804702659 + ,17.0684504812639 + ,-8.01753595387331 + ,6.35569466769995 + ,-2.60032119619318 + ,-7.97322121819088 + ,0.474914236016726 + ,0.881907301984615 + ,-3.17822346734117 + ,0.0444015505138453 + ,-0.540800309532536 + ,-9.78873069920956 + ,-2.03414207641092 + ,-3.06230973047706 + ,1.50425532351333 + ,-4.85529042287201 + ,0.410034725742567 + ,1.66019566755176 + ,-3.96305513607336 + ,0.684387977516154 + ,-1.89847459363048 + ,0.0480867514532602 + ,2.5494509464279 + ,-4.64881767537846 + ,-8.74299996921144 + ,-1.79546410603364 + ,1.96397842738251 + ,2.12213572236232 + ,-1.53142399696699 + ,1.82461479057761 + ,9.22054040777745 + ,-6.97196023637948 + ,-0.774714717935781 + ,2.06191262103993 + ,-4.21434318421486 + ,3.62507335058573 + ,-7.78130386499525 + ,0.779658698747533 + ,1.09923671121991 + ,2.71556553075362 + ,0.592949276738042 + ,-0.147668502564047 + ,-0.51098067347573 + ,4.78740197523357 + ,-0.381795477377891 + ,6.12264276013345 + ,0.0170407340533487 + ,-5.84232347860017 + ,5.9362928585775 + ,1.63438579549496 + ,1.32958953797671 + ,5.98629928303477 + ,-0.86080284025978 + ,7.91773301835756 + ,-1.50650994380422 + ,1.29959077672656 + ,-1.40137718962717 + ,2.91895923162218 + ,-3.56386703823271 + ,2.81294345886286 + ,-1.37234582856103 + ,2.28507528223062 + ,0.368778799663316 + ,3.53899458912853 + ,-0.672881818102939 + ,-2.04949207670170 + ,4.7582411547085 + ,3.40759883241841 + ,-1.75357047183921 + ,6.43650469111989 + ,-4.87735079298824 + ,2.67160071852449 + ,-0.790168693416824 + ,-1.89184726380014 + ,5.23405134930391 + ,8.4187956558534 + ,1.47772197731134 + ,8.32038412831031 + ,-1.36243758125905 + ,0.61540688364956 + ,2.99391088629770 + ,1.31120486231028 + ,2.22017623479942 + ,4.05311913748494 + ,2.42598548606079 + ,-10.0442327608798 + ,2.50581147941573 + ,-4.15917403511201 + ,-2.57550667346294 + ,-0.943220932829732 + ,-4.50486844476557 + ,-2.31125584451949 + ,0.384996477730022 + ,-2.93184928296967 + ,-8.6451572501056 + ,6.04335104830048 + ,-3.26868508868629 + ,-4.1897972420899 + ,-5.38225615279879 + ,1.24055235166111 + ,1.67798526998128 + ,3.85443348518267 + ,-3.87633363880561 + ,2.57755673784785 + ,-2.02029315704359 + ,1.16159156094252 + ,3.40473201376932 + ,5.64166702268803 + ,0.484316119705291 + ,2.00590897726214 + ,-6.29360619012433 + ,0.97518977263484 + ,4.43495858251776 + ,-8.0316928697743 + ,-6.050081628142 + ,0.93920996998291 + ,-6.9290958098531 + ,-8.75697368715944 + ,5.22042345042553 + ,3.92195734145979 + ,-9.2752221487274 + ,-9.39272630932555 + ,3.27926888042508 + ,-2.16395195315198 + ,-0.114846011299871 + ,9.40782044021352 + ,4.61855066194025 + ,1.60014482089509 + ,-0.166039750338521 + ,8.63074847797661 + ,-6.99745247734323 + ,-1.10653747669793 + ,-0.660719142002788 + ,-7.97440068590113 + ,-5.39456284617709 + ,2.27166283556795 + ,-2.56373762068679 + ,2.59060662299397 + ,0.196259421316251 + ,-2.84028127388367 + ,-1.50087391103705 + ,-6.18844424174075 + ,2.06223064022139 + ,12.0650962486017 + ,-8.24295482447917 + ,-13.4234258875712 + ,-4.54195165660215 + ,5.11211162969743 + ,5.38833892797317 + ,-2.23669560410401 + ,3.24756194895173 + ,-1.58983540879586 + ,2.58401266545852 + ,1.76719274974144 + ,-1.79943603194111 + ,3.18297816315878 + ,-8.57087340678738 + ,-0.88672731763819 + ,-4.59633070019224 + ,-0.810395983388456 + ,1.30677966147441 + ,-0.349034718004156 + ,-3.28524758189259 + ,7.73818597276439 + ,-3.75692382905960 + ,1.6802375430731 + ,-1.84955407436030 + ,9.90226639896084 + ,-6.52546967455787 + ,-0.167349480887034 + ,-5.2101737890817 + ,-4.29278113587457 + ,6.6692695246446 + ,5.48546608659903 + ,3.9754670428778 + ,0.347764122921001 + ,6.30586024629087 + ,-4.31422114034296 + ,0.603491075416112 + ,-1.12276733977906 + ,16.5297269906827 + ,0.7909194825211 + ,-17.3994748701379 + ,-0.243858488379779 + ,0.948849321602561 + ,8.5736867852737 + ,5.71993211904582 + ,-4.05267342242049 + ,2.81852966627613 + ,2.20676380363063 + ,7.23191421396672 + ,-18.6661221293120 + ,0.679228558768283 + ,8.98545184658175 + ,12.5119051563469 + ,-2.41037447249911 + ,1.96673803275929 + ,1.31198971491305 + ,-1.39028666831682 + ,6.87614948403731 + ,-1.95721006029587 + ,3.23335075614169 + ,4.19955058192922 + ,-7.19772310922096 + ,0.181904717332093 + ,-4.39447492235539 + ,-5.67945087210234 + ,4.81563528246852 + ,3.67392630322103 + ,2.60572175708290 + ,1.36054243616509 + ,-5.79870104743685 + ,0.926111434920953 + ,2.91830262097045 + ,-5.72748271256833 + ,0.0576334433091275 + ,6.54456131713305 + ,11.2180456538184 + ,-4.45388768408582 + ,-5.82552730669147 + ,-8.03037991196091 + ,1.26540551070925 + ,4.64092998451394 + ,-0.562529063429161 + ,0.443534846641978 + ,-0.145637272098788 + ,0.244599218899082 + ,-12.1584630964438 + ,3.78501345141859 + ,-7.20776253455389 + ,-0.508351359377894 + ,3.30890852172676 + ,-1.76747312067812 + ,0.491076716628438 + ,0.484897273950985 + ,4.82589228080054 + ,4.96845543824267 + ,-0.116097163086708 + ,-6.75408234776817 + ,-8.32537009136505 + ,0.730429682170979 + ,-17.5661435290056 + ,-2.84428124121664 + ,-4.78767494637244 + ,8.60629559348975 + ,-6.03515939814025 + ,-1.66840601473531 + ,3.38031954019794 + ,-8.33089717209478 + ,-6.96020352149993 + ,8.7667308447824 + ,1.05185235473861 + ,-14.3034922181559 + ,9.48827315971597 + ,5.07603762425064 + ,9.55314429551763 + ,-2.07496632642002 + ,0.52022359831456 + ,-0.928136133947836 + ,2.26055304648137 + ,-8.67855216925862 + ,16.2988700985574 + ,-7.27735529169542 + ,-5.95123724758947 + ,-1.07870057092587 + ,2.61044790019159 + ,14.3768116417951 + ,9.31036004060947 + ,5.20942793618914 + ,9.94795595795348 + ,10.6529978513669 + ,-0.783751499049832 + ,-2.30555788621834 + ,2.98534099464617 + ,-5.51900894497926 + ,-2.41335128475402 + ,-7.01512801024298 + ,-5.12147170034238 + ,4.98731238533973) > par3 = '0.005' > 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) Loading required package: survival Loading required package: splines Attaching package: 'Hmisc' The following object(s) are masked from package:survival : untangle.specials 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/1832m1260626903.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] 197 > 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/2ju291260626903.tab") > > system("convert tmp/1832m1260626903.ps tmp/1832m1260626903.png") > > > proc.time() user system elapsed 1.872 0.244 1.969