search for: emclust

Displaying 10 results from an estimated 10 matches for "emclust".

Did you mean: mclust
2005 Oct 22
0
package mclust: cdens, EMclust?
...ap, and I am trying to use the "mclust" package to solve this problem. I need to decompose that mixture into its 2 components which will need to be plotted. ----------- QUESTIONS ----------- What I don't know how to do is: (1) restrict the number of components to 2 in the "EMclust" function (2) obtain and plot a component gaussian density Regarding (1), I think this should be the 'G' value but the documentation is somewhat cryptic as regards the format. Regarding (2), I think I need to use the "cdens" function, but again the documentation is somewh...
2002 Feb 14
1
Subsets in mclust
...ing phase will use a sample of size `k' of the data in the initial hierarchical clustering phase. The default is to use the entire data set. m2 is a matrix (not very large, but the same happens always when specifying k): > dim(m2) [1] 242 2 > em321 <- emclust(m2,k=50) Error in mhtree.VVV(data[smpl, ]) : data must be a matrix with at least 2 columns > em321 <- emclust(m2) (No problems) > em321 <- emclust(m2,modelid="VEV",k=50) Error in mhtree.VVV(data[smpl, ]) : data must be a matrix with at least 2 columns > library(help=mclus...
2005 Oct 21
1
finite mixture model (2-component gaussian): plotting component gaussian components?
...e the components have a large overlap, and I am trying to use the "mclust" package to solve this problem. I need to decompose that mixture into its 2 components which will need to be plotted. What I don't know how to do is: (1) restrict the number of components to 2 in the "EMclust" function (2) obtain and plot a component gaussian density Regarding (1), I think this should be the 'G' value but the documentation is somewhat cryptic as regards the format. Regarding (2), I think I need to use the "cdens" function, but again the documentation is somewha...
2004 Oct 04
3
Help with normal distributions
Hi I have two questions, the first perhaps dumber than the second. Firstly, I have a data set, and when I plot a histogram it looks like a normal distribution. So I want to overlay a bell-shaped normal distribution on top of it, to demonstrate how similar it is to the normal distribution. I have read the help on dnorm(), rnorm(), pnorm() etc but still can't figure out how to plot a normal
2010 Jan 06
1
positive log likelihood and BIC values from mCLUST analysis
...0.1, wascores = TRUE, expand = TRUE, trace = FALSE, plot = FALSE, old.wa = FALSE) ######################### BEGIN EM ANALYSIS ######################### #Use the points determined by MDS to perform EM clustering. #Allow only the unconstrained models. Sometimes, constrained models mess things up! EMclusters <- mclustBIC(mds$points, G=Clusterrange, modelNames= c("VII", "VVI", "VVV"), prior=NULL, control=emControl(), initialization=list(hcPairs=NULL, subset=NULL, noise=NULL), Vinv=NULL, warn=FALSE, x=NULL) The input data are in the form of an N...
2004 May 06
0
Problem with mclust surfacePlot function
...oods maples. I am using R 1.8.1 with the code below. The same code works fine in S-Plus 6.2 Am I missing something or is this a bug? Thanks Mike White library(mclust) data(lansing) # R only maples<-lansing[as.character(lansing[,"species"]) == "maple", -3] maplesBIC <- EMclust(maples) maplesModel <- summary(maplesBIC, maples) plotMaples2 <- function(type, what, transformation) { out <- do.call("surfacePlot", c(maplesModel, list(data=maples, type=type, what=what,transformation=transformation))) invisible() } par(pty="s", mfrow=c(1,2)) plotMa...
2017 Aug 27
1
Fwd: Find maxima of a function
I have not followed the history of this thread, but I am quite flummoxed as to why the OP is rewriting code to estimate parameters from an univariate Gaussian mixture model when alternatives such as EMCluster (which generally appears to handle initialization better than MClust) exist. Or perhaps there is more to it in which case I apologize. But I thought that I would make the OP aware of the package (of which, in full disclosure, I am co-author). Best wishes, Ranjan On Sun, 27 Aug 2017 16:01:59 +0...
2003 Jul 17
3
univariate normal mixtures
Hello, I have a concrete statistical question: I have a sample of an univariate mixture of an unknown number (k) of normal distributions, each time with an unknown mean `m_i' and a standard deviation `k * m_i', where k is known factor constant for all the normal distributions. (The `i' is a subscript.) Is there a function in R that can estimate the number of normal distributions k
2017 Aug 27
0
Fwd: Find maxima of a function
Dear Niharika, As I said before, the problem is basically an optimization issue. You should isolate the problematic part from the rest of your study. Sometimes, more information does not help to solution. All the answers from us (Ulrik, David, me) are more or less are correct to find a maximum point. Newton?s method is also correct. But after answers, you only say, it didn?t find the right
2017 Aug 27
2
Fwd: Find maxima of a function
---------- Forwarded message ---------- From: niharika singhal <niharikasinghal1990 at gmail.com> Date: Sun, Aug 27, 2017 at 11:57 AM Subject: Re: Find maxima of a function To: "David Winsemius [via R]" <ml+s789695n4745009h56 at n4.nabble.com>, "Ismail SEZEN [via R]" <ml+s789695n4744993h60 at n4.nabble.com>, Ulrik Stervbo <ulrik.stervbo at gmail.com>