
確率論と統計学において、マルコフ連鎖またはマルコフ過程とは、各事象の確率が前の事象で達成された状態のみに依存する一連の事象を記述する確率過程である。非公式には、「次に何が起こるかは、現在の状況のみに依存する」と考えることができる。連鎖が離散的な時間ステップで状態を遷移する可算無限のシーケンスは、離散時間マルコフ連鎖(DTMC)となる。連続時間過程は、連続時間マルコフ連鎖(CTMC)と呼ばれる。マルコフ過程は、ロシアの数学者アンドレイ・マルコフにちなんで名付けられた。
マルコフ連鎖は、現実世界のプロセスの統計モデルとして多くの応用例があります。 [ 1 ]これらは、複雑な確率分布からのサンプリングをシミュレートするために使用されるマルコフ連鎖モンテカルロ法として知られる一般的な確率シミュレーション手法の基礎を提供し、ベイズ統計、生物学、化学、経済学、金融、情報理論、物理学、信号処理、音声処理などの分野で応用されています。[ 1 ] [ 2 ] [ 3 ]
形容詞のマルコフ的およびマルコフは、マルコフ過程に関連するものを説明するために使用されます。[ 4 ]

マルコフ過程は、マルコフ性(時には「記憶性がない」とも呼ばれる)を満たす確率過程です。簡単に言うと、現在の状態のみに基づいて将来の結果を予測できる過程であり、最も重要なのは、そのような予測は過程の完全な履歴を知っている場合と同じくらい正確であるということです。[ 5 ]言い換えれば、システムの現在の状態を条件として、その将来の状態と過去の状態は独立しています。
マルコフ連鎖は、離散状態空間または離散インデックス集合(多くの場合、時間を表す)を持つマルコフ過程の一種ですが、マルコフ連鎖の正確な定義は様々です。[ 6 ]例えば、マルコフ連鎖は、可算状態空間を持つ離散時間または連続時間のマルコフ過程として定義されるのが一般的です(したがって、時間の性質に関係なく)[ 7 ] [ 8 ] [ 9 ] [ 10 ]が、マルコフ連鎖は、可算状態空間または連続状態空間の離散時間を持つものとして定義されるのも一般的です(したがって、状態空間に関係なく)。[ 6 ]
システムの状態空間と時間パラメータのインデックスを指定する必要があります。以下の表は、離散時間と連続時間の両方について、状態空間の一般性の異なるレベルにおけるマルコフ過程のさまざまなインスタンスの概要を示しています。
マルコフ過程の特殊なケースを示す用語の使用法については、文献で明確な合意が得られていないことに注意してください。通常、「マルコフ連鎖」という用語は、離散的な時間の集合を持つプロセス、つまり離散時間マルコフ連鎖(DTMC)[ 11 ]に予約されていますが、一部の著者は、明示的な言及なしに連続時間マルコフ連鎖(CTMC)を指すために「マルコフ過程」という用語を使用しています。[ 12 ] [ 13 ] [ 14 ]さらに、マルコフ過程の他の拡張もそのように呼ばれていますが、必ずしもこれら 4 つのカテゴリのいずれにも属しません(マルコフ モデルを参照)。さらに、時間インデックスは必ずしも実数値である必要はありません。状態空間と同様に、他の数学的構成要素を持つインデックス集合を通過するプロセスも考えられます。一般的な状態空間連続時間マルコフ連鎖は、非常に一般的であるため、特定の用語がないことに注意してください。
時間パラメータは通常離散的ですが、マルコフ連鎖の状態空間には一般的に合意された制約はありません。この用語は任意の状態空間上のプロセスを指す場合があります。[ 15 ]ただし、マルコフ連鎖の多くのアプリケーションでは、より単純な統計分析が可能な有限または可算無限の状態空間が使用されます。時間インデックスと状態空間パラメータの他に、多くのバリエーション、拡張、および一般化があります(バリエーションを参照)。簡潔にするため、特に断りのない限り、この記事の大部分は離散時間、離散状態空間の場合に焦点を当てています。
システムの状態変化は遷移と呼ばれます。様々な状態変化に伴う確率は遷移確率と呼ばれます。プロセスは、状態空間、特定の遷移の確率を表す遷移行列、および状態空間全体にわたる初期状態(または初期分布)によって特徴付けられます。慣例として、プロセスの定義にはすべての可能な状態と遷移が含まれていると仮定するため、常に次の状態が存在し、プロセスは終了しません。
離散時間ランダムプロセスとは、各ステップで特定の状態にあるシステムであり、その状態がステップ間でランダムに変化するプロセスです。ステップはしばしば時間の瞬間として考えられますが、物理的な距離やその他の離散的な測定値を指す場合もあります。形式的には、ステップは整数または自然数であり、ランダムプロセスはこれらを状態にマッピングしたものです。マルコフ性とは、次のステップ(そして実際には将来のすべてのステップ)におけるシステムの条件付き確率分布は、システムの現在の状態のみに依存し、前のステップにおけるシステムの状態には依存しないという性質です。
システムはランダムに変化するため、マルコフ連鎖の将来のある時点における状態を確実に予測することは一般的に不可能です。しかし、システムの将来に関する統計的特性は予測可能です。多くの応用分野では、これらの統計的特性が重要となります。
アンドレイ・マルコフは20世紀初頭にマルコフ過程を研究し、1906年にこのテーマに関する最初の論文を発表しました。[ 16 ] [ 17 ] [ 18 ]連続時間におけるマルコフ過程は、20世紀初頭の彼の研究よりもずっと前に、ポアソン過程の形で発見されていました。[ 19 ] [ 20 ] [ 21 ]マルコフは、独立ランダムシーケンスの拡張の研究に興味を持ち、独立性が大数の弱法則が成り立つために必要だと主張したパベル・ネクラソフとの意見の相違がきっかけとなりました。 [ 22 ] 1906年に発表されたマルコフ連鎖に関する最初の論文で、マルコフは、特定の条件下ではマルコフ連鎖の平均結果が固定値のベクトルに収束することを示し、独立性の仮定なしに大数の弱法則を証明しました。[ 16 ] [ 17 ] [ 18 ]これは、そのような数学法則が成り立つための要件として一般的に考えられていたものです。[ 18 ]マルコフは後にマルコフ連鎖を用いてアレクサンドル・プーシキン作の『エヴゲーニー・オネーギン』における母音の分布を研究し、そのような連鎖の中心極限定理を証明した。[ 16 ]
1912年、アンリ・ポアンカレはカードシャッフルの研究を目的として有限群上のマルコフ連鎖を研究した。マルコフ連鎖のその他の初期の用途としては、 1907年にポールとタチアナ・エーレンフェストによって導入された拡散モデルや、マルコフの研究に先立つ1873年にフランシス・ゴルトンとヘンリー・ウィリアム・ワトソンによって導入された分岐過程などがある。 [ 16 ] [ 17 ]ゴルトンとワトソンの研究の後、彼らの分岐過程は、約30年前にイレーネ=ジュール・ビエネメによって独立して発見され研究されていたことが後に明らかになった。[ 23 ] 1928年からモーリス・フレシェはマルコフ連鎖に興味を持ち始め、最終的に1938年にマルコフ連鎖に関する詳細な研究を発表するに至った。[ 16 ] [ 24 ]
アンドレイ・コルモゴロフは1931年の論文で、連続時間マルコフ過程の初期理論の大部分を開発した。[ 25 ] [ 26 ]コルモゴロフは、ルイ・バシュリエの1900年の株式市場の変動に関する研究と、ノーバート・ウィーナーのアインシュタインのブラウン運動モデルに関する研究に部分的に影響を受けていた。[ 25 ] [ 27 ]彼は拡散過程として知られる特定のマルコフ過程を導入して研究し、その過程を記述する一連の微分方程式を導出した。[ 25 ] [ 28 ]コルモゴロフの研究とは独立して、シドニー・チャップマンは1928年の論文で、ブラウン運動を研究しながら、コルモゴロフよりも数学的に厳密ではない方法で、現在チャップマン・コルモゴロフ方程式と呼ばれる方程式を導出した。 [ 29 ]これらの微分方程式は現在、コルモゴロフ方程式[ 30 ]またはコルモゴロフ・チャップマン方程式[ 31 ]と呼ばれています。マルコフ過程の基礎に大きく貢献した他の数学者には、1930 年代からウィリアム・フェラー、そして後に1950 年代からユージン・ディンキン[ 26 ]がいます。
25セント相当の硬貨(25セント硬貨)が5枚、10セント相当の硬貨(10セント硬貨)が5枚、5セント相当の硬貨(5セント硬貨)が5枚入った小銭入れがあるとします。小銭入れから硬貨を1枚ずつランダムに取り出し、テーブルの上に置きます。n回の抽選後にテーブルに置かれたコインの合計値を表し、すると、次のシーケンスマルコフ過程ではない。
なぜそうなるのかを理解するために、最初の6回の抽選で、5枚の5セント硬貨と1枚の25セント硬貨がすべて引き当てられたと仮定してみましょう。したがって. もし私たちがただ知っているのではなくしかし、以前の値も考慮すれば、どのコインが引き出されたかを判断でき、次のコインはニッケルではないことがわかっているので、確率1で。しかし、以前の値がわからない場合は、値のみに基づいて4枚の10セント硬貨と2枚の5セント硬貨を引いたと推測できますが、その場合、次にもう1枚の5セント硬貨を引くことは確かに可能です。したがって、以前の価値観に関する知識によって影響を受ける。
しかし、このシナリオをマルコフ過程としてモデル化することは可能です。定義する代わりにテーブル上のコインの合計価値を表すために、次のように定義できます。テーブル上のさまざまなコインの種類の数を表す。たとえば、は、6回の1枚ずつの抽選の後、テーブル上に25セント硬貨が1枚、10セント硬貨が0枚、5セント硬貨が5枚ある状態を表すものとして定義できる。この新しいモデルは次のように表すことができる。考えられる状態。各状態は、テーブル上にある各種類のコインの枚数(0~5枚)を表します。(これらの状態すべてが6回の抽選で到達できるわけではありません。)
最初の抽選で状態が達成する確率今は例えば、州不可能です。2回目の抽選後、3回目の抽選はこれまでにどのコインが引かれたかに依存しますが、最初の状態のために引かれたコインだけに依存するわけではありません(確率的に重要な情報がその後シナリオに追加されたため)。このようにして、州は、州。
離散時間マルコフ連鎖とは、マルコフ性を持つ確率変数X 1、X 2、X 3、 ... のシーケンスであり、次の状態へ移行する確率は現在の状態のみに依存し、前の状態には依存しないという性質を持つ。
X iの取りうる値は、チェーンの状態空間と呼ばれる可算集合Sを形成する。
状態空間が有限の場合、遷移確率分布は遷移行列と呼ばれる行列で表すことができ、Pの ( i , j ) 番目の要素は、
Pの各行の合計が1であり、すべての要素が非負であるため、 Pは右確率行列である。
定常分布πは、(行)ベクトルであり、その要素は非負で合計が1であり、遷移行列Pの作用によって変化しないため、次のように定義されます。
この定義を固有ベクトルの定義と比較すると、2つの概念は関連しており、
正規化された()遷移行列Pの固有値が 1 である左固有ベクトルeの倍数。単位固有ベクトルが複数ある場合は、対応する定常状態の重み付き和も定常状態になります。しかし、マルコフ連鎖の場合、通常は、ある初期分布に対する分布のシーケンスの極限である定常状態の方が関心があります。
定常分布の値Pの状態空間に関連付けられており、その固有ベクトルは相対的な比率が保持されています。πの成分は正であり、それらの合計が1であるという制約は次のように書き換えることができます。πと成分がすべて1であるベクトルの内積が1であること、そしてπが単体上にあることがわかります。
マルコフ連鎖が時間的に均質である場合、遷移行列P は各ステップ後も同じであるため、kステップ遷移確率は遷移行列P kのk乗として計算できます。
マルコフ連鎖が既約かつ非周期的である場合、一意の定常分布πが存在する。[ 40 ]さらに、この場合、P kは各行が定常分布πであるランク1行列に収束する。
ここで、1はすべての要素が 1 に等しい列ベクトルです。これはペロン・フロベニウスの定理で述べられています。どのような手段であれ、が見つかれば、後述するように、任意の初期分布に対して、問題のマルコフ連鎖の定常分布を容易に決定できます。
いくつかの確率行列Pに対して、極限定常分布は存在するが、定常分布は存在しない。この例がそれを示している。
(この例は周期的なマルコフ連鎖を示しています。)
考慮すべき特殊なケースが多数あるため、この極限が存在する場合、それを見つけるプロセスは長い作業になる可能性があります。しかし、この極限を見つけるのに役立つ多くの手法があります。Pをn × n行列とし、次のように定義します。
常に真実である
両辺からQを引いて因数分解すると、
ここで、I nはnサイズの単位行列、0 n , nはn × nサイズの零行列です。確率行列を掛け合わせると必ず別の確率行列が得られるため、Q は確率行列でなければなりません(上記の定義を参照)。Qを解くには、上記の行列方程式とQが確率行列であるという事実を使用するだけで十分な場合もあります。P の各行の合計が 1 であるという事実を含めると、n 個の未知数を決定するための方程式がn+1 個あるため、一方ではQの 1 つの行を選択してその各要素を 1 で置き換え、他方では対応する要素 (同じ列の要素) をベクトル0に置き換え、次にこの後者のベクトルを変換された前の行列の逆行列で左から乗算してQ を見つける方が計算が容易です。
これを行う方法の 1 つは次のとおりです。まず、行列Aの右端の列をすべて 1 に置き換えた行列を返す関数f ( A ) を定義します。 [ f ( P − I n )] −1が存在する場合、 [ 41 ] [ 40 ]
注目すべき点の1つは、Pの主対角線上に1に等しい要素P i , iがあり、 i番目の行または列がそれ以外はすべて0で埋められている場合、その行または列は後続のすべてのべき乗P kで変更されないということです。したがって、Qのi番目の行または列には、Pと同じ位置に1と0が配置されます。
前述の通り、方程式から(存在する場合)定常(または安定状態)分布πは、行確率行列Pの左固有ベクトルである。次に、Pが対角化可能であるか、または同等にPがn個の線形独立な固有ベクトルを持つと仮定すると、収束速度は次のように詳述される。(非対角化可能、つまり欠陥のある行列の場合、 Pのジョルダン標準形から始めて、同様の方法でもう少し複雑な一連の議論を進めることができる。[ 42 ])
U を固有ベクトルの行列 (各固有ベクトルは L2 ノルムが 1 に正規化されている) とし、各列はPの左固有ベクトルとする。また、 Σ をPの左固有値の対角行列、すなわちΣ = diag( λ 1 , λ 2 , λ 3 ,..., λ n ) とする。すると、固有値分解により、
固有値を以下のように列挙する。
Pは行確率行列であるため、その最大の左固有値は1です。一意の定常分布が存在する場合、最大の固有値と対応する固有ベクトルも一意になります(上記の定常分布方程式を満たすπは他に存在しないため)。u iをU行列のi番目の列とします。つまり、u iはλ iに対応するPの左固有ベクトルです。また、xを有効な確率分布を表す長さnの行ベクトルとします。固有ベクトルu iは私たちは書くことができます
xにPを右から掛け、この操作を結果に対して続けると、最終的に定常分布πが得られます。言い換えると、π = a 1 u 1 ← xPP ... P = xP k as k → ∞ となります。つまり、
πはu 1(L2ノルムで正規化)と平行であり、π(k)は確率ベクトルであるため、 k →∞のとき、 π(k)はλ2 / λ1のオーダーの指数関数的な速度で1u 1 = πに近づきます。これは次の理由によります。したがって、λ 2 / λ 1が支配的な項です。比率が小さいほど、収束は速くなります。[ 43 ]状態分布πのランダムノイズも、定常分布へのこの収束を速めることができます。[ 44 ]
連続時間マルコフ連鎖は、有限または可算な状態空間S、状態空間と同じ次元を持つ遷移率行列Q、および状態空間上で定義された初期確率分布によって定義されます。 i ≠ jの場合、要素q ijは非負であり、状態iから状態jへのプロセスの遷移率を表します。要素q iiは、遷移率行列の各行の合計がゼロになるように選択されますが、(離散) マルコフ連鎖における確率遷移行列の行の合計はすべて 1 になります。
このプロセスには3つの同等の定義がある。[ 45 ]

