Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
59.8 60.7 59.7 60.2 61.3 59.8 61.2 59.3 59.4 63.1 68 69.4 70.2 72.6 72.1 69.7 71.5 75.7 76 76.4 83.8 86.2 88.5 95.9 103.1 113.5 115.7 113.1 112.7 121.9 120.3 108.7 102.8 83.4 79.4 77.8 85.7 83.2 82 86.9 95.7 97.9 89.3 91.5 86.8 91 93.8 96.8 95.7 91.4 88.7 88.2 87.7 89.5 95.6 100.5 106.3 112 117.7 125 132.4 138.1 134.7 136.7 134.3 131.6 129.8 131.9 129.8 119.4 116.7 112.8 116 117.5 118.8 118.7 116.3 115.2 131.7 133.7 132.5 126.9 122.2 120.2
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