β分布
P
r
(
λ
)
=
Γ
[
α
+
β
]
Γ
[
α
]
Γ
[
β
]
λ
α
−
1
(
1
−
λ
)
β
−
1
Pr(\lambda) = \frac{\Gamma[\alpha+\beta]}{\Gamma[\alpha]\Gamma[\beta]}\lambda^{\alpha-1}(1-\lambda)^{\beta-1}
Pr(λ)=Γ[α]Γ[β]Γ[α+β]λα−1(1−λ)β−1
α
>
0
,
β
>
0
,
0
<
=
λ
<
=
1
\alpha > 0, \beta > 0, 0 <= \lambda <= 1
α>0,β>0,0<=λ<=1
import numpy as np
import math
import matplotlib.pyplot as plt
alpha = 0.1
beta = 0.1
alpha_ = math.gamma(alpha)
beta_ = math.gamma(beta)
ab_ = math.gamma(alpha+beta)
print(alpha_,beta_,ab_)
# 模拟1000个点
pr = []
for i in range(1,1000):
lambda_ = i/1000
pr.append(ab_/(alpha_+beta_) * (math.pow(lambda_,(alpha-1))) * (math.pow((1-lambda_),(beta-1))))
x = [i for i in range(1,1000)]
plt.plot(x, pr)
plt.show
<function matplotlib.pyplot.show(*args, **kw)>
alpha = 2
beta = 2
alpha_ = math.gamma(alpha)
beta_ = math.gamma(beta)
ab_ = math.gamma(alpha+beta)
print(alpha_,beta_,ab_)
1.0 1.0 6.0
# 模拟1000个点
pr = []
for i in range(1,1000):
lambda_ = i/1000
pr.append(ab_/(alpha_+beta_) * (math.pow(lambda_,(alpha-1))) * (math.pow((1-lambda_),(beta-1))))
x = [i for i in range(1,1000)]
plt.plot(x, pr)
plt.show
<function matplotlib.pyplot.show(*args, **kw)>