R mahalanobis 马哈拉诺比斯距离

R语言 mahalanobis 位于 stats 包(package)。

说明

返回 x 中所有行和向量 \mu = center 相对于 \Sigma = cov 的马哈拉诺比斯距离平方。 (对于向量 x )定义为

用法

mahalanobis(x, center, cov, inverted = FALSE, ...)

参数

`x`	数据向量或矩阵，例如 `p` 列。
`center`	分布的平均向量或长度为 `p` 或可回收到该长度的第二数据向量。如果设置为 `FALSE` ，则跳过居中步骤。
`cov`	分布的协方差矩阵 (`p \times p`)。
`inverted`	合乎逻辑的。如果 `TRUE` ， `cov` 应该包含协方差矩阵的逆矩阵。
`...`	传递给`solve`以计算协方差矩阵的逆(如果`inverted`为假)。

例子

require(graphics)

ma <- cbind(1:6, 1:3)
(S <-  var(ma))
mahalanobis(c(0, 0), 1:2, S)

x <- matrix(rnorm(100*3), ncol = 3)
stopifnot(mahalanobis(x, 0, diag(ncol(x))) == rowSums(x*x))
        ##- Here, D^2 = usual squared Euclidean distances

Sx <- cov(x)
D2 <- mahalanobis(x, colMeans(x), Sx)
plot(density(D2, bw = 0.5),
     main="Squared Mahalanobis distances, n=100, p=3") ; rug(D2)
qqplot(qchisq(ppoints(100), df = 3), D2,
       main = expression("Q-Q plot of Mahalanobis" * ~D^2 *
                         " vs. quantiles of" * ~ chi[3]^2))
abline(0, 1, col = 'gray')

也可以看看

cov , var

相关用法

注：本文由纯净天空筛选整理自R-devel大神的英文原创作品 Mahalanobis Distance。非经特殊声明，原始代码版权归原作者所有，本译文未经允许或授权，请勿转载或复制。