Displaying 4 results from an estimated 4 matches for "mclust1dplot".
2008 Oct 20
1
Mclust problem with mclust1Dplot: Error in to - from : non-numeric argument to binary operator
...ar list members,
I am using Mclust in order to deconvolute a distribution that I
believe is a sum of two gaussians.
First I can make a model:
> my.data.model = Mclust(my.data, modelNames=c("E"), warn=T, G=1:3)
But then, when I try to plot the result, I get the following error:
> mclust1Dplot(my.data.model, parameters = my.data.model$parameters, what = "density")
Error in to - from : non-numeric argument to binary operator
Also, I'd like to allow for each gaussian to have a different variance
(modelNmaes=c("V")) , but then I get another error message:
> my.d...
2005 Oct 21
1
finite mixture model (2-component gaussian): plotting component gaussian components?
...############
# (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...
2007 Mar 23
0
plotting dnorm() issued from mclust models
Dear all
I have a problem in fitting lines() of the normal distributions
identified with Mclust on a histogram or a mclust1Dplot. Here is some
sample code to explain :
set.seed(22)
foo <- c(rnorm(400, 10, 2), rnorm(500, 17, 4))
mcl <- Mclust(foo, G=2)
mcl.sd <- sqrt(mcl$parameters$variance$sigmasq)
mcl.size <- c(length(mcl$classification[mcl$classification==2]),
length(mcl$classification[mcl$classi...
2005 Oct 22
0
package mclust: cdens, EMclust?
...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...