0

我很好奇如何将其转换为 Rcpp 格式？由于内存问题，R 的 as.matrix 方法不起作用

sparse.cor4 <- function(x){
    n <- nrow(x)
    cMeans <- colMeans(x)
    covmat <- (as.matrix(crossprod(x)) - n*tcrossprod(cMeans))/(n-1)
    sdvec <- sqrt(diag(covmat)) 
    cormat <- covmat/tcrossprod(sdvec)
    list(cov=covmat,cor=cormat)
}

该函数来自此链接：Running cor() (or any variant) over a sparse matrix in R

附加信息：我能够使用 rcpp 创建一个 500kx500k 的稀疏矩阵，但我需要关联由于内存大小而无法使用 R 的稀疏矩阵，cor()这就是为什么我问我是否可以将上述函数转换为 Rcpp 模式得到稀疏矩阵的相关性

谢谢

4

0 回答 0