从贝叶斯方法谈到贝叶斯网络

转自zdy0_2004的专栏

1 贝叶斯方法

  长久以来，人们对一件事情发生或不发生的概率θ$\theta$，只有固定的0和1，即要么发生，要么不发生，从来不会去考虑某件事情发生的概率有多大，不发生的概率又是多大。而且概率θ$\theta$虽然未知，但最起码是一个确定的值。比如如果问那时的人们一个问题：“有一个袋子，里面装着若干个白球和黑球，请问从袋子中取得白球的概率θ$\theta$是多少？”他们会想都不用想，会立马告诉你，取出白球的概率θ$\theta$就是1/2，要么取到白球，要么取不到白球，即θ$\theta$只能有一个值，不是1/2，就是0，而且不论你取了多少次，取得白球的概率θ$\theta$始终都是1/2，即不随观察结果X$\mathbf{X}$的变化而变化。
  这种频率派的观点长期统治着人们的观念，直到后来一个名叫Thomas Bayes的人物出现。

1.1 贝叶斯方法的提出

  托马斯$$贝叶斯（Thomas Bayes, 1702-1763）在世时，并不为当时的人们所熟知，很少发表论文或出版著作，与当时学术界的人沟通交流也很少，用现在的话来说，贝叶斯就是活生生一民间学术“屌丝”，可这个“屌丝”最终发表了一篇名为“An essay towards solving a problem in the doctrine of chances”，翻译过来则是：机遇理论中一个问题的解。你可能觉得我要说：这篇论文的发表随机产生轰动效应，从而奠定贝叶斯在学术史上的地位。

  事实上，上篇论文发表后，在当时并未产生多少影响，在20世纪后，这篇论文才逐渐被人们所重视。对此，与梵高何其类似，画的画生前一文不值，死后价值连城。
  回到上面的例子：“有一个袋子，里面装着若干个白球和黑球，请问从袋子中取得白球的概率θ$\theta$是多少？”贝叶斯认为取得白球的概率是个不确定的值，因为其中含有机遇的成分。比如，一个朋友创业，你明明知道创业的结果就两种，即要么成功要么失败，但你依然会忍不住去估计他创业成功的几率有多大？你如果对他为人比较了解，而且有方法、思路清晰、有毅力、且能团结周围的人，你会不由自主的估计他创业成功的几率可能在80%以上。这种不同于最开始的“非黑即白非0即1”的思考方式，便是贝叶斯式的思考方式。
  继续深入讲解贝叶斯方法之前，先简单总结下频率派与贝叶斯派各自不同的思考方式：

• 频率派把需要推断的参数θ$\theta$看做是固定的未知常数，即概率虽然是未知的，但最起码是确定的一个值，同时，样本X$\mathbf{X}$ 是随机的，所以频率派重点研究样本空间，大部分的概率计算都是针对样本X$\mathbf{X}$ 的分布；
• 而贝叶斯派的观点则截然相反，他们认为参数θ$\theta$是随机变量，而样本X$\mathbf{X}$是固定的，由于样本是固定的，所以他们重点研究的是参数θ$\theta$的分布。
    相对来说，频率派的观点容易理解，所以下文重点阐述贝叶斯派的观点。
    贝叶斯派既然把看做是一个随机变量，所以要计算θ$\theta$的分布，便得事先知道θ$\theta$的无条件分布，即在有样本之前（或观察到X$\mathbf{X}$之前），θ$\theta$有着怎样的分布呢？
    比如往台球桌上扔一个球，这个球落会落在何处呢？如果是不偏不倚的把球抛出去，那么此球落在台球桌上的任一位置都有着相同的机会，即球落在台球桌上某一位置的概率θ$\theta$服从均匀分布。这种在实验之前定下的属于基本前提性质的分布称为先验分布，或的θ$\theta$无条件分布。
    至此，贝叶斯及贝叶斯派提出了一个思考问题的固定模式：
• 先验分布π(θ)+$\pi (\theta )+$样本信息X⇒$\mathbf{X}\Rightarrow$ 后验分布π(θ|x)$\pi (\theta |x)$

  上述思考模式意味着，新观察到的样本信息将修正人们以前对事物的认知。换言之，在得到新的样本信息之前，人们对的认知是先验分布π(θ)$\pi (\theta )$，在得到新的样本信息后X$\mathbf{X}$，人们对的认知为π(θ|x)$\pi (\theta |x)$。
  其中，先验信息一般来源于经验跟历史资料。比如林丹跟某选手对决，解说一般会根据林丹历次比赛的成绩对此次比赛的胜负做个大致的判断，再比如，某工厂每天都要对产品进行质检，以评估产品的不合格率θ$\theta$，经过一段时间后便会积累大量的历史资料，这些历史资料便是先验知识，有了这些先验知识，便在决定对一个产品是否需要每天质检时便有了依据，如果以往的历史资料显示，某产品的不合格率只有0.01%，便可视为信得过产品或免检产品，只每月抽检一两次，从而省去大量的人力物力。
  而后验分布π(θ|x)$\pi (\theta |x)$一般也认为是在给定样本X$\mathbf{X}$的情况下θ$\theta$的条件分布，而使π(θ|x)$\pi (\theta |x)$达到最大的值θMD${\theta }_{MD}$称为最大后验估计，类似于经典统计学中的极大似然估计。
  综合起来看，则好比是人类刚开始时对大自然只有少得可怜的先验知识，但随着不断是观察、实验获得更多的样本、结果，使得人们对自然界的规律摸得越来越透彻。所以，贝叶斯方法既符合人们日常生活的思考方式，也符合人们认识自然的规律，经过不断的发展，最终占据统计学领域的半壁江山，与经典统计学分庭抗礼。
  此外，贝叶斯除了提出上述思考模式之外，还特别提出了举世闻名的贝叶斯定理。

