12.07.2015 Views

Télécharger le texte intégral - ISPED-Enseignement à distance

Télécharger le texte intégral - ISPED-Enseignement à distance

Télécharger le texte intégral - ISPED-Enseignement à distance

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Modè<strong>le</strong> nonlinéaire <strong>à</strong> classes latentes 2 1185.2.2 Troncature <strong>à</strong> gaucheIl est très fréquent dans une étude de cohorte sur <strong>le</strong> vieillissement que <strong>le</strong>s sujetssoient inclus <strong>à</strong> la condition qu’ils n’aient pas subi l’événement avant l’âge d’entréedans l’étude. Pour tenir compte de cette troncature <strong>à</strong> gauche, <strong>le</strong>s paramètres doiventêtre estimés par maximisation de la vraisemblance conditionnel<strong>le</strong> <strong>à</strong> la troncature.La contribution individuel<strong>le</strong> <strong>à</strong> la vraisemblance L i doit donc être divisée par laprobabilité de survie jusqu’<strong>à</strong> l’inclusion dans l’étude T 0i :donc :L (2)i (θ) = L i(θ)S(T 0i ; θ)où S(T 0i ; θ) est la fonction de survie margina<strong>le</strong> en T 0i définie comme suit :S(T 0i ; θ) =G∑P (c ig = 1)f(T i ≥ T 0i | c ig = 1; θ) =g=1G∑π ig S 0g (T 0i ; θ) exp(X 4iδ g)g=1(5.10)(5.11)En cas de données tronquées <strong>à</strong> gauche, la log-vraisemblance <strong>à</strong> maximiser devientL (2) (θ) = L(θ) −(N∑ G)∑ln π ig e −(ζ 1gT 0i ) ζ 2g eX 4i δgi=1g=1(5.12)5.2.3 Classification et probabilité a posteriori de l’événementA partir des estimations ˆθ, nous pouvons calcu<strong>le</strong>r <strong>le</strong>s probabilités a posteriorid’appartenance aux classes latentes sachant <strong>le</strong>s données :- la probabilité ˆπ y,Tigd’appartenir <strong>à</strong> chaque classe latente g sachant <strong>à</strong> la fois <strong>le</strong>smesures répétées aux tests y i , <strong>le</strong>s variab<strong>le</strong>s explicatives x i et <strong>le</strong> délai jusqu’<strong>à</strong>l’événement ou jusqu’<strong>à</strong> la sortie d’étude T i :ˆπ y,Tig= P (c ig = 1 | y i , (T i , E i ), x i ; ˆθ)= P (c ig = 1 | x i ; ˆθ)f(y i , (T i , E i ) | c ig = 1, x i ; ˆθ)∑ Gl=1 P (c il = 1 | x i ; ˆθ)f(y i , (T i , E i ) | c il = 1, x i ; ˆθ)(5.13)

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

Saved successfully!

Ooh no, something went wrong!