23.12.2013 Views

Mélanges de GLMs et nombre de composantes : application ... - Scor

Mélanges de GLMs et nombre de composantes : application ... - Scor

Mélanges de GLMs et nombre de composantes : application ... - Scor

SHOW MORE
SHOW LESS

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Annexe D. Espace <strong>de</strong>s paramètres <strong>de</strong>s <strong>GLMs</strong><br />

Calcul <strong>de</strong> la limite :<br />

lim log L cc(ψ G ; y j )<br />

β i →−∞<br />

I<strong>de</strong>m que le cas <strong>de</strong> lim log L cc(ψ G ; y j ), sauf que les calculs montrent que les résultats sont<br />

β i →+∞<br />

⎧<br />

⎪⎨ si ν i X j y j > 0 alors lim A i = +∞,<br />

β<br />

inversés par rapport au signe <strong>de</strong> ν i X j y j . Ainsi,<br />

i →−∞<br />

⎪⎩ si ν i X j y j < 0 alors lim A i = 0.<br />

β i →−∞<br />

⎧<br />

⎪⎨ si ν k X j y j > 0 alors lim b i = 1,<br />

β<br />

Pour l’étu<strong>de</strong> <strong>de</strong> B i , nous avons :<br />

i →−∞<br />

⎪⎩ si ν k X j y j < 0 alors lim b i = +∞.<br />

β i →−∞<br />

⇒ Finalement,<br />

lim log L cc(ψ G ; y j ) = +∞.<br />

β i →−∞<br />

D.5 Mélange d’Inverses Gaussiennes<br />

D.5.1<br />

Décomposition <strong>de</strong> la log-vraisemblance L cc<br />

La log-vraisemblance classifiante conditionnelle log L cc (ψ G ; y j ) peut s’écrire pour une observation<br />

y j comme la somme <strong>de</strong> <strong>de</strong>ux termes :<br />

{ }} {<br />

A i<br />

{ }} √ {<br />

( G∑ π i<br />

log L cc (ψ G ; y j ) = log √ exp ( − 1 (y j − 1<br />

X j β i<br />

) 2<br />

) )<br />

i=1 2πσi 2y3 2<br />

1<br />

j<br />

X j β i<br />

σi 2y +<br />

j<br />

G∑<br />

i=1<br />

∑ G<br />

k=1<br />

π i<br />

√<br />

2πσ 2 i y3 j<br />

√<br />

π k<br />

⎛<br />

⎜<br />

exp ⎝− 1 2<br />

⎛<br />

⎜<br />

exp ⎝− 1 2<br />

(<br />

y j −<br />

√<br />

1<br />

A<br />

X j β i<br />

) 2<br />

1<br />

X j β i<br />

σ 2 i y j<br />

(<br />

y j −<br />

√<br />

1<br />

⎞<br />

⎟<br />

⎠<br />

X j β k<br />

) 2<br />

2πσk 2y3 1<br />

j<br />

X j β k<br />

σk 2y j<br />

} {{ }<br />

b i<br />

⎛<br />

⎞ log<br />

⎟ ⎜<br />

⎠ ⎝<br />

∑ G<br />

k=1<br />

π i<br />

√<br />

2πσ 2 i y3 j<br />

π k<br />

√<br />

2πσ 2 k y3 j<br />

⎛<br />

⎜<br />

exp ⎝− 1 2<br />

⎛<br />

⎜<br />

exp ⎝− 1 2<br />

(<br />

y j −<br />

√<br />

1<br />

X j β i<br />

) 2<br />

1<br />

X j β i<br />

σ 2 i y j<br />

(<br />

y j −<br />

√<br />

1<br />

⎞<br />

⎟<br />

⎠<br />

X j β k<br />

) 2<br />

1<br />

X j β k<br />

σ 2 k y j<br />

⎞<br />

⎞<br />

.<br />

⎟⎟<br />

⎠⎠<br />

} {{ }<br />

}<br />

B i =b i log b i<br />

{{ }<br />

B= P i B i<br />

D.5.2<br />

Calcul <strong>de</strong>s limites <strong>et</strong> <strong>de</strong>s cas critiques<br />

Calcul <strong>de</strong> la limite :<br />

lim log L cc(ψ G ; y j )<br />

β i →+∞<br />

→ Nous remarquons par un simple calcul que lim A i = 0, donc lim A = K.<br />

β i →+∞ β i →+∞<br />

→ De plus,<br />

232<br />

lim B i =<br />

β i →+∞<br />

lim b i log b i =<br />

β i →+∞<br />

lim<br />

β i →+∞<br />

A i<br />

A log A i<br />

A = 0 car A i → 0 <strong>et</strong> A → K.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!