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>