#topicpath

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

トップ   編集 差分 添付 複製 名前変更 リロード   新規 単語検索 最終更新   ヘルプ   最終更新のRSS