- 追加された行はこの色です。
- 削除された行はこの色です。
#topicpath
** 一つのクラスに同じ誕生日の人がいる確率 [#z715b3fe]
- 概要~
一つのクラスに同じ誕生日の人が一組以上いる確率を求めてみる。~
一つのクラスに同じ誕生日の人が一組以上いる確率を[[R>http://cran.r-project.org/]]で数値的に求めてみる。~
~
- プログラム~
~
############################################
#
# クラスに同じ誕生日の人がいる確率.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);
&ref(https://robo.mydns.jp/Lecture/DATA/DAMPEN/SameBirthDay/Result.png,,335x335);