r - 如果有机会获得财产，那么 3 人中有 1 人的机会有财产吗？（如何在 R 中计算？）

Question

我在一篇论文中找到了这句话：

在北欧血统的个体中，多达 8% 的男性和 0.5% 的女性经历了常见的红绿色盲。如果提交的稿件碰巧投给了三名北欧血统的男性审稿人，那么至少有 22% 的人是色盲的。

Wong, B.：观点：色盲，Nat Methods, 8(6), 441–441, doi:10.1038/nmeth.1618, 2011。

我怎么能在R中计算这个？提前致谢。

标签： rstatisticsbinomial-coefficients

Answer 1

另一种方法是凭经验而不是乔回答的分析方式：

n = 100000 #Arbitrary large number
at.least.1 = numeric() #Empty vector to store "1" when at least one was colorblind, "0" if not
for(j in 1:n){
  at.least.1[j] = sum(sample(0:1, size=3, prob=c(0.92,0.08), replace=TRUE))>=1}
sum(at.least.1)/n

对于这个问题，分析起来很容易，但有时这会更好。