Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
31.1 31.8 32.5 34.4 35.5 35.5 36.6 37.1 37.9 38.1 39 41.5 41.8 41.9 44.6 46.1 46.4 47.2 47.7 49.2 49.3 49.3 49.5 50.1 51.9 52.6 53.2 53.5 53.7 53.7 53.9 54.1 54.8 55.4 55.9 56.8 58.4 59.3 60.3 60.5 60.8 61 61.1 61.3 61.4 61.5 63.9 63.9 64 64.1 64.5 64.5 65.9 66.8 68.7 69.2 69.6 70.2 70.6 70.7 70.7 71 72.1 73.7 77.4 79.7 91.6 93.6 94.3 97.3 101.7 103 103.1 104.6 107.2 107.7 108.3 108.8 113.1 113.8 113.8 116.5 116.9 117.6
Chart options
R Code
library(Hmisc) m <- mean(x) e <- median(x) bitmap(file='test1.png') op <- par(mfrow=c(2,1)) mydensity1 <- density(x,kernel='gaussian',na.rm=TRUE) plot(mydensity1,main='Density Plot - Gaussian Kernel',xlab='Median (0 -> full line) | Mean (0 -> dashed line)',ylab='density') abline(v=e,lty=1) abline(v=m,lty=5) grid() myseq <- seq(0.01, 0.99, 0.01) hd <- hdquantile(x, probs = myseq, se = TRUE, na.rm = FALSE, names = TRUE, weights=FALSE) plot(myseq,hd,col=2,main='Harrell-Davis Quantiles',xlab='quantiles',ylab='Median (0 -> full) | Mean (0 -> dashed)') abline(h=m,lty=5) abline(h=e,lty=1) grid() par(op) dev.off() load(file='createtable') a<-table.start() a<-table.row.start(a) a<-table.element(a,'Median versus Mean',2,TRUE) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'mean',header=TRUE) a<-table.element(a,mean(x)) a<-table.row.end(a) a<-table.row.start(a) a<-table.element(a,'median',header=TRUE) a<-table.element(a,median(x)) 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