ベイズ法のソースを表示

{{Expand English|Bayesian inference in phylogeny|date=2024年5月}}
'''ベイズ法'''（ベイズほう、[[英語|英]]：'''Bayesian inference method'''）は、[[生物]]の系統進化を示す[[系統樹]]を推定する手法の一つ。[[ベイズの定理]]に基づいて[[尤度]]を通してデータを加味した[[事後確率]]分布を目的関数にとり、[[マルコフ連鎖モンテカルロ法]]を適用して事後確率分布を[[ベイズ推定|推定]]し、その[[期待値]]としての最良の樹形を選択する。ベイズ法を利用した系統推定ソフトウェアではMrBayesが代表的である<ref name=三中2009>{{Cite journal|和書|author=[[三中信宏]] |pages=192-196 |journal=植物防疫 |publisher=[[日本植物防疫協会]] |volume=63 |issue=3 |year=2009 |url=http://jppa.or.jp/archive/pdf/63_03_66.pdf |title=分子系統学：最近の進歩と今後の展望}}</ref>。

== 理論 ==
=== 基本的理論 ===
ベイズ法では、[[最尤推定|最尤法]]が[[ベイズ推定]]の枠組みの中で再構築されている。最尤法ではある特定の仮定（進化モデル）の下で与えられた配列群の[[尤度]]を最大化する系統樹が最適樹として選択される。すなわち、与えられた配列に対する系統樹の枝長と進化モデルの尤度が最尤法では目的関数として取られている。一方で、ベイズ法の目的関数は、与えられた配列に対する系統樹の樹形と枝長と進化モデルの事後確率となる。最適樹は事後確率分布の[[最頻値]]を最大事後確率推定した系統樹として選択される<ref name=松井2021>{{Cite journal|和書|title=分子系統解析の最前線 |author=松井求 |journal=JSBi Bioinformatics Review |volume=2 |issue=1 |pages=30-57 |year=2021 |url= https://doi.org/10.11234/jsbibr.2021.7 |doi=10.11234/jsbibr.2021.7}}{{オープンアクセス}}</ref>。

ただし、特に樹形の事前分布を一様分布とした場合、最尤法による系統樹とベイズ法による系統樹はほぼ一致し、ベイズ法はむしろ最尤法で必要のない事後確率を計算する回りくどい手法となる。そこでベイズ法では単一の最適樹を求めるのではなく、複数の最適樹を事後確率と共に提示するか、あるいはそれらを統合して単一の合意樹を作成することを主目的とする。これにより、[[遺伝子の水平伝播]]などにより系統樹が一意に定まらない場合や、複数の系統仮説が求められる場合に、ベイズ法は効果を発揮することになる<ref name=松井2021/>。

=== アルゴリズム ===
系統樹<math>\tau_i</math>の事前確率分布を<math>f(\tau_i)</math>とする。配列データXが与えられた場合の<math>\tau_i</math>の事後確率分布<math>P(\tau_i|X)</math>は、<math>f(X|\tau_i)</math>を最尤法により求められる尤度関数とすると、[[ベイズの定理]]に従って以下のように表される<ref name=松井2021/>。

<math>P(\tau_i|X)=\frac {f(X|\tau_i)f(\tau_i)} {\sum^{|\tau|}_{j=i} f(X|\tau_i)f(\tau_i)}</math>

なお、系統樹<math>\tau_i</math>の枝長と、その進化モデルの持つ全てのパラメータも確率分布として予測でき、実際にベイズ法で求める確率分布は分母・分子ともに枝長とパラメータの[[多重積分|二重積分]]の形で表現できる<ref name=松井2021/>。ただしいずれにせよ、分母は非常に小規模な系統樹でなければ最尤法による計算が不可能であるため、ベイズ法では[[マルコフ連鎖モンテカルロ法]]（MCMC）を利用して近似的に計算し、同時に最適な樹形の探索も行っている<ref name=松井2021/>。

MCMCの手順は次の通りである<ref name=松井2021/>。
#[[近隣結合法]]や[[最大節約法]]により構築された系統樹か、あるいはランダム系統樹を初期系統樹とする。初期系統樹の枝長と進化モデルのパラメータも設定し、初期系統樹の尤度を求める。
#網羅的探索により新たな系統樹を作成し、枝長と進化モデルのパラメータにも変更を加え、新たな系統樹の尤度を求める。
#新旧系統樹の尤度を比較して、新系統樹の尤度が大である場合に、新系統樹を採用する。新系統樹の尤度が小であるならば低い確率で新系統樹を採用するが、多くの場合は棄却し、別の新たな系統樹を作成する。
#系統樹とその樹形およびパラメータを記録し、2へ戻る。このサイクルを延々と繰り返す。

手順を繰り返すうちに系統樹は定常状態に達する。定常状態に達するまでの記録は無視し、定常状態に達している記録から一定の回数ごとに標本を抽出すると、その標本は系統樹の分布に従うものになっているため、標本系統樹群から上式の分母を近似的に計算できる。こうして事後確率分布を求められる<ref name=松井2021/>。また、ここでの樹形の出現頻度が系統樹の事後確率になる<ref name=仲田2006/>。

その後は、事後確率が最大となる単一の系統樹を選択するか、事後確率の高い順に複数の系統樹を選択して合意樹を形成する<ref name=松井2021/>。

== MrBayes ==
ベイズ法は1996年頃からYangにより系統推定に導入され、2001年にはベイズ法を利用した系統推定ソフトウェアMrBayesが登場した<ref name=松井2021/>。MrBayesはベイズ法を用いる系統推定ソフトウウェアの中では最も広く使用されているものである<ref name=三中2009/>。MrBayesの登場により、より効率的なアルゴリズムや現実的な進化モデルが実装されたこともあり、ベイズ法は広く一般に普及した<ref name=松井2021/>。

MrBayesでは、事後確率分布がある程度高い場所（局所的な最適解）でマルコフ連鎖が定常状態に達しないよう、MCMCを改良したMetropolis Coupled MCMC（MCMCMC）が採用されている。MCMCMCでは複数本のMCMCを走らせており、そのうち1本をcold chain、残りをhot chainとする。cold chainは元々のMCMCであり、局所的な最適解を脱出しにくい一方で、膨大な時間をかければ正しい事後確率分布に収束する。hot chain では目標分布が平坦化されており、系統樹が撹乱を受けるため、正しい事後確率分布への収束の保証を破棄する一方で局所的な最適解からの脱出率を高めている。また、一定回数ごとに異なるchain同士の変数を交換し、脱出を容易にしつつ収束の効率を高めてもいる<ref name=松井2021/><ref name=仲田2006>{{Cite web|和書|url= http://www2.tba.t-com.ne.jp/nakada/takashi/bayes/idea.html |title= Bayes 法（ベイズ法）の原理 |author=仲田崇志 |date=2006-06-04 |accessdate=2021-10-29}}</ref>。

== 特徴 ==
ベイズ法は複雑な進化モデルを仮定しても最尤法よりも計算時間を必要としないという長所がある<ref name=三中2009/>。ただし、近隣結合法や最大節約法と比較すれば時間を要する。また、事後確率が過大評価される点も懸念されている<ref>{{Cite journal|和書|title=微生物の系統樹，どう描くの？|author1=飯野隆夫 |author2=伊藤隆 |year=2013 |volume=91 |issue=10 |pages=576-581 |journal=生物工学会誌 |publisher=[[日本生物工学会]] |url=https://www.sbj.or.jp/wp-content/uploads/file/sbj/9110/9110_yomoyama.pdf}}</ref>。

== 出典 ==
{{Reflist}}

== 関連項目 ==
*[[距離行列法]]
**[[非加重結合法]]（UPGMA）
**[[近隣結合法]]（NJ法）
*[[最大節約法]]（MP法）
*[[最尤推定|最尤法]]（ML法）

== 外部リンク ==
* [http://www2.tba.t-com.ne.jp/nakada/takashi/bayes/idea.html Bayes 法（ベイズ法）の原理] - 仲田崇志

{{確率論}}
{{デフォルトソート:へいすほう}}
[[Category:系統学]]
[[Category:分類学 (生物学)]]
[[Category:進化]]
[[Category:ベイズ統計]]