矩阵中的R替代NA

Question

RI中的数据帧具有一些丢失的值，因此read.table()函数使用NA而不是空白单元格。

我这样写：

a <- sample(1000:50000000, size=120, replace=TRUE)
values <- matrix(a, nrow=6, ncol=20)
mtx <- cbind.data.frame(values, c(rep(NA),6))
mtx <- apply(mtx, 2, function(x){
    if (x==NA) sample(100:500, replace=TRUE, size=nrow(mtx)) else (x)})

但是我有这个错误：

Error in if (x == NA) sample(100:500, replace = TRUE, size = nrow(mtx)) else (x) : 
  missing value where TRUE/FALSE needed
In addition: Warning message:
In if (x == NA) sample(100:500, replace = TRUE, size = nrow(mtx)) else (x) :
  the condition has length > 1 and only the first element will be used

有任何想法吗？

最佳里卡多

Answer 1

您不能使用比较运算符来测试NA ，原因是该值是NA或缺失。 is.na()是用于识别NA形式的缺失的适当函数。

这是在矩阵values中替换NA的示例。 这里的关键是以向量化的方式工作，只需确定哪些元素是NA然后使用索引将所有NA替换为所需的值。

> set.seed(2)
> values <- matrix(sample(1000:50000000, size=120, replace=TRUE),
+                  nrow=6, ncol=20)
> ## add some NA to simulate
> values[sample(120, 20)] <- NA
> 
> ## how many NA
> (tot <- sum(is.na(values)))
[1] 20
> 
> ## replace the NA
> values[is.na(values)] <- sample(100:500, tot, replace=TRUE)
> 
> ## now how many NA
> (sum(is.na(values)))
[1] 0