Usando sim() com lmer()

Question

Executei duas regressões logísticas multinível usando os mesmos preditores, mas em duas respostas diferentes:

fruitMLM <- lmer(InsuffFruit ~ Income + HDI + Income:HDI + (1 + Income | Country),family=binomial(link="logit"))  
fuelMLM <- lmer(Pollution ~ Income + HDI + Income:HDI + (1 + Income | Country),family=binomial(link="logit"))

Income é discreto com valores c(-2,-1,0,1,2), HDI é contínuo entre 0 e 1, Country é categórico e as respostas são ambas 1/0.

Para traçar faixas de confiança, executo uma simulação usando a função sim() do pacote arm:

sim(fruitMLM,100)  
sim(fuelMLM,100)

O primeiro calcula bem.O segundo retorna o seguinte erro:

Error in mvnorm(n.sims, bhat[j,], V.beta) :  
  'Sigma' is not positive definite

Na verdade, estou fazendo isso com 8 respostas diferentes.Seis deles funcionaram bem e dois deles retornaram este erro.

Alguém sabe como corrigir isso?

Answer 1

Você pode ver se a mensagem de erro está correta, ou seja, verificar "manualmente" se a matriz de covariância do modelo com falha é realmente positiva semidefinida?Não tenho certeza, mas imagino que existam algumas verificações que você possa realizar - talvez haja algo aqui: http://en.wikipedia.org/wiki/Positive-definite_matrix