Hi,
I am trying to understand how to use the flexmix package, I have read the Leisch paper but am very unclear what is needed for the M-step driver. I am just fitting a simple linear regression model. The documentation is far from clear what the FLXmclust function does, but, it in principle could do all I need, however, I do not get sentible results as if I try the following the result is poor:
>x<-c()
> for(i in 0:99){x$y[2*i]=(0+i);x$x[2*i]=i;x$x[2*i+1]=i;x$y[2*i+1]=i+1000;x$g[2*i]=1;x$g[2*i+1]=2}
> m1<-flexmix(y~x ,data=x,k=2)
> table(x$g,m1@cluster)
1 2
1 25 74
2 67 33
It all depends on the randomised starting values. So I think I need a better "driver", but, I cannot find a spec for what I have to do in the driver.
Where is FLXmclust documented? can anyone assist?
[[alternative HTML version deleted]]