Hola,
Soy nuevo manejando R y no tengo mucha experiencia. Estoy intentando modelar
una funci�n que me relacione el n� de cebas (n� de presas que los padres traen
a los pollos) con el tama�o de parche de un bosque (factor categ�rico; 2
niveles= grande y peque�o). Al ser un conteo (n� de cebas) he pensado utilizar
familia= poisson con link= logar�tmico. He construido un GLMM con: N� de cebas
(variable respuesta), Duraci�n del video (esfuerzo de muestreo: 60, 50 y 45
minutos) como una offset en el modelo (transform�ndola previamente por el
logaritmo). Tama�o de parche, factor explicativo, N� de pollos y fecha de
puesta covariables. Adem�s he incluido como random el a�o de estudio (3 a�os),
el ID del nido anidado dentro del ID de cada parche (tengo varios parches,
grandes y peque�os) y la identidad del padre y de la madre tambi�n como r�ndom,
puesto que se pueden repetir entre a�os. La estructura de mi modelo es la
siguiente:
#paquete(lme4), mi versi�n de R es la 3.1.1.
m1<- glmer(ncebas~offset(LogDur_video)+ Tama�oParche*Ldate*Csize+(1|YEAR) +
(1|Bosque/nideo)+ (1|Hembra) + (1|Macho), data=Cebas, family=poisson(link=log))
#con todas las variables e interacciones posibles para que sirva como
referencia a la hora de la selecci�n del modelo final.
El modelo me funciona, sin embargo me aparecen unos warnings que no s� a que se
deben ni la importancia real en el c�culo de mi modelo. Mi pregunta es si
alguien en el foro sabe porqu� aparecen y como solucionarlos. Los warnings son:
Warning messages:
1: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model failed to converge with max|grad| = 0.00431561 (tol = 0.001, component
7)
2: In checkConv(attr(opt, "derivs"), opt$par, ctrl = control$checkConv, :
Model is nearly unidentifiable: very large eigenvalue
- Rescale variables?;Model is nearly unidentifiable: large eigenvalue ratio
- Rescale variables?
Muchas gracias, atentamente,
Javier
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
[email protected]
https://stat.ethz.ch/mailman/listinfo/r-help-es
