Similaridad de 11 especies pertenecientes a las familias Falconidae y Strigidae a partir de sus vocalizaciones

 

Similarity of 11 species belonging to the families Falconidae and Strigidae from their vocalization

 

Javier Armando Azabache Requena1

 

1  Universidad Ricardo Palma, Av. Alfredo Benavides 5440, Santiago de Surco 15039. Lima, Perú.  

 

* Autor corresponsal:  javierazabache7@outlook.com (J. A. Azabache Requena).

 

ID ORCID de los autores

J. A. Azabache Requena:  http://orcid.org/0000-0001-9889-074X

 

 

 

RESUMEN

 

Conocer aspectos biológicos de las familias Falconidae y Strigidae como la vocalización es relevante, ya que podríamos utilizar este aspecto como base para agrupar especies y evaluar su similaridad. En ese sentido se caracterizaron las vocalizaciones de 11 especies pertenecientes a las familias Falconidae y Strigidae, a través de la descarga de 53 audios por especie de la base de datos de acceso público Xeno Canto. Las variables consideradas en este estudio fueron: frecuencia, amplitud, duración, números de picos de frecuencia, factor cresta, mediana de envolvente de amplitud e índice de complejidad acústica. Una vez caracterizadas las vocalizaciones se procedió a evaluar la similaridad de especies, utilizando para ello el agrupamiento jerárquico (índice de Conectividad = 9,415, índice de Dunn = 0,819,). Guardando similaridad a partir de las características de sus vocalizaciones Falco peregrinus y F. columbarius, F. rufigularis y F. femoralis; mientras que F. sparverius tuvo más similaridad con las 2 primeras especies. Todas estas especies se diferencian de las especies de la familia Strigidae teniendo aquí Glaucidium brasilianum, G. jardini y G. peruanum más similaridad; asimismo Glaucidium hardyi y Aegolius harrisii guardan más concordancia entre ellos y finalmente Megascops roboratus presenta similaridad con todas las especies del género Glaucidium (coeficiente de aglomeración = 0,67).

 

Palabras clave: Análisis por conglomerados; vocalización animal; aves.

 

 

ABSTRACT

 

Knowing biological aspects of the families Falconidae and Strigidae as vocalization is relevant, since we could use this aspect as a basis to group species and evaluate their similarity. In this sense, the vocalizations of 11 species belonging to the Falconidae and Strigidae families were characterized, through the download of 53 audios per species from the Xeno Canto public access database. The variables considered in this study were: frequency, amplitude, duration, number of frequency peaks, crest factor, mean amplitude envelope and index of acoustic complexity. Once the vocalizations were characterized, the species similarity was evaluated, using hierarchical grouping (connectivity index = 9.415, Dunn index = 0.819). Keeping similarity from the characteristics of their vocalizations Falco peregrinus and F. columbarius, F. rufigularis and F. femoralis; while F. sparverius had more similarity with the first 2 species. All these species differ from the species of the family Strigidae having here Glaucidium brasilianum, G. jardini and G. peruanum more similarity; Likewise, Glaucidium hardyi and Aegolius harrisii have more concordance between them and finally Megascops roboratus presents similarity with all species of the genus Glaucidium (agglomeration coefficient = 0.67).

 

Keywords:  Cluster analysis; animal vocalization; birds.

 

 

Recibido: 31-08-2023.

Aceptado: 11-12-2023.


  
Esta obra está publicada bajo la licencia CC BY 4.0
INTRODUCCIÓN

 


En general los animales utilizan señales acústicas para comunicarse, lo que implica que las propiedades de estas señales pueden estar sujetas a una fuerte selección. Se conoce también que la adaptación acústica predice que las especies en hábitats densos emiten sonidos de menor frecuencia que las de áreas abiertas porque los sonidos de baja frecuencia se propagan más lejos en la vegetación densa que los sonidos de alta frecuencia. La frecuencia de la señal también puede estar bajo selección sexual porque se correlaciona con el tamaño del cuerpo (Mikula et al., 2020).

En ese sentido para garantizar una comunicación acústica eficaz, las señales deben llegar a los receptores de la forma menos distorsionada posible. Los animales utilizan diversas estrategias a corto y largo plazo para evitar la degradación y el enmascaramiento de la señal (Budka, 2023).

Aunque muchas formas de comunicación animal no dependen del orden en que se combinan los componentes de las señales para que sean efectivas, existe evidencia de que el orden sí importa para algunos sistemas de comunicación (Campbell et al., 2023).

En cuanto a las aves emiten sonidos a partir de la siringe los cuales se llaman vocalizaciones, estas se clasifican en llamadas y cantos. En general, los llamados son sonidos de corta duración formados por una o pocas sílabas simples o medianas en complejidad. Los llamados son producidos tanto por machos como por hembras a lo largo del año, generalmente son innatos en todos los grupos de aves y normalmente pueden relacionarse con funciones específicas como vuelo, amenaza, alarma, forrajeo o contacto con parientes u otros miembros de la misma especie, hasta ahora se sabe que los cantos son aprendidos solamente en tres grupos de aves: aves canoras(oscinos), pericos y colibríes. Para poder estudiar los cantos en las aves y cualquier otra señal acústica, es necesario medirlos, esto se hace a partir de los espectro-gramas, que son representaciones visuales del sonido. Las características o medidas más usadas son la frecuencia mínima, frecuencia máxima, frecuencia pico y ancho debanda, las cuales miden la variación en el tono de los cantos, las características temporales más usadas son la duración del canto, el tiempo entre cantos o entre sílabas y tasa de canto. Las propiedades del sonido de los cantos de las aves nos ayudan a identificar a las especies, incluso ciertas propiedades pueden corresponder a determinados grupos de aves. (Moran & Gonzales, 2022).

Las vocalizaciones que producen las aves son una herramienta que ha sido ampliamente usada, por su papel en la delimitación e identificación de especies. En ese sentido representan un carácter diagnóstico que facilita la discriminación de entidades taxonómicas, esto ha sido de gran utilidad en la confirmación y registro de especies que dados sus comportamientos y requerimientos de hábitat son de difícil observación (Fernández et al., 2023).

De acuerdo con Griffiths & Aaronson (2023), la variación en la vocalización de especies del género Falco puede ser sólo el resultado de la variación en la estructura anatómica, siendo el tamaño el factor principal. Los autores encontraron que cuanto más grande es el ave, menor es la frecuencia máxima.

Por otro lado, las técnicas de minería de datos provienen de la inteligencia artificial y de la estadística, dichas técnicas son plasmadas en algoritmos, que después se aplican sobre un conjunto de datos para obtener resultados (Peralta et al., 2022). Dentro de la minería de datos se usa el análisis clúster el cual engloba una serie de técnicas que pueden ser usadas para crear agrupamientos, más concretamente un método clúster es un proce-dimiento estadístico multivariante que comienza con un conjunto de observaciones y los reorganiza en grupos o clúster con las siguientes cara-cterísticas: a) cada grupo o clúster es homogéneo o compacto con respecto a ciertas características, es decir las observaciones en cada grupo son similares entre sí,b) cada clúster ha de ser diferente de los otros con respecto a las mismas características, es decir las observaciones de un grupo serán diferentes de las observaciones de los otros grupos y estas diferencias serán mayores cuantos más alejados estén los clústers entre sí. Parte de los anális clúster son los llamados métodos jerár-quicos; los cuales se caracterizan fundamen-talmente porque tratan de agrupar varios clústers para formar uno nuevo o bien para dar origen a otros dos partiendo de uno ya existente, de tal forma que, si sucesivamente se va realizando este proceso de aglomeración o división, se minimicen las distancias o se maximicen similitudes entre los grupos (Broncano de Pablo, 2018).

