Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
14724.00 14404.00 14058.00 13427.00 19946.00 19631.00 14724.00 11462.00 11778.00 11778.00 12093.00 12760.00 13742.00 13427.00 11462.00 11778.00 20929.00 22893.00 17666.00 14724.00 15387.00 15707.00 17351.00 18964.00 19315.00 16022.00 16369.00 12093.00 24222.00 27800.00 19631.00 17004.00 18649.00 20613.00 23555.00 27164.00 27164.00 24853.00 23871.00 17982.00 27800.00 32391.00 28462.00 24222.00 24853.00 27164.00 30426.00 34355.00 31724.00 30111.00 30111.00 24853.00 32391.00 37297.00 33373.00 29129.00 30426.00 35653.00 37964.00 41222.00 38595.00 34355.00 33373.00 25520.00 30742.00 36315.00 30111.00 26502.00 30111.00 33689.00 35653.00 40906.00 38280.00 31724.00 32391.00 26186.00 31409.00 36000.00 30742.00 27164.00 30426.00 34355.00 33689.00 41542.00 40244.00 35017.00 35333.00 28462.00 32706.00 39262.00 34355.00 31409.00 36315.00 39262.00 36982.00 47431.00 44835.00 38946.00 37297.00 29764.00 34035.00 37964.00 33053.00 33053.00 38595.00 41542.00 39924.00 51355.00 48413.00 42871.00 40560.00 32391.00 35333.00 40560.00 36631.00 35653.00 40244.00 44168.00 39924.00 50057.00
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
par3 <- '0.01' par2 <- '0.99' par1 <- '0.01' 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')
Compute
Summary of computational transaction
Raw Input
view raw input (R code)
Raw Output
view raw output of R engine
Computing time
0 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation