一つのクラスに同じ誕生日の人がいる確率 †
- 概要
一つのクラスに同じ誕生日の人が一組以上いる確率を求めてみる。
############################################
#
# クラスに同じ誕生日の人がいる確率.R
#
############################################
ProbabilityOfTheSameBirthday <- function(n, NumberOfTheDaysInAYear=365){
ProbabilityNG <- NULL
ProbabilityNG[1] <- 1
for( i in 2:max(n) ){
ProbabilityNG[i] <- ProbabilityNG[i-1] * (NumberOfTheDaysInAYear - i + 1) / NumberOfTheDaysInAYear
}
return (1 - ProbabilityNG)
}
n <- 1:45
plot(n, ProbabilityOfTheSameBirthday(n), type='l', ylim=c(0,1))