similar to: Help with mclust package

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, 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

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 regarding parameters in 2 functions in the "mclust" package. -------------------- PROBLEM DESCRIPTION -------------------- I have a finite mixture modeling problem -- for example, a 2-component gaussian

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

Whenever I submit the following command >result <- em("E", m1, mu=c(25, 50), sigmasq=10, pro=c(0.49,0.51)) in MCLUST , I experience a problem with my Windows XP. Has anyone had this type of problem? I receive the general "send report" dialogue box and my program is no longer able to run. Should I change any of my system settings? Let me know if you think of anything

Hi! Why do I still get an error message on the try function?: > options(show.error.messages = FALSE) > > trab <-matrix(nrow=6,ncol=6) > trab[,] <-0.0001 > rownames(trab) <- c("A","B","C","D","E","I") > colnames(trab) <-

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,

I'm running Mclust 3.0 in R-2.7.0 and have run into a situation where the optimal number of clusters, 12, seems to be greater than the number of colors available to Mclust. The code below, which demonstrates errors, does not if the max no of clusters is set at 10. sessionInfo, after a restart of R and loading of packages, follows the output. Is this a known problem with Mclust? Can I

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,

Dear all I have a problem in fitting lines() of the normal distributions identified with Mclust on a histogram or a mclust1Dplot. Here is some sample code to explain : set.seed(22) foo <- c(rnorm(400, 10, 2), rnorm(500, 17, 4)) mcl <- Mclust(foo, G=2) mcl.sd <- sqrt(mcl$parameters$variance$sigmasq) mcl.size <- c(length(mcl$classification[mcl$classification==2]),

Dear all, I am attempting to model some one-dimensional data using Gaussian mixture model with mclust. Generally, the data that I have have 3 overlapping populations (with one of them being the majority, and the other two combining to less than 15%) and for some reason, mclust consistently ignores the smaller peaks, giving me strange values for the means (completely anti-intuitive in terms of

Hello, I'm trying out flexmix and mclust for the first time on some univariate data which is typically best described as lognormal, but can sometimes be gamma distributed as well. I first tried using EM on mclust assuming the data was lognormally distributed and could only get it to work in "E" mode, i.e. the equal variance mode. I could never get it to work on "V" mode [

> No HTML this time. Sorry Dear all, I am attempting to model some one-dimensional data using Gaussian mixture model with mclust.? Generally, the data that I have have 3 overlapping populations (with one of them being the majority, and the other two combining to less than 15%) and for some reason, mclust consistently ignores the smaller peaks, giving me strange values for the means

Hallo together, I am trying to estimate parameters by means of QMLE using the nlminb optimizer for a tree-structured GARCH model. I face two problems. First, the optimizer returns error[8] false convergence if I estimate the functions below. I have estimated the model at first with nlm without any problems, but then I needed to add some constraints so i choose nlminb.

Hi, My questions concern the function 'mclustBIC' which compute BIC for a range of clusters of several models on the given data and the other function 'mclustModel' which choose the best model and the best number of cluster accordind to the results of the previous cited function. 1) When trying the following example (see ?mclustModel), I get negative BIC computed by

Dear R help list, I have a plot with two different vertical scales that I want to display on either side of the plot. It's quite similar to the Fahrenheit-Centigrade example in the examples section of the documentation for axis.default. The right-side axis is clipped off, though, and I haven't been able to figure out anything with viewport() and clipping or trellis.par.set to fix

I'm applying a function (Cov.f) defined below to each element of a distance matrix. When I run the code below, I get a warning message (below) and elements of returned matrix [2,3] and [3,2] are not zero as I would expect. Clearly, there is an error... What am I doing wrong? Thanks. --Dale Warning message: In if (h <= phi) { : the condition has length > 1 and only the first element

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

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

Dear people, I've a problem in doing a power calculation. In Fryer and Nicholson (1993), ICES J. mar. Sci. 50: 161-168 page 164 an example is given with the following characteristics T=5, points in time R=5, replicates Var.within=0.1 q=10, a 10% increase per year The degrees of freedom for the test are calculated as Vl=T*R-2=23 and the non-centrality parameter Dl=4.54. Using this they get a