- 392 名前：１３２人目の素数さん mailto:sage [2020/11/24(火) 22:00:09.38 ID:PhesMNbr.net]
- >>378
ベイズの公式なんぞ使わずに、数を数えて計算するシミュレーションをしてみた。
sim <- function(p){ # p:Aの袋に赤が2個含まれる確率
f <- function(){# 1:赤 0:黒
A=c(1,rbinom(1,1,p)) # A:袋の中の玉
a=sample(A,2,replace=TRUE) # a:取り出した玉
c(sum(A),sum(a)) # 袋の中の赤の数、取り出した赤の数
}
k=1e4 # 試行回数
re=t(replicate(k,f()))
a2=re[re[,2]==2,] # a==2の場合の数
nrow(a2[a2[,1]==2,])/nrow(a2) # a==2&A==2の場合の数/a==2の場合の数
}
あとはpに一様乱数を与えて結果をヒ
