[pymc3][python]pymc3安装后测试代码2

测试环境:

pymc3==3.11.2

代码:

复制代码
import numpy as np
import pymc3 as pm
import matplotlib.pyplot as plt


if __name__ == '__main__':

    # 生成随机数据
    np.random.seed(123)
    x = np.linspace(0, 1, 100)
    y = 0.5 * x + np.random.normal(0, 0.1, size=100)

    # 定义概率模型
    with pm.Model() as model:
        # 定义先验分布
        alpha = pm.Normal('alpha', mu=0, sd=1)
        beta = pm.Normal('beta', mu=0, sd=1)
        sigma = pm.HalfNormal('sigma', sd=1)

        # 定义线性关系
        mu = alpha + beta * x

        # 定义似然函数
        likelihood = pm.Normal('y', mu=mu, sd=sigma, observed=y)

        # 进行贝叶斯推断
        trace = pm.sample(1000, tune=1000)

    # 可视化结果
    pm.traceplot(trace)
    plt.show()

结果:

相关推荐
木头左6 分钟前
Python实现ETF网格自动化交易集成动量阈值判断
开发语言·自动化
CodeCraft Studio10 分钟前
全球知名的Java Web开发平台Vaadin上线慧都网
java·开发语言·前端·vaadin·java开发框架·java全栈开发·java ui 框架
静水流深-刘申11 分钟前
算法继续刷起-2025年09月26日
开发语言·c++·算法
我是华为OD~HR~栗栗呀21 分钟前
前端面经-高级开发(华为od)
java·前端·后端·python·华为od·华为·面试
平平无奇。。。31 分钟前
C++11之异常
开发语言·c++·visual studio
木头左34 分钟前
跨周期共振效应在ETF网格参数适配中的应用技巧
开发语言·python·算法
爱蹦跶的精灵37 分钟前
降级版本Pillow解决freetypefont has no attribute getsize问题
python·pillow
almighty271 小时前
C# WPF实现ComboBox实时搜索与数据绑定
开发语言·c#·wpf·combobox
菜鸟小九1 小时前
SSM(MybatisPlus)
java·开发语言·spring boot·后端
一人の梅雨1 小时前
亚马逊 MWS 关键字 API 实战:关键字搜索商品列表接口深度解析与优化方案
python·spring