[R-es] random effects

gianlu.85 en virgilio.it
Thu, 20 Dec 2012 12:25:17 -0000

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]]