折りたたまれた正規分布からのサンプリングは、0で切り捨てられた正規分布からのサンプリングと同等ですか?


9

通常の密度(たとえば、mean = 1、sd = 1)からシミュレーションしたいのですが、正の値のみが必要です。

1つの方法は、法線からシミュレーションし、絶対値を取ることです。これは普通の折りたたみだと思います。

Rには、切り捨てられたランダム変数を生成するための関数があることがわかります。打ち切られた法線(0での打ち切り)からシミュレーションすると、これは折り畳みアプローチと同じですか?

回答:


10

はい、アプローチはゼロ平均正規分布に対して同じ結果を与えます。

ΦΦ((a,b])(a,b]0ab

Φtruncated((a,b])=Φ((a,b])/Φ([0,])=2Φ((a,b])

Φ([0,])=1/2

Φfolded((a,b])=Φ((a,b])+Φ([b,a))=2Φ((a,b])

Φ0

00

3つの分布

このグラフは、正規(1,1)分布(黄色)、折りたたまれた正規(1,1)分布(赤)、および切り捨てられた正規(1,1)分布(青)の確率密度関数を示しています。折りたたまれた分布が他の2つの分布と特徴的なベルカーブ形状を共有していないことに注意してください。青い曲線(切り捨てられた分布)は黄色の曲線の正の部分であり、単位面積を持つように拡大されています。一方、赤い曲線(折りたたまれた分布)は、黄色の曲線の正の部分とその負の尾の合計です(周囲に反映されています)。 y軸)。


1
私は絵が好きです。
Karl、

5

XN(μ=1,SD=1)X|X>0|X|

Rでの簡単なテスト:

x <- rnorm(10000, 1, 1)
par(mfrow=c(2,1))
hist(abs(x), breaks=100)
hist(x[x > 0], breaks=100)

これは以下を与えます。 シミュレーションヒストグラム

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