Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
228768.00 227916.00 227052.00 225264.00 242952.00 242016.00 228768.00 219960.00 220812.00 220812.00 221760.00 223464.00 226116.00 226116.00 224412.00 219960.00 242952.00 246456.00 241164.00 228768.00 234072.00 226116.00 229704.00 231420.00 233208.00 228768.00 229704.00 223464.00 242952.00 249108.00 243816.00 234072.00 244668.00 233208.00 243816.00 242952.00 245604.00 235860.00 246456.00 245604.00 261504.00 257916.00 243816.00 236712.00 246456.00 233208.00 242952.00 244668.00 248256.00 240312.00 244668.00 247320.00 257064.00 249108.00 238512.00 227052.00 237660.00 208500.00 222612.00 230556.00 238512.00 227052.00 227052.00 227052.00 233208.00 224412.00 212868.00 203208.00 210216.00 182856.00 199620.00 209364.00 211152.00 201408.00 202260.00 199620.00 208500.00 202260.00 189960.00 181068.00 196104.00 163452.00 184656.00 194316.00 194316.00 182856.00 172260.00 171408.00 181068.00 172260.00 155508.00 143964.00 156360.00 127212.00 153708.00 167808.00 172260.00 162516.00 150204.00 159012.00 162516.00 159864.00 133356.00 121056.00 129852.00 103356.00 130716.00 140460.00 148404.00 135156.00 122760.00 129852.00 133356.00 126348.00 99852.00 88308.00 98904.00 69756.00 101556.00 121056.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
1 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation