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