替换default默认解决方法（）ginv

Question

我想在mahanalobis函数中用ginv()替换solve()。

有什么办法可以强制R中的任何功能使用ginv()而不是solve()？

Answer 1

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

从?mahalanobis：

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