#|echo: false
data <- read.csv("ParametrosMolossus_revision.csv", sep = ";")
#str(data)
data$Especie <- as.factor(data$Especie)
data$Tipodellamado <- as.factor(data$Tipodellamado)
data$Autor <- as.factor(data$Autor)Sintesis_Molossus
Introducción
El objetivo de este análisis es graficar y visualizar la variación de los parámetros acústicos reportados en la literatura para las diferentes especies del género Molossus.
Cargar datos
Los datos provienen de la revisión bibliográfica hecha por Geraldine Buitrago y Orlando Hernández. Para cada referencia, se tabuló el promedio y la desviación estándar de cada parámetro acústico.
Datos en formato largo
# A tibble: 6 × 6
Especie Tipodellamado Autor Parametro Media SD
<fct> <fct> <fct> <chr> <dbl> <dbl>
1 Molossus_bondae Low Loureiroetal2020 BW NaN NaN
2 Molossus_bondae Low Loureiroetal2020 CD 12.1 1.8
3 Molossus_bondae Low Loureiroetal2020 EndF NaN NaN
4 Molossus_bondae Low Loureiroetal2020 MaxF 28.5 2.6
5 Molossus_bondae Low Loureiroetal2020 MinF 25.4 2.5
6 Molossus_bondae Low Loureiroetal2020 PI NaN NaN # A tibble: 6 × 6
Especie Tipodellamado Autor Parametro Media SD
<fct> <fct> <fct> <chr> <dbl> <dbl>
1 Molossus_sinaloe Middle Jungetal2014 EndF 34.0 0.9
2 Molossus_sinaloe Middle Jungetal2014 MaxF NaN NaN
3 Molossus_sinaloe Middle Jungetal2014 MinF NaN NaN
4 Molossus_sinaloe Middle Jungetal2014 PI 110. 13
5 Molossus_sinaloe Middle Jungetal2014 PeakF NaN NaN
6 Molossus_sinaloe Middle Jungetal2014 StartF 34.1 0.6Cálculo del consenso entre referencias
# A tibble: 6 × 6
Especie Tipodellamado Autor Parametro Media SD
<fct> <fct> <chr> <chr> <dbl> <dbl>
1 Molossus_bondae Low Loureiroetal2020 BW NaN NaN
2 Molossus_bondae Low Loureiroetal2020 CD 12.1 1.8
3 Molossus_bondae Low Loureiroetal2020 EndF NaN NaN
4 Molossus_bondae Low Loureiroetal2020 MaxF 28.5 2.6
5 Molossus_bondae Low Loureiroetal2020 MinF 25.4 2.5
6 Molossus_bondae Low Loureiroetal2020 PI NaN NaN # A tibble: 6 × 6
Especie Tipodellamado Autor Parametro Media SD
<fct> <fct> <chr> <chr> <dbl> <dbl>
1 Molossus_sinaloe Middle Consenso_refs EndF 34.0 0.9
2 Molossus_sinaloe Middle Consenso_refs MaxF NaN NaN
3 Molossus_sinaloe Middle Consenso_refs MinF NaN NaN
4 Molossus_sinaloe Middle Consenso_refs PI 110. 13
5 Molossus_sinaloe Middle Consenso_refs PeakF NaN NaN
6 Molossus_sinaloe Middle Consenso_refs StartF 34.1 0.6Función para graficar por especie
#|echo: false
graficar_especie <- function(nombre_especie, titulo){
datos_especie <- data_plot %>%
filter(Especie == nombre_especie)
lineas_consenso <- datos_especie %>%
filter(Autor == "Consenso_refs")
ggplot(datos_especie,
aes(x = Autor, y = Media)) +
geom_point(aes(color = color_error),
size = 2.5) +
geom_errorbar(aes(ymin = Media - SD,
ymax = Media + SD,
color = color_error),
width = 0.1) +
geom_hline(data = lineas_consenso,
aes(yintercept = Media),
color = "red",
linetype = "dashed",
linewidth = 0.6) +
scale_color_identity() +
facet_grid(~Parametro ~ Tipodellamado,
scales = "free_y") +
theme_bw(base_size = 12) +
labs(title = titulo,
y = "Media",
x = "Autor") +
theme(
axis.text.x = element_text(angle = 45, hjust = 1),
axis.text.y = element_text(color = "black"),
axis.ticks = element_line(color = "black"),
plot.title = element_text(hjust = 0.5, face = "bold")
)
}Gráficos por especie
graficar_especie("Molossus_molossus", "Molossus molossus")Warning: Removed 37 rows containing missing values or values outside the scale range
(`geom_point()`).
graficar_especie("Molossus_rufus", "Molossus rufus")Warning: Removed 28 rows containing missing values or values outside the scale range
(`geom_point()`).Warning: Removed 2 rows containing missing values or values outside the scale range
(`geom_hline()`).
graficar_especie("Molossus_coibensis", "Molossus coibensis")Warning: Removed 18 rows containing missing values or values outside the scale range
(`geom_point()`).
graficar_especie("Molossus_currentium", "Molossus currentium")Warning: Removed 18 rows containing missing values or values outside the scale range
(`geom_point()`).Warning: Removed 9 rows containing missing values or values outside the scale range
(`geom_hline()`).
graficar_especie("Molossus_bondae", "Molossus bondae")Warning: Removed 8 rows containing missing values or values outside the scale range
(`geom_point()`).Warning: Removed 4 rows containing missing values or values outside the scale range
(`geom_hline()`).
graficar_especie("Molossus_sinaloe", "Molossus sinaloe")Warning: Removed 28 rows containing missing values or values outside the scale range
(`geom_point()`).Warning: Removed 3 rows containing missing values or values outside the scale range
(`geom_hline()`).
graficar_especie("Molossus_pretiousus", "Molossus pretiousus")Warning: Removed 8 rows containing missing values or values outside the scale range
(`geom_point()`).Warning: Removed 4 rows containing missing values or values outside the scale range
(`geom_hline()`).
Síntesis general
#|echo: false
data_sintesis <- data_plot %>%
filter(Autor == "Consenso_refs")
ggplot(data_sintesis,
aes(x = Especie, y = Media)) +
geom_point(size = 2) +
geom_errorbar(aes(ymin = Media - SD,
ymax = Media + SD),
width = 0.1,
color = "red") +
facet_grid(~Parametro ~ Tipodellamado,
scales = "free_y") +
theme_bw(base_size = 14) +
labs(title = "Síntesis del género Molossus",
y = "Media",
x = "Especie") +
theme(
axis.text.x = element_text(angle = 45, hjust = 1),
plot.title = element_text(hjust = 0.5, face = "bold")
)Warning: Removed 22 rows containing missing values or values outside the scale range
(`geom_point()`).
data_wide <- data_sintesis %>%
pivot_wider(
names_from = Parametro,
values_from = c(Media, SD),
names_glue = "{Parametro}_{.value}"
)
data_sintesis_mol <- data_wide[ , - c(3,4)]
data_sintesis_mol# A tibble: 16 × 18
Especie Tipodellamado BW_Media CD_Media EndF_Media MaxF_Media MinF_Media
<fct> <fct> <dbl> <dbl> <dbl> <dbl> <dbl>
1 Molossus_bo… Low NaN 12.1 NaN 28.5 25.4
2 Molossus_co… High 4.7 5.35 35 40.5 35.2
3 Molossus_co… Low 5.6 7.06 29.8 36.7 30.9
4 Molossus_cu… High 3.26 14.5 30.4 NaN NaN
5 Molossus_cu… Low 4.36 14.0 24.5 NaN NaN
6 Molossus_cu… Middle 4.46 14.2 28.3 NaN NaN
7 Molossus_mo… High 3.30 9.68 38.1 41.4 37.1
8 Molossus_mo… Low 5.38 10.4 29.0 38.2 32.4
9 Molossus_mo… Middle 4.17 10.8 36.1 39.7 34
10 Molossus_pr… Low NaN 13.7 NaN 33.6 24.5
11 Molossus_ru… High 7.68 15.6 25.8 37.8 30.4
12 Molossus_ru… Low 4.33 13.2 24.8 31.7 26.4
13 Molossus_ru… Middle 6.5 16.4 NaN 30.5 24
14 Molossus_si… High 1.58 6.82 43.7 38.3 35.5
15 Molossus_si… Low 1.48 8.57 38.0 32.6 27.0
16 Molossus_si… Middle 0.26 7.66 34.0 NaN NaN
# ℹ 11 more variables: PI_Media <dbl>, PeakF_Media <dbl>, StartF_Media <dbl>,
# BW_SD <dbl>, CD_SD <dbl>, EndF_SD <dbl>, MaxF_SD <dbl>, MinF_SD <dbl>,
# PI_SD <dbl>, PeakF_SD <dbl>, StartF_SD <dbl>