Question

我想在矩阵的每个元素中设置NA，其中列中的值大于或等于给定向量的值。例如，我可以创建一个矩阵：

set.seed(1)
zz <- matrix(data = round(10L * runif(12)), nrow = 4, ncol = 3)

给出了zz：

     [,1] [,2] [,3]
[1,]    8    5    7
[2,]    6    5    1
[3,]    5   10    3
[4,]    9    1    9

和比较矢量（例如）：

xx <- round(10L * runif(4))

其中xx是：

[1] 6 3 8 2

如果我执行此操作：

apply(zz,2,function(x) x >= xx)

我明白了：

      [,1]  [,2]  [,3]
[1,]  TRUE FALSE  TRUE
[2,]  TRUE  TRUE FALSE
[3,] FALSE  TRUE FALSE
[4,]  TRUE FALSE  TRUE

我想要的是无处不在我有一个TRUE元素我想要一个NA而且每个地方我都有一个FALSE我得到zz矩阵中的数字（例如，手动......）：

NA  5  NA
NA  NA 1
5   NA 3
NA  1  NA

我可以拼凑一些“for”循环来做我想要的，但有没有基于矢量的方法来做到这一点？

感谢您的任何提示。

Answer 1

这是获得预期输出的一个选项。我们得到一个逻辑矩阵（zz >= xx)，使用NA^返回NA表示TRUE值，1表示FALSE，然后将其乘以原始矩阵＆＃39; zz＆＃39;以便NA保持不变因此，当1变为＆＃39; zz＆＃39;中的相应值时。

NA^(zz >= xx)*zz
#     [,1] [,2] [,3]
#[1,]   NA    5   NA
#[2,]   NA   NA    1
#[3,]    5   NA    3
#[4,]   NA    1   NA

或另一个选项是ifelse

ifelse(zz >= xx, NA, zz)

zz <- structure(c(8, 6, 5, 9, 5, 5, 10, 1, 7, 1, 3, 9), .Dim = c(4L, 3L))
xx <- c(6, 3, 8, 2)

Answer 2

你可以这样做：

zz[zz>=xx] <- NA

#     [,1] [,2] [,3]
#[1,]   NA    5   NA
#[2,]   NA   NA    1
#[3,]    5   NA    3
#[4,]   NA    1   NA