1.2 贝叶斯定理

  在引出贝叶斯定理之前，先学习几个定义：

• 条件概率就是事件A$A$在另外一个事件B$B$已经发生条件下的发生概率。条件概率表示为P(A|B)$P(A|B)$，读作“在B$B$条件下A$A$发生的概率”。
• 联合概率表示两个事件共同发生的概率。A$A$与B$B$的联合概率表示为P(A∩B)$P(A\cap B)$或者P(A,B)$P(A,B)$。
• 边缘概率（又称先验概率）是某个事件发生的概率。边缘概率是这样得到的：在联合概率中，把最终结果中那些不需要的事件通过合并成它们的全概率，而消去它们（对离散随机变量用求和得全概率，对连续随机变量用积分得全概率），这称为边缘化（marginalization），比如A$A$的边缘概率表示为P(A)$P(A)$，B$B$的边缘概率表示为P(B)$P(B)$。

  接着，考虑一个问题：P(A|B)$P(A|B)$是在B$B$发生的情况下A$A$发生的可能性。

1. 首先，事件B$B$发生之前，我们对事件A$A$的发生有一个基本的概率判断，称为A$A$的先验概率，用P(A)$P(A)$表示；
2. 其次，事件B$B$发生之后，我们对事件A$A$的发生概率重新评估，称为A$A$的后验概率，用P(A|B)$P(A|B)$表示；
3. 类似的，事件A$A$发生之前，我们对事件B$B$的发生有一个基本的概率判断，称为B$B$的先验概率，用P(B)$P(B)$表示；
4. 同样，事件A$A$发生之后，我们对事件B的发生概率重新评估，称为B的后验概率，用P(B|A)$P(B|A)$表示。

  贝叶斯定理便是基于下述贝叶斯公式：

P(A|B)=P(B|A)P(A)P(B)$$P(A|B)=\frac{P(B|A)P(A)}{P(B)}$$

  上述公式的推导其实非常简单，就是从条件概率推出。   根据条件概率的定义，在事件B发生的条件下事件A发生的概率是：

P(A|B)=P(A∩B)P(B)$$P(A|B)=\frac{P(A\cap B)}{P(B)}$$

  同样地，在事件A发生的条件下事件B发生的概率：

P(B|A)=P(A∩B)P(A)$$P(B|A)=\frac{P(A\cap B)}{P(A)}$$

  整理与合并上述两个方程式，便可以得到：

P(A|B)P(B)=P(A∩B)=P(B|A)P(A)$$P(A|B)P(B)=P(A\cap B)=P(B|A)P(A)$$

  接着，上式两边同除以P(B)$P(B)$，若P(B)$P(B)$是非零的，我们便可以得到贝叶斯定理的公式表达式：

P(A|B)=P(B|A)P(A)P(B)$$P(A|B)=\frac{P(B|A)P(A)}{P(B)}$$

