Distribuição da razão entre duas variáveis ​​aleatórias uniformes independentes


17

Supppse X e Y são padrão distribuídos uniformemente em [0,1] e são independentes, qual é o PDF de Z=Y/X ?

A resposta de algum livro de teoria da probabilidade é

fZ(z)={1/2,if 0z11/(2z2),if z>10,otherwise.

Estou pensando, por simetria, não deveria ? Este não é o caso, de acordo com o PDF acima.fZ(1/2)=fZ(2)


Qual é o domínio de X e Y ?
Sobi 08/12/2015


2
Por que você espera que isso seja verdade? A função de densidade diz-lhe como bem embalado a probabilidade é na vizinhança de um ponto, e é claramente mais difícil para Z para estar perto de 2 de 1/2 (considere, por exemplo, que Z pode ser sempre 1/2 não importa o X é, mas Z<2 quando X>1/2 ).
dsaxton


3
Não acho que seja uma cópia duplicada, essa pergunta está buscando o PDF, aqui estou o PDF, estou apenas questionando sua correção (talvez um tanto ingênua).
Qd #

Respostas:


19

A lógica correta é que, com X,YU(0,1) , Z=YX e Z1=XY tem a mesmadistribuiçãoe, portanto, para0<z<1

P{YXz}=P{XYz}=P{YX1z}FZ(z)=1FZ(1z)
onde a equação com CDFs usa o fato de queYX é uma variável aleatória contínua e, portanto,P{Za}=P{Z>a}=1FZ(a). Portanto, o pdf deZsatisfaz Assim, f Z ( 1
fZ(z)=z2fZ(z1),0<z<1.
e nãofZ(1fZ(12)=4fZ(2) como você pensou que deveria ser.fZ(12)=fZ(2)

14

Essa distribuição é simétrica - se você a observar da maneira certa.

A simetria que você observou (corretamente) é que e X / Y = 1 / ( Y / X ) devem ser distribuídos de forma idêntica. Ao trabalhar com relações e potências, você está realmente trabalhando no grupo multiplicativo dos números reais positivos. O análogo da medida invariante da localização d λ = d x nos números reais aditivos R é a medida invariante da escala d μ = d x / xY/XX/Y=1/(Y/X)dλ=dxR dμ=dx/x do grupo multiplicativo de números reais positivos. Possui estas propriedades desejáveis:R

  1. é invariável sob a transformação x a x para qualquer constante positiva a : d μ ( a x ) = d ( a x )dμxaxa

    dμ(ax)=d(ax)ax=dxx=dμ.
  2. é covariante sob a transformação x x b para números diferentes de zero b : d μ ( x b ) =dμxxbb

    dμ(xb)=d(xb)xb=bxb1dxxb=bdxx=bdμ.
  3. é transformado em d λ através do exponencial: d μ ( e x ) =dμdλ Da mesma forma,dλé transformado novamente emdμatravés do logaritmo.

    dμ(ex)=dexex=exdxex=dx=dλ.
    dλdμ

(3) estabelece um isomorfismo entre os grupos medidos e ( R , , d μ ) . A reflexão x - x no espaço aditivo corresponde à inversão x 1 / x no espaço multiplicativo, porque e - x =(R,+,dλ)(R,,dμ)xxx1/x .ex=1/ex

Vamos aplicar essas observações escrevendo o elemento de probabilidade de em termos de d μ (entendendo implicitamente que z > 0 ) em vez de d λ :Z=Y/Xdμz>0dλ

fZ(z)dz=gZ(z)dμ=12{1dz=zdμ,if 0z11z2dz=1zdμ,if z>1.

Ou seja, o PDF com relação à medida invariante dμ é , proporcional a z quando 0 < z 1 e a 1 / z quando 1 z , próximo ao que você esperava.gZ(z)z0<z11/z1z


Este não é um mero truque pontual. Compreender o papel de faz com que muitas fórmulas pareçam mais simples e mais naturais. Por exemplo, o elemento de probabilidade da função Gamma com o parâmetro k , x k - 1 e xdμk se torna x k e x d μ . É mais fácil trabalhar com d μ do que com d λ ao transformar x , redimensionando, obtendo poderes ou exponenciando.xk1exdxxkexdμdμdλx

A idéia de uma medida invariável em um grupo também é muito mais geral e tem aplicações nessa área de estatística em que os problemas exibem alguma invariância em grupos de transformações (como alterações de unidades de medida, rotações em dimensões mais altas etc.) )


3
Parece uma resposta muito perspicaz. É uma pena que eu não o entenda no momento. Voltarei mais tarde.
qed

4

Se você pensa geometricamente ...

No plano - Y , curvas da constante Z = Y / X são linhas através da origem. ( Y / X é a inclinação.) É possível ler o valor de Z de uma linha através da origem, encontrando sua interseção com a linha X = 1 . (Se você já estudou o espaço projetivo: aqui X é a variável de homogeneização, observar os valores na fatia X = 1 é uma coisa relativamente natural a se fazer.)XYZ=Y/XY/XZX=1XX=1

Considere um pequeno intervalo de s, ( a , b ) . Esse intervalo também pode ser discutido na linha X = 1 como o segmento de linha de ( 1 , a ) a ( 1 , b ) . O conjunto de linhas através da origem que passa por esse intervalo forma um triângulo sólido no quadrado ( X , Y ) U = [ 0 , 1 ] × [ 0 , 1 ]Z(a,b)X=1(1,a)(1,b)(X,Y)U=[0,1]×[0,1], qual é a região na qual realmente estamos interessados. Se , a área do triângulo é 10a<b1, portanto, mantendo o comprimento do intervalo constante e deslizando-o para cima e para baixo na linhaX=1(mas não após0ou1), a área é a mesma, portanto, a probabilidade de escolher um(X,Y)no triângulo é constante; portanto, a probabilidade de escolher umZno intervalo é constante.12(10)(ba)X=101(X,Y)Z

No entanto, para , o limite da região U se afasta da linha X = 1 e o triângulo é truncado. Se 1 a < b , as projeções abaixo das linhas desde a origem de ( 1 , a ) e ( 1 , b ) até o limite superior de U estão nos pontos ( 1 / a , 1 ) e ( 1 / b , 1 )b>1UX=11a<b(1,a)(1,b)U(1/a,1)(1/b,1). A área resultante do triângulo é . A partir disso, vemos que a área não é uniforme e, à medida que deslizamos(a,b)cada vez mais para a direita, a probabilidade de selecionar um ponto no triângulo diminui para zero.12(1a1b)(10)(a,b)

Então a mesma álgebra demonstrada em outras respostas termina o problema. Em particular, retornando a última pergunta do OP, corresponde a uma linha que atinge X = 1 , mas f Z ( 2 ) não faz, então a simetria desejada não se sustenta.fZ(1/2)X=1fZ(2)


3

Só para constar, minha intuição estava totalmente errada. Estamos falando de densidade , não de probabilidade . A lógica correta é verificar se

1kfZ(z)dz=1/k1fZ(z)=12(11k)
,

and this is indeed the case.


1

Yea the link Distribution of a ratio of uniforms: What is wrong? provides CDF of Z=Y/X. The PDF here is just derivative of the CDF. So the formula is correct. I think your problem lies in the assumption that you think Z is "symmetric" around 1. However this is not true. Intuitively Z should be a skewed distribution, for example it is useful to think when Y is a fixed number between (0,1) and X is a number close to 0, thus the ratio would be going to infinity. So the symmetry of distribution is not true. I hope this help a bit.

Ao utilizar nosso site, você reconhece que leu e compreendeu nossa Política de Cookies e nossa Política de Privacidade.
Licensed under cc by-sa 3.0 with attribution required.