Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
62.9 64.92 64 65.01 47 61.35 61 61.71 59.8 57.55 40 45.01 38 40.5 44.5 40.98 36.51 37.49 34.95 38.21 33.63 33.49 28 28 26.75 26.62 29.5 31 34 35 35.9 32.8 33 33.5 34.35 32.6 28.6 27.45 26.51 26.9 24 19.61 20 19.2 20.44 20.3 20.4 23.15 23.78 22 23.4 23.75 24.78 25.28 29.57 22.05 22.75 20.85 20.35 16.8 12.99 9.66 11 12.51 12.7 12.82 12.9 14.5 14.13 14 13.59 13.87
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