1.3 应用：拼写检查

  经常在网上搜索东西的朋友知道，当你不小心输入一个不存在的单词时，搜索引擎会提示你是不是要输入某一个正确的单词，比如当你在Google中输入“Julw”时，系统会提示你是不是要搜索“July”，如下图所示：    这叫做拼写检查。根据谷歌一员工写的[文章](http://norvig.com/spell-correct.html)显示，Google的拼写检查基于贝叶斯方法。下面我们就来看看，怎么利用贝叶斯方法，实现”拼写检查”的功能。   用户输入一个单词时，可能拼写正确，也可能拼写错误。如果把拼写正确的情况记做c$c$（代表correct），拼写错误的情况记做w$w$（代表wrong），那么”拼写检查”要做的事情就是：在发生w$w$的情况下，试图推断出c$c$。换言之：已知w$w$，然后在若干个备选方案中，找出可能性最大的那个c$c$，也就是求的最大值。   而根据贝叶斯定理，有：

P(c|w)=P(w|c)P(c)/P(w)$$P(c|w)=P(w|c)P(c)/P(w)$$

  由于对于所有备选的c$c$来说，对应的都是同一个w$w$，所以它们的P(w)$P(w)$是相同的，因此我们只要最大化P(w|c)P(c)$P(w|c)P(c)$即可。其中： - P(c)$P(c)$表示某个正确的词的出现”概率”，它可以用”频率”代替。如果我们有一个足够大的文本库，那么这个文本库中每个单词的出现频率，就相当于它的发生概率。某个词的出现频率越高，P(c)$P(c)$就越大。 - P(w|c)$P(w|c)$表示在试图拼写c$c$的情况下，出现拼写错误w$w$的概率。为了简化问题，假定两个单词在字形上越接近，就有越可能拼错，P(w|c)$P(w|c)$就越大。举例来说，相差一个字母的拼法，就比相差两个字母的拼法，发生概率更高。你想拼写单词July，那么错误拼成Julw（相差一个字母）的可能性，就比拼成Jullw高（相差两个字母）。   所以，我们只要找到与输入单词在字形上最相近的那些词，再在其中挑出出现频率最高的一个，就能实现P(w|c)P(c)$P(w|c)P(c)$的最大值。

2 贝叶斯网络

  贝叶斯网络(Bayesian network)，又称信念网络(Belief Network)，或有向无环图模型(directed acyclic graphical model)，是一种概率图模型，于1985年由Judea Pearl首先提出。它是一种模拟人类推理过程中因果关系的不确定性处理模型，其网络拓朴结构是一个有向无环图(DAG)。
  贝叶斯网络的有向无环图中的结点表示随机变量{x1,x2,⋯,xn}$\{x_1,x_2,\cdots ,x_n\}$，它们可以是可观察到的变量，或隐变量、未知参数等。认为有因果关系（或非条件独立）的变量或命题则用箭头来连接（换言之，连接两个节点的箭头代表此两个随机变量是具有因果关系，或非条件独立）。若两个节点间以一个单箭头连接在一起，表示其中一个节点是“因(parents)”，另一个是“果(children)”，两节点就会产生一个条件概率值。
  例如，假设结点E$E$直接影响到结点H$H$，即E→H$E\to H$，则用从E$E$指向H$H$的箭头建立结点E$E$到结点H$H$的有向弧(E,H)$(E,H)$，权值(即连接强度)用条件概率P(H|E)$P(H|E)$来表示，如下图所示：

  简言之，把某个研究系统中涉及的随机变量，根据是否条件独立绘制在一个有向图中，就形成了贝叶斯网络。其主要用来描述随机变量之间的条件依赖，用圈表示随机变量(random variables)，用箭头表示条件依赖(conditional dependencies)。

2.1 贝叶斯网络的定义

  令G=(I,E)$G=(I,E)$表示一个有向无环图(DAG)，其中I$I$代表图形中所有的结点的集合，而E$E$代表有向连接线段的集合，且令x=(xi),xi∈I$x=(x_i),x_i\in I$为其有向无环图中的某一结点xi$x_i$所代表的随机变量，若结点x$x$的联合概率可以表示成：

p(x)=∏xi∈Ip(xi|xpa(i))$$p(x)=\prod _{x_i\in I}p(x_i|x_{pa(i)})$$

  则称x$x$为相对于一有向无环图G$G$的贝叶斯网络，其中，表示结点xi$x_i$之“因”，或称xpa(i)$x_{pa(i)}$是xi$x_i$的parents（父母）。
  此外，对于任意的随机变量，其联合概率可由各自的局部条件概率分布相乘而得出：

p(x1,s,xK)=p(xK|x1,s,xK−1)sp(x2|x1)p(x1)$$p(x_1,s,x_K)=p(x_K|x_1,s,x_{K-1})sp(x_2|x_1)p(x_1)$$

  如下图所示，便是一个简单的贝叶斯网络：

  因为a$a$导致b$b$，a$a$和b$b$导致c$c$，所以有：

p(a,b,c)=p(c|a,b)p(b|a)p(a)$$p(a,b,c)=p(c|a,b)p(b|a)p(a)$$

2.2 贝叶斯网络的实例

  给定如下图所示的贝叶斯网络：

  其中，各个单词、表达式表示的含义如下：

• Smoking表示吸烟，其概率用P(S)$P(S)$表示，lung Cancer表示的肺癌，一个人在吸烟的情况下得肺癌的概率用P(C|S)$P(C|S)$表示，X-ray表示需要照医学上的X光，肺癌可能会导致需要照X光，吸烟也有可能会导致需要照X光（所以smoking也是X-ray的一个因），所以，因吸烟且得肺癌而需要照X光的概率用P(X|C,S)$P(X|C,S)$表示。
• Bronchitis表示支气管炎，一个人在吸烟的情况下得支气管炎的概率用P(B|S)$P(B|S)$，dyspnoea表示呼吸困难，支气管炎可能会导致呼吸困难，肺癌也有可能会导致呼吸困难（所以lung Cancer也是dyspnoea的一个因），因吸烟且得了支气管炎导致呼吸困难的概率用P(D|C,B)$P(D|C,B)$表示。

  lung Cancer简记为C$C$，Bronchitis简记为B$B$，dyspnoea简记为D$D$，且C=0$C=0$表示lung Cancer不发生，C=1$C=1$表示lung Cancer发生。B=0$B=0$表示B$B$不发生，B=1$B=1$表示B发生，也类似于C$C$。同样的，D=1$D=1$表示D$D$发生，D=0$D=0$表示D$D$不发生。便可得到dyspnoea的一张概率表，如上图的最右下角所示。

2.3 贝叶斯网络的3种结构形式

  给定如下图所示的一个贝叶斯网络，

  从图上可以比较直观的看出：

1. x1,x2,⋯,x7$x_1,x_2,\cdots ,x_7$的联合分布为：

p(x1)p(x2)p(x3)p(x4|x1,x2,x3)p(x5|x1,x3)p(x6|x4)p(x7|x4,x5)$$p(x_1)p(x_2)p(x_3)p(x_4|x_1,x_2,x_3)p(x_5|x_1,x_3)p(x_6|x_4)p(x_7|x_4,x_5)$$

1. x1$x_1$和x2$x_2$独立（对应head-to-head）；
2. x6$x_6$和x7$x_7$在x4$x_4$给定的条件下独立（对应tail-to-tail）。

  根据上图，第1点可能很容易理解，但第2、3点中所述的条件独立是啥意思呢？其实第2、3点是贝叶斯网络中3种结构形式中的其中二种。为了说清楚这个问题，需要引入D-Separation（D-分离）这个概念。
  D-Separation是一种用来判断变量是否条件独立的图形化方法。换言之，对于一个DAG(有向无环图)E$E$，D-Separation方法可以快速的判断出两个节点之间是否是条件独立的。

2.3.1 形式1：head-to-head

  贝叶斯网络的第一种结构形式如下图所示：

  所以有：P(a,b,c)=P(a)P(b)P(c|a,b)$P(a,b,c)=P(a)P(b)P(c|a,b)$成立，化简后可得：

∑cP(a,b,c)=∑cP(a)P(b)P(c|a,b)⇒P(a,b)=P(a)P(b)$$\begin{gathered} \sum _{c}P(a,b,c)=\sum _{c}P(a)P(b)P(c|a,b) \\ \Rightarrow P(a,b)=P(a)P(b) \end{gathered}$$

  即在c$c$未知的条件下，a$a$、b$b$被阻断(blocked)，是独立的，称之为head-to-head条件独立，对应本节中最开始那张图中的“x1$x_1$、x2$x_2$独立”。

2.3.2 形式2：tail-to-tail

  贝叶斯网络的第二种结构形式如下图所示：

  有P(a,b,c)=P(c)P(a|c)P(b|c)$P(a,b,c)=P(c)P(a|c)P(b|c)$，则：P(a,b|c)=P(a,b,c)/P(c)$P(a,b|c)=P(a,b,c)/P(c)$，然后将P(a,b,c)=P(c)P(a|c)P(b|c)$P(a,b,c)=P(c)P(a|c)P(b|c)$带入上式，得到：P(a,b|c)=P(a|c)P(b|c)$P(a,b|c)=P(a|c)P(b|c)$。
  即在c$c$给定的条件下，a$a$，b$b$被阻断(blocked)，是独立的，称之为tail-to-tail条件独立，对应本节中最开始那张图中的“x6$x_6$和x7$x_7$在x4$x_4$给定的条件下独立”。

2.3.3 形式3：head-to-tail

  贝叶斯网络的第三种结构形式如下图所示：

  有：P(a,b,c)=P(a)P(c|a)P(b|c)$P(a,b,c)=P(a)P(c|a)P(b|c)$。

  化简后可得：

= = = = P(a,b|c)P(a,b,c) / P(c)P(a)P(c|a)P(b|c) / P(c)P(a,c)P(b|c) / P(c)P(a|c)P(b|c)$$\begin{array}{rl}& P(a,b|c)\\ =& P(a,b,c)/P(c)\\ =& P(a)P(c|a)P(b|c)/P(c)\\ =& P(a,c)P(b|c)/P(c)\\ =& P(a|c)P(b|c)\end{array}$$

  即在c$c$给定的条件下，a$a$，b$b$被阻断(blocked)，是独立的，称之为head-to-tail条件独立。

  插一句：这个head-to-tail其实就是一个链式网络，如下图所示：

  在xi$x_i$给定的条件下，xi+1$x_{i+1}$的分布和x1,x2,⋯,xi−1$x_1,x_2,\cdots ,x_{i-1}$条件独立。即：xi+1$x_{i+1}$的分布状态只和xi$x_i$有关，和其他变量条件独立，这种顺次演变的随机过程，就叫做马尔科夫链（Markov chain）。且有：

P(Xn+1)=x|X0,X1,X2,⋯,Xn)=P(Xn+1=x|Xn)$$P(X_{n+1})=x|X_0,X_1,X_2,\cdots ,X_n)=P(X_{n+1}=x|X_n)$$

  OK，今天在总结贝叶斯网络中的上述3种结构时，发现跟河流关系比较相像，比如：

