Hola, gracias por el mail precedente. Me sirvió para resolver mi duda. Aquí tengo un problema más. Escribo esta función: fit01<-lmer(cr.ur~type.nest+hour+age+sit+type.nest:age+type.nest:sit+(1|years)+(1|nest), data=ceba1, family=poisson,na.action=na.omit) pero en el resultado no me puede calcular los valor de los efecto random. Eso es lo que sale de la parte relativa a los random effects: Random effects: Groups Name Variance Std.Dev. nest (Intercept) 0 0 years (Intercept) 0 0 Number of obs: 51, groups: nest, 21; years, 4 En mi base de dato tengo 4 años (pero en cada año el número de nido en cada año no es lo mismo. Ej 2008 7 nidos, 2009 40, etc etc). Es esto el problema? También he intentado con el anidado, pero sale exactamente igual. Que tengo que hacer? PD: aqui pongo la str de i datos por si hace falta. data.frame'': 90 obs. of 10 variables: $ years : int 2008 2008 2008 2008 2008 2009 2009 2009 2009 2009 ... $ hour : int 1 1 1 1 1 1 2 2 1 1 ... $ nest : Factor w/ 36 levels "21","27","3",..: 25 17 17 17 18 35 35 26 26 7 ... $ type.nest: Factor w/ 2 levels "m","s": 1 1 1 1 1 1 1 1 1 1 ... $ age : int 34 22 25 NA NA NA NA NA NA 28 ... $ sit : int 3 1 1 2 1 1 1 1 1 1 ... $ cr.ur : int 0 0 1 5 1 0 0 0 0 0 ... $ ur.cr : int 15 2 3 2 2 0 0 4 2 1 ... $ pur.urr : int 0 0 0 0 0 0 0 0 0 0 ... $ urr.pur : int 23 1 3 0 0 4 0 4 3 0 ... Gracias, Gianlu [[alternative HTML version deleted]]