Bom dia Pessoal, Estou com problemas no ajuste de um modelo GLM de Poisson, trata-se de uma relação entre número de insetos (n.insetos) com a temperatura (temp). O ajuste usando o modelo GLM de Poisson (mod.1) não parece ser o mais adequado e demonstra sobredispersão, além que a interpretação da curva do ponto de vista biológico não é realistica, pois eu teria que ter uma curva que aumenta com a temperatura, chega a um ponto ótimo e diminui, como uma equação de 2º grau, porém o modelo de Poisson que utilizei não faz isto, por isso pergunto, existe alguma sugestão de como eu poderia ajustar dados dessa natureza com mais semelhança biológica com o fenômeno observado? Segue CRM como exemplo: ##Dados artificias # Número de indivíduos n.insetos<-c(rpois(10,lambda=5),rpois(10,lambda=8),rpois(10,lambda=20),rpois(10,lambda=30),rpois(10,lambda=7),rpois(10,lambda=5)) #Temperatura temp<- c(rnorm(10,10,0.1),rnorm(10,15,0.1),rnorm(10,22,0.1),rnorm(10,24,0.1),rnorm(10,28,0.1),rnorm(10,32,0.1)) # #Visualiza o comportamentos das variáveis plot(temp,n.insetos) #Ajuste por GLM de Poisson mod.1<-glm(n.insetos~temp, family="poisson") summary(mod.1) anova(mod.1, test="Chi")## Anova #Representação do modelo a<-mod.1$coeff[1]###Coeficiente a b<-mod.1$coeff[2] ###Coeficiente b x<-seq(5,38,1) curve(exp(a+b*x),lty=1, add=T) # Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com