1

私は2つの対数正規分布からランダムな値のペアを生成しようとしています-キャッチは、それらの一方が他方よりも小さくなければならないということです。例えば：

a1 <- log(47.31)
b1 <- sqrt(2*log(50.84/47.31))
a2 <- log(47.31)
b2 <- sqrt(2*log(59.34/47.31))

x1 <- rlnorm(1,a1,b1)
x2 <- rlnorm(1,a2,b2)

x1<x2であることを確認する方法が必要です。これを行うための巧妙な方法はありますか？

4

1 に答える 1

1

ええ、はい、いいえ。最も簡単な方法は、条件が満たされているかどうかを確認し、満たされていない場合はランダムを再生成することです。しかし、この結果、変数は、開始した統計分布によって特徴付けられなくなります。つまり、フィルタリングプロセスはx1を低くバイアスし、x2を高くバイアスします。しかし、これで問題がなければ、目的の条件が満たされるまでループするだけです...理論的には、これには無限の反復が必要になる可能性がありますが、それほど不運ではないと思います:)。

2つの分布が同じである場合は、より簡単です。x1> x2の場合は、それらを交換するだけです（等しくないと思います！）

于 2012-06-21T20:52:09.143 に答える