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

