similar to: A Few MCLUST Questions

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

I keep on receiving the message below after submitting the following line using the mclust package. m2 is a 99 X 1 column vector. * em(modelName = "E", m2, mu = c(25, 50), sigmasq=10, pro = c(0.4, 0.6)) Error in as.double.default(data) : (list) object cannot be coerced to double. Why do I receive this error? Thank, Brian C. Newquist Research Statistician

Hello All... I gave a task to my students that involved using mclust to look for clusters in some bivariate data of isotopes vs various mining locations. They discovered something I didn?t expect; the data (called tur) is appended below. p <- qplot(x = dD, y = dCu65, data = tur, color = mine) print(p) # simple bivariate plot of the data; looks fine mod1 <- Mclust(tur[,2:3]) mod1$G mod2

Hi, I'm trying to use the library mclust for gaussian mixture on a numeric vector. The function Mclust(data,G=3) is working fine but the fitting is not optimal and is using modelNames="E". When I'm trying Mclust(data,G=3,modelName="V") I have the following message: Error in if (Sumry$G > 1) ans[c(orderedNames, "z")] else ans[orderedNames] : argument is

Hi, I'm trying to use the library mclust for gaussian mixture on a numeric vector. The function Mclust(data,G=3) is working fine but the fitting is not optimal and is using modelNames="E". When I'm trying Mclust(data,G=3,modelName="V") I have the following message: Error in if (Sumry$G > 1) ans[c(orderedNames, "z")] else ans[orderedNames] : argument is

Hi All, I need to fit a mutlivariate normal model to a dataset in order to obtain the mean and covariance parameters. I see that the MVN function in the MCLUST package can do this, however when I try to run even the simplest example provided in the documentation, as below, I get the following error: n <- 1000 set.seed(0) x <- rnorm(n, mean = -1, sd = 2) mvn(modelName = "X", x)

Dear list members, I am using Mclust in order to deconvolute a distribution that I believe is a sum of two gaussians. First I can make a model: > my.data.model = Mclust(my.data, modelNames=c("E"), warn=T, G=1:3) But then, when I try to plot the result, I get the following error: > mclust1Dplot(my.data.model, parameters = my.data.model$parameters, what = "density")

Hi All, I am using em() function to estimate a poisson-gaussian process from a univariate one dimension time series, but not sure how to do. In the help manual, it specify that in "pro" of the argument "parameter", if the model includes a Poisson term for noise, there should be one more mixing proportion than the number of Gaussian components. But in the example, the parameter

Hi all, I have to ask this and I know that the reason is that I am a newbie with R programming. So apologize if it is too obvious but I didn't find an answer after googling and reading "An introduction to R". So i have return data from 30 instruments and I am fitting a mixture of normal distributions for the asymmetric marginal distributions and then simulating from those

Hello all (especially MCLUS users). I'm trying to make use of the MCLUST package by C. Fraley and A. Raftery. My problem is trying to figure out how the (model) identifier (e.g, EII, VII, VVI, etc.) relates to the covariance matrix. The parameterization of the covariance matrix makes use of the method of decomposition in Banfield and Rraftery (1993) and Fraley and Raftery (2002) where

Hi, I am trying to understand how the functions em() and me() from the mclust package work. I cannot make sense of what the algorithm returns. Here is a basic, simple example: ######################################### # two bivariate normals, centered at (-5,0) and (5,0), with Id covariance x1 = cbind(rep(-5, 100), rep(0, 100)) + matrix(rnorm(100*2), 100, 2) x2 = cbind(rep(5, 100),

Dear R users, I use mclust to fit a mixture of normal distributions to many datasets. Usually the Mclust function finds 1 or two normal distributions, rarely, 3. But I hit a strange case today. my.data <- c(57.96920, 51.79415, 51.20538, 55.53637, 51.64291, 56.61476, 51.28855, 55.56169, 51.85113, 54.03330, 51.37370, 49.48561, 52.41580, 53.51176, 60.49293, 55.77012, 51.59270, 56.29660,

My question is with respect to mCLUST and the values of BIC and log likelihood. The relevant part of my R script is: ######################### BEGIN MDS ANALYSIS ######################### #load data data <- read.table("Ecoli33_Barry.dis", header = TRUE, row.names = 1) #perform MDS Scaling mds <- metaMDS(data, k = Dimensions, trymax = 20, autotransform =TRUE, noshare = 0.1,

Hello, I'm trying to use Mclust to fit a Gaussian Mixture Model to a mulitdimensional data set. Because of the specific source of my data, I know that all components have the same variance and that the covariance between dimensions is zero (modelname=VII). Furthermore, I have a reliable estimate of the variance of the components. I want to to use this estimate as a prior in mclust, hoping

Hi All, I have a n x m matrix. The n rows are individuals, the m columns are variables. The matrix is in itself a collection of 1s (if a variable is observed for an individual), and 0s (is there is no observation). Something like: [,1] [,2] [,3] [,4] [,5] [,6] [1,] 1 0 1 1 0 0 [2,] 1 0 1 1 0 0 [3,] 1 0 1 1 0 0 [4,] 0 1 0

Dear Knowledgeable R Community Members, Please excuse my ignorance, I apologize in advance if this is an easy question, but I am a bit stumped and could use a little guidance. I have a finite mixture modeling problem -- for example, a 2-component gaussian mixture -- where the components have a large overlap, and I am trying to use the "mclust" package to solve this problem. I need

Hi, all; I got an unusual error while using mclust library 2.1-7 on R 2.0.1. When I tried to run Mclust(rnorm(100)), I got segmentation fault error. Does anyone know what causes this problem? I came across the same problem even when I tried to run the example shown in Mclust help using iris data. Thanks in advance, Tae-Hoon Chung -------------------------------------------------- Tae-Hoon Chung

I've been trying to install some R packages such as mclust and flexmix on linux but have had the following error messages. > I've been trying to install mclust on my notebook which has linpus linux lite os and I have installed R as well as some packages all right. However, when I tried to install mclust, it gave me the following messages. Any suggestions? Tim ++++++++++++++++

Hi All, I am trying to cluster a one-dimensional data (see the file attached) using Mclust() but got an error message like: >Mclust(x) Error in rep(1, n) : Object "n" not found When I do a simulation sometimes it works sometimes doesn't. >Mclust(c(rnorm(50),rnorm(56,-0.5))) Error in rep(1, n) : Object "n" not found >Mclust(c(rnorm(56),rnorm(56,-0.5))) best

Hi, All; I tried to use library mclust in 64-bit compiled R 2.0.1 but failed. Installation went smoothly without any warning or error. However, when I tried to use them with the following simple code, it crashed. Library(mclust) Dat <- c(rnorm(20, mean=0, sd=0.2), rnorm(30, mean=1, sd=0.2)) Ind <- Mclust(dat, 1, 5)$classification cbind(Dat, Ind) The error message was: