モンテカルロ分析に必要なシミュレーション数


10

私の質問は、モンテカルロ分析法に必要なシミュレーション数についてです。私は、任意の許容されるパーセンテージエラーのためにシミュレーションに必要な数を見る限り、(例えば、5)は 、N = { 100 Z CSTD X E

n={100zcstd(x)Emean(x)}2,

ここで、は結果のサンプリングの標準偏差、z cは信頼水準係数です(たとえば、95%の場合は1.96)。したがって、この方法で、nシミュレーションの結果の平均と標準偏差が実際の平均と95%の信頼水準の標準偏差を表すことを確認できます。std(x)zcn

私の場合、シミュレーションを7500回実行し、7500シミュレーションから100サンプリングの各セットの移動平均と標準偏差を計算します。私が取得する必要なシミュレーションの数は常に100未満ですが、結果全体の平均と標準と比較した平均と標準誤差の%は常に5%未満ではありません。ほとんどの場合、平均誤差率は5%未満ですが、標準誤差は最大30%になります。

実際の平均と標準を知らずに必要なシミュレーションの数を決定する最良の方法は何ですか(私の場合、シミュレーションの対象となる結果は正規分布です)?

助けてくれてありがとう。


シミュレーションが無限に実行されたときにシミュレーション結果の分布がどのようになるかを理解するために、n回のシミュレーション後に結果の平均と分散を使用する代わりに、結果の分布のフィット関数を見つけることにしました。ただし、ここでは、nは%エラーをフルフィルする必要があります。そのようにすると、たとえば97.5%に関連する累積分布関数のより正確な結果を見つけることができると思います。400と7000のシミュレーションの結果を比較すると、両方のサンプリングの分布のフィット関数は互いに似ているため、2番目の曲線の曲線のみが滑らかになります。また、したがって、MATLAB / Simulinkのモデルは非線形ですが、生成された入力パラメーターは正規分布ですが、シミュレーションの結果のヒストグラムは通常ではないため、「一般化極値分布」を使用しました。これは、MATLABでは 'gev'という名前です。しかし、それでも、この方法論についてはよくわかりません。事前のコマンドに感謝します


シミュレーションの結果が任意の合格基準で評価されるときに見る限り、任意の信頼水準で必要なシミュレーション数を見つけることは可能ですが、私の場合、特定の信頼度で結果全体の平均と分散を見つけたいと思います任意の有限反復回数のレベル。したがって、n個のサンプルでは、​​分散を使用して平均の間隔を定義しますが、実際には、CPDFが0.975となる可能性のある値を見つけるためにも分散が必要です。コメントありがとう
maxwell

回答:


4

私は通常、収束スタディを実施し、必要なシミュレーションの数を決定して、その後のシミュレーションでこの数を使用します。また、エラーが選択した数で示されるよりも大きい場合は、警告を出します。

σ^N2σ^NN

または、各シミュレーションのエラーを計算し、特定のしきい値を超えるか、パスの最大数に達したときに停止することもできます。この数は、収束スタディによって決定されました。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.