Home » date » 2009 » Oct » 27 »

WS 3 Prediction future

*The author of this computation has been verified*

R Software Module: rwasp_meanversusmedian.wasp (opens new window with default values)

Title produced by software: Mean versus Median

Date of computation: Tue, 27 Oct 2009 09:58:13 -0600

Cite this page as follows:

Statistical Computations at FreeStatistics.org, Office for Research Development and Education, URL http://www.freestatistics.org/blog/date/2009/Oct/27/t12566591582urc501j1qtir90.htm/, Retrieved Fri, 24 May 2013 08:08:54 +0000

Original text written by user:

IsPrivate?

No (this computation is public)

User-defined keywords:

System-generated keywords (parent):

(pk = 0)

Estimated Impact

Dataseries X:

» Textfile « » CSV « » Stem and Leaf « » Histogram « » Kernel Density « » Harrell-Davis Quantiles « » Central Tendency « » Variability «

Output produced by software:

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	'Gwilym Jenkins' @ 72.249.127.135

Median versus Mean
mean	77.5381666666667
median	76.775

Charts produced by software:

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566591582urc501j1qtir90/1gnmw1256659090.png (opens in new window)

http://www.freestatistics.org/blog/date/2009/Oct/27/t12566591582urc501j1qtir90/1gnmw1256659090.ps (opens in new window)

Click here to open pdf file.

Parameters (Session):

Parameters (R input):

R code (references can be found in the software module):

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