R birthday 巧合的概率

说明

计算广义生日悖论问题的答案。 pbirthday 计算重合概率，qbirthday 计算至少具有指定重合概率所需的最小观测值数量。

用法

qbirthday(prob = 0.5, classes = 365, coincident = 2)
pbirthday(n, classes = 365, coincident = 2)

参数

classes	人们可以分为多少个不同的类别
prob	期望的重合概率
n	人的数量
coincident	属于同一类别的人数

细节

生日悖论是，极少数人(23 岁)足以有 50-50 的机会其中两个或更多人拥有相同的生日。此函数将计算推广到 0.5 以外的概率、2 以外的同时事件数以及 365 以外的类数。

使用的公式是 coincident > 2 的近似公式。该近似对于 prob 的中等值非常好，但对于非常小的概率则不太好。

值

qbirthday	至少需要 prob 且 k 或更多人在 classes 等概率标签中具有相同标签所需的最少人数。
pbirthday	指定巧合的概率。

例子

require(graphics)

## the standard version
qbirthday() # 23
## probability of > 2 people with the same birthday
pbirthday(23, coincident = 3)

## examples from Diaconis & Mosteller p. 858.
## 'coincidence' is that husband, wife, daughter all born on the 16th
qbirthday(classes = 30, coincident = 3) # approximately 18
qbirthday(coincident = 4)  # exact value 187
qbirthday(coincident = 10) # exact value 1181

## same 4-digit PIN number
qbirthday(classes = 10^4)

## 0.9 probability of three or more coincident birthdays
qbirthday(coincident = 3, prob = 0.9)

## Chance of 4 or more coincident birthdays in 150 people
pbirthday(150, coincident = 4)

## 100 or more coincident birthdays in 1000 people: very rare
pbirthday(1000, coincident = 100)

参考

Diaconis, P. and Mosteller F. (1989). Methods for studying coincidences. Journal of the American Statistical Association, 84, 853-861. doi:10.1080/01621459.1989.10478847.