• ①两条小河流入一条大河，叫head-to-head，由P(a,b,c)=P(c|a,b)P(b)P(a)$P(a,b,c)=P(c|a,b)P(b)P(a)$，可得：P(a,b)=P(a)P(b)$P(a,b)=P(a)P(b)$，即c$c$未知的条件下，a$a$、b$b$被阻断(blocked)，是独立的。同时，也谓之汇连，且汇连是条件依赖的（C$C$依赖于A$A$、B$B$的联合分布），汇连这种情况也称为一个v$v$结构；

• ②一条大河到某处分叉成两条支流，称之为tail-to-tail，由P(a,b,c)=P(b|c)P(a|c)P(c)$P(a,b,c)=P(b|c)P(a|c)P(c)$，可得：P(a,b|c)=P(a|c)P(b|c)$P(a,b|c)=P(a|c)P(b|c)$，即在c$c$给定的条件下，a$a$，b$b$被阻断(blocked)，是独立的。同时，也谓之分连；

• ③一条大河流到底，中间不分叉不汇入其它河流，但可能其中的某段叫什么江，另一段叫什么江，称之为head-to-tail，由P(a,b,c)=P(b|c)P(c|a)P(a)$P(a,b,c)=P(b|c)P(c|a)P(a)$，化简可得：P(a,b,c)=P(a)P(c|a)P(b|c)$P(a,b,c)=P(a)P(c|a)P(b|c)$，即在c$c$给定的条件下，a$a$，b$b$被阻断(blocked)，是独立的。同时，也谓之顺连。

  不知道读者对这个河流的比喻怎么看？^_^
  接着，将上述结点推广到结点集，则是：对于任意的结点集A$A$，B$B$，C$C$，考察所有通过A$A$中任意结点到B$B$中任意结点的路径，若要求A$A$，B$B$条件独立，则需要所有的路径都被阻断(blocked)，即满足下列两个前提之一：

