[R] finite mixture model (2-component gaussian): plotting component gaussian components?

From: <Karen.Green_at_sanofi-aventis.com>
Date: Fri 21 Oct 2005 - 11:42:18 EST


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 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 somewhat cryptic.

Here is a little test script to illustrate. (Note: my real dataset will not have peaks this well separated, but I needed to find a small example.)

##################

# (1) get some data
##################

data(faithful)

##################

# (2) get model
##################

library(mclust)
MyMixtureModel<-summary(EMclust(faithful$eruptions),faithful$eruptions)

##################

# (3) plot mixture model
##################

attach(MyMixtureModel)
mclust1Dplot(data=faithful$eruptions,z=z,mu=mu,sigmasq=sigmasq,pro=pro,ask=FALSE,type=c("density")) do.call("mclust1Dplot",c(list(data=faithful$eruptions,ask=FALSE,type=c("density")),MyMixtureModel))

##################

# (4) plot components
##################

???

Any information you might be able to shed on this would be very much appreciated.

With appreciation for your help,

Karen

---
Karen M. Green, Ph.D.
Research Investigator
Drug Design Group
Sanofi Aventis Pharmaceuticals
Tucson, AZ
USA


	[[alternative HTML version deleted]]

______________________________________________
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide! http://www.R-project.org/posting-guide.html
Received on Fri Oct 21 12:00:16 2005

This archive was generated by hypermail 2.1.8 : Fri 03 Mar 2006 - 03:40:46 EST