modele mélange de deux gaussiennes

[iam1899]

u1<-1
u2<-2
sig1<-1
sig2<-4
n<-100
x<-c(rnorm(n,u1,sig1))
z<-c(rnorm(n,u2,sig2))
w<-c(rbinom(n,1,pi))
for(i in 1:length(w)){
if w[i]==1 y[i]<-x[i]
else y[i]<-z[i]}

voilà on m'as demandé de faire une simulation pour un mélange de deux gaussiennes mais j y arrive pas , ça bloque dans la boucle for , svp si quelqu'un peut m'aider .. merci
-----

[minushabens]

Bonjour,

C'est bien quand on fournit un bout de code de préciser le langage. Il me semble reconnaître du S (le langage du logiciel R).

Si tu veux un mélange de n gaussiennes N(mu1,sigma1) et N(mu2,sigma2) avec les probabiités p et (1-p), une façon de faire sous R est:

melange<-function(n,p,mu1,sigma1,mu2,si gma2)
{
X1<-rnorm(n,mean=mu1,sd=sigma1)
X2<-rnorm(n,mean=mu2,sd=sigma2)
choix<-rbinom(n,size=1,prob=p)

return(choix*X1+(1-choix)*X2)
}

[iam1899]

merci beaucoup et navré j'aurai du en effet préciser le langage.. j'ai réussi a appliquer ceci p<-0.8
u1<-1
u2<-1
sig1<-1
sig2<-0.01
n<-1000
y<-rep(0,n)
x<-c(rnorm(n,u1,sig1))
z<-c(rnorm(n,u2,sig2))
w<-c(rbinom(n,1,p))
for(i in 1:length(w))
{
if (w[i]==1)
{
y[i]<-x[i]
}
else
{
y[i]<-z[i]}
}
plot(density(y),col=1)

maintenant je suis entrain d'essayer de réaliser l'algorithme em sur un mélange de gaussiennes .. si vous pouvez m'aider je vous serai très reconnaissant