
/////////////////////////////////////////////////////////////
//TEST du RAPPORT de VRAISEMBLANCE DANS UN MODELE DE MELANGE
//Auteur: F. Richard
//Date: 24/04/2006
/////////////////////////////////////////////////////////////


//Cette fonction calcule la classe de chq individu echantillone.

//Modele statistique:
// Y=(Yk)k variables aleatoires iid selon la loi de densite 
// f(y) =  p f0(y) + (1-p) f1(y)
// avec:
// - p et (1-p) probabilites que l'individu appartienne aux classes 0 et 1, respectivement.
// - f0 densite de loi sur la classe 0 (gaussienne de moyenne nulle et de variance sig0)
// - f1 densite de loi sur la classe 1 (gaussienne de moyenne m1 et de variance sig1)



//PARAMETRES D'ENTREE
///////////////////////

//Y: echantillon dont la loi de densite est f.
//sig0, m1, sig1, p: parametres du modele

//EN SORTIE
////////////

//classes: estimation des classes de chq individu



function [classes]=test_rapport_vraisemblance(Y,sig0,moy1,sig1,p),

a=0.5*(1/sig0-1/sig1);
b=1/sig1;

seuil=moy1^2/sig1*0.5-log(p*sqrt(sig0)/((1-p)*sqrt(sig1)))

classes=zeros(Y);
classes(a*Y.^2+b*Y*moy1>seuil)=1;

endfunction;

