r - prcomp() 缺失值插补的最佳选择

Question

我有近 1,000,000 个基因座（行）的大约 200 个个体基因组（列）的基因型数据集。由于测序数据不佳，大多数行包含 1-2 个缺失的基因型。

如果我使用

df_new = na.omit(df)

我的新数据框只包含几千行，这导致数据丢失比通过每行输入一两个缺失值得到的数据要大得多。我一直在网上寻找如何使用与 na.option 和 prcomp() 关联的插补选项，但找不到示例。我想从最简单的方法开始，例如用中值或类似的东西替换 NA。

有人可以指导我在 prcomp 的上下文中如何执行此操作的示例吗？

Answer 1

现在我理解了您的问题，请参见下面的示例：

library(plyr)

   
     ddply(df_new, ~ my_groups, transform,
         missing value column = ifelse(is.na(missing value column), 
                      median(missing value column, na.rm = TRUE), 
                                 missing value column))

  #missing value column is the column that consist the missing value
  #my_groups could be the first column of df_new

我希望这行得通。