Como provar essa desigualdade da mistura gaussiana? (Montagem / sobreajuste)


8

Seja f [x] um pdf de mistura gaussiana com n termos de peso uniforme, meios e variações correspondentes :{μ1,...,μn}{σ1,...,σn}

f(x)≡1n∑i=1n12πσi2e−(x−μi)22σi2

Parece intuitivo que a probabilidade de log amostrada nos n centros gaussianos seja maior que (ou igual a) a probabilidade de log média:

1n∑j=1nln(f(μj))≥∫f(x)ln(f(x))dx

Isto é obviamente verdadeiro para pequenas variações (cada está no topo de um gaussiano estreito) e para variações muito grandes (todas as estão no topo de um amplo gaussiano juntos), e isso é verdade para todos os conjuntos de e eu e , mas não consigo descobrir como provar que isso sempre é verdade. Socorro?μiμiμiσEu


Você provavelmente está perdendo uma expectativa no LHS?
— lacerbi

@lacerbi Não, eu não sou. Nada está faltando. Nas LHS, o é avaliada no indexada 'sx if(x)xi
— Jerry Guern

Sim, desculpe - eu estava com muito sono e interpretei mal a definição.
— 181116 lacerbi

Respostas:


2

Este é um comentário mais extenso, portanto, tome-o como tal. Defina: (estou usando a notação padrão para distribuições gaussianas).

f(x)≡1n∑i=1nN(x|xi,σi2)

Você quer provar que: que é {1

1n∑i=1nlog⁡f(xi)−∫f(x)log⁡f(x)dx≥0
{1n∑i=1nlog⁡f(xi)}+H[f]≥0.

Devido à desigualdade de Jensen (ver, por exemplo, Huber et al., On Entropy Approximation for Gaussian Mixture Random Vectors, 2008 ), comgi(x)≡1

H[f]≥−1n∑i=1nlog⁡∫f(x)N(x|xi,σi2)dx=−1n∑i=1nlog⁡gi(xi)
, que vem da convolução de duas densidades gaussianas. Então temos: {1gEu(x)≡1n∑j=1nN(x|xj,σEu2+σj2) Curiosamente, ogiainda são misturas de Gaussianas com meios de componentes iguais aos daf, mas cada componente degitem uma variância estritamente maior do que o seu correspondente componente emf. Você pode fazer alguma coisa com isso?
{1n∑Eu=1nregistro⁡f(xEu)}+H[f]≥1n∑Eu=1nregistro⁡f(xEu)gEu(xEu).
gEufgEuf

Obrigado. Parece que eu teria que provar que o RHS final é> = 0, o que também parece intuitivo, mas difícil de provar, mas esse é realmente um passo na direção certa. Eu já vi esse jornal antes.
— precisa

É tentador pensar que o RHS final é sempre positivo, mas também não posso provar isso.
— Jerry Guern

0

Eu acho que entendi. Leva apenas etapas elementares, embora você precise combiná-las corretamente.

fEuEu12πσEu2e(x-μEu)22σEu2

g(x)=xeuog(x)f(x)registro⁡(f(x))≤1n∑Eu=1nfEu(x)registro⁡(fEu(x))

∫f(x)registro⁡(f(x))dx≤1n∑Eu=1n∫fEu(x)registro⁡(fEu(x))dx

Euf≥fEu

euog(f(μEu))≥euog(fEu(μEu))
1n∑Eu=1neuog(f(μEu))≥1n∑Eu=1neuog(fEu(μEu))
1n∑Eu=1neuog(fEu(μEu))≥1n∑Eu=1nfEu(x)registro⁡(fEu(x))
euog(fEu(μEu))=∫fEu(x)euog(fEu(μEu))dx≥∫fEu(x)euog(fEu(x))dx
Eun

Estou confuso. Você definiu ag (x), mas nunca o usou, e eu não sei o que seu f_i significa.
— Jerry Guern

fEugg(1n∑Eu=1nfEu(x))≤1n∑Eu=1ng(fEu(x))

f> =fEu1/nfEu

1
1n∑Eu=1neuog(f(μEu))≥1n∑Eu=1neuog(fEu(μEu))

1
Sim, eu percebi isso ontem. Parece que essa desigualdade é bastante difícil, deixarei minha resposta de qualquer maneira com uma edição.
— Sjm.majewski
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.