From: Karen Green <karengreenphd_at_yahoo.com>

Date: Sat 22 Oct 2005 - 18:33:59 EST

**PROBLEM DESCRIPTION
**

**QUESTIONS
**

**SAMPLE CODE
**

data(faithful)

library(mclust)

Date: Sat 22 Oct 2005 - 18:33:59 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 regarding parameters in 2 functions in the "mclust" package.

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.)

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

data(faithful)

library(mclust)

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

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

# 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 ______________________________________________ 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 Sat Oct 22 18:40:55 2005

*
This archive was generated by hypermail 2.1.8
: Sun 23 Oct 2005 - 19:17:49 EST
*