没有 NA 的样本

Question

这个问题与stackoverflow下一页中提出的问题有关：Ignoring values or NAs in the sample function

问题：我有一个矩阵 K（比如 1528）行和许多（比如 1000）列。许多条目中都有 NA。从每一列中，我想按以下方式采样 40 个元素：如果列中此类元素的总数不少于 40，我从非 NA 元素中随机抽样。如果它们小于 40（比如=k），则将是所有 k 个元素和 40-k 个 NA。我试过以下代码：

mysample <- function(x){
  if(sum(is.na(x))>1488){
    sum(x[!is.na(x)])
    return(c(x[!is.na(x)],rep(NA,40-sum(x[!is.na(x)]))))
  }
  return(sample(x[!is.na(x)],40))
}

J=apply(K, 2, mysample)

在第四行，它在重复中显示无效的“次”参数。任何人都可以使代码可行（我想包括 NA，因为我想从中生成一个 40X1000 矩阵）

Answer 1

在您的rep 函数中，时间参数应该是40 - length_of_nona，它应该是40 - sum(!is.na(x))。没有测试，但我认为这会解决问题：

mysample <- function(x){
  if(sum(is.na(x))>1488){
    # sum(x[!is.na(x)])
    return(c(x[!is.na(x)],rep(NA,40-sum(!is.na(x)))))
                                        ^        ^
  }
  return(sample(x[!is.na(x)],40))
}

J=apply(K, 2, mysample)