Wessa.net

Send output to:


Data:
123.46 123.24 123.86 124.28 124.78 125.19 125.46 127.60 127.80 126.63 127.06 126.77 127.05 128.23 128.60 128.97 129.34 129.71 130.08 130.45 128.82 132.19 131.56 131.93 132.30 130.67 133.05 132.42 133.79 134.16 134.53 136.90 135.27 136.64 136.01 136.38 136.75 138.12 137.50 137.87 138.24 138.61 138.98 140.35 139.72 143.09 140.46 141.83 143.20 140.57 141.95 141.32 142.69 143.06 144.43 143.80 144.17 144.54 146.91 145.28 144.65 145.02 144.40 146.77 146.14 147.51 148.88 148.25 147.62 150.99 148.36 149.73 150.10

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')

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