En ese sentido el objetivo de este trabajo fue conocer que especies de las familias Falconidae y Strigidae guardan similaridad entre si al pertenecer a un mismo grupo después de aplicar algoritmos de aprendizaje no supervisado (análisis clúster), tomando como base las características cuantitativas de sus vocalizaciones.


 

MATERIAL Y MÉTODOS

 


Se descargaron grabaciones de la base de datos de acceso abierto Xeno canto de las especies: Falco sparverius, Falco rufigularis,  Falco femoralis, Falco peregrinus, Falco columbarius, (Familia Falconidae), Glaucidium peruanum, Glaucidium jardinii, Glaucidium hardyi, Megascops roboratus, Aegolius harrisii y Glaucidium brasilianum (Familia Strigidae).Debido a que hubo diferente cantidad de grabaciones de cada una de las especies y para poder analizar igual cantidad de grabaciones se realizó un muestreo no probabilístico por conveniencia de 53 audios, se eligió esta cantidad ya que la especie Megascops roboratus es la que tiene la menor cantidad de audios los cuales son 53.

Una vez descargados las grabaciones, desde el 3 de enero al 3 de diciembre del 2022 con ayuda del programa Audacity se limpiaron los audios de cualquier ruido de fondo con la opción efecto EQ y filtros y ecualizador gráfico. Una vez limpios los audios se obtuvieron datos de frecuencia mínima, media y máxima, duración mínima, media y máxima, amplitud mínima, media y máxima, detección de picos de frecuencia, factor cresta, índice de complejidad acústica y mediana de envolvente de amplitud. Estos se calcularon con la ayuda de la librería Seewave (Sueur et al., 2008) del programa R y el programa Raven pro Lite y se colectaron en una matriz de datos. En esta matriz de datos las columnas estuvieron formadas por las variables: frecuencia mínima, media y máxima, duración mínima, media y máxima, amplitud mínima, media y máxima, picos de frecuencia, factor cresta, índice de complejidad acústica y mediana de envolvente de amplitud. Y las filas por las especies estudiadas.

Con la matriz de datos, en el programa R se realizó un arreglo para que el nombre de las especies el programa no lo lea como una variable cualitativa. Para lo cual se reordenó el nombre de las especies

en una columna aparte, para que de esa manera solo queden las variables cuantitativas. Esto se realizó con la ayuda de la librería textshape (Rinker, 2021). Luego se prepararon los datos a través de un escalamiento. Asimismo, se midió la correlación entre las variables con la librería PerformanceAnalytics (Peterson & Carl, 2020) de R. Posteriormente se evaluó la tendencia de los datos al agrupamiento de acuerdo con el esta-dístico de Hopkins, con la ayuda del paquete Hopkins (Wright, 2022) del programa R, el cual establece que valores entre 0 - 0,3 indican datos espaciados regularmente, datos alrededor de 0,5 indican datos aleatorios y valores entre 0,7 - 1 indican datos agrupados. Luego con los algoritmos: kmeans, PAM, CLARA y clúster jerárquico, y con la ayuda del paquete clValid (Brock et al., 2008) de R se eligió el más adecuado, utilizando para este fin los índices de conectividad, de Dunn y la silueta teniendo en cuenta que el algoritmo más adecuado para aplicar a los datos debe minimizar el primer índice y maximizar los últimos.

Con las librerías clúster (Maechler et al., 2022) y facto extra (Kassambara & Mundt, 2020) se procedió a aplicar el algoritmo más idóneo, que para estos datos fue: el clúster jerárquico, específicamente el aglomerativo y se comparó cual, entre los métodos single, complete y Ward era el más idóneo, usando para ello el coeficiente de aglomeración, teniendo en cuenta que el coeficiente más cercano a 1 es el que mejor agrupa a las observaciones. Una vez obtenido el mejor método de clúster jerárquico aglomerativo, se procedió a graficarlos. Todo esto se realizó con las funciones hclust y agnes de la librería clúster.


 

 

