Erro ao rodar modelo em GLMM

[Bruna]

Boa tarde a todos,

Estou trabalhando em GLMM com distribuição de poisson. Ao rodar o modelo, dá a seguinte mensagem de erro:

#modelo para vi2

> vi2modelo1<-glmer(cbind(vi)~sp+sexo+(1|dia)+(1|id),poisson(link="log"),data=abf)   ### este deu certo

> vi2modelo2<-glmer(cbind(vi)~sp*sexo+(1|dia)+(1|id),poisson(link="log"),data=abf) 

Error in pwrssUpdate(pp, resp, tol = tolPwrss, GQmat = GQmat, compDev = compDev,  : 

  Downdated VtV is not positive definite

> vi2modelo3<-glmer(cbind(vi)~sp+sexo*(1|dia)+(1|id),poisson(link="log"),data=abf) ### deu certo

Uma das características de meus dados é que tenho muitos zeros nas variáveis respostas. Isso pode ser a causa do problema? Mas porque dá certo para alguns argumentos e outros não?

Estou usando o pacote lme4.

Não sei se fui clara na explicação, mas se precisarem de mais informações a respeito é só me dizer. Espero conseguir ajuda para resolver este problema, pois preciso dar continuidade às análises o quanto antes.

Obrigada.

[Augusto Ribas]

Ola, tem uma discussão sobre isso no github onde o pacote lme4 é desenvolvido.

https://github.com/lme4/lme4/issues/173

Talvez você possa sugerir os valores iniciais para os parametros, para facilitar para o otimizador, pelo menos o que entendi é que o algoritimo que otimizador encalha.

Tente produzir um exemplo reproduzivel, sem ter ideia de como são as variaveis é complicado sugerir codigo.

--
Sempre envie um código mínimo reproduzível.
Envie ou responda perguntas para o e-mail (r-eco...@googlegroups.com)
---
Você recebeu essa mensagem porque está inscrito no grupo "R-Ecology" dos Grupos do Google.
Para cancelar inscrição nesse grupo e parar de receber e-mails dele, envie um e-mail para r-ecologia+...@googlegroups.com.
Para postar nesse grupo, envie um e-mail para r-eco...@googlegroups.com.
Acesse esse grupo em http://groups.google.com/group/r-ecologia.
Para ver essa discussão na Web, acesse https://groups.google.com/d/msgid/r-ecologia/ac85e8a2-79f4-4189-9a21-d07924c4088a%40googlegroups.com.
Para mais opções, acesse https://groups.google.com/d/optout.

--

Grato
Augusto C. A. Ribas

 

Site Pessoal: http://recologia.com.br/

Github: https://github.com/Squiercg

Lattes: http://lattes.cnpq.br/7355685961127056

[Nicholas Marino]

Oi Bruna!

Os seus dados estao balanceados para todas as combinações de tratamentos?

Enviada do meu iPhone

[Bruna Ramos]

Oi Nicholas,

Desculpa, eu to começando a estudar GLMM e R recentemente, então você poderia me explicar o que você quer dizer com balanceados?

Obrigada

Abs, Bruna

Você recebeu essa mensagem porque está inscrito em um tópico no grupo "R-Ecology" dos Grupos do Google.
Para cancelar inscrição nesse tópico, acesse https://groups.google.com/d/topic/r-ecologia/gF0q4rmJgz4/unsubscribe.
Para cancelar inscrição nesse grupo e todos os seus tópicos, envie um e-mail para r-ecologia+...@googlegroups.com.

Para postar nesse grupo, envie um e-mail para r-eco...@googlegroups.com.
Acesse esse grupo em http://groups.google.com/group/r-ecologia.

Para ver essa discussão na Web, acesse https://groups.google.com/d/msgid/r-ecologia/38ABAF47-EBB1-4E9C-B984-2AD6E7656491%40gmail.com.

Para mais opções, acesse https://groups.google.com/d/optout.

--

Bruna Menezes Ramos

Doutoranda em Ecologia - UFRJ

Lab. de Interação Inseto-Planta

Av. Carlos Chagas Filho, 373- Cidade Universitária (Ilha do Fundão), RJ, Brasil 

Telefone: 3938. 6311

[Bruna Ramos]

Oi Nicholas,

sim, meus dados são balanceados (sem dados faltantes).

Obrigada.

Abs, Bruna