📊 対象級:準1級 ・ 1級 | 重要度:B(標準)
確率過程(マルコフ連鎖・ポアソン過程)
要点(BLUF)
確率過程とは、時間で添字付けられた確率変数の族 のこと。要するに「時間とともに値が確率的に動いていくもの」です。試験で問われる代表は次の2つです。
要するに、マルコフ連鎖は「未来は現在だけで決まり過去は忘れる(無記憶)離散ステップの過程」、ポアソン過程は「ランダムに到着する事象を数える連続時間の過程」です。準1級は推移確率行列と定常分布の計算・ポアソン過程の基本性質、1級は収束定理(エルゴード定理)と状態分類の理論まで問われます。
1. 確率過程とは
確率過程(stochastic process) は、添字集合 の各 に確率変数 を対応させた集まり
です。 が のような離散集合なら離散時間、 のような連続集合なら連続時間の過程です。各 が取りうる値の集合を状態空間 と呼びます。
要するに「1個の確率変数を時間方向に並べたもの」です。 を固定すれば普通の確率変数 、(試行結果)を固定すれば時間の関数 が1本得られ、これを標本路(サンプルパス) と呼びます。
| 過程 | 時間 | 状態空間 | 本ノートでの扱い |
|---|---|---|---|
| マルコフ連鎖 | 離散 | 離散(有限・可算) | 推移確率行列・定常分布 |
| ポアソン過程 | 連続 | 離散(計数 ) | 到着数・到着間隔 |
時系列モデル(AR・MA など)も離散時間連続状態の確率過程の一種で、こちらは 時系列解析(定常性・ACF/PACF・AR・MA・ARMA・ARIMA) で扱います。本ノートは確率過程の基礎理論に焦点を当て、マルコフ連鎖の MCMC への応用は MCMC(マルコフ連鎖モンテカルロ) に譲ります。
2. マルコフ性(無記憶性)
マルコフ性とは、未来の状態が「現在の状態だけ」で決まり、それ以前の履歴に依存しない性質です。離散時間で書くと
要するに「今どこにいるかが分かれば、そこに至るまでの経路は次の一歩に一切影響しない」。過去を条件に加えても確率が変わらない、という条件付き独立の言い換えです(条件付き確率の一般論は 条件付き確率・独立性・全確率の定理)。
この右辺を推移確率と呼び、 から への1ステップ確率を
と書きます。さらに、この が (時刻)によらず一定であるとき時間的に一様(斉時的, homogeneous) といい、本ノートではこの一様な連鎖だけを扱います。要するに「ルールが時間で変わらない」ということです。
3. 離散時間マルコフ連鎖
3.1 推移確率行列
状態空間が有限 のとき、推移確率 をすべて並べた行列
を推移確率行列(遷移行列) と呼びます。要するに「 行 列が、状態 から状態 へ1歩で移る確率」を集めた表です。
この行列には絶対的な制約があります。各成分は確率なので 、かつ各行の和が必ず1です。
要するに「今 にいるなら、次にどこかの状態へは必ず行く(行き先の確率を全部足せば1)」。この性質を持つ行列を確率行列(stochastic matrix) と呼びます。⚠️ 列の和ではなく行の和が1です(行=出発点、列=行き先)。
具体例として、天気が「晴 (S)・雨 (R)」の2状態で動くモデルを考えます。
stateDiagram-v2
direction LR
S : 晴 S
R : 雨 R
S --> S : 0.8
S --> R : 0.2
R --> S : 0.4
R --> R : 0.6
「晴の翌日が晴は 0.8、雨は 0.2(行和 1)」「雨の翌日が晴は 0.4、雨は 0.6(行和 1)」という読み方です。
3.2 n段推移とチャップマン–コルモゴロフ方程式
「 ステップ後」の推移確率を考えます。 から 歩で に到達する確率を
と書きます。これがどう計算できるかを与えるのがチャップマン–コルモゴロフ方程式です。 ステップと ステップに分けて、途中の状態 で場合分けします。
要するに「 から 歩で中継地 へ行き、そこから 歩で へ行く。中継地 をすべて足し合わせれば 歩の確率になる」(全確率の定理をマルコフ性のもとで適用したもの)。マルコフ性のおかげで「 に到達した経路」を忘れてよいので、 と の単純な積になります。
この式の右辺は、ちょうど行列の積の定義です。だから
すなわち「 段推移確率行列は、推移確率行列 の 乗そのもの」になります。要するに「 ステップ後の確率を知りたければ、 を 回掛ければよい」。これがマルコフ連鎖の計算が線形代数に帰着する理由です。
時刻 0 の状態分布を行ベクトル とすると、 時刻後の分布は
要するに「初期分布(行ベクトル)に を右から 回掛ければ、 時刻後にどの状態にいる確率がいくつか分かる」。
3.3 定常分布
連鎖を長く回したとき、状態分布が落ち着く先が定常分布です。行ベクトル が
を満たすとき、 を定常分布(不変分布) と呼びます。要するに「 を1回掛けても分布が変わらない(自分自身に戻る)分布」です。一度この分布になれば、その後どれだけ時刻が進んでも分布は のまま固定されます。
これは線形代数の言葉では「固有値 1 に対応する左固有ベクトルを、和が1になるよう規格化したもの」です(確率行列は必ず固有値1を持つ。なぜなら全成分1の列ベクトル に対し 、行和が1だから。右固有ベクトルが なら左にも固有値1の固有ベクトルが存在する)。
定常分布の求め方(完全に展開)
を成分で書くと、各 について
要するに「状態 にいる定常確率は、各状態 から へ流れ込む確率 をすべて足したもの(流入=定常)」。これは状態数 本の連立方程式ですが、 は と同値で、 はランク落ち(行和0)なので1本は冗長です。そこで規格化条件 を1本足して解きます。
先の天気の例 で実際に解きます。 として は
第1式を整理すると つまり (第2式も同じ式に帰着=冗長)。規格化 に代入して
要するに「長期的には、3日に2日が晴・1日が雨に落ち着く」。検算:。確かに不変です。
3.4 状態の分類(1級の山)
定常分布がいつ存在し一意になるかは、状態の性質で決まります。1級ではこの分類の理論が問われます。
-
到達可能・既約: から へ有限ステップで到達できる( となる が存在する)とき「 は から到達可能」。すべての状態が互いに到達可能なら連鎖は既約(irreducible) という。要するに「どの状態からでも、どの状態へもいつかは行ける(孤立した塊がない)」。
-
再帰的・過渡的:状態 から出発していつか必ず に戻る確率を とする。 なら再帰的(recurrent, 再帰)、 なら過渡的(transient, 過渡)。要するに「再帰=必ず戻ってくる状態、過渡=二度と戻らない可能性がある(やがて見捨てられる)状態」。有限既約な連鎖は全状態が再帰的。
-
周期性:状態 に戻れるステップ数の集合の最大公約数 を周期と呼ぶ。 なら非周期的(aperiodic)、 なら周期的。要するに「戻れるのが必ず の倍数ステップ後だけなら周期 、特定の周期に縛られず戻れるなら非周期」。例えば2状態を と必ず往復する連鎖は周期2です。
graph LR
A["既約<br/>どの状態にも到達可"] --> D["エルゴード的<br/>= 既約 + 非周期"]
B["非周期<br/>周期 d = 1"] --> D
D --> E["定常分布が一意<br/>かつ極限分布に収束"]
3.5 収束定理(エルゴード定理)

