
/////////////////////////////////////////////////////////////
//ESTIMATEUR DU MAXIMUM DE VRAISEMBLANCE
//Auteur: F. Richard
//Date: 24/04/2006
/////////////////////////////////////////////////////////////


//Cette fonction permet de calculer le maximum de vraisemblance des parametres du modele suivant.

// 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
//classes: classes des individus.


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

//moy0, sig0, moy1, sig1, p: estimations des parametres du modele



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

taille=size(Y,"c")*size(Y,"r");
//estimation des moyennes sur chq classe
moy0=mean(Y(classes==0));
moy1=mean(Y(classes==1));
//estimation des variances sur chq classe
sig0=st_deviation(Y(classes==0))^2;
sig1=st_deviation(Y(classes==1))^2;
//taille (absolue et relative) de la classe 0 dans l'image
taille0=sum(classes);
p=taille0/(taille)^2;

endfunction;

