起源と基本原理
モンテカルロ法は1946年、ロスアラモス研究所の数学者スタニスワフ・ウラムが療養中にソリティアの成功確率を計算しようとして思いついたアイデアに始まります。直接計算するのが難しければ、何度も試して成功の割合を数えればよい、というわけです。ジョン・フォン・ノイマンがこれを中性子拡散の計算に応用し、ニコラス・メトロポリスがウラムの叔父がよく通ったカジノにちなんで「モンテカルロ」と名付けました。最初の論文はメトロポリスとウラム(1949年)です。原理を示す古典的な例が円周率の推定で、一辺1の正方形にランダムに点を打ち、四分円の中に入った割合に4を掛けるとπに近づきます。標本が増えるほど平均が真の値に近づくのは大数の法則のおかげです。
精度と分散減少
中心極限定理により、モンテカルロ推定の誤差はおおよそ標本数の平方根に反比例して小さくなります。精度を10倍にするには標本が100倍必要ですが、この速さは問題の次元にほとんど依存しないため、変数の多い高次元積分では格子に基づく数値解析より有利です。効率を上げるため、重点サンプリング、制御変量、層化サンプリング(調査の確率標本抽出と同じ原理)、準モンテカルロ法などの分散減少法が使われます。結果の質は疑似乱数生成器にも左右され、松本眞と西村拓士(1998年)のメルセンヌ・ツイスタが広く使われています。ハマーズレーとハンズコム(1964年)、ロバートとカセラ(2004年)、オーウェン(2013年)が代表的な教科書です。
応用分野
メトロポリスら(1953年)のアルゴリズムとヘイスティングス(1970年)の一般化は、マルコフ連鎖を使って複雑な分布から標本を得るMCMCへと発展し、ゲルファンドとスミス(1990年)以降はベイズ統計の標準的な計算手段になりました。エフロン(1979年)のブートストラップも、データを再抽出して推定値の不確実性を評価するモンテカルロ的手法です。金融ではボイル(1977年)がオプション価格の計算に導入し、リスク管理のシナリオ分析にも使われます。物理学の粒子輸送や統計力学、コンピュータグラフィックスのパストレーシングにも用いられます。人工知能では、コチシュとセペシュヴァーリ(2006年)のUCTアルゴリズムに始まるモンテカルロ木探索が、シルバーら(2016年)のAlphaGoの中核となりました。これは人工知能の探索アルゴリズムの重要な系統です。
よくある質問
- モンテカルロ・シミュレーションは常に正確ですか?
- いいえ。結果は確率的な近似値であり、標本数、乱数の質、モデルの仮定によって誤差が生じます。そのため通常は信頼区間や標準誤差を併せて示します。
- モンテカルロ法とMCMCはどう違いますか?
- 通常のモンテカルロ法は互いに独立な標本を直接抽出しますが、MCMCは直接抽出が難しい分布からマルコフ連鎖を使って互いに依存する標本を順に生成します。
出典確認日:2026-10-26
参考文献
- Metropolis, N., & Ulam, S. (1949). The Monte Carlo method. Journal of the American Statistical Association, 44(247), 335–341.
- Metropolis, N., Rosenbluth, A. W., Rosenbluth, M. N., Teller, A. H., & Teller, E. (1953). Equation of state calculations by fast computing machines. Journal of Chemical Physics, 21(6), 1087–1092.
- Hastings, W. K. (1970). Monte Carlo sampling methods using Markov chains and their applications. Biometrika, 57(1), 97–109.
- Eckhardt, R. (1987). Stan Ulam, John von Neumann, and the Monte Carlo method. Los Alamos Science, 15, 131–137.
- Hammersley, J. M., & Handscomb, D. C. (1964). Monte Carlo Methods. London: Methuen.
- Boyle, P. P. (1977). Options: A Monte Carlo approach. Journal of Financial Economics, 4(3), 323–338.
- Efron, B. (1979). Bootstrap methods: Another look at the jackknife. Annals of Statistics, 7(1), 1–26.
- Gelfand, A. E., & Smith, A. F. M. (1990). Sampling-based approaches to calculating marginal densities. Journal of the American Statistical Association, 85(410), 398–409.
- Kocsis, L., & Szepesvári, C. (2006). Bandit based Monte-Carlo planning. In Machine Learning: ECML 2006 (pp. 282–293). Berlin: Springer.
- Silver, D., et al. (2016). Mastering the game of Go with deep neural networks and tree search. Nature, 529(7587), 484–489.
- Robert, C. P., & Casella, G. (2004). Monte Carlo Statistical Methods (2nd ed.). New York: Springer.
- Owen, A. B. (2013). Monte Carlo Theory, Methods and Examples. Stanford, CA: Stanford University.
- Matsumoto, M., & Nishimura, T. (1998). Mersenne twister: A 623-dimensionally equidistributed uniform pseudo-random number generator. ACM Transactions on Modeling and Computer Simulation, 8(1), 3–30.
