[glmmADMB] Correta modelagem para dados de contagem

Estou tentando analisar um conjunto de dados de contagem e fui cair na pagina do pacote glmmADMB. No entanto, achei pouca informação que pudesse me ajudar. Meus dados estão assim: Esp Mat Time Rep Germ Inf Sucupira M02 1 1 0 0 Sucupira M02 2 1 1 0 Sucupira M02 3 1 1 8 Sucupira M02 4 1 1 8 Sucupira M02 5 1 1 9 Sucupira M02 6 1 1 10 Sucupira M02 7 1 1 15 Sucupira M02 1 2 0 0 (...) Sucupira M02 7 3 1 10 Sucupira M04 1 1 0 0 Sucupira M04 2 1 0 1 (...) Vinhatico M15CV15 7 3 5 8 Esp = Espécie Mat = Matriz (três matrizes de locais diferentes) Time = Dias em que foi analisado Rep = Repetição (gerbox com 20 sementes cada) Germ = Quantidade de sementes que germinaram Inf = Quantidade de sementes com micro-organismos Trata-se de medidas repetidas no tempo, ou seja, no dia 1 contaram-se quantas sementes germinaram e quantas apresentavam micro-organismos. No dia 2 contaram-se quantas sementes germinaram e quantas apresentavam micro-organismos. Assim até o dia 7. Tentei o seguinte código: library(glmmADMB) fit = glmm.admb(Germ~Esp*Mat+Rep,random=~Rep,group="Time", zeroInflation=TRUE,data=ka.dat,family="poisson") Mas, além do erro: Maximum number of iterations exceeded in dvector eigen(_CONST dmatrix& m) Error in glmm.admb(Germ ~ Esp * Mat + Rep, random = ~Rep, : The function maximizer failed In addition: Warning message: running command './nbmm -maxfn 500 ' had status 1 não tenho certeza se a modelagem proposta por mim esta correta. Por exemplo, o 'group' deveria ser mesmo o tempo? E o efeito aleatorio seria mesmo so a repeticao? Eu tambem gostaria de ver se ha alguma relacao entre germinacao e infestacao. Se voce tiver experiencia com esse pacote, poderia dar algumas dicas? Obrigado -- Marcelo Brazil Linux user number 487797

Marcelo, Não conheço o pacote mas o seu dado tem mais cara de Binomial que Poisson, uma vez que se trata do número de sementes germinadas num total de 20. Pra ser mais criterioso, como você observa resultados indexados no tempo, é um binomial em que o tamanho (parâmetro n) diminui. No dia 0 é um binomial com n=20 e p1, se germinam 2 sementes, no dia 2 é um binomial com n=18 e p2. De qualquer forma, não tem função que modele dessa forma. Para modelar da forma tradicional, você pode usar a lme4::glmer(), declarando a parcela (gerbox com 20 sementes) como fator de efeito aleatório. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Marcelo Luiz de Laia
-
Walmes Zeviani