如何在R中隨機生成缺失數據（MNAR）？

Question

我有Y是以概率0.1和N = 100二項分佈的。自變量x的正態分佈均值爲0和0.5。我想在Y上生成MNAR機制。

Answer 1

要生成MNAR數據，您可以使用數據生成過程，其中 缺失機制取決於未觀察到的數據。

# Generate the true data 
y1 <- rbinom(100, size=1, prob=0.1) 
# Generate the missing process. Depends on the "true" observed value 
r <- rbinom(length(y1), size=1, prob=c(.25, .1)[y1+1]) 
y <- y1 
y[r==1] <- NA