# パッケージがインストール済なら次の行は不要
# import Pkg; Pkg.add("StatsBase")
import Random  # 乱数を扱うためのモジュールをインポート
using StatsBase
n = 100000  # 試行回数
x = zeros(n)  # 感染していないことを0で表し，感染していることを1で表す．
y = zeros(n)  # 陰性を0で表し，陽性を1で表す．
Random.seed!(1)  # 乱数のシード値を固定
for i in 1:n
    x[i] = sample([0, 1], ProbabilityWeights([0.99, 0.01]))    
    if x[i] == 0  # 感染していない場合
        y[i] = sample([0, 1], ProbabilityWeights([0.98, 0.02]))
    else  # 感染している場合
        y[i] = sample([0, 1], ProbabilityWeights([0.01, 0.99]))
    end
end
# 確率を表示 (※モンテカルロ法で求めた値，理論値の順)
println(count((x .== 0) .& (y .== 1)) / count(y .== 1))
println(2/3)