Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
37.3 39.5 40.6 41.4 41.3 43.5 44 44.9 46.4 47.4 48.7 49.7 51.1 53.2 56.2 58.1 60.6 64.1 67.4 68 70.9 72.8 74.9 76.1 77 78.1 80 79.7 82.7 84.3 83.5 85.9 87 88.6 90.6 91.3 91.6 93.2 95 95.2 97.4 98.6 99.6 100.6 101.3 102.8 103.2 103 105.4 104.7 105.2 105.2 102.8 100.3 99.8 99.4 100.6 100.2 100.4 98.8 96.9 96.3 96.1 93.5 92.1 91.7 87.9 86.4 84.9 81.7 82.6 83.1
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