分散減少法

OR事典

分散減少法

読み方：ぶんさんげんしょうほう
【英】：variance reduction method

単純なモンテカルロ法では,実験を繰り返して, 同一分布からの独立標本と考えられる数値を得て, これらの算術平均によってシステムの特性値(前記の同一分布の平均値)を推定する. この推定値に含まれる誤差の標準偏差は, 標本数の平方根に反比例するので, 高精度の推定値を得るのは困難である. この点を改善し, 反比例の係数を小さくするための方法が分散減少法である.

「OR事典」の他の用語

シミュレーション：

事象  事象リスト  共通乱数法  分散減少法  時間制御方式  確定的シミュレーション  確率的シミュレーション

ウィキペディア

分散減少法

出典: フリー百科事典『ウィキペディア（Wikipedia）』 (2019/12/02 01:12 UTC 版)

単位正方形内にランダムに生成される点の分散は層化抽出法によって減少させることができる。

数学、特にモンテカルロ法の理論における分散減少法（ぶんさんげんしょうほう、英: variance reduction）は推定の精度を改善するのに用いられる手法であり、与えられたシミュレーション、計算量（computational effort）に応じて適用し得る。シミュレーションの出力値となる確率変数は、その結果の精度を左右する量である分散と結び付いている。シミュレーションを統計上効果的に、つまり、注目している確率変数の出力がより高い精度・より狭い信頼区間となるようにするために、分散減少法が利用できる場合がある。代表的なものに共通乱数法、対称変量法、制御変量法、重点サンプリング法、層化抽出法がある。ブラックボックスモデルを使ったシミュレーションに対しては、部分空間シミュレーション法やラインサンプリング法が用いられることもある。これらの項目の下位区分に、様々な特化型の技法が存在する。例えば、粒子輸送シミュレーションでは広範にわたって「ウェイト・ウインドウ法（weight windows）」や「セルインポータンス法（splitting/Russian roulette）」の技法が用いられるが、これらは重点サンプリング法の一形式である。

共通乱数法

共通乱数法（Common Random Numbers, CRN）は、（系の単一の設定（configuration）に注目するのではなく）複数の設定間の比較をする際に適用できる、頻用され有用な手法である。correlated sampling、matched streams、matched pairs とも呼ばれる。

共通乱数法では乱数列を同期（synchronization）させる。同一の乱数列を全ての設定に対して共通に用い、また発生する各乱数を全ての設定で同一の用途に用いるということである。例えば待ち行列理論において銀行の窓口係の2通りの配置を比較する場合であれば、各配置での N 番目の顧客のランダムな到着時刻を、共通の乱数列から生じる乱数を用いてそれぞれ発生させる。

共通乱数法の原理

$X_{1j}$ と $X_{2j}$ とを、第1の設定と第2の設定とで独立に反復した j 回目の出力結果とし、今

$\xi =E(X_{1j})-E(X_{2j})=\mu _{1}-\mu _{2}$

を推定したいとする。各設定について出力を n 回実行し

$Z_{j}=X_{1j}-X_{2j}\quad {\mbox{for }}j=1,2,\ldots ,n$

とおくと $E(Z_{j})=\xi$ であり、$Z(n)={\frac {\sum _{j=1,\ldots ,n}Z_{j}}{n}}$ は $\xi$ の不偏推定量である。

確率変数列 $Z_{j}$ は独立同分布なので、

$\operatorname {Var} [Z(n)]={\frac {\operatorname {Var} (Z_{j})}{n}}={\frac {\operatorname {Var} [X_{1j}]+\operatorname {Var} [X_{2j}]-2\operatorname {Cov} [X_{1j},X_{2j}]}{n}}$

乱数列を共有せず、それぞれに別個のものを用いる場合は Cov(X_1j, X_2j) = 0 だが、X₁ と X₂ との間に正の相関をもたらす（Cov(X_1j, X_2j) > 0）ような要素をうまく導入することができれば、上式の通り分散を減少させられる。

共通の乱数列によって相関が負になる（Cov(X_1j, X_2j) < 0）場合、この手法は逆効果、つまり意に反して分散を増大させてしまい得ることもわかる。

関連項目

• Explained variance (回帰分析等における、"説明される変動")

参考文献

• Hammersley, J. M.; Handscomb, D. C. (1964). Monte Carlo Methods. London: Methuen. ISBN 0-416-52340-4. 
• Kahn, H.; Marshall, A. W. (1953). “Methods of Reducing Sample Size in Monte Carlo Computations”. Journal of the Operational Research Society 1 (5): 263–271. doi:10.1287/opre.1.5.263. 
• MCNP — A General Monte Carlo N-Particle Transport Code, Version 5 Los Alamos Report LA-UR-03-1987