#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)) - 結果~ 横軸:クラスの人数 / 縦軸:同じ誕生日の人の組がある確率~ ~ &ref(http://robo.mydns.jp/Lecture/DATA/DAMPEN/SameBirthDay/Result.png,,335x335);