[R-br] GLM de Poisson para dados com ponto máximo

[ASANTOS]

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 em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================


---
Este email foi escaneado pelo Avast antivírus.
http://www.avast.com