分散減少法のソースを表示

[[File:StratifiedPoints.gif|thumbnail|right|単位正方形内にランダムに生成される点の分散は層化抽出法によって減少させることができる。]]

[[数学]]、特に[[モンテカルロ法]]の理論における'''分散減少法'''（ぶんさんげんしょうほう、{{Lang-en-short|variance reduction}}）は推定の精度を改善するのに用いられる手法であり、与えられたシミュレーション、計算量（computational effort）に応じて適用し得る<ref name="varred17">{{cite journal|last1=Botev|first1=Z.|last2=Ridder|first2=A.|title=Variance Reduction|journal= Wiley StatsRef: Statistics Reference Online|date=2017|pages=1–6|doi=10.1002/9781118445112.stat07975}}</ref>。シミュレーションの出力値となる[[確率変数]]は、その結果の精度を左右する量である[[分散 (確率論)|分散]]と結び付いている。シミュレーションを統計上効果的に、つまり、注目している確率変数の出力がより高い精度・より狭い[[信頼区間]]となるようにするために、分散減少法が利用できる場合がある。代表的なものに共通乱数法、{{仮リンク|対称変量法|en|antithetic variates}}、{{仮リンク|制御変量法|en|control variate}}、{{仮リンク|重点サンプリング法|en|importance sampling}}、[[層化抽出法]]がある。[[ブラックボックス]]モデルを使ったシミュレーションに対しては、{{仮リンク|部分空間シミュレーション法|en|subset simulation}}や{{仮リンク|ラインサンプリング法|en|line sampling}}が用いられることもある。これらの項目の下位区分に、様々な特化型の技法が存在する。例えば、粒子輸送シミュレーションでは広範にわたって「ウェイト・ウインドウ法（weight windows）」や「セルインポータンス法（splitting/Russian roulette）」の技法が用いられるが、これらは重点サンプリング法の一形式である。

== 共通乱数法 ==

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

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

==共通乱数法の原理==

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

:<math>\xi= E(X_{1j})-E(X_{2j})=\mu_1-\mu_2 </math>

を推定したいとする。各設定について出力を ''n'' 回実行し
:<math>Z_j=X_{1j}-X_{2j} \quad\mbox{for } j=1,2,\ldots, n</math>
とおくと <math>E(Z_j)=\xi</math> であり、<math>Z(n) = \frac{\sum_{j=1,\ldots,n} Z_j}{n}</math> は <math>\xi</math> の不偏推定量である。

確率変数列 <math>Z_j</math> は[[独立同分布]]なので、
:<math>\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} </math>

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

共通の乱数列によって相関が負になる（Cov(''X''<sub>1''j''</sub>, ''X''<sub>2''j''</sub>) < 0）場合、この手法は逆効果、つまり意に反して分散を増大させてしまい得ることもわかる<ref name="wolfram">{{cite web|last1=Hamrick|first1=Jeff|title=The Method of Common Random Numbers: An Example|url=http://demonstrations.wolfram.com/TheMethodOfCommonRandomNumbersAnExample/|website=Wolfram Demonstrations Project|accessdate=29 March 2016}}</ref>。

==関連項目==
*[[Explained variance]] ([[回帰分析]]等における、"説明される変動")

==脚注==
{{Reflist}}

==参考文献==
*{{cite book |last=Hammersley |first=J. M. |last2=Handscomb |first2=D. C. |year=1964 |title=Monte Carlo Methods |publisher=Methuen |location=London |isbn=0-416-52340-4 }}
*{{cite journal |last=Kahn |first=H. |last2=Marshall |first2=A. W. |year=1953 |title=Methods of Reducing Sample Size in Monte Carlo Computations |journal=[[Journal of the Operational Research Society]] |volume=1 |issue=5 |pages=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

{{DEFAULTSORT:ふんさんけんしようほう}}
[[Category:統計的偏差と分散]]
[[Category:モンテカルロ法]]
[[Category:乱数]]
[[Category:数値解析]]
[[Category:数学に関する記事]]