﻿rm(list = ls())  # ワークスペースをクリア
n <- 100000  # 試行回数
x <- integer(n)  # 感染していないことを0で表し，感染していることを1で表す．
y <- integer(n)  # 陰性を0で表し，陽性を1で表す．
set.seed(1)  # 乱数のシード値を固定
for (i in 1:n) {
  x[i] <- sample(c(0,1), 1, replace=TRUE, prob = c(99,1))
  if (x[i] == 0) {  # 感染していない場合
    y[i] <- sample(c(0,1), 1, replace=TRUE, prob = c(98,2))
  } else {  # 感染している場合
    y[i] <- sample(c(0,1), 1, replace=TRUE, prob = c(1,99))
  }
}
# 確率を表示
sum(x == 0 & y == 1) / sum(y == 1)  # モンテカルロ法
2/3  # 参考：理論値