estoy intentando realizar un estudio sobre invertebrados en diferentes tipos de olivares. Mi duda es la siguiente: quiero saber si alguna especie de invertebrado se asocia significativamente a alguno de los olivares (son tres tipos: superintensivo, convencional y ecológico). Para ello, he intentado elaborar un GLM anidado, con Abundancia de cada especie (nº de individuos) como variable respuesta, y el Muestreo (se hicieron 3: abril, mayo y junio), Manejo (tipo de olivar), y Especies (nombre de cada especie) como factores. Además, he intentado que Especies vaya anidada dentro de Muestreo, anidado a su vez dentro de Manejo. El resultado que he obtenido lo he traducido como que, donde los valores de p son menores de 0,05, significa que esa especie está asociada a ese olivar.
Los datos analizados siguen el siguiente formato:
Especies | Manejo | Muestreo | Abundancia |
---|---|---|---|
Especie 1 | Superin | Muestreo1 | 70 |
Especie 1 | Superin | Muestreo2 | 50 |
Especie 1 | Superin | Muestreo3 | 5 |
Especie 1 | Convenc | Muestreo1 | 40 |
Especie 1 | Convenc | Muestreo2 | 34 |
Especie 1 | Convenc | Muestreo3 | 12 |
Especie 1 | Ecologi | Muestreo1 | 20 |
Especie 1 | Ecologi | Muestreo2 | 6 |
Especie 1 | Ecologi | Muestreo3 | 4 |
Especie 2... |
El script de R que he utilizado es el siguiente:
summary(glm(Abundancia ~ Muestreo+Manejo+ï..Especies%in%Muestreo%in%Manejo,data=bioin,family=poisson(link="log")))
Me gustaría saber si el modelo que he hecho sería correcto para responder a la pregunta que me hago (cuáles son las especies que se asocian a los olivares).
Muchas gracias.