From: <Karen.Green_at_sanofi-aventis.com>

Date: Fri 21 Oct 2005 - 11:42:18 EST

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.htmlReceived 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
*