貝氏網路

貝氏網路（Bayesian network），又稱信念網路（belief network）或是有向無環圖模型（directed acyclic graphical model），是一種機率圖型模型，藉由有向無環圖（directed acyclic graphs, or DAGs）中得知一組隨機變數及其 n 組條件機率分配（conditional probability distributions, or CPDs）的性質。舉例而言，貝氏網路可用來表示疾病和其相關症狀間的機率關係；倘若已知某種症狀下，貝氏網路就可用來計算各種可能罹患疾病之發生機率。

一般而言，貝氏網路的有向無環圖中的節點表示隨機變數，它們可以是可觀察到的變數，抑或是潛在變數、未知引數等。連線兩個節點的箭頭代表此兩個隨機變數是具有因果關係或是非條件獨立的；而兩個節點間若沒有箭頭相互連線一起的情況就稱其隨機變數彼此間為條件獨立。若兩個節點間以一個單箭頭連線在一起，表示其中一個節點是「因（parents）」，另一個是「果（descendants or children）」，兩節點就會產生一個條件機率值。比方說，我們以表示第 i 個節點，而的「因」以表示，的「果」以表示；圖一就是一種典型的貝氏網路結構圖，依照先前的定義，我們就可以輕易的從圖一可以得知：

，以及

大部分的情況下，貝氏網路適用在節點的性質是屬於離散型的情況下，且依照此條件機率寫出，此條件機率表的每一行（row）列出所有可能發生的，每一列（column）列出所有可能發生的，且任一行的機率總和必為 1。寫出條件機率表後就很容易將事情給條理化，且輕易地得知此貝氏網路結構圖中各節點間之因果關係；但是條件機率表也有其缺點：若是節點是由很多的「因」所造成的「果」，如此條件機率表就會變得在計算上既複雜又使用不便。下圖為圖一貝氏網路中某部分結構圖之條件機率表。

數學定義

令 G = (I,E) 表示一個有向無環圖（DAG），其中 I 代表圖中所有的節點的集合，而 E 代表有向連線線段的集合，且令 X = (Xi) i∈I 為其有向無環圖中的某一節點 i 所代表之隨機變數，若節點 X 的聯合機率分配可以表示成：

則稱 X 為相對於一有向無環圖 G 的貝氏網路，其中表示節點 i 之「因」。

對任意的隨機變數，其聯合分配可由各自的區域性條件機率分配相乘而得出：

依照上式，我們可以將一貝氏網路的聯合機率分配寫成：

（對每個相對於 Xi 的「因」變數 Xj 而言)

上面兩個表示式之差別在於條件機率的部分，在貝氏網路中，若已知其「因」變數下，某些節點會與其「因」變數條件獨立，只有與「因」變數有關的節點才會有條件機率的存在。

如果聯合分配的相依數目很稀少時，使用貝氏函式的方法可以節省相當大的記憶體容量。舉例而言，若想將 10 個變數其值皆為 0 或 1 儲存成一條件機率表型式，一個直觀的想法可知我們總共必須要計算個值；但若這 10 個變數中無任何變數之相關「因」變數是超過三個以上的話，則貝氏網路的條件機率表最多隻需計算個值即可。另一個貝式網路優點在於：對人類而言，它更能輕易地得知各變數間是否條件獨立或相依與其區域性分配（local distribution）的型態來求得所有隨機變數之聯合分配。

馬可夫毯（Markov blanket）

定義一個節點之馬可夫毯為此節點的因節點、果節點與果節點的因節點所成之集合。一旦給定其馬可夫毯的值後，若網路內之任一節點 X 皆會與其他的節點條件獨立的話，就稱 X 為相對於一有向無環圖 G 的貝氏網路。

求解方法

以上例子是一個很簡單的貝氏網路模型，但是如果當模型很複雜時，這時使用列舉式的方法來求解機率就會變得非常複雜且難以計算，因此必須使用其他的替代方法。一般來說，貝氏機率有以下幾種求法：

4.1 精確推論

列舉推理法（如上述例子）
變數消元演演算法（variable elimination）

4.2 隨機推論（蒙地卡羅方法）

直接取樣演演算法
拒絕取樣演演算法
概似加權演演算法
馬可夫鏈蒙地卡羅演演算法（Markov chain Monte Carlo algorithm）

在此，以馬可夫鏈蒙地卡羅演演算法為例，又馬可夫鏈蒙地卡羅演演算法的型別很多，故在這裡只說明其中一種吉布斯取樣的操作步驟：

