I probably don't understand problem. I'd assume that variance$sigmasq are
the three estimated component variances (probably estimated by maximum a
posteriori, but consult the mclust documentation).
What's wrong with that?
(The values you submit as scale in "prior" are not fixed variances, but
parameters of the prior distribtion - your problem may be that you
believe that they are meant to be variances fixed by you!?)
Christian
On Tue, 6 Sep 2011, Nico902 wrote:
Hi,
Thanks a lot for your answer. I effectively was able to get rid of this
message by doing:
resClust <-
Mclust(data,G=3,modelName="V",prior=priorControl(scale=c(1.44,0.81,0.49)));
However, I would like to be able to retrieve the variances I defined in the
result. I found:
resClust$parameters
$Vinv
NULL
$pro
[1] 0.5502496 0.1986852 0.2510652
$mean
1 2 3
-2.8390006980 -0.0003267873 3.1072574619
$variance
$variance$modelName
[1] "V"
$variance$d
[1] 1
$variance$G
[1] 3
$variance$sigmasq
[1] 0.840267666 0.009466821 1.510263146
$variance$scale
[1] 0.840267666 0.009466821 1.510263146
I do not manage to get where the sigmasq is coming from. I tried to sqrt or
square the sigmasq but it does not correspond to what I defined. I found
nothing in the manual. If I am missing something obvious or if somebody has
the solution it will help me a lot. I want to retrieve those values
automatically to plot the different curves of the fitting and to be sure
this is doing what I want.
Thank you very much again.
--
View this message in context:
http://r.789695.n4.nabble.com/mclust-modelName-E-vs-modelName-V-tp3789167p3793697.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
*** --- ***
Christian Hennig
University College London, Department of Statistical Science
Gower St., London WC1E 6BT, phone +44 207 679 1698
chr...@stats.ucl.ac.uk, www.homepages.ucl.ac.uk/~ucakche
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
