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。非經特殊聲明，原始代碼版權歸原作者所有，本譯文未經允許或授權，請勿轉載或複製。