18/07/01 10:18:46.34 F3KIhVUK.net
>>380
バグをしてしてみ！
ド底辺シリツ医大卒はスクリプトすら読めないんじゃねぇの？
簡単にデバッグできた。
date=1:366
p=c(97/400,rep(1,365))
same.birth <- function(n,lwr=2,upr=1e6){
x=sample(date,n,replace=TRUE,prob=p)
di=max(table(x))
lwr<=di & di<=upr
}
birth <- function(n,lwr=2,upr=1e6,k=1e4){
mean(replicate(k,same.birth(n,lwr,upr)))
}
Vectorize(birth)(100,3:5)
100人中同じ誕生月日の人が3,4,5人以上いる確率も算出できるな。