させてを時刻tにおけるプロセスの状態を表す確率変数とし、時刻tにおいてプロセスが状態iにあると仮定する。すると、、以前の値とは無関係である、そしてすべてのjおよびすべてのtに対してh → 0となる。 どこは、小文字のoを用いたクロネッカーデルタです。これは、 iからjへの遷移がどれだけ速く起こるかを測定するものと見なすことができる。
プロセスのn番目のジャンプを記述する離散時間マルコフ連鎖Y n と、各状態における保持時間を記述する変数S 1、S 2、S 3 、 ... を定義します。ここでS iはレートパラメータ − q Y i Y iの指数分布に従います。
任意の値n = 0, 1, 2, 3, ... およびこのnの値までのインデックス付き時刻t 0 , t 1 , t 2 , ... およびこれらの時刻で記録されたすべての状態i 0 , i 1 , i 2 , i 3 , ... に対して、以下が成り立つ。
ここでp ijは順方向方程式(1 階微分方程式)の解である。
初期条件P(0)は単位行列である。
「局所的に相互作用するマルコフ連鎖」とは、他のマルコフ連鎖の状態を考慮した進化を持つマルコフ連鎖のことです。これは、状態空間が(デカルト)積形式を持つ状況に対応します。相互作用する粒子システムと確率的セルオートマトン(確率的セルオートマトン)を参照してください。例えば、マルコフ過程の相互作用[ 46 ] または[ 47 ]を参照してください。
有限状態空間を持つ離散時間マルコフ連鎖に関する多くの結果は、ハリス連鎖を通して非可算状態空間を持つ連鎖に一般化することができる。
マルコフ連鎖モンテカルロ法におけるマルコフ連鎖の使用は、プロセスが連続状態空間に従う場合を対象としています。
一般状態空間における連続時間マルコフ過程の定義は、上記よりも技術的な内容となる。
連続時間マルコフ過程は、濾過に適した確率過程である。局所的にコンパクトなポーランド空間における値(例えば、後者は基本的に、それらは規則的であり、簡単に言えば「行儀よく」振る舞うということです。は、マルコフ性を満たす場合、マルコフ過程と呼ばれます。すなわち、すべての に対して が成り立ちます。そして[ 5 ]
さらに、すべてのに対して弱マルコフ性を満たす場合、それは時間的に均質であると呼ばれる。:
機能は、いわゆる遷移関数である。そしてプロセスの遷移半群。遷移関数は、有限状態空間の設定で使用される遷移行列の一般化である。
より抽象的な方法では、マルコフ過程は逆の方向にも定義または構築できます。遷移半群である、すなわち、
どこディラック測度は、 そして。 それから自然濾過に関して均質なマルコフ過程であるすべての場合、基礎となる確率尺度満たす
または、確率尺度がない場合が指定されており、上記の式は尺度を定義します。の上そのプロセスの下で開始これは構成上マルコフ過程である。
言い換えれば、マルコフ過程は確率過程として定義できる。フィルタリングされた確率空間上、または間接的に遷移半群(すなわち、プロセスの遷移確率)の観点から、確率空間を誘導し、マルコフ性を持つ。
2つの状態は、正の確率を持つ遷移のシーケンスによって互いに到達可能な場合、互いに通信可能であると言われます。これは同値関係であり、通信可能なクラスの集合を生成します。クラスから出る確率がゼロの場合、そのクラスは閉じていると言われます。マルコフ連鎖は、通信可能なクラスが1つだけ存在する場合、すなわち状態空間が存在する場合に、既約であると言われます。
状態i の周期がkであるとは、k が、状態iから出発して状態iに到達できる遷移回数の最大公約数である場合をいう。すなわち、
状態が周期的であるのは、; さもないとそしてその状態は非周期的である。
状態iは、 iから開始して、連鎖がiに決して戻らない確率がゼロでない場合、一時的であると言われます。そうでない場合は、再帰的(または持続的)と呼ばれます。 [ 48 ]再帰的状態iの場合、平均到達時間は次のように定義されます。
状態iが正の再帰的である場合それ以外の場合は有限でヌル再帰である。周期性、一時性、再帰性、正の再帰性、ヌル再帰性はクラス特性である。つまり、ある状態がその特性を持つ場合、その状態と通信するクラスのすべての状態がその特性を持つ。[ 49 ]
状態iから外部への遷移がない場合、その状態は吸収状態と呼ばれます。
周期性はクラス特性であるため、マルコフ連鎖が既約であれば、そのすべての状態は同じ周期を持つ。特に、1つの状態が非周期であれば、マルコフ連鎖全体が非周期となる。[ 50 ]
有限マルコフ連鎖が既約である場合、すべての状態は正の再帰性を持ち、次の式で与えられる一意の定常分布を持つ。。
状態iは、非周期的かつ正の再帰性を持つ場合、エルゴード的であると言われる。言い換えれば、状態iは、再帰性があり、周期が 1 であり、平均再帰時間が有限である場合にエルゴード的である。
既約マルコフ連鎖のすべての状態がエルゴード的である場合、その連鎖はエルゴード的であると言われる。言い換えれば、ある整数が存在する。すべてのエントリ肯定的である。
有限状態の既約マルコフ連鎖は、非周期状態を持つ場合、エルゴード的であることが示される。
複数の状態を持ち、各状態から出る遷移が1つだけのマルコフ連鎖は、既約でも非周期的でもないため、エルゴード的ではない。
一部の著者は、周期的なものも含め、既約で正の再帰的なマルコフ連鎖をエルゴード的と呼ぶ。[ 51 ]実際、既約なマルコフ連鎖だけが、エルゴード理論に従って定義されるエルゴード過程に対応する。[ 52 ]
ある整数が存在する場合、一部の著者は行列をプリミティブと呼ぶ。すべてのエントリ正の値です。[ 53 ]一部の著者はこれを規則的と呼んでいます。[ 54 ]
正則行列の原始性指数、または指数は、最小値である。すべてのエントリは正である。指数は純粋にグラフ理論的な性質であり、各エントリがはゼロまたは正であり、したがって有向グラフ上で見つけることができる。隣接行列として。
有限個の状態がある場合、指数に関する組み合わせ論的な結果がいくつかあります。状態の数を とすると、[ 55 ]
マルコフ連鎖が定常分布を持つ場合、それは測度保存力学系に変換できる。確率空間を次のように定義する。、 どこはマルコフ連鎖のすべての状態の集合です。確率空間上のシグマ代数は円筒集合によって生成されるものとします。確率測度は定常分布とマルコフ連鎖遷移によって生成されるものとします。シフトオペレーターになる:同様に、次のような動的システムを構築できます。代わりに。[ 57 ]
有限状態空間を持つ既約マルコフ連鎖は一意の定常分布を持つため、上記の構成は既約マルコフ連鎖に対しては曖昧さがない。
エルゴード理論では、測度保存力学系は、任意の可測部分集合がそのため暗示するまたは(空集合を除いて)
用語に一貫性がない。すべての状態において厳密に正である定常分布を持つマルコフ連鎖が与えられた場合、対応する測度保存力学系がエルゴード的であれば、マルコフ連鎖は既約である。[ 52 ]
場合によっては、一見非マルコフ過程であっても、「現在」と「未来」の状態の概念を拡張することで、マルコフ的な表現を持つことがあります。例えば、X を非マルコフ過程とします。このとき、 Yの各状態がXの状態の時間間隔を表すような過程Yを定義します。数学的には、これは次の形式になります。
Yがマルコフ性を持つ場合、YはXのマルコフ表現である。
マルコフ表現を持つ非マルコフ過程の例として、次数が1より大きい自己回帰時系列が挙げられる。 [ 58 ]
到達時間とは、与えられた状態群から開始して、連鎖が特定の状態または状態群に到達するまでの時間のことである。このような時間間隔の分布は位相型分布となる。最も単純な分布は、単一の指数分布遷移の分布である。
状態のサブセットA ⊆ Sに対して、ヒット時間のベクトルk A (要素 は、状態iから始めて、チェーンがセットA ) のいずれかの状態に入るときの期待値を表します。は[ 59 ]の最小非負解です。
一般的なマルコフ過程の場合連続時間(CTMCまたは一般状態空間を持つプロセス)では、逆プロセス一定の時間からは再びマルコフ過程です。これはマルコフ性から直接導かれます。非公式に言えば、現在が与えられた場合、未来と過去は独立しています。時間反転では、それらの役割が入れ替わるだけです。ただし、逆過程は一般に時間的に均質ではありません。あるランダムな時間に対して、(必ずしも停止時間ではない)停止したプロセス時間的に均質なマルコフ過程である場合、逆の過程再び時間的に均質である。[ 60 ]
もしはCTMCであり、ケリーの補題により順方向プロセスと同じ定常分布を持つ。
マルコフ連鎖は、逆方向の過程が順方向の過程(分布において)と同じである場合に可逆であると言われます。コルモゴロフの基準によれば、マルコフ連鎖が可逆であるための必要十分条件は、閉ループ内の遷移率の積が両方向で同じであることとされています。
エルゴード的な連続時間マルコフ連鎖Qの定常確率分布πを求める方法の一つは、まずその埋め込みマルコフ連鎖 (EMC) を求めることです。厳密に言えば、EMC は離散時間マルコフ連鎖であり、ジャンプ過程と呼ばれることもあります。EMC の 1 ステップ遷移確率行列Sの各要素はs ijで表され、状態iから状態jへの遷移の条件付き確率を表します。これらの条件付き確率は次のように求めることができます。
このことから、Sは次のように書ける。
ここで、Iは単位行列であり、diag( Q ) は行列Qから主対角線を選択し、他のすべての要素をゼロに設定することによって形成される対角行列です。
定常確率分布ベクトルを見つけるには、次にそのため
と行ベクトルであり、すべての要素が0より大きく、= 1. これから、π は次のように求められます。
(Qが周期的でなくても、Sは周期的である可能性がある。πが求められたら、それを単位ベクトルに正規化する必要がある。)
連続時間マルコフ連鎖から派生するもう1つの離散時間プロセスは、δスケルトンです。これは、δ単位の時間間隔でX ( t )を観測することによって形成される(離散時間)マルコフ連鎖です。確率変数X (0)、X (δ)、X (2δ)、...は、δスケルトンが訪れる状態のシーケンスを表します。
マルコフモデルは、変化するシステムをモデル化するために使用されます。マルコフ連鎖を一般化した主なモデルは4種類あり、各連続状態が観測可能かどうか、また観測に基づいてシステムを調整するかどうかによって分類されます。
ベルヌーイスキームは、遷移確率行列の行がすべて同一であるマルコフ連鎖の特殊なケースであり、これは次の状態が現在の状態から独立しているだけでなく、過去の状態からも独立していることを意味します。可能な状態が2つしかないベルヌーイスキームは、ベルヌーイ過程として知られています。
Note, however, by the Ornstein isomorphism theorem, that every aperiodic and irreducible Markov chain is isomorphic to a Bernoulli scheme;[61] thus, one might equally claim that Markov chains are a "special case" of Bernoulli schemes. The isomorphism generally requires a complicated recoding. The isomorphism theorem is even a bit stronger: it states that anystationary stochastic process is isomorphic to a Bernoulli scheme; the Markov chain is just one such example.
When the Markov matrix is replaced by the adjacency matrix of a finite graph, the resulting shift is termed a topological Markov chain or a subshift of finite type.[61] A Markov matrix that is compatible with the adjacency matrix can then provide a measure on the subshift. Many chaotic dynamical systems are isomorphic to topological Markov chains; examples include diffeomorphisms of closed manifolds, the Prouhet–Thue–Morse system, the Chacon system, sofic systems, context-free systems and block-coding systems.[61]
Markov chains have been employed in a wide range of topics across the natural and social sciences, and in technological applications.
Markovian systems appear extensively in thermodynamics and statistical mechanics, whenever probabilities are used to represent unknown or unmodelled details of the system, if it can be assumed that the dynamics are time-invariant, and that no relevant history need be considered which is not already included in the state description.[62][63] For example, a thermodynamic state operates under a probability distribution that is difficult or expensive to acquire. Therefore, Markov Chain Monte Carlo method can be used to draw samples randomly from a black-box to approximate the probability distribution of attributes over a range of objects.[63]
Markov chains are used in lattice QCD simulations.[64]
反応ネットワークは、複数の反応と化学種を含む化学システムです。このようなネットワークの最も単純な確率モデルでは、システムを連続時間マルコフ連鎖として扱い、状態は各種の分子の数であり、反応は連鎖の可能な遷移としてモデル化されます。[ 65 ]マルコフ連鎖と連続時間マルコフ過程は、物理システムがマルコフ特性に非常に近い場合に化学で役立ちます。たとえば、溶液中に多数のn個の分子が状態 A にあり、それぞれが一定の平均速度で状態 B に化学反応を起こすことができると想像してください。おそらくその分子は酵素であり、状態はそれがどのように折り畳まれているかを示しています。任意の単一の酵素の状態はマルコフ連鎖に従い、分子は基本的に互いに独立しているため、ある時点で状態 A または B にある分子の数は、特定の分子がその状態にある確率のn倍になります。
酵素活性の古典的なモデルであるミカエリス・メンテン速度論は、各時間ステップで反応が何らかの方向に進行するマルコフ連鎖として見なすことができる。ミカエリス・メンテンは比較的単純であるが、はるかに複雑な反応ネットワークもマルコフ連鎖でモデル化することができる。[ 66 ]
マルコフ連鎖に基づくアルゴリズムも、医薬品や天然物などの目的の化合物クラスに向けて、インシリコでの化学物質のフラグメントベースの成長を集中させるために使用されました。 [ 67 ]分子が成長するにつれて、新生分子からフラグメントが「現在の」状態として選択されます。それは過去を認識していません(つまり、すでに何に結合しているかを認識していません)。その後、フラグメントが結合すると、次の状態に遷移します。遷移確率は、実際の化合物クラスのデータベースでトレーニングされます。[ 68 ]
また、共重合体の成長(および組成)は、マルコフ連鎖を用いてモデル化できる。成長中のポリマー鎖を構成するモノマーの反応性比に基づいて、鎖の組成(例えば、モノマーが交互に付加される傾向があるか、あるいは同じモノマーが連続して付加される傾向があるかなど)を計算できる。立体効果により、一部のポリマー鎖の成長には二次マルコフ効果も関与する可能性がある。
同様に、一部のエピタキシャル超格子酸化物材料の結晶化と成長はマルコフ連鎖によって正確に記述できることが示唆されている。 [ 69 ]
マルコフ連鎖は生物学の様々な分野で用いられています。代表的な例としては以下のようなものがあります。
マルコフ連鎖は情報処理のあらゆる場面で使用されています。クロード・シャノンの有名な1948年の論文「通信の数学的理論」は、情報理論の分野を一挙に創始したもので、自然言語(英語など)のテキストをエルゴード的マルコフ過程によって生成されるものとしてモデル化することでエントロピーの概念を導入することから始まります。この過程では、各文字が前の文字に統計的に依存する可能性があります。 [ 72 ]このような理想化されたモデルは、システムの統計的規則性の多くを捉えることができます。システムの完全な構造を完全に記述しなくても、このような信号モデルは、算術符号化などのエントロピー符号化技術によって非常に効果的なデータ圧縮を可能にします。また、効果的な状態推定とパターン認識も可能にします。マルコフ連鎖は強化学習でも重要な役割を果たします。
マルコフ連鎖は隠れマルコフモデルの基礎でもあり、電話ネットワーク(誤り訂正にビタビアルゴリズムを使用)、音声認識、バイオインフォマティクス(再配置検出など[ 73 ])といった多様な分野で重要なツールとなっています。
LZMAロスレスデータ圧縮アルゴリズムは、マルコフ連鎖とレンペル・ジブ圧縮を組み合わせることで、非常に高い圧縮率を実現します。
マルコフ連鎖は、待ち行列の解析的処理(待ち行列理論)の基礎となるものです。アグナー・クラルップ・エルランが1917年にこの分野を創始しました。[ 74 ]このため、マルコフ連鎖は、メッセージが限られたリソース(帯域幅など)をめぐって競合することが多い電気通信ネットワークのパフォーマンスを最適化する上で非常に重要となります。[ 75 ]
数多くの待ち行列モデルは、連続時間マルコフ連鎖を使用します。例えば、M/M/1待ち行列は、非負整数上の連続時間マルコフ連鎖であり、iからi + 1への上方遷移はポアソン過程に従ってレートλで発生し、ジョブの到着を表します。一方、 iからi - 1への遷移( i > 1の場合)はレートμで発生し(ジョブのサービス時間は指数分布に従います)、待ち行列からの完了したサービス(出発)を表します。

Googleが使用するウェブページのPageRankは、マルコフ連鎖によって定義されます。[ 76 ] [ 77 ] [ 78 ]これは、ページが表示される確率です。すべての(既知の)ウェブページ上の以下のマルコフ連鎖の定常分布において。は既知のウェブページの数であり、ページはもっているそこから発信リンクがあると、遷移確率がリンクされているすべてのページとリンクされていないすべてのページについて。パラメーター約0.85とみなされている。[ 79 ]
マルコフモデルは、ユーザーのウェブナビゲーション行動の分析にも用いられています。特定のウェブサイトにおけるユーザーのリンク遷移は、一次または二次マルコフモデルを用いてモデル化することができ、将来のナビゲーション行動を予測したり、個々のユーザーに合わせてウェブページをパーソナライズしたりするために活用できます。
マルコフ連鎖法は、マルコフ連鎖モンテカルロ法(MCMC)と呼ばれる手法を用いて、非常に複雑な確率分布を正確に反映する乱数列を生成する上で、非常に重要な役割を果たすようになりました。近年、この手法はベイズ推論法の実現可能性を大きく向上させ、幅広い事後分布をシミュレーションし、そのパラメータを数値的に求めることが可能になりました。
マルコフ連鎖は、金融や経済学において、所得の分布、企業の規模分布、資産価格、市場暴落など、さまざまな現象をモデル化するために使用されています。DGチャンパーノウンは1953 年に所得分布のマルコフ連鎖モデルを構築しました。 [ 80 ]ハーバート A. サイモンと共著者のチャールズ ボニーニは、マルコフ連鎖モデルを使用して、企業の規模の定常ユール分布を導出しました。[ 81 ]ルイ バシェリエは、株価がランダム ウォークに従うことを最初に観察しました。[ 82 ]ランダム ウォークは後に効率的市場仮説を支持する証拠と見なされ、ランダム ウォーク モデルは 1960 年代の文献で人気がありました。[ 83 ]景気循環のレジーム スイッチング モデルは、マルコフ連鎖を使用して高 GDP 成長期間と低 GDP 成長期間 (または、経済拡大と景気後退) の切り替えをモデル化したジェームズ D. ハミルトン(1989) によって普及しました。[ 84 ]より最近の例としては、ローラン・E・カルヴェとアドレー・J・フィッシャーによるマルコフスイッチング多重フラクタルモデルがあり、これは以前のレジームスイッチングモデルの利便性を基に構築されています。[ 85 ] [ 86 ]これは、任意に大きなマルコフ連鎖を使用して、資産収益の変動レベルを制御します。
動学的マクロ経済学では、マルコフ連鎖が多用される。例えば、一般均衡設定において、株式価格を外生的にモデル化するためにマルコフ連鎖が用いられる。 [ 87 ]
マルコフ連鎖は一般的に、現在の構造的構成が将来の結果を左右する経路依存的な議論を記述する際に用いられる。例としては、カール・マルクスの『資本論』に由来する、経済発展と資本主義の台頭を結びつける考え方の再定式化が挙げられる。現在の研究では、マルコフ連鎖を用いて、ある国が特定の経済発展レベルに達すると、中間層の規模、都市部と農村部の居住比率、政治動員率などの構造的要因の構成が、権威主義体制から民主主義体制への移行確率を高めることをモデル化するのが一般的である。[ 89 ]
マルコフ連鎖は、アルゴリズムによる音楽作曲、特にCsound、Max、SuperColliderなどのソフトウェアで利用されています。一次連鎖では、システムの各状態が音符またはピッチ値となり、各音符の確率ベクトルが構築され、遷移確率行列が完成します(下記参照)。遷移行列の重みに基づいて出力音符値を生成するアルゴリズムが構築され、これはMIDI音符値、周波数(Hz)、またはその他の任意の望ましい指標となる可能性があります。[ 90 ]
2 番目の表に示すように、現在の状態と前の状態の両方を考慮することで、2 次マルコフ連鎖を導入できます。より高次のn次連鎖は、特定の音符をまとめて「グループ化」する傾向があり、時折他のパターンやシーケンスに「分岐」します。これらの高次連鎖は、1 次システムによって生成される「目的のない彷徨」ではなく、フレーズ構造の感覚を持つ結果を生成する傾向があります。 [ 91 ]
マルコフ連鎖は、クセナキスの『アナロジケ A』と『B』のように構造的に使用できます。[ 92 ]マルコフ連鎖は、マルコフモデルを使用して音楽入力にインタラクティブに反応するシステムでも使用されます。[ 93 ]
通常、音楽システムは生成する有限長のシーケンスに特定の制御制約を適用する必要がありますが、制御制約はマルコフモデルと互換性がありません。なぜなら、制御制約はマルコフの限定記憶仮説に違反する長距離依存性を引き起こすからです。この制限を克服するために、新しいアプローチが提案されています。[ 94 ]
マルコフ連鎖は、多くのギャンブルゲームをモデル化するために使用できます。例えば、子供向けのゲームである「ヘビとはしご」や「ハイホー!チェリーオー」は、マルコフ連鎖によって正確に表現されます。各ターンで、プレイヤーは特定の状態(特定のマス)からスタートし、そこから特定の他の状態(マス)に移動できる確率が固定されています。
マルコフ連鎖モデルは1960年以来、高度な野球分析に用いられてきたが、その使用は依然として稀である。野球の試合の各イニングは、走者数とアウト数を考慮すると、マルコフ連鎖の状態に適合する。どの打席においても、アウト数と走者の位置の組み合わせは24通りある。マーク・パンキンは、マルコフ連鎖モデルが個々の選手とチームの両方で生み出された得点を評価するために使用できることを示している。[ 95 ]また、彼はさまざまな戦略とプレー条件についても論じている。マルコフ連鎖モデルがバントや盗塁 などの試合状況の統計を分析するためにどのように使用されてきたか、また天然芝と人工芝でのプレーの違いについても論じている。[ 96 ]
マルコフ過程は、サンプル文書から表面上は本物そっくりのテキストを生成するためにも使用できます。マルコフ過程は、さまざまな娯楽用「パロディ生成」ソフトウェアで使用されています( dissociated press、Jeff Harrison、[ 97 ] Mark V. Shaney、[ 98 ] [ 99 ] 、およびAcademias Neutroniumを参照)。マルコフ連鎖を使用したオープンソースのテキスト生成ライブラリもいくつか存在します。
{{cite book}}ISBN /日付の不一致(ヘルプ){{citation}}: CS1メンテナンス: ISBNを使用した作業パラメータ (リンク)