
確率論と統計学において、マルコフ連鎖またはマルコフ過程とは、各事象の確率が前の事象で達成された状態のみに依存する一連の事象を記述する確率過程である。非公式には、「次に何が起こるかは、現在の状況のみに依存する」と考えることができる。連鎖が離散的な時間ステップで状態を遷移する可算無限のシーケンスは、離散時間マルコフ連鎖(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)は単位行列である。
"Locally interacting Markov chains" are Markov chains with an evolution that takes into account the state of other Markov chains. This corresponds to the situation when the state space has a (Cartesian-) product form. See interacting particle system and stochastic cellular automata (probabilistic cellular automata). See for instance Interaction of Markov Processes[46] or.[47]
Many results for discrete-time Markov chains with finite state space can be generalized to chains with uncountable state space through Harris chains.
The use of Markov chains in Markov chain Monte Carlo methods covers cases where the process follows a continuous state space.
The definition of Markov processes in continuous time with general state space is more technical than the above.
A continuous-time Markov process is a stochastic process adapted to a filtration with values in a locally compactPolish space (e.g., ). The latter essentially ensures that the conditional expectations of are regular, which, in simple terms, means that they behave "nicely". Then is called a Markov process, if it satisfies the Markov property, i.e., for all and [5]
Moreover, is called time-homogeneous, if it satisfies the weak Markov property for all :
The function is the so-called transition function of and the transition semigroup of the process. Transition functions are generalizations of the transition matrices used in the setting with finite state space.
In a more abstract way, Markov processes can also be defined or constructed the other way around: Let be a transition semigroup, i.e.,
where is the Dirac-measure in , and . Then is a homogeneous Markov process w.r.t. the natural filtration , if for all , the underlying probability measure satisfies
Or, if no probability measure has been specified, the above equation defines a measure on under which the process started in is a Markov process by construction.
In other words, Markov processes can be defined either as stochastic processes on a filtered probability space, or indirectly in terms of a transition semigroup (i.e., the transition probabilities of the process), which induces a probability space under which has the Markov property.
Two states are said to communicate with each other if both are reachable from one another by a sequence of transitions that have positive probability. This is an equivalence relation which yields a set of communicating classes. A class is closed if the probability of leaving the class is zero. A Markov chain is irreducible if there is one communicating class, the state space.
A state i has period k if k is the greatest common divisor of the number of transitions by which i can be reached, starting from i. That is:
The state is periodic if ; otherwise and the state is aperiodic.
A state i is said to be transient if, starting from i, there is a non-zero probability that the chain will never return to i. It is called recurrent (or persistent) otherwise.[48] For a recurrent state i, the mean hitting time is defined as:
State i is positive recurrent if is finite and null recurrent otherwise. Periodicity, transience, recurrence and positive and null recurrence are class properties — that is, if one state has the property then all states in its communicating class have the property.[49]
A state i is called absorbing if there are no outgoing transitions from the state.
Since periodicity is a class property, if a Markov chain is irreducible, then all its states have the same period. In particular, if one state is aperiodic, then the whole Markov chain is aperiodic.[50]
If a finite Markov chain is irreducible, then all states are positive recurrent, and it has a unique stationary distribution given by .
A state i is said to be ergodic if it is aperiodic and positive recurrent. In other words, a state i is ergodic if it is recurrent, has a period of 1, and has finite mean recurrence time.
If all states in an irreducible Markov chain are ergodic, then the chain is said to be ergodic. Equivalently, there exists some integer such that all entries of are positive.
有限状態既約マルコフ連鎖は、非周期状態を持つ場合、エルゴード的であることが示される。
複数の状態を持ち、各状態から出る遷移が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種類あり、各連続状態が観測可能かどうか、また観測に基づいてシステムを調整するかどうかによって分類されます。
A Bernoulli scheme is a special case of a Markov chain where the transition probability matrix has identical rows, which means that the next state is independent of even the current state (in addition to being independent of the past states). A Bernoulli scheme with only two possible states is known as a Bernoulli process.
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]
A reaction network is a chemical system involving multiple reactions and chemical species. The simplest stochastic models of such networks treat the system as a continuous time Markov chain with the state being the number of molecules of each species and with reactions modeled as possible transitions of the chain.[65] Markov chains and continuous-time Markov processes are useful in chemistry when physical systems closely approximate the Markov property. For example, imagine a large number n of molecules in solution in state A, each of which can undergo a chemical reaction to state B with a certain average rate. Perhaps the molecule is an enzyme, and the states refer to how it is folded. The state of any single enzyme follows a Markov chain, and since the molecules are essentially independent of each other, the number of molecules in state A or B at a time is n times the probability a given molecule is in that state.
The classical model of enzyme activity, Michaelis–Menten kinetics, can be viewed as a Markov chain, where at each time step the reaction proceeds in some direction. While Michaelis-Menten is fairly straightforward, far more complicated reaction networks can also be modeled with Markov chains.[66]
An algorithm based on a Markov chain was also used to focus the fragment-based growth of chemicals in silico towards a desired class of compounds such as drugs or natural products.[67] As a molecule is grown, a fragment is selected from the nascent molecule as the "current" state. It is not aware of its past (that is, it is not aware of what is already bonded to it). It then transitions to the next state when a fragment is attached to it. The transition probabilities are trained on databases of authentic classes of compounds.[68]
Also, the growth (and composition) of copolymers may be modeled using Markov chains. Based on the reactivity ratios of the monomers that make up the growing polymer chain, the chain's composition may be calculated (for example, whether monomers tend to add in alternating fashion or in long runs of the same monomer). Due to steric effects, second-order Markov effects may also play a role in the growth of some polymer chains.
Similarly, it has been suggested that the crystallization and growth of some epitaxial superlattice oxide materials can be accurately described by Markov chains.[69]
Markov chains are used in various areas of biology. Notable examples include:
Markov chains are used throughout information processing. Claude Shannon's famous 1948 paper A Mathematical Theory of Communication, which in a single step created the field of information theory, opens by introducing the concept of entropy by modeling texts in a natural language (such as English) as generated by an ergodic Markov process, where each letter may depend statistically on previous letters.[72] Such idealized models can capture many of the statistical regularities of systems. Even without describing the full structure of the system perfectly, such signal models can make possible very effective data compression through entropy encoding techniques such as arithmetic coding. They also allow effective state estimation and pattern recognition. Markov chains also play an important role in reinforcement learning.
Markov chains are also the basis for hidden Markov models, which are an important tool in such diverse fields as telephone networks (which use the Viterbi algorithm for error correction), speech recognition and bioinformatics (such as in rearrangements detection[73]).
The LZMA lossless data compression algorithm combines Markov chains with Lempel-Ziv compression to achieve very high compression ratios.
マルコフ連鎖は、待ち行列の解析的処理(待ち行列理論)の基礎となるものです。アグナー・クラルップ・エルランが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 /日付の不一致(ヘルプ)