x <- c(93.2 ,96 ,95.2 ,77.1 ,70.9 ,64.8 ,70.1 ,77.3 ,79.5 ,100.6 ,100.7 ,107.1 ,95.9 ,82.8 ,83.3 ,80 ,80.4 ,67.5 ,75.7 ,71.1 ,89.3 ,101.1 ,105.2 ,114.1 ,96.3 ,84.4 ,91.2 ,81.9 ,80.5 ,70.4 ,74.8 ,75.9 ,86.3 ,98.7 ,100.9 ,113.8 ,89.8 ,84.4 ,87.2 ,85.6 ,72 ,69.2 ,77.5 ,78.1 ,94.3 ,97.7 ,100.2 ,116.4 ,97.1 ,93 ,96 ,80.5 ,76.1 ,69.9 ,73.6 ,92.6 ,94.2 ,93.5 ,108.5 ,109.4 ,105.1 ,92.5 ,97.1 ,81.4 ,79.1 ,72.1 ,78.7 ,87.1 ,91.4 ,109.9 ,116.3 ,113 ,100 ,84.8 ,94.3 ,87.1 ,90.3 ,72.4 ,84.9 ,92.7 ,92.2 ,114.9 ,112.5 ,118.3 ,106 ,91.2 ,96.6 ,96.3 ,88.2 ,70.2 ,86.5 ,88.2 ,102.8 ,119.1 ,119.2 ,125.1 ,106.1 ,102.1 ,105.2 ,101 ,84.3 ,87.5 ,92.7 ,94.4 ,113 ,113.9 ,122.9 ,132.7 ,106.9 ,96.6 ,127.3 ,98.2 ,100.2 ,89.4 ,95.3 ,104.2 ,106.4 ,116.2 ,135.9 ,134 ,104.6 ,107.1 ,123.5 ,98.8 ,98.6 ,90.6 ,89.1 ,105.2 ,114 ,122.1 ,138 ,142.2 ,116.4 ,112.6 ,123.8 ,103.6 ,113.9 ,98.6 ,95 ,116 ,113.9 ,127.5 ,131.4 ,145.9 ,131.5 ,131 ,130.5 ,118.9 ,114.3 ,85.7 ,104.6 ,105.1 ,117.3 ,142.5 ,140 ,159.8 ,131.2 ,125.4 ,126.5 ,119.4 ,113.5 ,98.7 ,114.5 ,113.8 ,133.1 ,143.4 ,137.3 ,165.2 ,126.9 ,124 ,135.7 ,130 ,109.4 ,117.8 ,120.3 ,121 ,132.3 ,142.9 ,147.4 ,175.9 ,132.6 ,123.7 ,153.3 ,134 ,119.6 ,116.2 ,118.6 ,130.7 ,129.3 ,144.4 ,163.2 ,179.4 ,128.1 ,138.4 ,152.7 ,120 ,140.5 ,116.2 ,121.4 ,127.8 ,143.6 ,157.6 ,166.2 ,182.3 ,153.1 ,147.6 ,157.7 ,137.2 ,151.5 ,98.7 ,145.8 ,151.7 ,129.4 ,174.1 ,197 ,193.9 ,164.1 ,142.8 ,157.9 ,159.2 ,162.2 ,123.1 ,130 ,150.1 ,169.4 ,179.7 ,182.1 ,194.3 ,161.4 ,169.4 ,168.8 ,158.1 ,158.5 ,135.3 ,149.3 ,143.4 ,142.2 ,188.4 ,166.2 ,199.2 ,182.7 ,145.2 ,182.1 ,158.7 ,141.6 ,132.6 ,139.6 ,147 ,166.6 ,157 ,180.4 ,210.2 ,159.8 ,157.8 ,168.2 ,158.4 ,152 ,142.2 ,137.2 ,152.6 ,166.8 ,165.6 ,198.6 ,201.5 ,170.7 ,164.4 ,179.7 ,157 ,168 ,139.3 ,138.6 ,153.4 ,138.9 ,172.1 ,198.4 ,217.8 ,173.7 ,153.8 ,175.6 ,147.1 ,160.3 ,135.2 ,148.8 ,151 ,148.2 ,182.2 ,189.2 ,183.1 ,170 ,158.4 ,176.1 ,156.2 ,153.2 ,117.9 ,149.8 ,156.6 ,166.7 ,156.8 ,158.6 ,210.8 ,203.6 ,175.2 ,168.7 ,155.9 ,147.3 ,137 ,141.1 ,167.4 ,160.2 ,191.9 ,174.4 ,208.2 ,159.4 ,161.1 ,172.1 ,158.4 ,114.6 ,159.6 ,159.7 ,159.4 ,160.7 ,165.5 ,205 ,205.2 ,141.6 ,148.1 ,184.9 ,132.5 ,137.3 ,135.5 ,121.7 ,166.1 ,146.8 ,162.8 ,186.8 ,185.5 ,151.5 ,158.1 ,143 ,151.2 ,147.6 ,130.7 ,137.5 ,146.1 ,133.6 ,167.9 ,181.9 ,202 ,166.5 ,151.3 ,146.2 ,148.3 ,144.7 ,123.6 ,151.6 ,133.9 ,137.4 ,181.6 ,182 ,190 ,161.2 ,155.5 ,141.9 ,164.6 ,136.2 ,126.8 ,152.5 ,126.6 ,150.1 ,186.3 ,147.5 ,200.4 ,177.2 ,127.4 ,177.1 ,154.4 ,135.2 ,126.4 ,147.3 ,140.6 ,152.3 ,151.2 ,172.2 ,215.3 ,154.1 ,159.3 ,160.4 ,151.9 ,148.4 ,139.6 ,148.2 ,153.5 ,145.1 ,183.7 ,210.5 ,203.3 ,153.3 ,144.3 ,169.6 ,143.7 ,160.1 ,135.6 ,141.8 ,159.9 ,145.7 ,183.5 ,198.2 ,186.8 ,172 ,150.6 ,163.3 ,153.7 ,152.9 ,135.5 ,148.5 ,148.4 ,133.6 ,194.1 ,208.6 ,197.3 ,164.4 ,148.1 ,152 ,144.1 ,155 ,124.5 ,153 ,146 ,138 ,190 ,192 ,192 ,147 ,133 ,163 ,150 ,129 ,131 ,145 ,137 ,138 ,168 ,176 ,188 ,139 ,143 ,150 ,154 ,137 ,129 ,128 ,140 ,143 ,151 ,177 ,184 ,151 ,134 ,164 ,126 ,131 ,125 ,127 ,143 ,143 ,160 ,190 ,182 ,138 ,136 ,152 ,127 ,151 ,130 ,119 ,153) 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) 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/1n4dj1268077315.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() #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) 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/2w5yh1268077315.tab") try(system("convert tmp/1n4dj1268077315.ps tmp/1n4dj1268077315.png",intern=TRUE))