β分布

β分布

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)>

在这里插入图片描述

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值