RESULTADOS Y DISCUSIÓN

 


De acuerdo con la Tabla 1, cada especie y cada familia presenta sus características específicas de frecuencia, amplitud y duración. De acuerdo con Moran & Gonzales (2022), las propiedades del sonido de los cantos de las aves nos ayudan a identificar las especies, incluso ciertas propiedades pueden corresponder a determinados grupos de aves. Por ejemplo, las palomas (familia Columbidae) emiten cantos de bajas frecuencias (graves), con sílabas largas, simples (sin modulaciones de frecuencia) y de baja intensidad. Los saltaparedes (familia Troglodytidae) emiten cantos muy diversos, pero algunas propiedades que los caracterizan son trinos y sílabas con modulaciones de frecuencia en un rango aproximado de 1 a 10 kHz, según sea la especie y generalmente son intensos. Otros trabajos también corroboran lo mencionado, por ejemplo, López (2015), quien encontró diversas frecuencias, amplitudes y duraciones en los cantos de los distintos taxa del género Bubo en América.

Esta diferencia de las características de las vocalizaciones no solo se aprecia entre las aves sino también en otros grupos como anfibios (Batallas & Brito, 2016) o en los quirópteros (Malo de la Molina, 2011). Las duraciones, amplitudes y frecuencias también están relacionados con el comportamiento de la especie a estudiar.

Orellana (2014) encontró diversas voces del zorzal andino (Turdus chiguanco) cada voz asociada a un comportamiento y a una frecuencia, amplitud y duración característica.


 

Tabla 1

Frecuencia(f), amplitud(A), duración(T), Picos de frecuencia(fpeak), factor cresta(C), índice de complejidad acústica (ACI) y mediana de envolvente de amplitud(M) de las vocalizaciones de las especies estudiadas

 

Especie

f. mínima

f.

media

f.

máxima

A.

mínima

A.

media

A.

máxima

T.

mínima

T.

media

T.

máxima

fpeak

C

ACI

M

F. sparverius

2,13

4,13

11,3

31,87

43,45

53,85

0,57

2,35

56,91

40

6,84

195,51

0,021

F. rufigularis

2

9,6

13,5

35,12

46,52

58,25

1,14

29

66,48

38

14,48

178,2

0,018

F. femoralis

1,44

6,44

11,93

36,09

48,47

64,64

1,59

8,96

42,09

42

13,4

197,26

0,007

F. peregrinus

1,53

4,82

8,79

33,76

45,07

54,65

2,3

4,67

43,78

40

8,67

158,48

0,015

F. columbarius

1,52

4,79

8,65

33,86

48,13

56,85

0,88

4,1

25,02

43

9,75

171,28

0,012

M. roboratus

0,51

0,79

1,86

22,18

32,84

51,56

0,18

2,15

4,35

54

4,8

180,07

0,03

G. jardinii

1,03

1,84

2,83

27,99

36,22

60,44

3

17,56

60,41

49

4,51

171,87

0,018

G. hardyi

0,88

1,32

1,84

29,72

40,25

63,57

9,3

28,2

47,36

58

3,36

174,16

0,039

G. peruanum

1

2,15

3,44

32,31

50,37

61,03

3,17

21,52

53,62

44

4,44

177,71

0,013

G. brasilianum

0,91

1,23

1,77

31,39

37,52

53,25

3,3

10

95

54

4,2

162,6

0,02

A. harrisii

0,78

1,02

1,45

25,75

33,71

45,18

6

28,5

75

55

5,74

187,59

0,015

 

Nota: Las especies pertenecientes a la familia Falconidae presentaron mayor frecuencia, factor cresta, índice de complejidad acústica y similar mediana de envolvente de amplitud en comparación con las especies de la familia Strigidae.

 


