独家 | ​PyMC3 介绍：用于概率编程的Python包

作者：Tung T. Nguye

翻译：王雨桐

校对：廖倩颖

本文约1900字，建议阅读8分钟

本文为你介绍PyMC3原理，并结合一个实际案例教你如何使用包实现计算。

介绍

我们经常从天气预报中听到：明天的降水率是80％。这意味着什么？我们很难直白地解释这种说法，尤其是从概率学派的角度：无限次（或非多次）地重复下雨/不下雨实验是不现实的。

 

贝叶斯方法可以解释这种说法。以下句子摘自《为黑客设计的概率规划与贝叶斯方法》一书，它完美地总结了贝叶斯学派的关键思想之一。

 

贝叶斯世界观将概率解释为事件可信度的量度，即我们对事件发生有多少信心。

 

这意味着在贝叶斯方法中，我们永远不能绝对确定自己的“信念”，但可以肯定表达我们对于相关事件发生有多少信心。此外随着收集到更多数据，我们可以对自己的信念更加信心。

 

作为一名科学家，我被训练着去相信数据，并且对所有事物都很谨慎。所以我认为贝叶斯推理是相当直观的。

 

但是使用贝叶斯推断在计算和概念上通常具有挑战性。完成工作经常需要大量耗时而复杂的数学计算。即使作为数学家，我有时也觉得这些计算很乏味；特别是要快速了解待解决的问题时。

 

幸运的是我的导师AustinRochford最近向我介绍了一个名为PyMC3的程序包，它使我们能够进行数值贝叶斯推理。本文将通过一个具体示例快速介绍PyMC3。

 

一个具体的例子

假设我们有一枚硬币，我们将其翻转三遍，结果是：

[0，1，1]

 

其中0表示硬币背面向上，1表示人头向上。我们有信心说这是一个公平的硬币吗？换句话说，如果让θ为人头向上的概率，那么证据是否足以支持θ= 0.5的说法？

 

由于除了上述实验的结果外，我们对硬币一无所知，因此很难确定地说什么。从概率学派的角度来看，θ的点估计为：