首先將已給定數值的變數固定，然後將未給定數值的其他變數隨意給定一個初始值，接著進入以下迭代步驟：

（1）隨意挑選其中一個未給定數值的變數

（2）從條件分配抽樣出新的的值，接著重新計算

當迭代結束後，刪除前面若干筆尚未穩定的數值，就可以求出的近似條件機率分配。馬可夫鏈蒙地卡羅演演算法的優點是在計算很大的網路時效率很好，但缺點是所抽取出的樣本並不具獨立性。

當貝氏網路上的結構跟引數皆已知時，我們可以透過以上方法來求得特定情況的機率，不過，如果當網路的結構或引數未知時，我們必須藉由所觀測到的資料去推估網路的結構或引數，一般而言，推估網路的結構會比推估節點上的引數來的困難。依照對貝氏網路結構的瞭解和觀測值的完整與否，我們可以分成下列四種情形：

以下就結構已知的部分，作進一步的說明。

4.2.1 1. 結構已知，觀測值完整：

此時我們可以用最大概似估計法（MLE）來求得引數。其對數概似函式為

其中代表的因變數，代表第個觀測值，N 代表觀測值資料的總數。

以圖二當例子，我們可以求出節點 U 的最大概似估計式為

由上式就可以藉由觀測值來估計出節點 U 的條件分配。如果當模型很複雜時，這時可能就要利用數值分析或其它最佳化技巧來求出引數。

4.2.2 2. 結構已知，觀測值不完整（有遺漏資料）：

如果有些節點觀測不到的話，可以使用 EM 演演算法（Expectation-Maximization algorithm）來決定出引數的區域最佳概似估計式。而 EM 演演算法的主要精神在於如果所有節點的值都已知下，在 M 階段就會很簡單，如同最大概似估計法。而 EM 演演算法的步驟如下：

（1）首先給定欲估計的引數一個起始值，然後利用此起始值和其他的觀測值，求出其他未觀測到節點的條件期望值，接著將所估計出的值視為觀測值，將此完整的觀測值帶入此模型的最大概似估計式中，如下所示（以圖二為例）：

其中代表在目前的估計引數下，事件 x 的條件機率期望值為

（2）最大化此最大概似估計式，求出此引數之最有可能值，如此重複步驟（1）與（2），直到引數收斂為止，即可得到最佳的引數估計值。

補充例子（列舉推理法）

讓我們考慮一個應用在醫藥上的機率推論例子，在此病人會被診斷出是否有呼吸困難的症狀。表一代表一個我們所觀測到的資料集合，包含 10 筆觀測值，S 代表的是吸菸與否（Smoker），C 代表是否為罹癌者（Cancer），B 代表是否罹患支氣管炎（bronchitis），D 代表是否有呼吸困難及咳嗽（dyspnea and asthma）的症狀。『1』和『0』分別代表『是』和『否』。此醫藥網路結構顯示於圖三。

表二代表的是整個網路的經驗聯合機率分配，是由所收集到的資料所建構而成，利用此表可建構出節點的聯合機率分配。見圖四。此貝氏公式可利用節點的邊際機率和聯合機率去計算節點的條件機率，待會會應用在建立條件機率表格（Conditional probability Table; CPT）上。見圖五。貝氏網路的聯合機率可由下列式子計算：

其值見表三。

使用整個網路經驗聯合機率分配所計算出來的值會與使用 CPT 所計算出來的值不同，其差異可由表二和表三得知。其中差異不只是值的不同，也出現了新事件的機率（原本所沒觀察到的事件）。

建立在觀測資料上的機率推論演演算法：

資料集合，其中（下標代表第幾個觀測值，上標代表第幾個變數），且一共有 n 個觀測值。每一個觀測值包含個變數，第 j 個變數有狀態。

此貝氏網路的結構 G 代表 N 個前代（predecessors）節點集合，也就是對第 j 個節點，為其親代節點的集合，j=1,2,…,N

範例點（instantiated node）為節點在已知狀態，即在此狀態的機率為 1。如果範例點為空集合，將使用古典機率推論

使用表一的觀測值和圖一的貝氏網路結構，並且已知範例點（instantiated node）為，也就是病人為非吸菸者和罹癌者：

問題：

病人患有支氣管炎的機率

病人會有呼吸困難的機率

解答：

故為 0.2

:

:

::

:

:

::

應用

貝氏網路目前應用在類比計算生物學與生物資訊學、蛋白質結構、基因表達分析、醫學、文件分類、資訊檢索、決策支援系統、工程學、、影像處理等。