1. A$A$和B$B$的“head-to-tail型”和“tail-to-tail型”路径都通过C$C$；
2. A$A$和B$B$的“head-to-head型”路径不通过C$C$以及C$C$的子孙。

  最后，举例说明上述D-Separation的3种情况，则是如下图所示：

  上图中左边部分是head-to-tail，右边部分的右上角是tail-to-tail，右边部分的右下角是head-to-head。

2.4 因子图

  回到2.2节中那个实例上，如下图所示：

  对于上图，在一个人已经呼吸困难的情况下，其抽样的概率是多少呢？即：

  咱们来一步步计算推导下：

  注：解释下上述式子的第二行到最后一行第三行的推导过程。最开始，所有变量都在∑d=1,b,x,c$\sum _{d=1,b,x,c}$的后面，但由于P(s)$P(s)$和“d=1,b,x,c$d=1,b,x,c$”都没关系，所以，可以提到式子的最前面。而且P(b|s)$P(b|s)$和x$x$、c$c$没关系，所以，也可以把它提出来，放到∑b$\sum _b$的后面，从而式子的右边剩下∑x$\sum _x$和∑c$\sum _c$。
  此外，Variable elimination表示的是变量消除的意思。为了更好的解决此类问题，咱们得引入因子图的概念。

2.4.1 因子图的定义

  wikipedia上是这样定义因子图的：将一个具有多变量的全局函数因子分解，得到几个局部函数的乘积，以此为基础得到的一个双向图叫做因子图（Factor Graph）。
  比如，假定对于函数g(X1,X2,⋯,Xn)$g(X_1,X_2,\cdots ,X_n)$，有下述式子成立：

g(X1,X2,⋯,Xn)=∏j=1mfj(Sj)$$g(X_1,X_2,\cdots ,X_n)=\prod _{j=1}^m{f}_j(S_j)$$

  其中Sj⊆{X1,X2,⋯,Xn}$S_j\subseteq \{X_1,X_2,\cdots ,X_n\}$，其对应的因子图G=(X,F,E)$G=(X,F,E)$包括：

1. 变量结点X={X1,X2,⋯,Xn}$X=\{X_1,X_2,\cdots ,X_n\}$；
2. 因子（函数）结点F={f1,f2,⋯,fm}$F=\{f_1,f_2,\cdots ,f_m\}$；
3. 边E$E$，边通过下列因式分解结果得到：在因子（函数）结点fj$f_j$和变量结点Xk$X_k$之间存在边的充要条件是Xk∈Sj$X_k\in S_j$存在。

  官方正式的定义果然晦涩！我相信你没看懂。通俗来讲，所谓因子图就是对函数进行因子分解得到的一种概率图。一般内含两种结点，变量节点和函数结点。我们知道，一个全局函数通过因式分解能够分解为多个局部函数的乘积，这些局部函数和对应的变量关系就体现在因子图上。举个例子，现在有一个全局函数，其因式分解方程为：

