如何将稀疏矩阵中的替换值设置为NA而不是0？

我想使用稀疏矩阵进行分析。稀疏矩阵中的每个单元格包含集合{0,1，NA}中的一个值。这里的NA代表缺失值。如何将稀疏矩阵中的替换值设置为NA而不是0？

例如，我可以使用下面的代码来创建一个稀疏矩阵：

library(Matrix); 
    toy <- Matrix(c(0,1,NA), nrow=3, ncol=3, sparse = TRUE)

，我得到下面的输出：

. . . 
1 1 1 
NA NA NA

是否有可能使用矩阵函数来创建数据点“遗漏”矩阵的函数是NA值而不是0？是否还有其他功能可以使用？

因此，从我的例子，我的期望输出为：

0 0 0 
1 1 1 
. . .

我已经探索了帮助页面和网站，但我还没有找到答案。我相信我需要设置'replValue'。

来源

2014-01-30 joel38237

有两个单独的问题，其实。第一个是如何显示零。这是很容易通过寻找被调度之后所使用的确切的方法来解决：

Matrix::printSpMatrix(toy, zero.print="0") 

[1,] 0 0 0 
[2,] 1 1 1 
[3,] NA NA NA

第二个问题是，是否NA输出可以与一些其它字符被抑制。那么，这是不可能的：没有适合的参数。

然而，修改源始终是一个选项。注意：这是一种黑客攻击，可能会导致不可预料的后果！

toy_print <- function (x, digits = NULL, maxp = getOption("max.print"), cld = getClassDef(class(x)), 
         zero.print = ".", col.names, note.dropping.colnames = TRUE, 
         col.trailer = "", align = c("fancy", "right")) 
{ 
    stopifnot(extends(cld, "sparseMatrix")) 
    x.orig <- x 
    cx <- formatSpMatrix(x, digits = digits, maxp = maxp, cld = cld, 
         zero.print = zero.print, col.names = col.names, note.dropping.colnames = note.dropping.colnames, 
         align = align) 
    if (col.trailer != "") 
     cx <- cbind(cx, col.trailer, deparse.level = 0) 
    # here's the NA hack 
    cx[cx=="NA"] <- "." 
    print(cx, quote = FALSE, right = TRUE, max = maxp) 
    invisible(x.orig) 
} 

toy_print(toy, zero.print="0") 

[1,] 0 0 0 
[2,] 1 1 1 
[3,] . . .

来源

2014-01-31 08:46:02 tonytonov

感谢您的帮助。我跑这个代码，我收到了不同的输出给你：[1，] 0 0 0 [2，] 1 1 1 [3，] 1 1 1 – joel38237

我刚刚仔细检查了代码，似乎很好。 – tonytonov

我已经清理了我的工作区，现在我得到和你一样的答案。再次感谢您的帮助和有趣的解决方案。 – joel38237

如何将稀疏矩阵中的替换值设置为NA而不是0？

回答

相关问题