En la Tabla 1 también podemos observar que unas especies tienen frecuencias más altas que otras. De acuerdo con Arévalo (2019), estos cambios se deben al hábitat donde se encuentra el ave, así aquellas que viven en zonas urbanas presentan una frecuencia más alta esto se pudo apreciar por ejemplo en Falco sparverius y Falco rufigularis quienes viven en zonas urbanas, además de las zonas donde viven estos animales también influye la masa corporal la cual está relacionada con las frecuencias mínimas de canto. Wheeldon et al. (2022), sin embargo, mencionan que los cantos de Laniarius atroflavus se propagan mejor en los bosques que en sitios perturbados cubiertos de arbustos y los tipos de canciones masculinas utilizadas para el canto transmitido se propagan más lejos que las canciones femeninas. Por otro lado, Miluka et al. (2020) encontraron en un estudio global de paseriformes que la frecuencia de las canciones disminuye con el aumento de la masa corporal y con el dimorfismo de tamaño sexual sesgado por los machos; mas no encontró relación prevista entre frecuencia y hábitat.

De acuerdo con la Tabla 1 existen más picos de frecuencia en la familia Strigidae, de acuerdo con Loayza (2021) esta métrica funciona bajo el supuesto de que a mayor cantidad de picos de frecuencia existe una mayor actividad acústica.

En la Figura 1 el factor cresta (C) se relaciona con las frecuencias mínimas (0,70), máximas y medias (0,90) positivamente, También podemos observar que las frecuencias, duraciones y amplitudes mínimas y medias tienen relaciones positivas entre ellas. El factor cresta mide cuanto fluctúa la amplitud y por lo tanto podemos mencionar que en los Falconiformes la amplitud de la señal varía mucho o lo que se traduciría de acuerdo Morán & Gonzales (2022) que el volumen de la señal varía. Por otro lado, este factor al guardar una correlación positiva con la frecuencia medía y máxima nos da la idea que estaría relacionado con tonos de las vocalizaciones.

 

El índice de complejidad acústica (ACI) es mayor en unas especies más que en otras y en la familia Falconidae más que en la Strigidae según Lozano et al. (2014) al analizar los ACI, pero a nivel de especies no de comunidad mencionan que el mayor ACI por especie puede deberse al repertorio acústico que es muy variado. Por otro lado, Pieretti et al. (2011) encontraron el ACI relacionado con la cantidad de notas de las vocalizaciones (correlación de Spearman = 0,94) más que de las frecuencias. En este estudio tampoco se encontraron correlaciones del ACI con las frecuencias.

 

En este trabajo el estadístico de Hopkins resultó ser de 0,90 lo que indica que esta data es clusterizable. La mayoría de los índices indican que el mejor algoritmo para los datos es la agrupación jerárquica (Tabla 2).

 

Tabla 2

Scores obtenidos en la validación del mejor algoritmo para los datos de las vocalizaciones

 

Índice

Score

Método

Conectividad

Dunn

Silueta

9,415

0,810                                   

0,295

Hierarchical

Hierarchical

K-means

Nota: 2 de los 3 índices señalan a clúster jerárquico como el mejor algoritmo para aplicar a nuestros datos.

 

Además, en el algoritmo de clúster jerárquico el método que presentó mayor coeficiente de aglomeración fue el Ward (0,67), mientras que el método single y completo sus coeficientes fueron de: 0,284 y 0,582 respectivamente. Por lo que se aplicó el primer método.


 

 

Figura 1. Correlación de variables.

Figura 2. Clúster jerárquico de las especies.

 


