我正在处理遗传数据。我有一个矩阵，其中包含每个SNP的种群名称（第一列）和频率。
我正在尝试使用函数从此矩阵中估算NA。这些功能在个体的每个SNP处使用最常见的基因型。
但是，即使我得到100670 NA，我仍然无法以某种方式使其无法运行，就像在没有NA时一样，在尝试运行该函数时出现错误。我已经在具有相同格式的其他数据集中使用了此功能，并且可以正常工作，所以我看不到为什么这里不起作用。 谢谢！
sum（is.na（as.matrix（glSA）））＃100670
str（Pop.alleles.SA）
1] 153844 num [1:10，1：31400] 0.0357 0 0.3333 0.1071 0.25 ... -attr（*，“假名”）=列表2 .. $：chr [1:10]“ GN”“ LE”“ LO”“ ME” ... .. $：chr [1：31400]“ vchr1_79537”“ vchr1_674851”“ vchr1_675774”“ vchr1_841100” ...
Pop.alleles.SA.F <-apply（SA，2，function（x）replace（x，is.na（x），> as.numeric（names（that.max（table（x）））） ））） sum（is.na（Pop.alleles.SA.F））＃没有NAs
x [列表]中的错误<-值：替换长度为零