j'ai peut être une solution. la manière la plus élégante de programmer ce modèle est de le voir comme un cas particulier de régression logistique. j'ai d'ailleurs trouver des sources pour Matlab qui vont dans ce sens. Un package trés complet est egalement disponible pour le logiciel libre R, avec une solide documentation à l'appui.