g(x1,x2,x3,x4,x5)=fA(x1)fB(x2)fC(x1,x2,x3)fD(x3,x4)fE(x3,x5)$$g(x_1,x_2,x_3,x_4,x_5)=f_A(x_1)f_B(x_2)f_C(x_1,x_2,x_3)f_D(x_3,x_4)f_E(x_3,x_5)$$

  其中fA,fB,fC,fD,fE$f_A,f_B,f_C,f_D,f_E$为各函数，表示变量之间的关系，可以是条件概率也可以是其他关系(如马尔可夫随机场Markov Random Fields中的势函数)。
  为了方便表示，可以写成：

g(X)=∏E⊆Qfε(Xε)$$g(X)=\prod _{E\subseteq Q}{f}_{\epsilon }(X_{\epsilon })$$

  其对应的因子图为：

  且上述因子图等价于：

  所以，在因子图中，所有顶点不是变量节点就是函数节点，边线表示它们之间的函数关系。
  但搞了半天，虽然知道了什么是因子图，但因子图到底是干嘛的呢？为何要引入因子图，其用途和意义何在？事实上，因子图跟贝叶斯网络和马尔可夫随机场（Markov Random Fields）一样，也是是概率图的一种。从上文中，我们可以看到，在概率图中，求某个变量的边缘分布是常见的问题。这问题有很多求解方法，其中之一就是可以把贝叶斯网络或马尔科夫随机场转换成因子图，然后用sum-product算法求解。换言之，基于因子图可以用sum-product 算法高效的求各个变量的边缘分布。
  先举个例子说明如何把贝叶斯网络（和马尔科夫随机场）转换成因子图。给定下图所示的贝叶斯网络或马尔科夫随机场：

  根据各个变量对应的关系，可得：

p(u,w,x,y,z)=p(u)p(w)p(x|u,w)p(y|x)p(z|x)$$p(u,w,x,y,z)=p(u)p(w)p(x|u,w)p(y|x)p(z|x)$$

  其对应的因子图为（以下两种因子图的表示方式皆可）：

  由上述例子总结出由贝叶斯网络构造因子图的方法：

• 贝叶斯网络中的一个因子对应因子图中的一个结点
• 贝叶斯网络中的每一个变量在因子图上对应边或者半边
• 结点g$g$和边x$x$相连当且仅当变量x$x$出现在因子g$g$中

  举几个例子。比如，对于如下的一个因子图：

  有：

f(x1,x2,x3,x4,x5)=fA(x1,x2,x3)fB(x3,x4,x5)fC(x4)$$f(x_1,x_2,x_3,x_4,x_5)=f_A(x_1,x_2,x_3)f_B(x_3,x_4,x_5)f_C(x_4)$$

  而对于下图所示的马尔科夫链：

  有：

pXYZ(x,y,z)=pX(x)pY|X(y|x)pZ|Y(z|y)$$p_{XYZ}(x,y,z)=p_X(x)p_{Y|X}(y|x)p_{Z|Y}(z|y)$$

  另对于如下图所示的隐马尔科夫模型：

  有：

p(x0,x1,x2,⋯,xn,y1,y2,⋯,yn)=p(x0)∏k=1np(xk|xk−1)p(yk|xk−1)$$p(x_0,x_1,x_2,\cdots ,x_n,y_1,y_2,\cdots ,y_n)=p(x_0)\prod _{k=1}^{n}p(x_k|x_{k-1})p(y_k|x_{k-1})$$

2.4.2 Sum-product算法

  下面，咱们来考虑一个问题：即如何由联合概率分布求边缘概率分布。
  首先回顾下联合概率和边缘概率的定义，如下：

• 联合概率表示两个事件共同发生的概率。A与B的联合概率表示为P(A∩B)$P(A\cap B)$或者P(A,B)$P(A,B)$。
• 边缘概率（又称先验概率）是某个事件发生的概率。边缘概率是这样得到的：在联合概率中，把最终结果中不需要的那些事件合并成其事件的全概率而消失（对离散随机变量用求和得全概率，对连续随机变量用积分得全概率）。这称为边缘化（marginalization）。A的边缘概率表示为P(A)，B的边缘概率表示为P(B)$P(B)$。

  事实上，某个随机变量fk$f_k$的边缘概率可由x1,x2,x3,⋯,xn$x_1,x_2,x_3,\cdots ,x_n$的联合概率求到，具体公式为：

f¯k(xk)≜∑x1,⋯,xnexceptxkf(x1,⋯,xn)$${\overline{f}}_k(x_k)\triangleq \sum _{x_1,\cdots ,x_n\mathrm{e}\mathrm{x}\mathrm{c}\mathrm{e}\mathrm{p}\mathrm{t}{x}_k}f(x_1,\cdots ,x_n)$$

  啊哈，啥原理呢？原理很简单，还是它：对x3$x_3$外的其它变量的概率求和，最终剩下x3$x_3$的概率！
  此外，换言之，如果有：

