#topicpath ** 一つのクラスに同じ誕生日の人がいる確率 [#z715b3fe] - 概要~ 一つのクラスに同じ誕生日の人が一組以上いる確率を求めてみる。~ ############################################ # # クラスに同じ誕生日の人がいる確率.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))