Distribuição de

Digamos, (com ) tem uma densidade . O que podemos dizer sobre a distribuição de $X \in \mathbb{R}^n$$n > 1$$f_X(x)$

$$Y = -\log f_X(X)?$$

O livro mencionado por Xi'an é de 2004. Refere-se a um artigo do ano de 1991 no qual o seguinte teorema aparece.

De: Troutt MD 1991 Um teorema da densidade das ordenadas de densidade e uma interpretação alternativa do método Box-Muller

Se uma variável aleatória X tiver uma densidade , , e se a variável aleatória tiver uma densidade , então onde é a medida de Lebesgue do conjunto$f(x)$$x \in \mathbb{R}^n$$v = f(x)$$g(v)$

$$g(v) = -vA^\prime(v),$$ $A(v)$ $$S(v) = \lbrace x: f(x) \geq v \rbrace$$

Intuitivamente e não formal:

$$\begin{array}\\ f_Z(z) dz = P(z<Z<z+dz) &= P(x(z)<X<x(z+dz)) \\ &= P(x(z)<X<x(z)+dz \frac{dx}{dz}) \\ &= f_X(X) \frac{dx}{dz} dz = z \frac{-dA(z)}{dz} dz \end{array}$$

De maneira semelhante, quando usamos uma variável transformada então:$Y = g(f_x(x))$

$$\begin{array}\\ f_Y(y) dy = P(y<Y<y+dy) &= P(x(y)<X<x(y+dy)) \\ &= P(x(y)<X<x(y)+dy \frac{dx}{dy}) \\ &= f_X(X) \frac{dx}{dy} dy = g^{-1}(y) \frac{-dA(y)}{dy} dy \end{array}$$

assim

$$f_Y(y) = -e^{-y} \frac{A(y)}{dy}$$

---

exemplo de distribuição normal padrão:

$$f_X(x) = \frac{1}{\sqrt{2\pi}} e^{-0.5 x^2}$$

$$y = \log(\sqrt{2\pi}) + 0.5 x^2$$

$$A(y) = C-\sqrt{8(y-\log(\sqrt{2\pi}))}$$

portanto

$$f_Y(y) = \frac{\sqrt{2} e^{-y}}{\sqrt{y-\frac{\log(2\pi)}{2}}}$$

---

exemplo uma distribuição normal multivariada:

$$f_X(x_1,x_2) = \frac{1}{2\pi} e^{-0.5 (x_1^2 + x_2^2)}$$

$$y = \log(2\pi) + 0.5 (x_1^2+x_2^2)$$

$$A(y) = C-2\pi(y-\log(2\pi))$$

portanto

$$f_Y(y) = 2\pi e^{-y} \qquad \qquad \text{for $y \geq log(2\pi)$}$$

verificação computacional:

# 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)