Este informe corresponde al Quiz #8 de la asignatura. El objetivo principal es realizar un análisis exploratorio de la abundancia y diversidad de morfotipos de artrópodos recolectados en distintas zonas de muestreo (GAM, Zona Fría y Zona Caliente) mediante trampas de plato.
A continuación se muestra una representación gráfica de la diversidad de artrópodos estudiada en biología de campo.
Cargamos los datos directamente desde el archivo de Excel
Hoja datos Diversidad Bioesta I (1).xlsx.
# Cargar la hoja principal de datos
ruta_excel <- "Hoja datos Diversidad Bioesta I (1).xlsx"
datos_raw <- read_excel(ruta_excel, sheet = "GAM")
# Limpieza básica de la base de datos
datos <- datos_raw %>%
filter(!is.na(ZONA)) %>%
mutate(across(starts_with("Morfotipo"), ~ replace_na(.x, 0)))
# Mostrar los primeros registros de la base de datos
kable(head(datos[, 1:6]), caption = "Primeros registros de la base de datos de diversidad")
| ZONA | Platos | MF | CANTIDAD | peso | Morfotipo 2: Escarabajos, abejones, mariquita (Coleópteros) |
|---|---|---|---|---|---|
| GAM | Plato día 1 | MT1 | 3 | NA | 0 |
| GAM | Plato día 2 | NA | 5 | NA | 1 |
| GAM | Plato Día 3 | NA | NA | NA | 0 |
| GAM | Plato Día 4 | NA | 1 | NA | 0 |
| GAM | Plato día 5 | NA | 4 | NA | 0 |
| GAM | Plato noche 1 | NA | NA | NA | 0 |
Calculamos la abundancia total por punto de muestreo y su resumen según la zona ecológica:
# Sumar la abundancia total de morfotipos por fila
columnas_morfotipo <- grep("Morfotipo|morfotipo", colnames(datos), value = TRUE)
datos$Abundancia_Total <- rowSums(datos[, columnas_morfotipo], na.rm = TRUE)
# Resumen estadístico agrupado por Zona
resumen_zona <- datos %>%
group_by(ZONA) %>%
summarise(
Muestreos = n(),
Promedio_por_Plato = round(mean(Abundancia_Total), 2),
Desviacion_Estandard = round(sd(Abundancia_Total), 2),
Abundancia_Total = sum(Abundancia_Total)
)
kable(resumen_zona, caption = "Resumen de abundancia total de artrópodos según zona")
| ZONA | Muestreos | Promedio_por_Plato | Desviacion_Estandard | Abundancia_Total |
|---|---|---|---|---|
| GAM | 20 | 2.3 | 3.11 | 46 |
| ZC | 10 | 0.0 | 0.00 | 0 |
| ZF | 10 | 3.1 | 3.07 | 31 |
Generamos un gráfico con ggplot2 para comparar la
distribución de abundancias según la zona de muestreo.
ggplot(datos, aes(x = ZONA, y = Abundancia_Total, fill = ZONA)) +
geom_boxplot(alpha = 0.7, outlier.color = "red", outlier.shape = 16) +
geom_jitter(width = 0.1, alpha = 0.5) +
scale_fill_brewer(palette = "Set2") +
theme_minimal() +
labs(
title = "Comparación de Abundancia de Artrópodos por Zona",
subtitle = "Muestreo realizado con trampas tipo plato",
x = "Zona Ecológica",
y = "Abundancia Total de Individuos",
fill = "Zona"
)
El gráfico permite comparar la distribución de la abundancia total de artrópodos entre las diferentes zonas de muestreo, identificando diferencias en la cantidad de individuos recolectados.
# Transformar datos a formato largo para graficar por morfotipo
datos_largos <- datos %>%
pivot_longer(
cols = all_of(columnas_morfotipo),
names_to = "Morfotipo",
values_to = "Conteo"
) %>%
mutate(Morfotipo_Corto = gsub(":.*", "", Morfotipo)) %>%
group_by(Morfotipo_Corto) %>%
summarise(Total = sum(Conteo)) %>%
filter(Total > 0)
ggplot(datos_largos, aes(x = reorder(Morfotipo_Corto, Total), y = Total, fill = Morfotipo_Corto)) +
geom_col(show.legend = FALSE) +
coord_flip() +
theme_minimal() +
labs(
title = "Abundancia Total por Morfotipo Registrado",
x = "Morfotipo",
y = "Número Total de Individuos"
)
Se observa la abundancia registrada para cada morfotipo, lo que permite identificar cuáles fueron los grupos más representativos durante el muestreo.