Hola a todos!

Tengo una duda respecto a la estructura del factor aleatorio que debo usar en un modelo lineal mixto.
Hice el siguiente experimento en nidos de aves:
Medí la intensidad de la respuesta agresiva de los dueños de los nidos frente a tres tipos intrusos (una especie control,  depredador 1 y depredador 2).
Cada nido/pareja fue expuesto a los tres intrusos el mismo día, uno a continuación del otro (con un intervalo de descanso entre presentaciones).  
Estos tratamientos se hicieron en dos etapas del ciclo de nidificación (Etapa 1 y Etapa 2). Los nidos/parejas no fueron utilizados en ambas etapas, son independientes.
Me interesa comparar la respuesta frente a los 3 modelos en cada etapa del ciclo y entre ambas etapas del ciclo.

El modelo que corrí es el siguiente: 




<-lme(ACPr~Intruso*Etapa+Date+Adultos+PV+MB,random=~1|ID,data=acp)
El efecto aleatorio que utilice  para contemplar las pseudoréplicas al medir al mismo individuo varias veces es la identidad del nido (ID) y tengo dudas respecto a si esto está bien así o me falta algo para obtener todas las comparaciones que necesito. 
Date, Adultos, PV, MB son co-variables.

Muchas graciasCarolina





	[[alternative HTML version deleted]]

