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
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.