Displaying 2 results from an estimated 2 matches for "mymixturemodel".
2005 Oct 21
1
finite mixture model (2-component gaussian): plotting component gaussian components?
...ptic.
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(da...
2005 Oct 22
0
package mclust: cdens, EMclust?
...the documentation is
somewhat cryptic.
--------------------
SAMPLE CODE
--------------------
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"...