Wessa.net

Send output to:


Data:
76 75 74 72 70 69 70 72 73 73 74 76 74 67 66 58 55 58 64 68 66 76 75 88 85 83 77 66 65 65 63 62 57 68 69 79 74 76 82 75 75 76 78 77 67 74 68 87 76 88 95 96 96 105 108 113 101 107 102 116 105 121 134 140 131 141 131 128 123 129 125 144 135 141 156 159 146 154 145 133 126 127 122 148


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