De acuerdo con la Figura 2, a partir de las características de sus vocalizaciones las 11 especies estudiadas se juntaron con miembros de su propia familia si lo comparamos con la lista de aves del Perú (Plenge, 2023). De acuerdo con Moreno (2006), las vocalizaciones permiten la comunicación de las aves; asimismo Vásquez & Vásquez (2018) señalan también que el canto permite la comunicación de individuos de la misma especie o entre especies distintas, buscando también modificar el compor-tamiento del individuo que la recibe, por ejemplo, Campbell et al. (2022) encontraron que individuos de menos un año de Poecile atricapillus son capaces de percibir y actuar según categorías de llamadas. Otro ejemplo de cómo el canto influye en el comportamiento lo señala Romani et al. (2023) quienes encontraron que los adultos de Ardea cinerea utilizan una amplia gama de tipos de llamadas, principalmente relacionadas con el cortejo y la fase reproductiva. Las vocalizaciones son emitidas especialmente por los machos y son un elemento clave durante la atracción de pareja, la construcción de nidos, en los eventos de reproducción y en el fortalecimiento de los vínculos dentro de las parejas. Por lo que agrupación resultante producto de la similaridad de sus vocalizaciones sería de alguna manera para asegurar su supervivencia en la naturaleza. Por otro lado, Vásquez & Vásquez (2018) también encontraron similitud en las vocalizaciones de las especies de diferentes familias las cuales fueron Mimus longicaudatus (Mimidae) y Campylorhynchus fasciatus (Troglodytidae) en especial en sus frecuencias máximas y mínimas.

En la Figura 2, las especies de la familia falconidae se separan de la Strigidae, esto se podría explicar  de acuerdo a Batallas & Brito(2016) a condiciones ambientales y de contexto social. Calderón (2022) concluyó que el factor morfológico es relevante tanto por la producción de nuevas vocalizaciones aumentando y disminuyendo los repertorios en cinco especies de la familia Rhinocryptidae, incluso estaría afectado a los rangos de actividad junto a los aspectos reproductivos de estas aves debido a los costos energéticos que demanda la reproducción que no se pueden sobreponer con las vocaliza-ciones. En 7 especies del genero Falco, Griffiths & Aaronson (2023) encontraron que  la variación en la vocalización de estas especies puede ser sólo el resultado de la variación en la estructura anató-mica, siendo el tamaño el factor principal. Los autores encontraron que cuanto más grande es el ave menor es la frecuencia fundamental y más lenta es la tasa de cambio de frecuencia (pendiente de frecuencia más pequeña) dentro de las notas, también encontraron una fuerte correlación entre la pendiente de frecuencia y  de las estructuras de la siringe. Vargas (2021) evaluó la comunidad de aves con relación al paisaje sonoro como una propuesta para el establecimiento de un Corredor Biológico Interurbano en la microcuenca del río Bermúdez en Heredia, Costa Rica, caracterizando la variación acústica de cada sitio mediante índices de paisaje sonoro que permitieron interpretar la interacción acústica de las biofonías, geofonías y antropofonías a lo largo de la microcuenca con los que se logró evidenciar la relación de las variables acústicas con las variables biológicas. En ese sentido, las vocalizaciones podrían sugerirse como un elemento para establecer semejanzas y dife-rencias entre especies o familias, concordonado con Fernández et al. (2022), para quienes los sonidos que producen los animales, son una herramienta que ha sido ampliamente usada en el estudio de las aves por su papel en la delimitación e identificación de especies; es decir  representan un carácter diagnóstico que facilita la discrimi-nación de entidades taxonómicas.


CONCLUSIONES

 


A partir de las características cuantitativas de las vocalizaciones después de aplicar el algoritmo de clúster jerárquico, las especies que guardaron similitud fueron: Falco peregrinus y F. columbarius, F. rufigularis y F. femoralis; mientras que F. sparverius tuvo más similaridad con las 2 primeras especies. Todas estas especies se diferencian de las especies de la familia Strigidae teniendo aquí Glaucidium brasilianum, G. jardini y G. peruanum más similaridad; asimismo Glaucidium hardyi y Aegolius harrisii guardan más concordancia entre ellos y finalmente Megascops roboratus presenta similaridad con todas las especies del género Glaucidium (coeficiente de aglomeración = 0,67).