初期分布が違っても(実線=状態1から/破線=状態3から)数ステップで同じ定常分布 π に収束する(エルゴード性)。図は simulations/markov_shusoku.py で生成。
有限状態で既約かつ非周期的なマルコフ連鎖(エルゴード的という)では、初期分布によらず分布が定常分布へ収束します。
要するに「十分時間が経つと、どこから出発したか()を忘れて、状態 にいる確率が定常確率 に近づく」。さらにこのとき定常分布 はただ1つに定まります(既約性が一意性を、非周期性が収束を保証する)。
この極限としての分布を極限分布と呼びます。エルゴード的なら「定常分布=極限分布」で一致しますが、両者は定義が別物であり、周期があると一致しない点が重要です(次章)。
なぜ非周期性が要るのか(直観)
周期2の往復連鎖 を考えます。 は を満たす立派な定常分布です。ところが状態1から出発すると分布は と永遠に振動し、どこにも収束しません。 は が偶数なら1、奇数なら0で、 に近づきません。要するに「周期があると、定常分布は存在しても極限としては到達できない」。非周期性は、この振動を消して収束を保証するための条件です。
4. 詳細釣り合い(可逆連鎖)
定常分布を「解かずに見抜く」便利な条件が詳細釣り合い(detailed balance) です。分布 がすべての状態対 について
を満たすとき、連鎖は について可逆(reversible) といいます。要するに「定常状態で、 の向きに流れる確率と、 の向きに流れる確率がちょうど釣り合う」。各リンクごとに行きと帰りの流量が等しい、という強い局所条件です。
詳細釣り合い → 定常分布(導出)
詳細釣り合いを満たす は、必ず定常方程式 を満たします。各 について両辺を で和を取ります。
左から順に、(1) 詳細釣り合い を代入、(2) を和の外に出す、(3) ( の行和は1)。結果 はまさに定常方程式の 成分です。
要するに「詳細釣り合い(各リンクの局所的な釣り合い)が成り立てば、各状態への流入=流出(大域的な釣り合い=定常)が自動的に従う」。この含意こそが詳細釣り合いの価値で、MCMC(MCMC(マルコフ連鎖モンテカルロ))では「狙った分布 を定常分布に持つ連鎖」を、解きにくい定常方程式の代わりに詳細釣り合いを満たすように推移確率を設計することで作ります(メトロポリス・ヘイスティングス法の土台)。
⚠️ 逆は成り立ちません。詳細釣り合いは定常分布であるための十分条件であって必要条件ではない。詳細釣り合いを満たさない(非可逆な)定常分布はいくらでもあります。1級ではこの「十分だが必要でない」が頻出の引っかけです(次章)。
5. ポアソン過程
ポアソン過程は、時刻 0 から始めて事象(到着・故障・来店)がランダムに起きるのを数える連続時間の確率過程です。時刻 までの累積到着数を とします()。
5.1 定義(3つの公理)
ポアソン過程(強度 )は次の3条件で特徴づけられます。
- 独立増分:互いに重ならない時間区間での到着数は独立。要するに「過去の到着が未来の到着に影響しない(区間どうしが独立)」。
- 定常増分:区間 の到着数の分布は、始点 によらず長さ だけで決まる。要するに「いつ測っても、同じ長さなら同じ分布(時間的に一様)」。
- 希少性(稀発性):微小時間 の間に到着が1回起きる確率は 、2回以上起きる確率は 。要するに「ごく短い時間には、たかだか1回しか起きない(同時多発しない)」。
これら(独立増分+定常増分+無記憶な発生)から、後述のとおり到着数がポアソン分布に、到着間隔が指数分布になることが導かれます。これがポアソン過程・ポアソン分布・指数分布が密結合する理由です。
5.2 到着数はポアソン分布
3公理から、時刻 までの到着数 は平均 のポアソン分布に従います。
要するに「強度 (単位時間あたり平均 回)の過程を時間 走らせると、到着回数は平均 のポアソン分布になる」。、(ポアソン分布の平均=分散、ポアソン分布)。 を強度(到着率, rate) と呼びます。
導出の概略:区間 を 等分し、各小区間で「到着1回 ≒ 確率 のベルヌーイ試行」とみなす(希少性より1区間に高々1回)。到着総数は で、独立増分よりこれらは独立。 で二項分布のポアソン極限( 一定)により に収束します(少数の法則は ポアソン分布)。
5.3 到着間隔は指数分布(完全導出)
ポアソン過程で「次の事象が起きるまでの待ち時間(到着間隔)」 の分布を求めます。ここが準1級頻出の核心です。
待ち時間 が より大きい、つまり「時刻 までまだ1回も到着していない」という事象は、 と同じです。
要するに「次が来るまで 以上待つ」=「最初の の間に到着0回」。 なので を代入すると
よって累積分布関数は 、これを微分して密度は
これはまさに指数分布 (指数分布・ガンマ分布・ベータ分布)です。要するに「ポアソン過程の到着間隔は、強度 の指数分布に従う」。平均待ち時間は (強度が高いほど待ち時間は短い)で直観に合います。
指数分布の無記憶性 は、ポアソン過程の独立増分・定常増分の裏返しです。要するに「すでに だけ待ったという事実は、これから先の待ち時間に一切影響しない(バスはあと何分かを、これまで待った時間と無関係に決める)」。
5.4 第 n 到着時刻はガンマ分布
回目の到着が起きる時刻 は、独立な到着間隔 (各 )の和です。
独立な指数分布の和はガンマ分布になる(指数分布・ガンマ分布・ベータ分布)ので
要するに「 回目の到着までの所要時間は、形状 ・強度 のガンマ分布」。(間隔の平均 が 個分)。ポアソン過程は「到着数=ポアソン分布」「到着間隔=指数分布」「第 到着時刻=ガンマ分布」の三位一体で覚えるのが定石です。
| 問いたいもの | 分布 | パラメータ |
|---|---|---|
| 時刻 までの到着数 | ポアソン分布 | 平均 |
| 次の到着までの待ち時間 | 指数分布 | 率 、平均 |
| 回目の到着時刻 | ガンマ分布 | 形状 ・率 、平均 |
6. 引っかけ・頻出論点
- ⚠️ 定常分布と極限分布は別物:定常分布は「 を満たす分布(存在すれば複数のこともある)」、極限分布は「 で実際に収束する先」。エルゴード的(既約+非周期)なら一致して一意ですが、周期があると定常分布は存在しても極限分布は存在しない(振動して収束しない)。「定常分布があるから収束する」は誤り。
- ⚠️ 周期性で極限が消える:周期 の連鎖は が に収束しない( の偶奇などで値が振れる)。収束を主張するには非周期性が必須。
- ⚠️ 詳細釣り合いは十分条件であって必要条件でない:詳細釣り合い 定常分布、は正しい。しかし定常分布 詳細釣り合い、は偽。非可逆な定常分布は普通に存在します(一方向に循環する連鎖など)。「定常分布なら詳細釣り合いを満たす」は誤り。
- ⚠️ 推移確率行列は行和が1(列和ではない):。行=出発状態、列=到着状態。定常分布は左固有ベクトル(、行ベクトルに右から )。(右固有ベクトル)と混同しない。
- ⚠️ ポアソン過程の3公理を1つでも欠くと別物:独立増分・定常増分・希少性のどれかが崩れると到着数はポアソン分布にならない(例:定常増分を外すと非斉時ポアソン過程で率が )。
- ⚠️ 「率 」と「平均 」の取り違え:指数分布のパラメータは率 (密度 )だが平均は 。「1時間に平均3件」なら (時間単位)で平均待ち時間は 時間。
試験での問われ方(級ごとの差)
確率過程は準1級・1級の出題範囲です(マルコフ連鎖は準1級ワークブック第14章相当、ポアソン過程は確率過程・時系列の単元)。級で深さが明確に違います(出題範囲・配点は改訂されうるため要最新確認)。
準1級レベル
ここで問われるのは「計算と基本性質」。推移確率行列から 段推移や定常分布を実際に計算できるか、ポアソン過程の到着数・到着間隔の分布を正しく当てはめられるか。
- 推移確率行列 が与えられ、 と を連立して定常分布を数値で求める(2〜3状態)。
- や で ステップ後の状態確率を計算する。
- ポアソン過程で「時刻 までに 回」「次まで 以上待つ」確率を、・到着間隔 から計算する。
- 有限マルコフ連鎖のパラメータ推定(観測した遷移回数から を最尤推定する)。
1級レベル
ここで問われるのは「理論と収束」。状態分類(既約・再帰/過渡・周期)を定義し、エルゴード定理(既約+非周期 定常分布が一意・極限に収束)を理由とともに説明できるか。
- 定常分布と極限分布の区別を論じ、周期性があると極限分布が存在しないことを示す。
- チャップマン–コルモゴロフ方程式から を導く。
- 詳細釣り合い 定常分布の証明、および「十分条件であって必要条件でない」ことの理解。
- ポアソン過程の3公理から到着数のポアソン分布・到着間隔の指数分布を導出する。連続時間マルコフ連鎖(生成行列)への接続。
確率の土台は 条件付き確率・独立性・全確率の定理、分布の前提は ポアソン分布 と 指数分布・ガンマ分布・ベータ分布、連続状態の確率過程としての時系列は 時系列解析(定常性・ACF/PACF・AR・MA・ARMA・ARIMA) が関連します。
よくある疑問(Q&A)
Q1. 定常分布と極限分布は、結局何が違うのですか?
定義が別物です。定常分布は方程式 を満たす分布で、「いったんこの分布になれば以後変わらない」という静的な性質。極限分布は初期分布から出発して で「実際に収束する先」という動的な性質です。エルゴード的(既約+非周期)な連鎖では両者は一致して一意になりますが、周期 の連鎖は定常分布を持っていても分布が振動して収束しないので極限分布を持ちません。「定常分布が存在する=収束する」ではない、ここが1級頻出の区別です。
Q2. なぜ非周期性がないと収束しないのですか?
周期があると分布が「決まった周期で振動」して落ち着かないからです。例えば は状態1と2を必ず交互に往復します。状態1から出発すれば分布は と振れ続け、定常分布 には近づきません。 は偶数ステップで1・奇数ステップで0となり、 に収束しないのです。非周期性はこの振動を消し、 を1つの値 に収束させるための条件です。
Q3. 詳細釣り合いを満たさないと定常分布は存在しないのですか?
いいえ。詳細釣り合いは定常分布であるための十分条件であって必要条件ではありません。詳細釣り合い(各リンクで行き帰りの流量が等しい)を満たせば必ず定常分布になりますが、逆に定常分布が詳細釣り合いを満たすとは限りません。例えば3状態を と一方向に循環する連鎖は、一様分布が定常分布になりますが、リンクは一方通行なので (行き帰りの釣り合い)は満たしません。詳細釣り合いはあくまで「定常分布を作りやすくする・MCMCで設計しやすくする」ための強い十分条件です。
Q4. ポアソン分布・指数分布・ガンマ分布は、なぜ全部ポアソン過程に出てくるのですか?
同じポアソン過程を「3つの角度から見ている」からです。ポアソン過程は「ランダムに到着する事象の列」で、(1)「ある時間 までの到着回数」を見るとポアソン分布 、(2)「次の到着までの待ち時間」を見ると指数分布 、(3)「 回目の到着までの総時間(=指数分布の和)」を見るとガンマ分布 になります。「回数を数える」か「時間を測る」かの視点の違いで、根は1つの過程です。だから「到着間隔が指数 到着数がポアソン」は同値な言い換えになります。
Q5. マルコフ性は「1ステップ前だけ」見ますが、2ステップ前を見たいときはどうするのですか?
状態の定義を拡張します。「過去2ステップを使うと未来が決まる」過程(2階のマルコフ)は、 というペアを新しい状態とみなせば、その は1階のマルコフ連鎖になります。要するに「必要な履歴を状態に畳み込めば、見かけ上『現在だけ』に依存する形に戻せる」。これがマルコフ連鎖が非常に汎用的な理由で、状態空間を広げる代償でマルコフ性を回復できます。だから「マルコフ性は本質的に限定的」というより「状態の設計次第で広く適用できる」枠組みです。
まとめ
- 確率過程は時間で添字付けられた確率変数の族 。試験の主役はマルコフ連鎖(離散時間)とポアソン過程(連続時間)。
- マルコフ性=未来は現在だけで決まる(過去を忘れる)。推移確率 を集めたのが推移確率行列(各行和1の確率行列)。 段推移はチャップマン–コルモゴロフ方程式より 。
- 定常分布 は固有値1の左固有ベクトル。連立方程式+規格化で解く。既約+非周期(エルゴード的)なら定常分布は一意かつ極限分布に一致(エルゴード定理)。
- 詳細釣り合い は定常分布の十分条件(必要条件ではない)。MCMC の設計原理。
- ポアソン過程(独立増分・定常増分・希少性):到着数 、到着間隔 、第 到着時刻 。
- 引っかけ:定常分布≠極限分布(周期で極限消失)/詳細釣り合いは十分のみ/行和1・左固有ベクトル。準1級は計算、1級は収束定理と状態分類の理論。
関連ノート
- ポアソン分布 ポアソン過程の到着数 が従う分布。少数の法則による導出の土台
- 指数分布・ガンマ分布・ベータ分布 到着間隔(指数)と第 到着時刻(ガンマ)の分布。無記憶性
- 条件付き確率・独立性・全確率の定理 マルコフ性(条件付き独立)とチャップマン–コルモゴロフの基礎
- MCMC(マルコフ連鎖モンテカルロ) 詳細釣り合いを定常分布に持つ連鎖を設計する応用(メトロポリス・ヘイスティングス)
- 時系列解析(定常性・ACF/PACF・AR・MA・ARMA・ARIMA) 連続状態の離散時間確率過程(AR・MA など)
- 確率変数の変換・モーメント母関数・積率 分布の導出・和の分布(ガンマ=指数の和)に使う道具
- 確率過程・時系列・応用(Phase 8 目次) 確率過程・時系列・応用ドメインの全体地図