機率布爾網路模型由Shmulevich等人於2002年提出,該模型是布爾網路模型的拓展。該模型中每個基因對應一個或多個布爾函式,每個疊代步驟中以一定機率選取其中一個函式預測其下一步的狀態。機率布爾網路保留了布爾網路的部分特性並且可以描述基因調控網路中的隨機現象。
機率布爾網為基因調控網路提供了一個更靈活的機率模型,以模擬生物系統中的不確定性。
基本介紹
- 中文名:機率布爾網、機率布爾網路
- 外文名:Probabilistic Boolean Network
- 縮寫:PBN
- 套用領域:系統生物學、布爾邏輯動力學系統
- 更新函式:布爾邏輯函式
- 相關概念:布爾網路
定義,例子,與布爾網路模型的異同,研究難點,
定義
機率布爾網路
中包含一個節點集合
和一個布爾函式集合序列
。集合
中的
個布爾函式作為節點
的預測函式,
![](/img/9/e9b/9208865cb0bed4e85531c2f9a473.jpg)
![](/img/d/c0d/aec65fff9f24c831b173ec41c840.jpg)
![](/img/f/d83/8ead4f280cf75af995b1792724e0.jpg)
![](/img/1/e81/bebdf249e0e7013f90315bf3c1c0.jpg)
![](/img/8/7cc/f929f613414313b9fd7437dbf00e.jpg)
![](/img/5/ec5/eb148c6838f58e01fbd19b58c947.jpg)
![](/img/7/a84/575b04ca00771cc848d362df90fb.jpg)
在每一步網路狀態轉移過程中,從集合
中以一定機率隨機選擇一個布爾函式作為節點
的預測函式,所有被選取的預測函式構成了一個函式向量
,該函式向量稱為
的一個實現(realization)。
的每個實現實際上都對應一個BN,因此可以看做是若干BN的機率組合。若每個節點預測函式的選擇是相互獨立的,則稱該為獨立的。獨立的包含的BN的個數為
![](/img/1/e81/bebdf249e0e7013f90315bf3c1c0.jpg)
![](/img/b/e3e/9d0e1eaa2ad3477fdb7efd9b7b0f.jpg)
![](/img/9/a6e/e2cb6bf3e54d6270aa9768dfb1d9.jpg)
![](/img/4/e47/fb7e18eb274f1e2321d2e9f00fd2.jpg)
![](/img/4/e47/fb7e18eb274f1e2321d2e9f00fd2.jpg)
![](/img/3/510/9e43b984464523c92a4f5c96384a.jpg)
![](/img/2/2f0/18cf7b84ed7a2053b8d00b6f3edc.jpg)
![](/img/d/6fd/89d0c83b45a1a7032b32a3bba4fb.jpg)
若對所有
都與
,則BN的個數
,此時
退化為BN。
![](/img/2/73c/bf3d338da894aac10878823c6951.jpg)
![](/img/e/397/091dfd4ed57f3ca96adf18027b13.jpg)
![](/img/3/cf2/02c865e1ad5a6f3a5a861a1ffbfc.jpg)
![](/img/4/e47/fb7e18eb274f1e2321d2e9f00fd2.jpg)
機率布爾網路的動態特性可用馬爾科夫鍊表述,其狀態轉移矩陣是一個
的矩陣,定義為
![](/img/6/845/a690aecd02db4b609f2b11db8e6f.jpg)
![](/img/8/fa3/f9ee12661992e7a831ba67695b12.jpg)
![](/img/4/e47/fb7e18eb274f1e2321d2e9f00fd2.jpg)
![](/img/6/b57/204f004494dd1665c9b399ef7e07.jpg)
![](/img/d/1e3/86d5e55a30f7ee718452ac07476a.jpg)
例子
假設
是一個包含三個節點的機率布爾網路,其中
,
,
,
,
,預測函式的真值表如表所示
![](/img/9/e9b/9208865cb0bed4e85531c2f9a473.jpg)
![](/img/a/bff/6ea3cfe759f4e1d3878b3110cd78.jpg)
![](/img/1/7ac/e84ec917eb7c50d488e3c2b7ea73.jpg)
![](/img/f/2e4/bcf5e9ba09919a74bfce78226cd3.jpg)
![](/img/2/6d9/b0b79f95af09bd66958df28a4130.jpg)
![](/img/b/32d/26bc57fc650c99d536d83d624627.jpg)
000 | 0 | 0 | 0 | 0 | 0 |
001 | 0 | 1 | 0 | 0 | 0 |
010 | 1 | 0 | 0 | 0 | 0 |
011 | 1 | 1 | 1 | 1 | 0 |
100 | 0 | 0 | 0 | 0 | 0 |
101 | 1 | 1 | 1 | 0 | 1 |
110 | 1 | 1 | 1 | 0 | 0 |
111 | 1 | 1 | 1 | 1 | 1 |
![]() | 1 | 0.6 | 0.4 | 0.6 | 0.4 |
該機率布爾網路包含八個狀態,其狀態之間的轉移關係如圖所示
![機率布爾網 機率布爾網](/img/2/4f0/nBnauIWMzI2MmF2MwAjNilzMlRmMmRmNzcDNlFWZ5MWNiJWNiNzMwIWO1MzLtVGdp9yYpB3LltWahJ2Lt92YuUHZpFmYuMmczdWbp9yL6MHc0RHa.jpg)
根據狀態轉移矩陣的定義可知,
![](/img/0/e1d/1449f58381df8ab7d38085fd5717.jpg)
與布爾網路模型的異同
- 機率布爾網路與布爾網路的狀態空間一致,大小均為2n;
- 布爾網路中任意狀態都唯一確定一個後續狀態,而機率布爾網路由於包含多個布爾網路,因此每個狀態可能對應多個後續狀態;
- 在布爾網路模型中,每個節點對應一個布爾函式,以該布爾函式作為其預測函式,布爾網路是一個決定性網路,網路當前狀態完全且唯一決定網路的下一步狀態。網路的狀態沿初始狀態在狀態空間中形成一條軌跡。由於狀態空間大小是有限的,這條軌跡最終將進入不動點或極限圈。而在機率布爾網路中,每個節點對應一個或多個預測函式,網路當前狀態可能對應多個後續狀態。
研究難點
機率布爾網路研究中的一個重要問題是求其馬爾科夫鏈相應的機率轉移矩陣,並計算該矩陣的穩態分布。機率轉移矩陣的穩態分布並不一定是存在的,為保證其穩態分布存在,可以在機率布爾網路中添加一個參數 。機率布爾網路中的節點狀態以機率 受控於其預測函式集,而以機率 隨機產生一個後續狀態。這樣,機率布爾網路的狀態是強連通的,從而保證機率轉移矩陣的穩態分布存在。
機率布爾網路的狀態轉移矩陣是一個
的矩陣,當網路節點個數力很大時,求解其穩態分布很困難,因此如何簡化機率布爾網路模型是一個有意義的問題。目前機率布爾網路的一個簡化方法是忽略其中的一些小機率的布爾網路,其它的簡化方法還有待探索。
![](/img/6/845/a690aecd02db4b609f2b11db8e6f.jpg)