0

この関数にモンテカルロ法を実装しようとすると、大きな問題が発生します。

D=log(T)

T測定された時間はどこにあるのでT>0、、そして明らかに、それは正規分布を持っています。

実験では10個の測定値があるTので、次のように計算します。

m_T (mean of T)  = 3.0 seconds
s_T (standard deviation of T)= 1.5 seconds

そして、このパラメーターを使用して、シミュレートTし、次にD：

T = Normal(m_T, s_T)
D=log(Normal(m_T, s_T)

しかしD、プログラムではエラーを返します。浄化すると、エラーはNormalにいくつかのNEGATIVE値があるためであることがわかり、(m_T, s_T)log（NEGATIVE）がクラッシュします！

ブロックされています。続行する方法がわかりません…何か提案はありますか？どうもありがとうございます！

4

2 に答える 2

1

定義により、正規分布は常に負の値に対して有限の確率をもたらします。次に、測定したもの (時間) は、厳密には正規分布ではありません。

切り捨てられた正規分布では、特定の境界に収まらないすべての値に 0 の確率が割り当てられますが、0 未満の値を無視することで、分布の平均と分散が変更されます。

于 2013-03-11T09:15:30.210 に答える
0

いくつかのコメント:

  1. obviously, it has a normal distribution.これは少しも明らかではありません。より良いのは、対数正規分布を使用することです。切り捨てられた法線はあなたが求めているものではないと強く思います。あなたのパラメータを使用して、数字が終わりになるようにします。では確率がかなり高いことに注意してくださいx=0。

  2. 代わりに、対数正規分布、指数分布、またはその他のより適切な分布を使用する必要があります。真の分布からのモーメントを観測値に一致させるか、それらの最尤推定量を使用できます。

ここに画像の説明を入力

于 2013-03-11T09:30:20.573 に答える