归因于中值加上抖动

Question

我想在每个单元格中有效地使用略有不同的值来计算缺失值。

例如：

df <- data_frame(x = rnorm(100), y = rnorm(100))
df[1:5,1] <- NA
df[1:5, 2] <- NA

df %<>% mutate_all(funs(ifelse(is.na(.), jitter(median(., na.rm = TRUE)), .))) 

然而，这在所有单元格中都被认为是相同的数字。如何向每个单元格添加不同的噪波？当然，我可以使用循环来完成这项工作，但我的数据帧很大，我希望高效地完成这项工作

Answer 1

我们可以在n()中使用rep

library(dplyr)
library(magrittr)
df %<>%
   mutate_all(list(~ case_when(is.na(.) ~ jitter(rep(median(., na.rm = TRUE), n())),
         TRUE ~ .)))