Hi,
I am trying to run th mclust package on a variable "Tuberculin indurations"
recorded as mm. The file has only one variable. When I run the package I get
NULL value for mu and sigma. Can anybody say why?
This is the program:
library("mclust")
mc<-Mclust(x.trab,G=1:9,warn=TRUE)
mc
mc$mu
sqrt(mc$sigmasq)
and the output I get is
> library("mclust")
> mc<-Mclust(x.trab,G=1:9,warn=TRUE)
There were 11 warnings (use warnings() to see them)
> mc
best model: unequal variance with 3 components
> mc$mu
NULL
> sqrt(mc$sigmasq)
Error in sqrt(mc$sigmasq) : Non-numeric argument to mathematical function
> warnings()
Warning messages:
1: In meV(data = data, z = z, prior = prior, control = control, ... :
sigma-squared falls below threshold
2: In mclustBIC(data = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, ... :
there are missing groups
3: In meV(data = data, z = z, prior = prior, control = control, ... :
sigma-squared falls below threshold
4: In mclustBIC(data = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, ... :
there are missing groups
5: In meV(data = data, z = z, prior = prior, control = control, ... :
sigma-squared falls below threshold
6: In mclustBIC(data = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, ... :
there are missing groups
7: In meV(data = data, z = z, prior = prior, control = control, ... :
sigma-squared falls below threshold
8: In mclustBIC(data = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, ... :
there are missing groups
9: In meV(data = data, z = z, prior = prior, control = control, ... :
sigma-squared falls below threshold
10: In mclustBIC(data = c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, ... :
there are missing groups
11: In meV(data = data, z = z, prior = prior, control = control, ... :
sigma-squared falls below threshold
Please help.
Basilea Watson
Biostatistician
Tuberculosis Research Centre, India
[[alternative HTML version deleted]]
______________________________________________
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.
