有人知道如何在R中的矩阵中随机删除或子集条目吗? 例如,我想在下面的代码中随机删除给出miss的条目。我跑的时候给了A [4,2]和A [3,3]。我真的想避免使用for循环,因为我正在使用的真实矩阵要大得多。
A=matrix(1:15, 5,3)
nmiss=2
obs=outer(1:nrow(A),1:ncol(A), FUN='paste', sep=',')
miss=sample(obs,nmiss, replace=TRUE)
不幸的是
A[miss]=NA
不起作用
答案 0 :(得分:3)
A[sample(nrow(A)*ncol(A), nmiss)] <- NA
您不必使用行索引和列索引随机地对矩阵进行子集化。