![[PukiWiki]](image/ROBO192.png "[PukiWiki]")
############################################
#
# クラスに同じ誕生日の人がいる確率.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))
