bonjour,
je cherche à faire une régression (ici logistique) avec des données génotypiques (peuvent prendre la valeur 0,1 ou 2 avec les probabilités P0, P1 et P2).
le modèle est log(p/1-p)=a+b1*G1+b2*G2+b3*G1*G2
En sachant que G1 et G2 ne sont pas toujours des observations directes mais sont données sous forme de distribution de proba.
Pour individu 1, par exemple, on a
P(G1=0)=0.5,P(G1=1)=0.3,P(G1=2)=0.2 etc ....
Existe-t-il une fonction permettant de boucler sur ces probabilités des variables indépendantes et de donner le bon test LRT ?
En vous remerciant,
Christian