O teorema de Bayes é válido para expectativas?

18

É verdade que, para duas variáveis ​​aleatórias e ,AB

E(AB)=E(BA)E(A)E(B)?
tomka
fonte
3
Hmm ... eu não acho que esses dois lados são equivalentes
Jon
6
Como apontado nas respostas, a pergunta é probabilisticamente sem sentido, devido à integração de variáveis ​​aleatórias de um lado, que são as variáveis ​​condicionantes do outro lado.
Xian

Respostas:

25

(1)E[AB]=?E[BA]E[A]E[B]
O resultado conjecturado é trivialmente verdadeiro para variáveis ​​aleatórias independentes e com médias diferentes de zero.(1)AB

Se E[B]=0 , o lado direito de (1) envolve uma divisão por 0 e, portanto, (1) não faz sentido. Observe que se A e B são independentes não é relevante.

Em geral , não se aplica a variáveis ​​aleatórias dependentes, mas exemplos específicos de e dependentes que satisfazem podem ser encontrados. Note que devemos continuar insistindo que , caso contrário, o lado direito de não tem sentido. Tenha em mente que é uma variável aleatória que é uma função da variável aleatória , digamos enquanto é uma variável aleatória que é uma função do variável aleatória , digamosA B ( 1 ) E [ B ] 0 ( 1 ) E [ A B ] B g ( B ) E [ B A ] A h ( A ) ( 1 )(1)AB(1)E[B]0(1)E[AB]Bg(B)E[BA]Ah(A) . Então, é semelhante a perguntar se(1)

g(B)h(A)

(2)g(B)=?h(A)E[A]E[B]
pode ser uma afirmação verdadeira, e obviamente a resposta é que não pode ser uma múltiplo de em geral.g(B)h(A)

Que eu saiba, existem apenas dois casos especiais em que pode ser mantido.(1)

  • Como notado acima, para independentes variáveis aleatórias e , e são degeneradas variáveis aleatórias (chamados constantes por pessoas estatisticamente-analfabetos) que igual e , respectivamente, e por isso, se , temos igualdade em .B g ( B ) h ( A ) E [ A ] E [ B ] E [ B ] 0 ( 1 )ABg(B)h(A)E[A]E[B]E[B]0(1)

  • No outro extremo do espectro da independência, suponha que onde é uma função invertível e, portanto, e sejam totalmente variáveis ​​aleatórias dependentes. Nesse caso, e então torna-se que vale exatamente quando onde pode ser qualquer número real diferente de zero. Assim, é válido sempre que é um múltiplo escalar de e, é claro,g ( ) A = g ( B ) B = g - 1 ( A ) E [ A B ] = g ( B ) ,A=g(B)g()A=g(B)B=g1(A)( 1 ) g ( B ) ? = B E [ A ]

    E[AB]=g(B),E[BA]=g1(A)=g1(g(B))=B
    (1) g(x)=αxα(1)ABE[B]
    g(B)=?BE[A]E[B]
    g(x)=αxα(1)ABE[B]deve ser diferente de zero (cf. resposta de Michael Hardy ). O desenvolvimento acima mostra que deve ser uma função linear e que não pode conter funções afins com . No entanto, observe que Alecos Papadopolous em sua resposta e em seus comentários posteriores afirma que se é uma variável aleatória normal com média diferente de zero, então para valores específicos de e que ele fornece, e satisfazemg(x)(1)β 0g(x)=αx+ββ0Bβ 0 A = α B + β B ( 1 )αβ0A=αB+βB(1) . Na minha opinião, o exemplo dele está incorreto.

Em um comentário sobre esta resposta, Huber sugeriu considerar a igualdade conjecturada simétrica quais dos naturalmente sempre é válido para as variáveis aleatórias independentes independentemente dos valores de e e para múltiplos escalares também. Naturalmente, mais trivialmente, vale para quaisquer variáveis ​​aleatórias com média zero e (independente ou dependente, múltiplo escalar ou não; isso não importa!): é suficiente pela igualdade em . Assim, pode não ser tão interessante quanto E[A]E[B]A=αB(3)

(3)E[AB]E[B]=?E[BA]E[A]
E[A]E[B]A=αB(3) ABE[A]=E[B]=0 (3)(3)(1) como um tópico para discussão.
Dilip Sarwate
fonte
9
+1. Para ser generosa, a questão poderia ser interpretada como perguntando se , onde a questão da divisão por zero desaparece. E(A|B)E(B)=E(B|A)E(A)
whuber
1
@whuber Obrigado. Minha edição aborda a questão mais geral sobre se é possível ter . E[AB]E[B]=E[BA]E[A]
Dilip Sarwate
11

O resultado é falso em geral, vamos ver isso em um exemplo simples. Seja uma distribuição binomial com os parâmetros e com a distribuição beta com os parâmetros , ou seja, um modelo bayesiano com conjugado anterior. Agora apenas calcule os dois lados da sua fórmula, o lado esquerdo é , enquanto o lado direito é e esses certamente não são iguais.XP=pn,pP(α,β)EXP=nP

E(PX)EXEP=α+Xn+α+βα/(α+β)nα/(α+β)
kjetil b halvorsen
fonte
2

