Bonjour,


j'aimerais utiliser l'algo EM de mclust sur le cas le plus simple (mélange de 4 gaussiennes en 4 dimensions) et obtenir les 4 moyennes, les 4 variances et les 4 proportions.
je suppose que la fonction "em" effectue l'algo sur un certain nombre d'itérations (combien ? ...). Mais lorsque je place le code suivant, il me dit que la parameters$pro est incorrect (j'ai mis rep(.25,4)...)

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
 
N=1000
w=sample(1:4, N, replace=TRUE,prob=c(.1,.5,.15,.15) )
echantillon=(w==1)*rnorm(N,-5,1)+(w==2)*rnorm(N,10,4)+(w==3)*rnorm(N,0,.5)+(w==4)*rnorm(N,5,10)
 
 
 ech=as.matrix(echantillon)
 par=data.frame(pro=rep(.25,4),mean=rep(1,4),variance=rep(1,4))
 em(modelName="EEE", data=ech,parameters=par)
 
 estep(modelName="EEE", data=ech,parameters=par)

Si quelqu'un peut m'expliquer ce qui se passe, ce serait super.
Merci par avance.