Wessa.net

Send output to:


Data:
92.76 93.12 93.6 93.24 93.4 93.32 93.13 93.19 93.84 94.01 93.78 93.47 93.6 92.85 92.91 92.29 92.5 93.1 92.86 93.19 93.73 93.88 93.85 93.45 93.43 93.59 95.28 94.95 94.49 94.45 94.35 95.52 96.89 97.54 97.65 97.35 98.2 99.46 100.35 99.72 99.69 99.62 99.77 100.19 100.82 100.36 101.08 100.73 101.51 102.12 102.88 103.47 103.53 103.67 103.68 103.76 103.67 103.01 103.39 103.43 103.4 104.8


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	0 seconds
R Server	Big Analytics Cloud Computing Center