デフォルトのsolve（）をginvに置き換えます。

Question

mahanalobis関数のでsolve()を置き換えます。

solve()の代わりにginv()を強制的に使用する方法はありますか？

Answer 1

mahalanobisの機能はかなりシンプルです。理由だけではなく、?mahalanobisからすなわち

、適切な置換で

mahalanobis_ginv <- function (x, center, cov, 
      inverted = FALSE, ...) { 
    x <- if (is.vector(x)) 
     matrix(x, ncol = length(x)) 
    else as.matrix(x) 
    if (!identical(center, FALSE)) 
     x <- sweep(x, 2L, center) 
    if (!inverted) 
     cov <- MASS::ginv(cov, ...) 
    setNames(rowSums(x %*% cov * x), rownames(x)) 
} 

を独自に定義しない：

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