En este estudio se analizaron 53 audios por especie, pero quizás en futuras investigaciones se podrían aumentar el número de audios analizados; pero teniendo en cuenta cantos colectados de museos de historia natural ya que sería una forma de asegurarnos de tener cantos estándar para un mejor análisis y además considerar separar las vocalizaciones de machos y hembras.


 

 

AGRADECIMIENTOS

 


Mi agradecimiento al comité editorial de la Revista Manglar por la revisión para la mejora de este trabajo y a la Licenciada en idiomas Karla González por su apoyo desinteresado.

 


 

REFERENCIAS BIBLIOGRÁFICAS

 


Arévalo, J. (2019). Variación morfológica y del canto en poblaciones de Zonotrichia capensis (Aves: Emberizidae) de ambientes urbanos y silvestres en Chile Central. [Tesis para optar el grado de Magister en Ciencias biológicas]. Universidad de Chile.

Batallas, D., & Brito J. (2016). Análisis bioacústico de las vocalizaciones de seis especies de anuros de la laguna Cormorán, complejo lacustre de Sardinayacu, Parque Nacional Sangay, Ecuador. Revista Mexicana de Biodiversidad, 87, 1292–1300.

Brock, G., Pihur, V., & Datta S. (2008). clValid: An R Package for Cluster Validation. Journal of Statistical Software, 25(4), 1–22.

Broncano de Pablo, G. (2018). Aplicación de análisis cluster jerárquicos al proceso de estimación de tráficos por el método de exposición cuasi-inducida. [Tesis para optar el título de Ingeniería Industrial]. Universidad Politécnica de Madrid.

Budka, M., Staniewicz, A., & Sokołowska, E. (2023). Interspecific avoidance of song overlap in tropical songbirds: species-specific responses to acoustically similar and different intruders. Behavioral Ecology and Sociobiology, 77, 1-11.

         https://doi.org/10.1007/s00265-023-03356-0

Calderón, S. (2022). Efecto de la estructura del hábitat en las vocalizaciones y estudio de repertorio y actividad vocal de cinco especies de la familia Rhinocryptidae. [Tesis para optar al grado de Magíster en Ciencias con mención en Zoología]. Universidad de Concepción.

 Campbell, K., Hoeschele, M., Mann, D., Congdon, J., Scully, E., Mischler, S., Montenegro, C., Service, W., & Sturdy, C. (2023).         Black-capped chickadees (Poecile atricapillus) discriminate between naturally-ordered and scramble-ordered chick-a-dee calls and individual preference is related to rate of learning, Behavioural Processes, 206, 1-9.

        https://doi.org/10.1016/j.beproc.2023.104842

Griffiths, C. & Aaronson, N. (2023). Analysis of vocal communication in the genus Falco. Scientific Reports,13(1),1-12.https://doi.org/10.1038/s41598-023-27716-y

Fernández, R., Ku, W., Botero, D., Niño, N., Laverde, O. Pantoja, H., Bravo, G., Álvarez, M., Marín, O., Duque, F. & Ocampo.N. (2023). La voz de nuestras aves: contribuciones de la bioacústica a la ornitología colombiana. Ornitología Colombiana, 23,3-30 https://doi.org/10.59517/oc.e555

 Kassambara, A. & Mundt, F. (2020). Factoextra: Extract and Visualize the Results of Multivariate Data Analyses. https://acortar.link/K0EsoA

López, B. (2015). Análisis comparativo de las vocalizaciones de distintos taxa del género bubo en américa. Hornero, 30(2), 69–88.

Loayza, J. (2021). Influencia del tipo de grabadora en los índices acústicos para análisis automático de cambios en ecosistemas. [Tesis para optar el título de Ingeniera electrónica]. Universidad de Antioquia.