f(x1,⋯,xn)=f1(x1)f2(x2)⋯fn(xn)$$f(x_1,\cdots ,x_n)=f_1(x_1)f_2(x_2)\cdots f_n(x_n)$$

  那么：

f¯k(xk)=(∑x1f1(x1))⋯(∑xk−1fk−1(xk−1))fk(xk)⋯(∑xnfn(xn))$${\overline{f}}_k(x_k)=\left(\sum _{x_1}{f}_1(x_1)\right)\cdots \left(\sum _{x_{k-1}}{f}_{k-1}(x_{k-1})\right)f_k(x_k)\cdots \left(\sum _{x_n}{f}_n(x_n)\right)$$

  上述式子如何进一步化简计算呢？考虑到我们小学所学到的乘法分配率，可知a×b+a×c=a×(b+c)$a\times b+a\times c=a\times (b+c)$，前者2次乘法1次加法，后者1次乘法，1次加法。我们这里的计算是否能借鉴到分配率呢？别急，且听下文慢慢道来。
  假定现在我们需要计算计算如下式子的结果：

f¯3(x3)=∑x1,⋯,x7exceptx3f(x1,⋯,x7)$${\overline{f}}_3(x_3)=\sum _{x_1,\cdots ,x_7\mathrm{e}\mathrm{x}\mathrm{c}\mathrm{e}\mathrm{p}\mathrm{t}{x}_3}f(x_1,\cdots ,x_7)$$

  同时，f$f$能被分解如下：

  借鉴分配率，我们可以提取公因子：

  因为变量的边缘概率等于所有与他相连的函数传递过来的消息的积，所以计算得到：

  仔细观察上述计算过程，可以发现，其中用到了类似“消息传递”的观点，且总共两个步骤。

  第一步、对于f$f$的分解图，根据蓝色虚线框、红色虚线框围住的两个box外面的消息传递：

  计算可得：

  第二步、根据蓝色虚线框、红色虚线框围住的两个box内部的消息传递：

  根据μ⃗ X1(x1)≜f1(x1),μ⃗ X2(x2)≜f2(x2)${\overrightarrow{\mu }}_{X_1}(x_1)\triangleq f_1(x_1),{\overrightarrow{\mu }}_{X_2}(x_2)\triangleq f_2(x_2)$，我们有：

  就这样，上述计算过程将一个概率分布写成两个因子的乘积，而这两个因子可以继续分解或者通过已知得到。这种利用消息传递的观念计算概率的方法便是sum-product算法。前面说过，基于因子图可以用sum-product算法可以高效的求各个变量的边缘分布。
  到底什么是sum-product算法呢？sum-product算法，也叫belief propagation，有两种消息：

• 一种是变量(Variable)到函数(Function)的消息：mx→f$m_{x\to f}$，如下图所示：

  此时，变量到函数的消息为mx→f=1$m_{x\to f}=1$。

• 另外一种是函数(Function)到变量(Variable)的消息：mf→x$m_{f\to x}$。如下图所示：

  此时，函数到变量的消息为：mf→x=f(x)$m_{f\to x}=f(x)$。
  以下是sum-product算法的总体框架：

• 1、给定如下图所示的因子图：

• 2、sum-product 算法的消息计算规则为：

μ⃗ X(x)=∑y1,⋯,yng(x,y1,⋯,yn)μ⃗ Y1(y1)⋯μ⃗ Yn(yn)$${\overrightarrow{\mu }}_X(x)=\sum _{y_1,\cdots ,y_n}g(x,y_1,\cdots ,y_n){\overrightarrow{\mu }}_{Y_1}(y_1)\cdots {\overrightarrow{\mu }}_{Y_n}(y_n)$$

• 3、根据sum-product定理，如果因子图中的函数f$f$没有周期，则有：

f¯X(x)=μ⃗ X(x)μ←X(x)$${\overline{f}}_X(x)={\overrightarrow{\mu }}_X(x){\overleftarrow{\mu }}_X(x)$$

  值得一提的是：如果因子图是无环的，则一定可以准确的求出任意一个变量的边缘分布，如果是有环的，则无法用sum-product算法准确求出来边缘分布。
  比如，下图所示的贝叶斯网络：

  其转换成因子图后，为：

  可以发现，若贝叶斯网络中存在“环”（无向），则因此构造的因子图会得到环。而使用消息传递的思想，这个消息将无限传输下去，不利于概率计算。
   解决方法有3个：

• 1、删除贝叶斯网络中的若干条边，使得它不含有无向环
    比如给定下图中左边部分所示的原贝叶斯网络，可以通过去掉C和E之间的边，使得它重新变成有向无环图，从而成为图中右边部分的近似树结构：
  
  
  
    具体变换的过程为最大权生成树算法MSWT，通过此算法，这课树的近似联合概率P′(x)$P^{\prime }(x)$和原贝叶斯网络的联合概率P(x)$P(x)$的相对熵（如果忘了什么叫相对熵，请参阅：最大熵模型中的数学推导）最小。
