Please read the FAQ and report this in the correct place: Bug reports on contributed packages should be sent first to the package maintainer, and only submitted to the R-bugs repository by package maintainers, mentioning the package in the subject line.
You are not the named maintainer. On Mon, 7 Mar 2005 [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) > > For bivariate data the following example works: > > # bivariate example > require(mclust) > Y <- as.matrix(cbind(c(rnorm(100, mean=1), rnorm(100, mean=5)), > c(rnorm(100, mean=1), rnorm(100, mean=5)))) > Ymap <- cbind(c(rep(1, 100), rep(0, 100)), c(rep(0, 100), rep(1, 100))) > Ym <- mstep(modelName="EEE", data=Y, z=Ymap) > # NO CRASH: > em(modelName=Ym$modelName, data=Y, mu=Ym$mu, sigmasq=Ym$Sigma, pro=Ym$pro) > > But when the variables are not columns (as written in the mclust manual) > but rows, R crashes again: > > # CRASH: > em(modelName=Ym$modelName, data=t(Y), mu=Ym$mu, sigmasq=Ym$Sigma, > pro=Ym$pro) > > Perhaps this is a hint to the origin of the problem. > I could not use the em function for univariate data successfully yet. > Perhaps anyone can provide a working example for univariate data? > > Thanks in advance. > > ______________________________________________ > R-devel@stat.math.ethz.ch mailing list > https://stat.ethz.ch/mailman/listinfo/r-devel > -- Brian D. Ripley, [EMAIL PROTECTED] Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595 ______________________________________________ R-devel@stat.math.ethz.ch mailing list https://stat.ethz.ch/mailman/listinfo/r-devel
