如何在R中生成Missing Not at random（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