Send output to:
Browser Blue - Charts White
Browser Black/White
CSV
Data:
1844.33 1845.39 1861.84 1862.85 1869.46 1870.8 1870.8 1871.52 1875.52 1880.38 1885.05 1886.42 1886.42 1891.65 1903.11 1905.29 1904.26 1905.37 1905.37 1905.12 1908.62 1915.08 1916.36 1916.68 1916.24 1922.05 1922.63 1922.47 1920.64 1920.66 1920.66 1921.19 1921.44 1921.73 1921.81 1921.81 1921.81 1921.48 1917.07 1912.64 1901.15 1898.12 1900.02 1900.02 1900.82 1901.9 1902.19 1901.84 1903.73 1889.7 1891.27 1894.48 1894.27 1893.98 1893.98 1895.62 1901.72 1905.4 1898.14 1898.09 1855.87 1868.53 1865.71 1872.59 1875.95 1875.95 1875.95 1878.08 1878.26 1876.39 1876.77 1876.88
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