Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
30.37 30.41 30.46 30.7 30.85 30.93 31 31.16 31.14 31.15 31.2 31.22 31.25 31.39 31.49 31.71 31.73 31.96 32.05 32.12 32.28 32.42 32.48 32.89 33.7 34.59 35.1 35.87 37.15 37.61 37.97 38.94 39.18 39.49 39.86 40.02 40.2 40.85 41.45 41.7 41.92 41.97 42.31 42.61 42.82 43.07 43.51 43.57 43.86 44.49 45.99 48.22 49.46 50.39 50.4 50.59 51.32 51.86 52.47 52.73 52.73 53.59 54.11 54.8 55.72 56.06 56.66 57.05 57.31 57.89 58.32 58.72 59.02 59.54 61.49 62.26 63.49 64.36 65.93 66.82 68.85 71.27 72.27 73.4
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
2 seconds
R Server
Big Analytics Cloud Computing Center
Click here to blog (archive) this computation