Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
58.1 57.9 57.3 55.9 55 55.9 56.6 57.3 56.2 57.7 56.8 57.9 58.3 58.2 56.9 57.1 56.7 54.2 54.2 52.1 51.5 51.8 53 52.4 52.41 52.36 52.94 52.34 51.84 51.42 50.85 50.66 51.53 51.59 52.32 51.98 51.17 50.57 49.84 50.12 49.08 48.57 47.22 46.78 46.04 45.05 44.42 44.09 44.46 44.34 43.04 42.87 42.32 42.49 41.94 41.6 41.42 41.12 41.28 40.21 39.69 39.16 38.8 38.44 37.02 36.75 35.95 36.29 36.35 36.07 36.6 36.5
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