______________________________________________
R-h...@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
Something like this ( I think it's correct ):
--%<---
DBhat <- function(X1,X2){
# define means
mX1 <- mean(X1)
mX2 <- mean(X2)
# define difference of means
mDiff <- mX1 - mX2
# define cov
cvX1 <- cov(X1)
cvX2 <- cov(X2)
# define halfsum of cv's
p <- (cvX1+cvX2)/2
# the equation
0.125 * t(mDiff) * p^(-1) * mDiff + 0.5 * log10( det(p) /
sqrt( det(cvX1) * det(cvX2) ))
--%<---
A related thread (in grass-stats):
http://www.mail-archive.com/grass...@lists.osgeo.org/msg00183.html
Regards, Nikos
Notice: This e-mail message, together with any attachme...{{dropped:10}}
http://en.wikipedia.org/wiki/Bhattacharyya_distance
There is also the Hellinger Distance and the Rao distance. For the Rao
distance, see:
http://www.scholarpedia.org/article/Fisher-Rao_metric
Jude
___________________________________________
Jude Ryan
Director, Client Analytical Services
Strategy & Business Development
UBS Financial Services Inc.
1200 Harbor Boulevard, 4th Floor
Weehawken, NJ 07086-6791
Tel. 201-352-1935
Fax 201-272-2914
Email: jude...@ubs.com