TRABAJO FINAL BIOESTADÍSTICA

Vista previa del material en texto

UNIVERSIDAD NACIONAL AGRARIA LA MOLINA
FACULTAD DE ECONOMÍA Y PLANIFICACIÓN
DEPARTAMENTO ACADÉMICO DE ESTADÍSTICA E INFORMÁTICA
CURSO: BIOESTADÍSTICA
Trabajo Final de Bioestadística
Artículo: Selección de hojas de Carludovica palmata como refugio por el
murciélago toldero Dermanura watsoni en la Estación Biológica Piro, Península
de Osa, Costa Rica.
Profesor: Febres Huamán, Grimaldo José
Integrantes:
- Custodio Jaimes, Rosa María 20181002
- Rojas Espinoza, José Miguel 20181022
- Huaranga Manrique, Alexander Martín Fabián 20181012
Ciclo: 2021-II
La Molina, 2021
1 INTRODUCCIÓN
El refugio es un recurso importante para los murciélagos, debido a que pasan la mayor
parte de su vida en ellos. Muchas especies de murciélagos usan cuevas, grietas, troncos de
árboles, puentes, construcciones humanas, entre otros (Altringham 1996, Kunz y Lumsden
2003). Existen actualmente 22 especies de murciélagos, sobre todo de la subfamilia
Stenodermatinae (familia Phyllostamidae), que se refugian en hojas modificadas conocidas
como «tiendas» (RodríguezHerrera et al. 2007). Estas tiendas le proporcionan a los
murciélagos protección contra factores ambientales como radiación solar, lluvia y viento,
que pueden afectar su termorregulación (Neuweiler 2000).
Dermanura watsoni se refugia bajo hojas
modificadas, siendo la especie que utiliza
la mayor diversidad de plantas como
refugio. Para varias de estas especies se
han descrito las características
morfológicas de las plantas que son
preferidas por D. watsoni. Sin embargo,
para C. palmata, la cual es una de las
plantas más usadas como refugio por esta
especie debido a su longevidad como
refugio (16 meses), no está muy claro qué
características de las hojas son relevantes
a la hora de elegirlas como refugio.
El objetivo de este artículo fue determinar
qué características físicas y estructurales
de las plantas y hojas y su entorno, influyen en la selección de hojas de C. palmata por el
murciélago constructor de tiendas D. watsoni.
Se hipotetiza que estas características son similares a las hojas de otras especies de
plantas que D. watsoni también prefiere.
2 METODOLOGÍA
La investigación se realizó en la estación Piro, en el Área de Conservación Osa, ubicada en
el suroeste de Costa Rica, durante el mes de julio de 2012. Se realizó un muestreo a lo
largo de los senderos de la estación Piro (5 km aproximadamente), buscando cualquier
planta de C. palmata con hojas con presencia de cortes que hacen los murciélagos con los
dientes y que tuvieran la forma tipo sombrilla que le dan a la hoja modificada (tienda). Una
vez encontrada la tienda se midieron las siguientes características: altura de la hoja con
respecto al suelo, radio de la hoja (circunferencia laminar), ángulo del peciolo, diámetro del
peciolo, distancia al árbol más cercano, y porcentaje de cobertura de dosel. Se usó el
método del vecino más cercano para tomar las medidas a las hojas sin modificar y
seleccionadas al azar como hojas control.
Todos los atributos estructurales de las hojas y las medidas del medio circundante se
analizaron tanto para hojas modificadas como no modificadas, a través de una regresión
logística, utilizando modelos lineales generalizados bajo la distribución binomial con el
paquete estadístico R-Project versión 2.8.1. Para escoger los modelos que mejor explican la
selección de hojas para la construcción de tiendas, se usó el Criterio de Información de
Akaike (AIC) con una corrección para muestras pequeñas (AICc) (Burnham y Anderson
2002). Además, se generó un modelo nulo para determinar si existe relación con las
características evaluadas en los otros modelos o fue solo azar. Una vez seleccionado el
mejor modelo se dividió en cuatro, para poder interpretar los coeficientes en términos de
probabilidad.
3 RESULTADOS
Del total de 46 hojas de C. palmata colectadas,
28 estaban modificadas como tiendas, y en la
comparación a priori de éstas con las no
modificadas se observó que las variables
ángulo de peciolo, distancia al árbol más
cercano y cobertura de dosel no variaron
significativamente.
Las variables que mejor explican la
selección de hojas por parte de los
murciélagos tolderos según el AICc
fue el modelo aditivo (radio
hoja+diámetro peciolo+altura hoja,
AICc=38.7) (Tabla 2), en donde el
radio de la hoja tuvo un coeficiente
de 0.006/4 (Tabla 3). El diámetro
del peciolo tuvo un coeficiente de
1.060/4=0.265, lo cual quiere decir
que por cada unidad que aumenta
el diámetro, aumenta 26.5% la probabilidad de ser seleccionada la hoja, siendo más factible
la selección de hojas cuando tienen un diámetro de 11 mm. Asimismo y de acuerdo con el
coeficiente arrojado para la altura (0.0426/4=0.01) se interpreta que por cada incremento de
la altura en una unidad (cm) se incrementa en 1% la probabilidad de ser seleccionadas
como tiendas, donde hojas de dos metros de altura son las óptimas para ser seleccionadas
4 OBTENCIÓN DE RESULTADOS SIMILARES
Se generó una tabla en Excel con todas las variables de estudio. Con el comando
ALEATORIO.ENTRE se generaron 46 valores aleatorios para cada una de las variables,
gracias a que contábamos con el promedio y la desviación estándar de estas, además
agregamos la variable dependiente “selección” (1 = hoja escogida; 0 = hoja no escogida).
Luego, trabajamos con esta tabla en R.
Fig. 1. Captura de pantalla de valores de las variables generados en Excel pasados a R.
Se generó un modelo de regresión con todas las variables y se fueron descartando una a
una las no significativas, hasta obtener el modelo final, el cual presentaba el menor AIC.
Fig. 2. Captura de pantalla del modelo final generado en R.
● Deviance Residual < Deviance Nulo; por lo tanto el modelo es significativo
● El AIC es el menor de entre todos los modelos que se generaron.
● Todas las variables tienen un p-valor < 0.01 → Variables significativas.
Ecuación del modelo de predicción que resulta de nuestro modelo:
𝑃𝑟(𝑌) = 1
1+𝑒−(−17.6845+0.0541*𝐴𝑙𝑡𝑢𝑟𝑎𝐷𝑒𝐻𝑜𝑗𝑎+0.0093*𝑅𝑎𝑑𝑖𝑜𝐷𝑒𝐻𝑜𝑗𝑎+1.1383*𝐷𝑖á𝑚𝑒𝑡𝑟𝑜𝐷𝑒𝑙𝑃𝑒𝑐𝑖𝑜𝑙𝑜)
5 OTRA METODOLOGÍA
Fig. 3. Clusters obtenidos con R Fig. 4. Número óptimo de clusters.
Se realizó un análisis particional de conglomerados en R en el que esperábamos obtener
dos grupos, uno que agrupe los valores de las variables de las hojas modificadas por el
murciélago, y otro grupo con los valores de las variables de las hojas no modificadas.
En efecto, obtuvimos que el número óptimo de clusters era 2 y agrupaban los valores que
esperábamos.
Adicionalmente realizamos un análisis jerárquico para 2 clusters, con el cual obtuvimos 28
grupos de valores en uno, y 18 en el otro. Precisamente los mismos valores de hojas
modificadas y no modificadas recogidas en campo, lo que refuerza la teoría de que
Dermanura watsoni posee un grado alto de selección de las hojas que usa para refugio.
28 valores 18 valores
Cluster 1 Cluster 2
Hojas modificadas en tiendas Hojas no modificadas
6 REFERENCIAS BIBLIOGRÁFICAS
● Gutiérrez-Sanabria, D. (2014) Selección de hojas de Carludovica palmata como
refugio por el murciélago toldero Dermanura watsoni en la Estación Biológica Piro,
Península de Osa, Costa Rica. En:
https://dialnet.unirioja.es/servlet/articulo?codigo=5168115
https://dialnet.unirioja.es/servlet/articulo?codigo=5168115