On Mon, 7 Mar 2005 09:35:39 +0100 (CET), [EMAIL PROTECTED] wrote : >Hi, > >I got the same problem like >http://tolstoy.newcastle.edu.au/R/devel/04/11/1204.html > >R crashes when I use the em function from the mclust package on >univariate data and on a special case on bivariate data (when the matrix >is not provided as written in the manual). >It seems as if the problem is the format of the data to be analyzed. > >Operating System: Windows XP (SP2) >R version: R-2.0.1 > >The following example causes a crash of R: > ># univariate example >require(mclust) >X <- c(rnorm(100, mean=1), rnorm(100, mean=5)) ># I also tried: X <- as.matrix(c(rnorm(100, mean=1), rnorm(100, mean=5))) >Xmap <- cbind(c(rep(1, 100), rep(0, 100)), c(rep(0, 100), rep(1, 100))) >Xm <- mstep(modelName="V", data=X, z=Xmap) ># CRASH: >em(modelName=Xm$modelName, data=X, mu=Xm$mu, sigmasq=Xm$sigmasq, pro=Xm$pro)
Which version of mclust? I just installed 2.1-8 and ran your code in R-2.0.1 with no problem. In any case, this looks more like an mclust problem than an R problem; if you're using the current version, you should follow up with the maintainer, Ron Wehrens <[EMAIL PROTECTED]>. Duncan Murdoch ______________________________________________ R-devel@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
