期中考試 9.27
生日概率
> pBirthday <- function(k)
+ 1 - exp(lfactorial(365)-lfactorial(365-k)-k*log(365))
> pBirthday(50)
The capture and recapture problem
極大似然估計(jì)
> Pn <-function(n) {
+ tmp <- choose (50,3)*choose(n-50,47)
+ tmp/choose(n,50)
+}
> plot(n, Pn(n),type='l')
Combinatorics
Permutation
Binomial identity #$\Sigma$(n,k)＝2n
-
n個(gè)球放入n個(gè)盒子，剛好有一個(gè)盒子是空的概率==$\frac{n(n-1)}{(2n-1,n) }$==
>prob <- function(n){n*(n-1)/choose(2*n-1,n)} >prob(10)
Conditional Probability
-
P(A|B) 給定事件B發(fā)生的情況下，A發(fā)生的概率
Q：家裡倆孩子惰许，至少一個(gè)是男孩。那麼兩個(gè)都是男孩的概率隆圆？
A：P(X=2|X$\ge$1)=$\frac{1/4}{1-1/4}$
Bayes' theorem
P(A1|B)=$\frac{P(A1B)}{P(B)}$=$\frac{p(P(B|A1)P(A1))}{P(B|A1)P(A1)+P(B|A2)P(A2)}$
-
The law of total probability
對(duì)于獨(dú)立事件A1和A2使得S=A1$\cup$A2，P(B)=P(B|A1)P(A1)+P(B|A2)P(A2)
-
先驗(yàn)分析 因果關(guān)系
Q：五個(gè)人輪流拿球，四白一紅。第二個(gè)人拿到紅球的概率缓窜？
A：P(B1)P(A2|B1)=4/5*1/4=1/5
無后效性
Q：郵件分為三類，A1谍咆，A2禾锤，A3，互相排斥摹察。B＝郵件包含“free”
P(A1)=0.7 P(A2)=0.2 P(A3)=0.1
P(B|A1)=0.9 P(B|A2)= 0.1 P(B|A3)= 0.1
那么含有“free”的A1郵件可能性多大恩掷？
A：by Bayes，0.955
Example 2.13 誤診