O valor condicional esperado de uma variável aleatória dado o evento que B = b é um número que depende de qual número b é. Então, chame-o de h ( b ) . Em seguida, o valor esperado condicional E ( A | B ) é H ( B ) , uma variável aleatória cujo valor é completamente determinada pelo valor da variável aleatória B . Assim, E ( A B ) é uma função de B e E (AB=bbh(b).E(AB)h(B),BE(AB)B é uma função da A .E(BA)A

O quociente é apenas um número.E(A)/E(B)

Portanto, um lado da igualdade proposta é determinado por e o outro por B , portanto, geralmente não podem ser iguais.AB

(Talvez eu deva acrescentar que eles podem ser iguais no caso trivial quando os valores de e B se determinam, como quando, por exemplo, A = α B , α 0 e E [ B ] 0 , quando E [ A B ] = α B = E [ B A ] α = E [ B A ] α E [ B ]ABA=αB,α0E[B]0Mas funções iguais entre si apenas em alguns pontosnãosãoiguais.)

E[AB]=αB=E[BA]α=E[BA]αE[B]E[B]=E[BA]E[A]E[B].
Michael Hardy
fonte
Você quer dizer que eles não são necessariamente iguais? Quero dizer que eles podem ser iguais?
BCLC
1
@BCLC: Eles são iguais apenas em casos triviais. E duas funções iguais entre si em alguns pontos e não em outros não são iguais.
Michael Hardy
2
"Mas somente nesse caso trivial eles podem ser iguais" (grifo nosso) não está correto. Considere independentes e B com E [ B ] 0 . Então, E [ A | B ] = E [ A ] , enquanto E [ B | A ] = E [ B ] e por isso E [ B | A ] E [ A ]ABE[B]0E[AB]=E[A]E[BA]=E[B]
E[BA]E[A]E[B]=E[B]E[A]E[B]=E[A]=E[AB].
Dilip Sarwate
@DilipSarwate Eu estava prestes a dizer isso haha!
BCLC 10/02
Editei sua resposta para adicionar alguns detalhes ao caso que você apontou. Reverta se não gosta das alterações.
precisa
-1

A expressão certamente não se aplica em geral. Por diversão, mostro abaixo que, se e B seguem conjuntamente uma distribuição normal bivariada e têm médias diferentes de zero, o resultado será mantido se as duas variáveis ​​forem funções lineares uma da outra e tiverem o mesmo coeficiente de variação ( a razão do desvio padrão sobre a média) em termos absolutos.AB

Para normais comuns, temos

E(AB)=μA+ρσAσB(BμB)

e queremos impor

μA+ρσAσB(BμB)=[μB+ρσBσA(AμA)]μAμB

μA+ρσAσB(BμB)=μA+ρσBσAμAμB(AμA)

Simplifique e depois ρ e reorganize para obterμAρ

B=μB+σB2σA2μAμB(AμA)

Portanto, esta é a relação linear que deve ser mantida entre as duas variáveis ​​(para que elas sejam certamente dependentes, com coeficiente de correlação igual à unidade em termos absolutos) para obter a igualdade desejada. O que isso implica?

Primeiro, também deve estar satisfeito que

E(B)μB=μB+σB2σA2μAμB(E(A)μA)μB=μB

portanto, nenhuma outra restrição é imposta à média de (ou de A ), exceto por serem diferentes de zero. Também uma relação para a variação deve ser satisfeita,BA

Var(B)σB2=(σB2σA2μAμB)2Var(A)

(σA2)2σB2=(σB2)2σA2(μAμB)2

(σAμA)2=(σBμB)2(cvA)2=(cvB)2

|cvA|=|cvB|

which was to be shown.

Note that equality of the coefficient of variation in absolute terms, allows the variables to have different variances, and also, one to have positive mean and the other negative.

Alecos Papadopoulos
fonte
1
Isn't this a convoluted way to A=αB where α is some scalar?
Matthew Gunn
1
@MatthewGunn Your comment is right on target. Normality has nothing to do with the matter. For random variables A and B such that A=αB, E[AB]=αB=A and similarly, E[BA]=B. Consequently, assuming that E[B]0,
E[AB]=αB=E[BA]α=E[BA]αE[B]E[B]=E[BA]E[A]E[B].
No normality, no |cvA|=|cvB| etc, and actually just a rehash of a comment in Michael Hardy's answer.
Dilip Sarwate
If you write \text{Var} instaed of \operatorname{Var} then you'll see aVarX and aVar(X) instead of aVarX and aVar(X). That's why the latter is standard usage.
Michael Hardy
@MatthewGun It seems to me that providing answers that contain specific examples is considered valuable content in this site. So yes, when a random variable is an affine function of another, and they are jointly normal with non-zero means, then one needs to have equal coefficients of variation, while, also there are no restrictions on the means of these rv's. On the other hand, when a random variable is just a linear function of another, the relation holds always. So no my answer is not a convoluted way to say A=aB. (cc:@DilipSarwate)
Alecos Papadopoulos
2
If B is a non-normal random variable with E[B]=μB0 and A=cB+d (and so B=Adc), then
E[AB]=cB+d=A,E[BA]=Adc=B.
Now, if we want to have E[AB]=cB+d to equal E[BA]μAμB=BμAμB, it must be that
cB+d=BμAμBd=0,c=μAμB
and so A=cB=μAμBB. So, for nonnormal B, the OP's conjectured result holds if A=cB but not if A=cB+d,d0.Of course, as you have proved, the result holds for normal random variables if A=cB+d,d0 .
Dilip Sarwate