• 2、重新构造没有环的贝叶斯网络
• 3、选择loopy belief propagation算法（你可以简单理解为sum-product 算法的递归版本），此算法一般选择环中的某个消息，随机赋个初值，然后用sum-product算法，迭代下去，因为有环，一定会到达刚才赋初值的那个消息，然后更新那个消息，继续迭代，直到没有消息再改变为止。唯一的缺点是不确保收敛，当然，此算法在绝大多数情况下是收敛的。

  此外，除了这个sum-product算法，还有一个max-product 算法。但只要弄懂了sum-product，也就弄懂了max-product 算法。因为max-product 算法就在上面sum-product 算法的基础上把求和符号换成求最大值max的符号即可！
  最后，sum-product 和 max-product 算法也能应用到隐马尔科夫模型hidden Markov models上，后面有机会的话可以介绍。

3 参考文献和推荐阅读

1. Thomas Bayes “An essay towards solving a Problem in the Doctrine of Chances”（贝叶斯定理原始论文）：http://www.sbs-bvs.be/bsn57/bsn57-6.pdf；
2. 《数理统计学简史 第三章 贝叶斯方法》；
3. 《贝叶斯统计 茆诗松著》；
4. 北京10月机器学习班第9次课，邹博讲贝叶斯网络的PPT：http://pan.baidu.com/s/1o69Lp1K；
5. 相关wikipedia，比如贝叶斯定理的wiki：http://zh.wikipedia.org/zh/%E8%B4%9D%E5%8F%B6%E6%96%AF%E5%AE%9A%E7%90%86，贝叶斯网络的wiki：http://zh.wikipedia.org/wiki/%E8%B2%9D%E6%B0%8F%E7%B6%B2%E8%B7%AF。因子图中文wiki：http://zh.wikipedia.org/zh/%E5%9B%A0%E5%AD%90%E5%9B%BE，英文wik：http://en.wikipedia.org/wiki/Factor_graph。
6. 《统计决策论及贝叶斯分析 James O.Berger著》；
7. 贝叶斯定理：http://www.guokr.com/question/547339/；
8. 贝叶斯推断及其互联网应用（一）：定理简介http://www.ruanyifeng.com/blog/2011/08/bayesian_inference_part_one.html；
9. 贝叶斯推断及其互联网应用（三）：拼写检查http://www.ruanyifeng.com/blog/2012/10/spelling_corrector.html；
10. Google研发总监Peter Norvig解释拼写检查的原理：http://norvig.com/spell-correct.html；
11. http://www.eng.yale.edu/pjk/eesrproj_02/luckenbill_html/node4.html(sum-product)；
12. Pattern Recognition and Machine Learning Chapter 8, M. Jordan, J. Kleinberg, ect, 2006；
13. D-Separation(D分离)-PRML-8.22-Graphical Model by 小军：http://www.zhujun.me/d-separation-separation-d.html；
14. 因子图介绍 by Hans-Andrea Loeliger：http://www.robots.ox.ac.uk/~parg/mlrg/papers/factorgraphs.pdf；
15. http://netclass.csu.edu.cn/jpkc2003/rengongzhineng/rengongzhineng/kejian/ai/ai/chapter4/442.htm；
16. 贝叶斯网的R实现（ Bayesian networks in R）（二）bnlearn(2)：http://site.douban.com/182577/widget/notes/12817482/note/283039795/；
17. 知乎上关于贝叶斯学派跟频率派的区别的讨论：http://www.zhihu.com/question/20587681；
18. factor graph，因子图，势函数potential function，Template models：http://www.cnblogs.com/549294286/archive/2013/06/06/3121454.html；
19. Online Bayesian Probit Regression介绍之Factor Graph：http://www.doingkong.com/?p=68；
20. An Introduction to Factor Graphs,Hans-Andrea Loeliger,MLSB 2008：http://people.binf.ku.dk/~thamelry/MLSB08/hal.pdf；
21. Factor graph and sum-product algorithm, Frank R. Kschischang, Brendan J.Frey, ect, 1998：http://filebox.vt.edu/~rmtaylor/Graphical_Modeling/Intro_and_tutorial/Kschischang_ffg_sumproduct.pdf；
22. A Tutorial on Inference and Learning in Bayesian Networks, Irina Rish：http://www.ee.columbia.edu/~vittorio/Lecture12.pdf；
23. Probabilistic Graphical Models Directed GMs: Bayesian Networks：http://www.cs.cmu.edu/~epxing/Class/10708/lectures/lecture2-BNrepresentation.pdf；
24. A Brief Introduction to Graphical Models and Bayesian Networks By Kevin Murphy, 1998：http://www.cs.ubc.ca/~murphyk/Bayes/bayes.html。

- - - - - - - - - - 转自zdy0_2004的专栏 - - - - - - - - - -