Wessa.net

Send output to:


Data:
41086 39690 43129 37863 35953 29133 24693 22205 21725 27192 21790 13253 37702 30364 32609 30212 29965 28352 25814 22414 20506 28806 22228 13971 36845 35338 35022 34777 26887 23970 22780 17351 21382 24561 17409 11514 31514 27071 29462 26105 22397 23843 21705 18089 20764 25316 17704 15548 28029 29383 36438 32034 22679 24319 18004 17537 20366 22782 19169 13807 29743 25591 29096 26482 22405 27044 17970 18730 19684 19785 18479 10698


Sample Range: (leave blank to include all observations)
From:
To:

lowest quantile


highest quantile


step size


Chart options


Title:
Label y-axis:
Label x-axis:

R Code
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) bitmap(file='test1.png') plot(myseq,hd,col=2,main=main,xlab=xlab,ylab=ylab) grid() dev.off() load(file='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='mytable.tab')

Summary of computational transaction
Raw Input	view raw input (R code)
Raw Output	view raw output of R engine
Computing time	1 seconds
R Server	Big Analytics Cloud Computing Center