Lozano, A., Farina, A., & Márquez, R. (2014). ACI (Acoustic Complexity Index): Nueva herramienta para el estudio del canto de anuros. Quehacer Científico en Chiapas, 9(2), 17-27.

Maechler, M., Rousseeuw, P., Struyf, A., Hubert, M., & Hornik, K. (2022). cluster: Cluster Analysis Basics and Extensions. https://CRAN.R-project.org/package=cluster.

Malo de la Molina, J., Velazco, S., Pacheco, V., & Robledo, J. (2011). Análisis de las vocalizaciones del murciélago longirrostro peruano Platalina genovensium Thomas, 1928 (Chiroptera: Phyllostomidae). Rev. Perú. biol., 18(3), 311 – 318.

Mikula, P., Valcu, M., Brumm, H., Forstmeier, W., Petruskova, T., Kempenaers, B., & Albrecht, T. (2020). A global analysis of song frequency in passerines provides no support for the acoustic adaptation hypothesis but suggests a role for sexual selection. Ecology Letters 24, 477–486. https://doi.org/10.1111/ele.13662

Morán, D., & González, C. (2022) El mejor lugar para cantar: hipótesis de adaptación acústica en aves. Elementos, 126, 81-88.

Moreno, G. (2006). El comportamiento de las aves como herramienta para su identificación. Acta Granatense, 4, 85-93.

Orellana, M. (2014). Estudio de la bioacústica del zorzal andino (Turdus chiguanco) con la aplicación del software Audacity, v.2.0. Apunt. cienc. soc, 4(2), 259-266.

Peralta, M., Merma, J., Chavez.E., Soto.C., & Jiménez, W. (2022). Aplicación del Algoritmo K-Means en estudiantes universitarios del área de sistemas e informática para caracterizar la salud mental, durante el aislamiento en COVID-19. https://www.iiis.org/DOI2022/CA001MJ/

Peterson, B., & Carl, P. (2020). Performance Analytics: Econometric Tools for Performance and Risk Analysis_ R package version 2.0.4.

Plenge, M. A. Versión [2023] Lista de las aves del Perú. Unión de Ornitólogos del Perú: https://sites.google.com/site/boletinunop/checklist

Pieretti, N., Farina, A., & Morri, D. (2011). A new methodology to infer the singing activity of an avian community: The Acoustic Complexity Index (ACI). Ecological Indicators, 11, 868-873.

Rinker, T. (2021). Textshape: Tools for Reshaping Text. version 1.7.3, https://github.com/trinker/textshape.

 Romani, F., Ramella, E., & Pavan, G. (2023). Repertoire and vocal behaviour of the grey heron (Ardea cinerea) during the breeding season: a case study in the western Po Valley, Ethology Ecology & Evolution.        https://doi.org/10.1080/03949370.2023.2248067

Sueur, J., Aubin, T., & Simonis, C. (2008). Seewave: a free modular tool for sound analysis and synthesis. Bioacoustics, 18, 213-226.

Vargas, R. (2021). Caracterización del paisaje sonoro asociado a la comunidad de aves como una propuesta para el establecimiento de un Corredor Biológico Interurbano en la microcuenca del río Bermúdez en Heredia, Costa Rica. [Tesis para optar al grado de Licenciatura en Biología con Énfasis en Manejo de Recursos Naturales].

Vásquez, R., & Vásquez, D. (2018). Patrón de vocalización de Mimus longicaudatus en ámbito urbano. [Tesis para obtener el título de biólogo]. Universidad Nacional de Trujillo.

Wheeldon, A., Kwiatkowska, K., Szymański, P., & Osiejuk, T. (2022) Male and female songs propagation in a duetting tropical bird species in its preferred and secondary habitat. PLOS ONE, 17(10). https://doi.org/10.1371/journal.pone.0275434

Wright, K. (2022). Will the Real Hopkins Statistic Please Stand Up?, The R Journal. https://doi.org/10.32614/RJ-2022-055