の分布


7

いう、 XRn (と n>1)密度があります fX(x)。の分布について何が言えるか

Y=logfX(X)?

4
まあそれは何に依存するだろう f され、それはないですか?
jbowman

2
1. mgf(またはより一般的にはcf)を検討することから始めて、そこから何が言えるかを理解するのは興味深いかもしれません。あるいは、漸近的な振る舞いに関心がある場合(特に、独立性を扱う場合)、次の漸近性について知られていることを検討することをお勧めします。2logL... 2.これは演習用ですか?
Glen_b-モニカを

4
Troutt et al。によるこれに捧げられた本全体があります。(1991)
西安

回答:


9

西安が書いた本は2004年のものです。1991年の記事で、次の定理が出ています。

From:Troutt MD 1991密度縦座標の密度に関する定理とBox-Muller法の別の解釈

確率変数Xに密度がある場合 f(x)xRn、そして確率変数が v=f(x) 密度があります g(v)、その後

g(v)=vA(v),
どこ A(v) セットのルベーグ測度です
S(v)={x:f(x)v}

直感的で非公式:

fZ(z)dz=P(z<Z<z+dz)=P(x(z)<X<x(z+dz))=P(x(z)<X<x(z)+dzdxdz)=fX(X)dxdzdz=zdA(z)dzdz

同様に、変換された変数を使用する場合 Y=g(fx(x)) 次に:

fY(y)dy=P(y<Y<y+dy)=P(x(y)<X<x(y+dy))=P(x(y)<X<x(y)+dydxdy)=fX(X)dxdydy=g1(y)dA(y)dydy

そう

fY(y)=eyA(y)dy

標準正規分布の例:

fX(x)=12πe0.5x2

y=log(2π)+0.5x2

A(y)=C8(ylog(2π))

したがって

fY(y)=2eyylog(2π)2

多変量正規分布の例:

fX(x1,x2)=12πe0.5(x12+x22)

y=log(2π)+0.5(x12+x22)

A(y)=C2π(ylog(2π))

したがって

fY(y)=2πeyfor ylog(2π)

計算チェック:

ここに画像の説明を入力してください

# random draws/simulation
x_1 = rnorm(100000,0,1)
x_2 = rnorm(100000,0,1)
y = -log(dnorm(x_1,0,1)*dnorm(x_2,0,1))

# display simulation along with theoretic curve
hist(y,breaks=c(0,log(2*pi)+c(0:(max(y+1)*5))/5),
     main = "computational check for distribution f_Y")
y_t <- seq(1,10,0.01)
lines(y_t,2*pi*exp(-y_t),col=2)

1
この視点の難しさは、変換 fX(X) に依存する X、 とは対照的に FX(X)(次元1)。
西安
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.