Presencia de cuatro especies de aves rapaces (Accipitridae, Falconidae y Pandionidae) en diferentes zonas de Costa Rica

Fiorella Gazo Bryan1, Pablo Gudiel Quesada1 y Jennifer Hernández Chavarría1
1 Escuela de Ciencias Biológicas, Universidad Nacional, Heredia, Costa Rica; , ,

Resumen: Los rapaces suelen encontrarse en sitios variados desde los bosques hasta las zonas urbanas, siendo estos últimos espacios altamente modificados y alterados por la actividad del ser humano. La humedad y la temperatura ambientales presentaron una relación directamente proporcional, mientras que la talla y la abundancia de las especies de rapaces evidencian una relación inversamente proporcional. Por ende, se pretende analizar las variables que determinan la presencia de aves rapaces en distintos sitios en Costa Rica, determinar si existen diferencias entre la abundancia de las especies según el sitio estudiado y, de forma adicional, comparar la tasa de crecimiento de polluelos de Falco peregrinus con dos tipos de alimentación. Sin embargo, la reinserción de rapaces al medio natural es lo suficientemente complicada como para tomar en cuenta solo la temperatura y la alimentación de los polluelos. La fragmentación y destrucción de los ambientes naturales ha provocado que las aves generalistas busquen adaptarse a las zonas urbanas; mientras que las especialistas luchan por sobrevivir a los cambios generados por el ser humano.

Palabras clave: Aves migratorias; Falco peregrinus; Rapaces diurnos; Repoblación de fauna silvestre; Zonas urbanas

INTRODUCCIÓN

Los rapaces son aves que generalmente presentan adaptaciones a su hábito alimenticio carnívoro y a diversos ambientes. Suelen categorizarse en generalistas y especialistas, donde las primeras tienen una dieta constituida por más de un tipo de alimento y las especialistas se restringen a uno solo (Vides-Hernández et al., 2017, p. 276). En el caso de Costa Rica, suelen encontrarse en sitios variados desde los bosques hasta las zonas urbanas. Mientras que en los bosques la vegetación es densa y con gran cantidad de árboles con alturas mínimas entre 2 y 5 metros en su madurez (Rojas & Iza, 2009, p. 264), los pastizales se caracterizan por un bajo número de estos y muchos arbustos. Por otra parte, los humedales son áreas en las que el agua es el factor principal que determina la vida vegetal y animal dentro del mismo. Las zonas urbanas son espacios modificados y alterados por la actividad del ser humano. En los sitios mencionados anteriormente se puede encontrar una variedad considerable de rapaces.

Buteo jamaicensis (Gavilán colirrojo). Rapaz generalista y solitario con tamaños entre los 48 y 60 centímetros (Facultad de Ciencias de la UNAM & Dirección General de Vida Silvestre, 2014). Normalmente, vuela entre los 1000 y 1500 metros, sobre parches de vegetación extensos y áreas abiertas. En el caso de las zonas extensamente cubiertas por bosque, es raro encontrarlo a menores alturas. En Costa Rica la población está constituida por individuos residentes y migratorios neárticos (provienen de Norteamérica).

Falco peregrinus (Halcón peregrino). Rapaz fuerte y de tamaño considerable (40-47 cm) que habita en ciudades, marismas y playas entre los 0 y 1200 msnm. Sigue una dieta intermedia que consiste en aves más pequeñas que captura en pleno vuelo a grandes velocidades. Presenta una distribución cosmopolita y cada vez es más frecuente avistarlos en edificios altos donde pueden anidar y capturar palomas (Santillán, Travaini & Fernández, 2010).

Falco sparverius (Cernícalo americano). Halcón pequeño (aproximadamente 30 cm) en comparación con las especies mencionadas anteriormente, con alas angostas y largas y una cola con forma cuadrada. Es un ave generalista, solitaria y sedentaria que sobrevuela áreas abiertas con vegetación corta y poca cantidad de árboles (Simonetti et al., 1982, p. 122). Es común observarlo a lo largo de carreteras en los cables telefónicos donde acostumbra perchar para después lanzarse sobre la presa.

Pandion haliaetus (Gavilán pescador). Rapaz especialista de gran tamaño que va desde 56 a 63 cm y que puede alcanzar una envergadura alar de hasta 180 cm (Muñoz et al., 2019). Sobrevuela cuerpos de agua en busca de peces. Suele anidar en superficies a gran altura ubicadas cerca del agua. En Costa Rica esta especie se puede encontrar en canales, estuarios y playas del Caribe (Garrigues & Dean, 2014, p. 92).

Figura 1. Ejemplares de las especies. a) Buteo jamaicensis; b) Falco peregrinus; c) Falco sparverius; d) Pandion haliaetus

Fotografías: a) y c) tomadas por Rio Dante Barrantes; b) y d) obtenidas de https://www.flickr.com/photos/archer10/2587959890/ y https://www.audubon.org/es/guia-de-aves/ave/aguila-pescadora, respectivamente.

Es importante mencionar que Falco peregrinus experimentó una seria afectación en su distribución y abundancia mundial entre los años 1960 y mediados de 1970, debido al uso de plaguicidas organoclorados que ocasionaban muerte embrionaria al reblandecer la cáscara de los huevos (McNutt et al., 1988, p. 247). Sumado a esto, el halcón peregrino es muy cotizado para cetrería, volviéndole altamente vulnerable a la captura de sus juveniles para tráfico ilegal (Ramos, 1986, p. 34). En la actualidad, gracias a los esfuerzos de conservación de los rapaces, se considera que la tendencia poblacional es estable (The IUCN Red List of Threatened Species, 2019). La reintroducción al medio natural de polluelos fue una de las medidas tomadas para frenar el descenso poblacional de Falco peregrinus, entre los múltiples y complejos procesos que conlleva, la alimentación es de suma importancia. Se requiere una dieta variada alta en proteínas, compuesta de animales recién sacrificados o vivos y de la que puedan obtener grasas para evitar enfermedades como la pododermatitis (Contreras & Ubilla, 2013, p. 7).

OBJETIVO GENERAL

Analizar las variables que determinan la presencia de aves rapaces en distintos sitios en Costa Rica.

OBJETIVOS ESPECÍFICOS

  1. Determinar si existen diferencias entre la abundancia de las especies según el sitio estudiado.

  2. Comparar la tasa de crecimiento de polluelos de Falco peregrinus con dos tipos de alimentación.

METODOLOGÍA

Se seleccionó como sitios de interés cuatro diferentes ambientes en territorio costarricense: bosque (BOS), humedal (HUM), pastizal (PAS) y zona urbana (ZU); en ellos se realizó la división en 10 transectos (A-J) de 1 km2 cada uno y se procedió con el muestreo de individuos de Buteo jamaicensis (BUT_JAM), Falco peregrinus (FAL_PER), Falco sparverius (FAL_SPA) y Pandion haliaetus (PAN_HAL) por un periodo de dos años. Para analizar la relación entre variables se corroboró la asimetría de los datos mediante la prueba de Shapiro-Wilk, por lo que la correlaciones entre la temperatura (°C) y la humedad (%) del ambiente, y la talla (cm) y la abundancia de las especies fueron obtenidas mediante el coeficiente de correlación de Spearman. Además, se planteó un modelo de regresión para la predicción del valor de la humedad a 30 °C y otro para la abundancia de individuos con una talla de 70 cm usando el método de Spearman.

De forma adicional, se obtuvo polluelos de Falco peregrinus en cautiverio, se registró su talla inicial (cm) y se dividieron en dos grupos. Uno de los grupos fue alimentado con granos, mientras que el otro fue alimentado con una dieta alta en proteínas. Cada grupo de alimentación se incubó a cuatro diferentes temperaturas de 10, 15, 25 y 35 °C. Finalmente, se registró la talla final de los polluelos para la comparación de la tasa de crecimiento (cm) con ambos tratamientos (tipo de alimentación y temperatura) mediante análisis de varianza (ANDEVA) de dos vías. También se realizó el análisis de varianza (dos vías) para la abundancia de las especies en cada sitio muestreado. Para ambos ANDEVA se aplicó interacción entre los factores y se comprobaron todos los supuestos estadísticos asociados (independencia, aleatoriedad, normalidad de residuos por pruebas de Shapiro-Wilk y homocedasticidad de varianza por pruebas de Fligner). No obstante, al realizar la prueba de homocedasticidad para ambos análisis de varianza se obtuvo que los datos presentan heterocedasticidad, razón por la que se transformó las variables numéricas (tasa de crecimiento y abundancia) con la aplicación de la raíz cuadrada y de la escala logarítmica, en diferentes intentos. Al persistir la heterocedasticidad, se optó por asumir la homocedasticidad de los datos. Finalmente, se realizaron pruebas de Kruskall-Wallis en la búsqueda de diferencias significativas entre tratamientos y posteriormente se demostró las diferencias existentes por grupo mediante pruebas a posteriori de Bonferroni. Todos los análisis de datos fueron hechos en el lenguaje de programación R, versión 1.3.1056 (Development Core Team, 2020).

RESULTADOS

La humedad y la temperatura presentaron una correlación estadísticamente significativa (rho<0.4; p<0.05), a pesar de ser una correlación débil y positiva. Por el contrario, la correlación realizada para la talla y la abundancia de las especies fue moderada negativa; aunque también, estadísticamente significativa (rho<0.4; p<0.05).

Figura 2. Correlación de Spearman para la temperatura y el porcentaje humedad en el ambiente

Figura 3. Correlación de Spearman para la talla y la abundancia de especies

En el caso de la predicción del valor de humedad a 30 °C se obtuvo un 84 % (F106.91; p<0.05) por medio del modelo de regresión lineal dado por la ecuación y = 55.85 + 0.94x (ecuación 1). Además, el modelo de predicción (y = 21.089334 + -0.257080x) para la correlación entre la talla y la abundancia arrojó que se puede encontrar, aproximadamente 3 individuos (F1118; p<0.05) con tallas de 70 cm. Tanto las correlaciones de Spearman como los modelos de regresión lineal fueron obtenidos y estimados con nivel de significancia del 99.9%.

Figura 4. Modelo de regresión lineal para la predicción de valores de la humedad en el ambiente a partir de la temperatura

Figura 5. Modelo de regresión lineal para la predicción de valores de la abundancia de las especies a partir de la talla

El ANDEVA realizado para la interacción de la tasa de crecimiento de los polluelos con el tipo de alimentación y la temperatura de incubación arrojó que existen diferencias significativas entre los tratamientos (KW=58.11; df=7;p<0.05), por lo que se procedió a corroborarlas por medio de un test de Bonferroni (datos balanceados y muestras mayores a 6 por tratamiento).

Figura 6. Crecimiento de los polluelos a diferentes temperaturas y con diferente tipo de alimentación
Cuadro 1. Diferencias en el promedio del crecimiento de los polluelos de acuerdo con los tratamientos aplicados

En el caso del ANDEVA para la abundancia de las especies por sitio se obtuvo diferencia significativa (KW=567.22, df= 15; p<0.05) que se confirmó por medio de la prueba de Bonferroni, comprobando que la abundancia de cada especie no presenta diferencias entre sitios, pero sí entre especies.

Figura 7. Abundancia de las cuatro especies en los diferentes sitios
Cuadro 2. Diferencias en el promedio de la abundancia de las especies de acuerdo con el sitio muestreado

De forma adicional, se formularon dos pruebas de hipótesis para la abundancia de Buteo jamaicensis. Una en busca de diferencias significativas entre los años 2018 y 2019, para la que se encontró que la abundancia no presentó mayores diferencias entre ambos años (p>0.05); y otra en la que se comparó la abundancia de la especie en las estaciones seca y lluviosa durante los dos años de muestreo, que tampoco presentó diferencia significativa (p>0.05). Ambas hipótesis de igualdad fueron aceptadas por medio de la prueba de Wilcoxon al 95% de confianza.

Figura 8. Abundancias de Buteo jamaicensis por años y estaciones (durante el 2018 y el 2019)

Además, se realizaron pruebas de hipótesis para una población (muestra de 30) para dos de las especies muestreadas, asumiendo normalidad de los datos y utilizando la prueba t-Student. En la primera la pregunta a responder fue si era posible que la abundancia de Falco peregrinus en el bosque, en estación seca del año 2018 fuera igual a 11, para la que se demostró que sí es posible con un nivel de confianza de 95% (p>0.05). La segunda pregunta consistió en demostrar que la abundancia de Falco sparverius en pastizales, en la estación lluviosa del año 2018 fue igual a 17 con un nivel de confianza de 95% (p>0.05). También, como ejercicio y con el objetivo de evaluar las posibilidades de éxito de las capturas, se obtuvo las probabilidades de encontrar un número específico de especies en diez muestreos; esto asumiendo la distribución normal de los datos y utilizando probabilidad binomial. Para dos de las cuatro especies requeridas la probabilidad fue de 28%; para al menos dos fue de 76%; y por último, para más de dos de las cuatro especies fue de un 47%. Se plantearon preguntas de carácter teórico asumiendo distribución normal de los datos. Primero, se encontró que la probabilidad de obtener una talla menor a 6 cm fue de 99%; mientras que la probabilidad de obtener una abundancia de especies menor a 5 fue de 73%. Tanto las probabilidades binomiales como las de distribución normal fueron obtenidas con un nivel de confianza de 95%.

DISCUSIÓN

En el caso de la correlación entre temperatura y humedad, las variables muestran una relación directamente proporcional que tiene sentido porque entre mayor sea la temperatura del ambiente mayor será el porcentaje de humedad, la explicación para esto radica en que entre mayor sea la temperatura ambiental, más fácil será que el agua empiece a evaporarse porque está más cerca de su punto de ebullición y cuanta más agua se evapore mayor será el porcentaje de humedad (Rodríguez, 2009, p. 190). También es importante tomar en cuenta que el aire caliente tiene una mayor capacidad de contener humedad que el aire frío, lo cual ayuda a explicar también esta relación (Rojas & Castillo, 2006, p. 68).

En la correlación entre talla y abundancia, las variables muestran una relación inversamente proporcional, lo que tiene sentido desde el punto de vista biológico debido a que las especies con menor talla suelen ser más abundantes, entre más pequeño es un individuo menos espacio y energía requiere; lo cual implica que para las especies con individuos más pequeños es más fácil tener poblaciones grandes. Se puede decir entonces que la abundancia presenta una relación negativa con respecto al tamaño del cuerpo (Granados, 2017, p. 40). Los datos obtenidos siguen este principio, ya que la especie con menor talla (18.64±6.46 cm) fue la más abundante (Falco sparverius), mientras que la de mayor talla (61.11±3.07 cm) es la segunda menos abundante (Pandion haliaetus). Es decir, los datos de este muestreo coinciden con las medias reportadas para cada especie en la literatura relacionada.

Por otra parte, el éxito de la reinserción de polluelos de aves rapaces está determinado por múltiples factores, entre ellos se destaca la alimentación y la temperatura de incubación. Como demuestran nuestros resultados, la dieta alta en proteínas favorece el crecimiento debido al hábito carnívoro de los rapaces y a que si se pretende conseguir la liberación del juvenil este debe haber entrenado sus habilidades de caza durante la crianza (Soler, 2002). Además, las aves rapaces deben consumir suficientes vísceras para obtener carotenoides y transformarlos en vitamina A (Redig, 1997, p.926). La temperatura de incubación en cautiverio ha sido un aspecto poco descrito en el material referente a conservación y reintegración de rapaces; no obstante, los rangos varían según las familias e incluso especies, así como de factores ambientales externos (Ruiz et al., 2016, p. 59). En nuestro caso, el mayor crecimiento de los polluelos se dio a los 25°C (ver figura 6), temperatura comúnmente utilizada en la incubación de aves.

En cuanto a la distribución, los resultados obtenidos indican que el promedio de abundancia de las cuatro especies de aves rapaces no mostró diferencias estadísticamente significativas por sitio (ver cuadro 2) pero si comparamos las abundancias entre especies, sí observamos diferencias. La explicación de por qué los distintos sitios no mostraron diferencias en sus valores promedio de abundancia, radica en que las aves rapaces se han adaptado muy bien a las zonas urbanas y esto implica que actualmente es posible encontrarlas, tanto en sus hábitats naturales como en las ciudades. Los rapaces se han visto forzados a adaptarse a ambientes altamente modificados debido la fragmentación y destrucción de su hábitat natural (Carrara et al., 2015, p. 118). Lo mencionado anteriormente, explica de manera adecuada las abundancias de los rapaces generalistas e intermedios, pero no del único especialista muestreado (Pandion haliaetus); ya que se esperaba que las abundancias del águila pescadora estuvieran agrupadas en los humedales principalmente, esto porque las aves especialistas se ven limitadas por sus hábitos y son las más afectadas por el desplazamiento obligatorio de sus dominios (Ferrer, 2015, p. 4).

CONCLUSIONES

La reinserción de rapaces al medio natural es lo suficientemente complicada como para tomar en cuenta solo la temperatura y la alimentación de los polluelos. Además, difícilmente se puede interpretar de forma confiable la estadística inferida sobre la tasa de crecimiento, debido a que la literatura sobre el papel de la temperatura de incubación es muy escasa. En cuanto a la distribución de las especies según el sitio, la fragmentación y destrucción de los ambientes naturales ha provocado que las aves generalistas busquen adaptarse a las zonas urbanas; mientras que las especialistas luchan por sobrevivir a los cambios generados por el ser humano. Por último, se recomienda investigaciones futuras acerca del impacto que tiene la urbanización descontrolada de los sitios naturales en los rapaces especialistas.

REFERENCIAS BIBLIOGRÁFICAS

Carrara, E., Arroyo, V., Vega, J., Schondube, J., de Freitas, S., & Fahrig, L. (2015). Impact of landscape composition and configuration on forest specialist and generalist bird species in the fragmented Lacandona rainforest, Mexico. Biological Conservation, 184, 118. https://doi.org/10.1016/j.biocon.2015.01.014

Contreras, P., & Ubilla, M. (2013). Evaluación del bienestar animal de aves rapaces en rehabilitación, descripción de técnicas que lo promuevan y mejoren su tasa de reintroducción. Avances en Ciencias Veterinarias, 28(2), 7. https://doi.org/10.5354/0719-5273.2014.30201

Development Core Team. (2020). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. Facultad de Ciencias de la UNAM & Dirección General de Vida Silvestre. (2014). Ficha de identificación Buteo jamaicensis. MGBE.

Garrigues, R., & Dean, R. (2014). Pandionidae. En The Birds of Costa Rica (II ed., p. 92). Amsterdam University Press.

Granados, C. L. (2017). Avances en ecología: Hacia un mejor conocimiento de la naturaleza (1.a ed). Universidad de Sevilla-Secretariado de Publicaciones.

McNutt, J., Ellis, D., Roundy, T., Vasina, W., & White, C. (1988). Distribution and status of the peregrine falcon in South America. The Peregrine Fund Inc, 237-249. http://globalraptors.org/grin/researchers/uploads/476/1988_pf_sa.pdf

Muñoz, A., Rau, J., Yáñez, J., Norambuena, H., & González, D. (2019). Aves rapaces de Chile (Segunda ed.). CEA Ediciones.

Ramos, M. (1986). Birds in Peril in Mexico: The Diurnal Raptors. Birds of Prey, 3, 34-36. http://www.raptors-international.org/book/birds_of_prey_1986/Ramos_1986_26-42.pdf

Redig, P. (1997). Raptors. Avian Medicine And Surgery, Primera, 918-928. http://globalraptors.org/grin/ResearcherResults.asp?lresID=266

Rodríguez, E. G. (2009), Influencia de la temperatura y de la humedad en la dinámica de la materia orgánica de los suelos de Galicia y su relación con el cambio climático. Universidad de Santiago de Compostela (España).

Rojas, D. C. & Castillo, J. M. (2006). La química al alcance de todos (1.a ed) Universidad Nacional de Colombia.

Rojas, G., & Iza, A. (2009). Derecho ambiental en Centroamérica: Vol. Tomo I. Diseño Editorial S.A.

Ruiz, N., Orrego, G., Reyes, M., & Silva, M. (2016). Aumento de la Temperatura de Incubación en Huevos de Gallina Araucana (Gallus inauris): Efecto sobre la Mortalidad Embrionaria, Tasa de Eclosión, Peso del Polluelo, Saco Vitelino y de Órganos Internos. International Journal of Morphology, 34(1), 57-62. https://doi.org/10.4067/S0717-95022016000100009

Santillán, M., Travaini A., & Fernández, J. (2010, 16 marzo). Dieta del halcón peregrino (Falco peregrinus) en la Ría Deseado, Patagonia austral Argentina. Cecara. http://www.cecara.com.ar/web/wp-content/uploads/2010-Santillan-et-al.-b.pdf

Simonetti, J., Nuñez, H., & Yañez, J. (1982b). Falco sparverius L.: Rapaz generalista en Chile central (aves: Falconidae). Boletín del museo nacional de historia natural, 119-124.

Soler, D. (2002). Medicina de aves rapaces diurnas (orden Falconiformes). Grupo de investigación en fauna silvestre-Asociación de veterinarios de vida silvestre. http://www.veterinariosvs.org/miembros/diegosoler42.Stringfield

Vides-Hernández, G., Velado-Cano, M., Pablo-Cea, J., & Carmona-Galindo, V. (2017). Patrones de riqueza y diversidad de aves en áreas verdes del centro urbano de San Salvador, El Salvador. Revista Mexicana de Ornitología, 18(2), 272-280. https://doi.org/10.28947/hrmo.2017.18.2.294

Y. Ferrer. (2015, julio). Variables que influyen en la distribución y abundancia de rapaces diurnas y en la ubicación de sus sitios de anidación en Cuba (TFG). Centro de Investigaciones Biológicas del Noroeste, S.C.

Anexos

 # # PROYECTO FINAL # #

# # ---------------------------------------------------------------------- # #

library(readxl)
BD_AvesRapaces <- read_excel("Data/BD_AvesRapaces.xlsx", 
                             sheet = "Muestreo", col_types = c("text", 
                                                               "text", "text", "text", "numeric", 
                                                               "text", "numeric", "numeric", "numeric", 
                                                               "numeric", "numeric"))

EXP <- read_excel("Data/EXP.xlsx", col_types = c("text", 
                                                 "text", "numeric"))

# # ---------------------------------------------------------------------- # #

                            # Prob. Binomial # 

# 1. Cual es la probabilidad de que en 10 muestreos encuentre exactamente 2 
# de las 4 especies de mi base de datos?
dbinom(2,10,0.25)                           # 0.2815676
## [1] 0.2815676
# 2. Cual es la probabilidad de que en 10 muestreos encuentre al menos 2 de 
# las 4 especies de mi base de datos?
1-pbinom(1,10,0.25)                         # 0.7559748
## [1] 0.7559748
# 3. Cual es la probabilidad de que en 10 muestreos encuentre más de 2 de 
# las 4 especies de mi base de datos?
1-pbinom(2,10,0.25)                         # 0.4744072
## [1] 0.4744072
# # ---------------------------------------------------------------------- # #
                            
                               # Dist. Normal # 

# SON PREGUNTAS TEORICAS, SE ASUME VARIANZA CONOCIDA. 

mean(BD_AvesRapaces$Talla)       # 45.14948 
## [1] 45.14948
sd(BD_AvesRapaces$Talla)          # 17.60134
## [1] 17.60134
mean(BD_AvesRapaces$Abundancia)    # 9.482292
## [1] 9.482292
sd(BD_AvesRapaces$Abundancia)       # 7.455737
## [1] 7.455737
library(TeachingDemos)

# 1. Cual es la probabilidad de que obtenga un valor de talla menor a 6 si 
# el promedio del valor de talla es 45.14948 y su desviacion estandar es de 
# 17.60134?
z.test (6, 45.14948, 17.60134, alternative = "g")  # p-value = 0.9869
## 
##  One Sample z-test
## 
## data:  6
## z = -2.2242, n = 1.000, Std. Dev. = 17.601, Std. Dev. of the sample
## mean = 17.601, p-value = 0.9869
## alternative hypothesis: true mean is greater than 45.14948
## 95 percent confidence interval:
##  -22.95163       Inf
## sample estimates:
## mean of 6 
##         6
# 2. Cual es la probabilidad de que obtenga un valor de abundancia menor a 5 
# si el promedio del valor de abundancia es 9.482292 y su desviacion estandar
# es de 7.455737?
z.test(5, 9.482292, 7.455737, alternative = "g")  # p-value = 0.7261
## 
##  One Sample z-test
## 
## data:  5
## z = -0.60119, n = 1.0000, Std. Dev. = 7.4557, Std. Dev. of the sample
## mean = 7.4557, p-value = 0.7261
## alternative hypothesis: true mean is greater than 9.482292
## 95 percent confidence interval:
##  -7.263596       Inf
## sample estimates:
## mean of 5 
##         5
# # ---------------------------------------------------------------------- # #

                          # PRUEBAS DE HIPOTESIS #

# Pruebas para una poblacion:
#1. Es posible concluir que la abundancia de la especie Falco peregrinus de la 
# estación seca, año 2018, sitio BOS y tamaño de muestra de 30 es igual a 
# 11 con un nivel de confianza del 95%?

Ab1<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Anio == "2018" & BD_AvesRapaces$Estacion == "SECA" & BD_AvesRapaces$Especies == "FAL_PER" & BD_AvesRapaces$Sitio=="BOS"]
mean(Ab1)  # 9.933333
## [1] 9.933333
sd(Ab1)  # 3.81407
## [1] 3.81407
shapiro.test(Ab1)  # W = 0.95419, p-value = 0.2186
## 
##  Shapiro-Wilk normality test
## 
## data:  Ab1
## W = 0.95419, p-value = 0.2186
library(TeachingDemos)
z.test(9.933333, 11, 3.81407) # p-value = 0.7797
## 
##  One Sample z-test
## 
## data:  9.933333
## z = -0.27967, n = 1.0000, Std. Dev. = 3.8141, Std. Dev. of the sample
## mean = 3.8141, p-value = 0.7797
## alternative hypothesis: true mean is not equal to 11
## 95 percent confidence interval:
##   2.457893 17.408773
## sample estimates:
## mean of 9.933333 
##         9.933333
# Se concluye que la abundancia de la especie Falco peregrinus de la estación 
# seca, año 2018, sitio BOS y tamaño de muestra de 30 puede ser igual a 11 
# con un nivel de confianza del 95%.


#2. Es posible concluir que la abundancia de Falco sparverius en la estacion 
# lluviosa del ano 2018 en pastizales es igual a 17 con un nivel de confianza
# del 95%?
Ab2<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Anio == "2018" & BD_AvesRapaces$Estacion == "LLUVIOSA" & BD_AvesRapaces$Especies == "FAL_SPA" & BD_AvesRapaces$Sitio=="PAS"]
mean(Ab2)   # 16.36667
## [1] 16.36667
sd(Ab2)   #  8.99227
## [1] 8.99227
shapiro.test(Ab2)  #  p-value = 0.02494
## 
##  Shapiro-Wilk normality test
## 
## data:  Ab2
## W = 0.91878, p-value = 0.02494
(16.36667-17)/(8.99227/sqrt(30))  # Zc: -0.3857637
## [1] -0.3857637
# Zc > Zt , ACEPTO NULA. 
# Se concluye que la abundancia de Falco sparverius en la estacion 
# lluviosa del ano 2018 en pastizales puede ser igual a 17 con un nivel de 
#confianza del 95%


# Prueba para dos poblaciones: 
# a. Anos:
#1. Existe una diferencia significativa entre las abundancias de la especie 
# Buteo jamaicensis en los anos 2018 y 2019 con un nivel de confianza del 95 %?
ab18<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Anio == "2018" & BD_AvesRapaces$Especies == "BUT_JAM"]
ab19<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Anio == "2019" & BD_AvesRapaces$Especies == "BUT_JAM"]
shapiro.test(ab18)  # W = 0.93147, p-value = 4.014e-09
## 
##  Shapiro-Wilk normality test
## 
## data:  ab18
## W = 0.93147, p-value = 4.014e-09
shapiro.test(ab19)  # W = 0.92422, p-value = 9.778e-10
## 
##  Shapiro-Wilk normality test
## 
## data:  ab19
## W = 0.92422, p-value = 9.778e-10
wilcox.test(ab18, ab19, paired = FALSE, conf.int = TRUE, exact = FALSE, conf.level = 0.95)
## 
##  Wilcoxon rank sum test with continuity correction
## 
## data:  ab18 and ab19
## W = 30196, p-value = 0.3558
## alternative hypothesis: true location shift is not equal to 0
## 95 percent confidence interval:
##  -3.775172e-05  9.999901e-01
## sample estimates:
## difference in location 
##           1.718474e-05
# p-value = 0.3558
# No existen diferencias 
boxplot(ab18, ab19, names = c("2018", "2019"), col = c("#F9EBEA", "#E8F8F5"), ylab = "Abundancia", xlab = "Año")


# b. Estaciones:
#2. Existe una diferencia significativa entre las tallas de la especie 
# Falco sparverius de la estación seca y del sitio BOS con las tallas de la 
# misma especie de la estación lluviosa y del mismo sitio con un nivel de 
# confianza del 95 % ?
t1<-BD_AvesRapaces$Talla[BD_AvesRapaces$Especies == "FAL_SPA" & BD_AvesRapaces$Estacion=="SECA" & BD_AvesRapaces$Sitio=="BOS"]
t2<-BD_AvesRapaces$Talla[BD_AvesRapaces$Especies == "FAL_SPA" & BD_AvesRapaces$Estacion=="LLUVIOSA" & BD_AvesRapaces$Sitio=="BOS"]
shapiro.test(t1)  # W = 0.9494, p-value = 0.01459
## 
##  Shapiro-Wilk normality test
## 
## data:  t1
## W = 0.9494, p-value = 0.01459
shapiro.test(t2)  # W = 0.93943, p-value = 0.005105
## 
##  Shapiro-Wilk normality test
## 
## data:  t2
## W = 0.93943, p-value = 0.005105
wilcox.test(t1,t2, paired = FALSE, conf.int = TRUE, exact = FALSE, conf.level = 0.95)
## 
##  Wilcoxon rank sum test with continuity correction
## 
## data:  t1 and t2
## W = 2004, p-value = 0.2848
## alternative hypothesis: true location shift is not equal to 0
## 95 percent confidence interval:
##  -1.000035  3.999980
## sample estimates:
## difference in location 
##               1.000037
# Se concluye que no existe una diferencia significativa entre las tallas de 
# la especie Falco sparverius  de la estación seca y del sitio BOS con las 
# tallas de la misma especie de la estación lluviosa

par(mfrow=c(1,2))
boxplot(ab18, ab19, names = c("2018", "2019"), col = c("#F9EBEA", "#E8F8F5"), ylab = "Abundancia", xlab = "Año")
boxplot(t1, t2, names = c("Seca", "Lluviosa"), col = c("#F7DC6F", "#85C1E9"), ylab = "Abundancia", xlab = "Estaciones")

# # ---------------------------------------------------------------------- # #

                         # CORRELACION Y RELACION #

TEMPHUM<-cor.test(BD_AvesRapaces$Temperatura, BD_AvesRapaces$Humedad, method = "spearman", conf.level = 0.999)
TEMPHUM   
## 
##  Spearman's rank correlation rho
## 
## data:  BD_AvesRapaces$Temperatura and BD_AvesRapaces$Humedad
## S = 906868130, p-value < 2.2e-16
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## 0.2312382
d1<-data.frame(BD_AvesRapaces$Temperatura, BD_AvesRapaces$Humedad)
names (d1) = c("Temperatura (°C)", "Humedad (%)")
names (d1)
## [1] "Temperatura (°C)" "Humedad (%)"
library(PerformanceAnalytics)
## Warning: package 'PerformanceAnalytics' was built under R version 4.0.3
## Loading required package: xts
## Warning: package 'xts' was built under R version 4.0.3
## Loading required package: zoo
## Warning: package 'zoo' was built under R version 4.0.3
## 
## Attaching package: 'zoo'
## The following objects are masked from 'package:base':
## 
##     as.Date, as.Date.numeric
## 
## Attaching package: 'PerformanceAnalytics'
## The following object is masked from 'package:graphics':
## 
##     legend

chart.Correlation(d1, method = "spearman")

TALABUN<-cor.test(BD_AvesRapaces$Talla, BD_AvesRapaces$Abundancia, method = "spearman", exact = F, conf.level = 0.999)
TALABUN 
## 
##  Spearman's rank correlation rho
## 
## data:  BD_AvesRapaces$Talla and BD_AvesRapaces$Abundancia
## S = 1.733e+09, p-value < 2.2e-16
## alternative hypothesis: true rho is not equal to 0
## sample estimates:
##       rho 
## -0.469111
library(PerformanceAnalytics)
d2<-data.frame(BD_AvesRapaces$Talla, BD_AvesRapaces$Abundancia) 
names (d2) = c("Talla (cm)", "Abundancia")
names (d2)
## [1] "Talla (cm)" "Abundancia"
chart.Correlation(d2, method = "spearman")
## Warning in cor.test.default(as.numeric(x), as.numeric(y), method = method):
## Cannot compute exact p-value with ties

# # ---------------------------------------------------------------------- # #

                          # MODELOS DE REGRESION #

ml1<-lm(BD_AvesRapaces$Humedad~BD_AvesRapaces$Temperatura, method = "spearman")
## Warning in lm(BD_AvesRapaces$Humedad ~ BD_AvesRapaces$Temperatura, method =
## "spearman"): method = 'spearman' is not supported. Using 'qr'
summary(ml1) 
## 
## Call:
## lm(formula = BD_AvesRapaces$Humedad ~ BD_AvesRapaces$Temperatura, 
##     method = "spearman")
## 
## Residuals:
##     Min      1Q  Median      3Q     Max 
## -34.466  -8.514   3.779   9.627  18.883 
## 
## Coefficients:
##                            Estimate Std. Error t value Pr(>|t|)    
## (Intercept)                 55.8463     2.2739   24.56   <2e-16 ***
## BD_AvesRapaces$Temperatura   0.9421     0.0911   10.34   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 11.82 on 1918 degrees of freedom
## Multiple R-squared:  0.05281,    Adjusted R-squared:  0.05231 
## F-statistic: 106.9 on 1 and 1918 DF,  p-value: < 2.2e-16
d3<-data.frame(BD_AvesRapaces$Temperatura, BD_AvesRapaces$Humedad)
names (d3) = c("Temperatura", "Humedad")
library(ggplot2)
ggplot(d3, aes(x=Temperatura, y=Humedad)) + geom_point(shape=1) + geom_smooth(method=lm) + labs(x="Temperatura (°C)",y="Humedad (%)")  
## `geom_smooth()` using formula 'y ~ x'

ml2<-lm(BD_AvesRapaces$Abundancia~BD_AvesRapaces$Talla, method = "spearman")
## Warning in lm(BD_AvesRapaces$Abundancia ~ BD_AvesRapaces$Talla, method =
## "spearman"): method = 'spearman' is not supported. Using 'qr'
summary(ml2)  
## 
## Call:
## lm(formula = BD_AvesRapaces$Abundancia ~ BD_AvesRapaces$Talla, 
##     method = "spearman")
## 
## Residuals:
##      Min       1Q   Median       3Q      Max 
## -17.5185  -3.7212  -0.3061   3.5713  19.1089 
## 
## Coefficients:
##                       Estimate Std. Error t value Pr(>|t|)    
## (Intercept)          21.089334   0.372497   56.62   <2e-16 ***
## BD_AvesRapaces$Talla -0.257080   0.007687  -33.44   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 5.927 on 1918 degrees of freedom
## Multiple R-squared:  0.3683, Adjusted R-squared:  0.368 
## F-statistic:  1118 on 1 and 1918 DF,  p-value: < 2.2e-16
d4<-data.frame(BD_AvesRapaces$Talla, BD_AvesRapaces$Abundancia) 
names (d4) = c("Talla", "Abundancia")
library(ggplot2)
ggplot(d4, aes(x=Talla, y=BD_AvesRapaces$Abundancia)) + geom_point(shape=1) + geom_smooth(method=lm) + labs(x="Talla (cm)",y="Abundancia")  
## `geom_smooth()` using formula 'y ~ x'

# # ---------------------------------------------------------------------- # #

                 # ANALISIS DE VARIANZA - HOJA EXPERIMENTO #

str(EXP)
## tibble [80 x 3] (S3: tbl_df/tbl/data.frame)
##  $ Alimentacion: chr [1:80] "GRANOS" "GRANOS" "GRANOS" "GRANOS" ...
##  $ Temperatura : chr [1:80] "10" "10" "10" "10" ...
##  $ Crecimiento : num [1:80] 0.447 0.18 0.628 1.01 0.13 ...
length(EXP$Alimentacion)==length(EXP$Crecimiento)
## [1] TRUE
EXP$Alimentacion<-as.factor(EXP$Alimentacion)
is.factor(EXP$Alimentacion)
## [1] TRUE
EXP$Temperatura<-as.factor(EXP$Temperatura)
is.factor(EXP$Temperatura)
## [1] TRUE
factores<-paste(EXP$Alimentacion, EXP$Temperatura, sep = "*")
is.factor(factores)
## [1] FALSE
prueba1<-aov(EXP$Crecimiento ~ factores)
summary(prueba1)
##             Df Sum Sq Mean Sq F value Pr(>F)    
## factores     7 189.91  27.129   36.22 <2e-16 ***
## Residuals   72  53.94   0.749                   
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
shapiro.test(prueba1$residuals)  #  # W = 0.98598, p-value = 0.5333
## 
##  Shapiro-Wilk normality test
## 
## data:  prueba1$residuals
## W = 0.98598, p-value = 0.5333
fligner.test(EXP$Crecimiento ~ factores)
## 
##  Fligner-Killeen test of homogeneity of variances
## 
## data:  EXP$Crecimiento by factores
## Fligner-Killeen:med chi-squared = 36.821, df = 7, p-value = 5.07e-06
y1<-sqrt(EXP$Crecimiento)
# Con log no porque hay 0 cm en algunos casos 
fligner.test(y1 ~ factores) 
## 
##  Fligner-Killeen test of homogeneity of variances
## 
## data:  y1 by factores
## Fligner-Killeen:med chi-squared = 27.89, df = 7, p-value = 0.0002302
# Fligner-Killeen:med chi-squared = 27.89, df = 7,p-value = 0.0002302

# Verificando con datos crudos y con transformacion, no homocedast
# ASUMIR HOMOCEDASTICIDAD DE LOS RESIDUOS
kruskal.test(EXP$Crecimiento ~ factores)  
## 
##  Kruskal-Wallis rank sum test
## 
## data:  EXP$Crecimiento by factores
## Kruskal-Wallis chi-squared = 58.112, df = 7, p-value = 3.592e-10
# Kruskal-Wallis chi-squared = 58.112, df = 7, p-value = 3.592e-10
# Hay diferencia signif entre los tratamientos
pairwise.wilcox.test(EXP$Crecimiento,factores,p.adj = "b", exact=F)
## 
##  Pairwise comparisons using Wilcoxon rank sum test with continuity correction 
## 
## data:  EXP$Crecimiento and factores 
## 
##              GRANOS*10 GRANOS*15 GRANOS*25 GRANOS*35 PROTEINAS*10 PROTEINAS*15
## GRANOS*15    1.0000    -         -         -         -            -           
## GRANOS*25    0.0092    0.0163    -         -         -            -           
## GRANOS*35    1.0000    1.0000    0.0163    -         -            -           
## PROTEINAS*10 1.0000    1.0000    1.0000    1.0000    -            -           
## PROTEINAS*15 0.0051    0.0051    0.3925    0.0051    1.0000       -           
## PROTEINAS*25 0.0051    0.0051    0.0051    0.0051    0.0163       0.1623      
## PROTEINAS*35 0.0051    0.0051    0.0051    0.0051    0.0617       0.8739      
##              PROTEINAS*25
## GRANOS*15    -           
## GRANOS*25    -           
## GRANOS*35    -           
## PROTEINAS*10 -           
## PROTEINAS*15 -           
## PROTEINAS*25 -           
## PROTEINAS*35 1.0000      
## 
## P value adjustment method: bonferroni
Crecimiento<-EXP$Crecimiento
Temperatura<-as.factor(EXP$Temperatura)
is.factor(Temperatura)
## [1] TRUE
Alimentación<-as.factor(EXP$Alimentacion)
is.factor(Alimentación)
## [1] TRUE
dafra1<-data.frame(Alimentación, Temperatura, Crecimiento)
ggplot(data=dafra1, aes(x=reorder(Alimentación,Crecimiento), y=Crecimiento, fill=Temperatura)) + geom_bar(stat="identity", position="dodge") + labs(x="Tipo de alimentación",y="Crecimiento (cm)") + scale_fill_manual("Temperatura (°C)", values = c("10" = "#FFCC80", "15" = "#FF9800", "25" = "#FF6D00", "35"="#E65100"))

# SE TOMA LA DECISION DE USAR OTRO GRAF XQ ESCALA NO COINCIDE

                 # ANALISIS DE VARIANZA - HOJA MUESTREO #
str(BD_AvesRapaces)
## tibble [1,920 x 11] (S3: tbl_df/tbl/data.frame)
##  $ Anio             : chr [1:1920] "2018" "2018" "2018" "2018" ...
##  $ Estacion         : chr [1:1920] "SECA" "SECA" "SECA" "SECA" ...
##  $ Sitio            : chr [1:1920] "BOS" "BOS" "BOS" "BOS" ...
##  $ Transecto        : chr [1:1920] "A" "A" "A" "B" ...
##  $ Area(KM2)        : num [1:1920] 1 1 1 1 1 1 1 1 1 1 ...
##  $ Especies         : chr [1:1920] "FAL_SPA" "FAL_SPA" "FAL_SPA" "FAL_SPA" ...
##  $ Talla            : num [1:1920] 10 13 25 23 17 13 9 24 27 22 ...
##  $ Abundancia       : num [1:1920] 31 14 27 23 25 20 31 6 33 8 ...
##  $ Temperatura      : num [1:1920] 18.6 17.6 18.5 18 18.2 ...
##  $ Humedad          : num [1:1920] 76.8 79.1 79.7 83.3 81 ...
##  $ Densidad(org/km2): num [1:1920] 31 14 27 23 25 20 31 6 33 8 ...
BD_AvesRapaces$Anio<-as.factor(BD_AvesRapaces$Anio)
BD_AvesRapaces$Sitio<-as.factor(BD_AvesRapaces$Sitio)
is.factor(BD_AvesRapaces$Sitio)
## [1] TRUE
BD_AvesRapaces$Especies<-as.factor(BD_AvesRapaces$Especies) 
is.factor(BD_AvesRapaces$Especies)
## [1] TRUE
length(BD_AvesRapaces$Sitio)==length(BD_AvesRapaces$Abundancia)
## [1] TRUE
length(BD_AvesRapaces$Sitio)==length(BD_AvesRapaces$Especies)
## [1] TRUE
fact<-paste(BD_AvesRapaces$Especies, BD_AvesRapaces$Sitio, sep = "*")
prueba2<-aov(BD_AvesRapaces$Abundancia ~ fact)
shapiro.test(prueba2$residuals)  # W = 0.98501, p-value = 2.694e-13
## 
##  Shapiro-Wilk normality test
## 
## data:  prueba2$residuals
## W = 0.98501, p-value = 2.694e-13
fligner.test(BD_AvesRapaces$Abundancia ~ fact) # p-value < 2.2e-16
## 
##  Fligner-Killeen test of homogeneity of variances
## 
## data:  BD_AvesRapaces$Abundancia by fact
## Fligner-Killeen:med chi-squared = 642.54, df = 15, p-value < 2.2e-16
y2<-sqrt(BD_AvesRapaces$Abundancia)
fligner.test(y2 ~ fact)  
## 
##  Fligner-Killeen test of homogeneity of variances
## 
## data:  y2 by fact
## Fligner-Killeen:med chi-squared = 232.77, df = 15, p-value < 2.2e-16
# Fligner-Killeen:med chi-squared = 232.77, df = 15, p-value < 2.2e-16


# SE ASUME HOMOCEDASTICIDAD DE LOS RESIDUOS #
kruskal.test(BD_AvesRapaces$Abundancia ~ fact)  
## 
##  Kruskal-Wallis rank sum test
## 
## data:  BD_AvesRapaces$Abundancia by fact
## Kruskal-Wallis chi-squared = 567.22, df = 15, p-value < 2.2e-16
# Kruskal-Wallis chi-squared = 567.22, df = 15, p-value < 2.2e-16
pairwise.wilcox.test(BD_AvesRapaces$Abundancia,fact,p.adj = "b", exact=F)
## 
##  Pairwise comparisons using Wilcoxon rank sum test with continuity correction 
## 
## data:  BD_AvesRapaces$Abundancia and fact 
## 
##             BUT_JAM*BOS BUT_JAM*HUM BUT_JAM*PAS BUT_JAM*ZU FAL_PER*BOS
## BUT_JAM*HUM 1.00000     -           -           -          -          
## BUT_JAM*PAS 1.00000     1.00000     -           -          -          
## BUT_JAM*ZU  1.00000     1.00000     1.00000     -          -          
## FAL_PER*BOS 3.3e-05     9.4e-07     1.8e-06     6.3e-08    -          
## FAL_PER*HUM 0.00016     7.2e-06     1.2e-05     8.9e-07    1.00000    
## FAL_PER*PAS 5.2e-05     2.3e-06     4.0e-06     2.1e-07    1.00000    
## FAL_PER*ZU  0.00074     7.9e-05     0.00011     1.2e-05    1.00000    
## FAL_SPA*BOS < 2e-16     < 2e-16     < 2e-16     < 2e-16    2.1e-11    
## FAL_SPA*HUM < 2e-16     < 2e-16     < 2e-16     < 2e-16    1.6e-14    
## FAL_SPA*PAS < 2e-16     < 2e-16     < 2e-16     < 2e-16    1.3e-11    
## FAL_SPA*ZU  2.1e-15     2.4e-16     3.5e-16     < 2e-16    1.6e-08    
## PAN_HAL*BOS 1.00000     1.00000     1.00000     1.00000    0.01062    
## PAN_HAL*HUM 0.05581     0.00246     0.00536     0.00028    1.00000    
## PAN_HAL*PAS 1.00000     1.00000     1.00000     0.56582    0.01290    
## PAN_HAL*ZU  0.19910     0.01273     0.02383     0.00173    0.60981    
##             FAL_PER*HUM FAL_PER*PAS FAL_PER*ZU FAL_SPA*BOS FAL_SPA*HUM
## BUT_JAM*HUM -           -           -          -           -          
## BUT_JAM*PAS -           -           -          -           -          
## BUT_JAM*ZU  -           -           -          -           -          
## FAL_PER*BOS -           -           -          -           -          
## FAL_PER*HUM -           -           -          -           -          
## FAL_PER*PAS 1.00000     -           -          -           -          
## FAL_PER*ZU  1.00000     1.00000     -          -           -          
## FAL_SPA*BOS 3.4e-11     9.9e-12     3.1e-11    -           -          
## FAL_SPA*HUM 4.5e-14     1.0e-14     6.1e-14    1.00000     -          
## FAL_SPA*PAS 2.9e-11     5.3e-12     2.1e-11    1.00000     1.00000    
## FAL_SPA*ZU  1.4e-08     7.9e-09     8.4e-09    1.00000     1.00000    
## PAN_HAL*BOS 0.01860     0.02714     0.07652    < 2e-16     < 2e-16    
## PAN_HAL*HUM 1.00000     1.00000     1.00000    9.9e-16     < 2e-16    
## PAN_HAL*PAS 0.02230     0.03443     0.09586    < 2e-16     < 2e-16    
## PAN_HAL*ZU  0.74586     1.00000     1.00000    3.5e-16     < 2e-16    
##             FAL_SPA*PAS FAL_SPA*ZU PAN_HAL*BOS PAN_HAL*HUM PAN_HAL*PAS
## BUT_JAM*HUM -           -          -           -           -          
## BUT_JAM*PAS -           -          -           -           -          
## BUT_JAM*ZU  -           -          -           -           -          
## FAL_PER*BOS -           -          -           -           -          
## FAL_PER*HUM -           -          -           -           -          
## FAL_PER*PAS -           -          -           -           -          
## FAL_PER*ZU  -           -          -           -           -          
## FAL_SPA*BOS -           -          -           -           -          
## FAL_SPA*HUM -           -          -           -           -          
## FAL_SPA*PAS -           -          -           -           -          
## FAL_SPA*ZU  1.00000     -          -           -           -          
## PAN_HAL*BOS < 2e-16     6.7e-14    -           -           -          
## PAN_HAL*HUM < 2e-16     9.4e-12    1.00000     -           -          
## PAN_HAL*PAS < 2e-16     1.0e-13    1.00000     1.00000     -          
## PAN_HAL*ZU  < 2e-16     4.3e-12    1.00000     1.00000     1.00000    
## 
## P value adjustment method: bonferroni
Abundancia<-BD_AvesRapaces$Abundancia
Especies<-BD_AvesRapaces$Especies
Sitios<-BD_AvesRapaces$Sitio
dataf<-data.frame(Especies, Sitios, Abundancia)
dataf
##      Especies Sitios Abundancia
## 1     FAL_SPA    BOS         31
## 2     FAL_SPA    BOS         14
## 3     FAL_SPA    BOS         27
## 4     FAL_SPA    BOS         23
## 5     FAL_SPA    BOS         25
## 6     FAL_SPA    BOS         20
## 7     FAL_SPA    BOS         31
## 8     FAL_SPA    BOS          6
## 9     FAL_SPA    BOS         33
## 10    FAL_SPA    BOS          8
## 11    FAL_SPA    BOS         25
## 12    FAL_SPA    BOS         11
## 13    FAL_SPA    BOS         27
## 14    FAL_SPA    BOS         33
## 15    FAL_SPA    BOS         12
## 16    FAL_SPA    BOS         13
## 17    FAL_SPA    BOS         11
## 18    FAL_SPA    BOS         32
## 19    FAL_SPA    BOS         11
## 20    FAL_SPA    BOS         17
## 21    FAL_SPA    BOS          7
## 22    FAL_SPA    BOS          6
## 23    FAL_SPA    BOS         29
## 24    FAL_SPA    BOS         15
## 25    FAL_SPA    BOS         13
## 26    FAL_SPA    BOS          6
## 27    FAL_SPA    BOS         23
## 28    FAL_SPA    BOS         31
## 29    FAL_SPA    BOS         21
## 30    FAL_SPA    BOS         21
## 31    FAL_PER    BOS         11
## 32    FAL_PER    BOS          7
## 33    FAL_PER    BOS          7
## 34    FAL_PER    BOS          8
## 35    FAL_PER    BOS         11
## 36    FAL_PER    BOS         16
## 37    FAL_PER    BOS          8
## 38    FAL_PER    BOS         10
## 39    FAL_PER    BOS          8
## 40    FAL_PER    BOS         14
## 41    FAL_PER    BOS          3
## 42    FAL_PER    BOS          6
## 43    FAL_PER    BOS          6
## 44    FAL_PER    BOS          9
## 45    FAL_PER    BOS          4
## 46    FAL_PER    BOS          9
## 47    FAL_PER    BOS         16
## 48    FAL_PER    BOS         13
## 49    FAL_PER    BOS         14
## 50    FAL_PER    BOS         13
## 51    FAL_PER    BOS          5
## 52    FAL_PER    BOS          5
## 53    FAL_PER    BOS         14
## 54    FAL_PER    BOS         16
## 55    FAL_PER    BOS          9
## 56    FAL_PER    BOS          8
## 57    FAL_PER    BOS          9
## 58    FAL_PER    BOS         11
## 59    FAL_PER    BOS         16
## 60    FAL_PER    BOS         12
## 61    PAN_HAL    BOS          1
## 62    PAN_HAL    BOS          4
## 63    PAN_HAL    BOS          2
## 64    PAN_HAL    BOS          1
## 65    PAN_HAL    BOS          4
## 66    PAN_HAL    BOS         10
## 67    PAN_HAL    BOS          5
## 68    PAN_HAL    BOS          8
## 69    PAN_HAL    BOS         10
## 70    PAN_HAL    BOS          2
## 71    PAN_HAL    BOS          3
## 72    PAN_HAL    BOS          4
## 73    PAN_HAL    BOS          5
## 74    PAN_HAL    BOS         11
## 75    PAN_HAL    BOS          1
## 76    PAN_HAL    BOS          3
## 77    PAN_HAL    BOS          8
## 78    PAN_HAL    BOS          5
## 79    PAN_HAL    BOS         12
## 80    PAN_HAL    BOS          4
## 81    PAN_HAL    BOS         12
## 82    PAN_HAL    BOS         12
## 83    PAN_HAL    BOS          5
## 84    PAN_HAL    BOS          6
## 85    PAN_HAL    BOS         12
## 86    PAN_HAL    BOS          8
## 87    PAN_HAL    BOS          3
## 88    PAN_HAL    BOS          6
## 89    PAN_HAL    BOS          9
## 90    PAN_HAL    BOS          1
## 91    BUT_JAM    BOS          8
## 92    BUT_JAM    BOS          5
## 93    BUT_JAM    BOS          5
## 94    BUT_JAM    BOS          9
## 95    BUT_JAM    BOS          2
## 96    BUT_JAM    BOS          9
## 97    BUT_JAM    BOS          8
## 98    BUT_JAM    BOS          5
## 99    BUT_JAM    BOS          9
## 100   BUT_JAM    BOS          1
## 101   BUT_JAM    BOS          3
## 102   BUT_JAM    BOS          6
## 103   BUT_JAM    BOS          7
## 104   BUT_JAM    BOS         10
## 105   BUT_JAM    BOS          3
## 106   BUT_JAM    BOS         10
## 107   BUT_JAM    BOS          1
## 108   BUT_JAM    BOS         10
## 109   BUT_JAM    BOS          1
## 110   BUT_JAM    BOS         10
## 111   BUT_JAM    BOS          7
## 112   BUT_JAM    BOS          6
## 113   BUT_JAM    BOS          8
## 114   BUT_JAM    BOS         10
## 115   BUT_JAM    BOS          5
## 116   BUT_JAM    BOS          3
## 117   BUT_JAM    BOS          6
## 118   BUT_JAM    BOS          7
## 119   BUT_JAM    BOS          6
## 120   BUT_JAM    BOS          4
## 121   FAL_SPA    BOS          8
## 122   FAL_SPA    BOS          2
## 123   FAL_SPA    BOS          8
## 124   FAL_SPA    BOS         32
## 125   FAL_SPA    BOS          6
## 126   FAL_SPA    BOS         27
## 127   FAL_SPA    BOS         18
## 128   FAL_SPA    BOS          7
## 129   FAL_SPA    BOS          3
## 130   FAL_SPA    BOS         23
## 131   FAL_SPA    BOS         33
## 132   FAL_SPA    BOS          6
## 133   FAL_SPA    BOS         12
## 134   FAL_SPA    BOS          8
## 135   FAL_SPA    BOS         18
## 136   FAL_SPA    BOS         27
## 137   FAL_SPA    BOS          6
## 138   FAL_SPA    BOS         26
## 139   FAL_SPA    BOS         18
## 140   FAL_SPA    BOS         14
## 141   FAL_SPA    BOS          8
## 142   FAL_SPA    BOS         32
## 143   FAL_SPA    BOS         33
## 144   FAL_SPA    BOS         16
## 145   FAL_SPA    BOS          6
## 146   FAL_SPA    BOS         14
## 147   FAL_SPA    BOS         20
## 148   FAL_SPA    BOS         13
## 149   FAL_SPA    BOS         13
## 150   FAL_SPA    BOS         28
## 151   FAL_PER    BOS          4
## 152   FAL_PER    BOS          4
## 153   FAL_PER    BOS          1
## 154   FAL_PER    BOS         11
## 155   FAL_PER    BOS         10
## 156   FAL_PER    BOS         16
## 157   FAL_PER    BOS         11
## 158   FAL_PER    BOS          2
## 159   FAL_PER    BOS         10
## 160   FAL_PER    BOS          8
## 161   FAL_PER    BOS         12
## 162   FAL_PER    BOS         10
## 163   FAL_PER    BOS          4
## 164   FAL_PER    BOS          1
## 165   FAL_PER    BOS          6
## 166   FAL_PER    BOS          5
## 167   FAL_PER    BOS          2
## 168   FAL_PER    BOS          4
## 169   FAL_PER    BOS          3
## 170   FAL_PER    BOS         11
## 171   FAL_PER    BOS          8
## 172   FAL_PER    BOS          4
## 173   FAL_PER    BOS         13
## 174   FAL_PER    BOS          2
## 175   FAL_PER    BOS          5
## 176   FAL_PER    BOS          1
## 177   FAL_PER    BOS          6
## 178   FAL_PER    BOS         11
## 179   FAL_PER    BOS          6
## 180   FAL_PER    BOS          2
## 181   PAN_HAL    BOS          2
## 182   PAN_HAL    BOS          3
## 183   PAN_HAL    BOS          4
## 184   PAN_HAL    BOS         12
## 185   PAN_HAL    BOS          2
## 186   PAN_HAL    BOS          1
## 187   PAN_HAL    BOS          8
## 188   PAN_HAL    BOS          1
## 189   PAN_HAL    BOS         12
## 190   PAN_HAL    BOS          1
## 191   PAN_HAL    BOS          9
## 192   PAN_HAL    BOS          6
## 193   PAN_HAL    BOS         11
## 194   PAN_HAL    BOS         12
## 195   PAN_HAL    BOS          5
## 196   PAN_HAL    BOS          7
## 197   PAN_HAL    BOS         10
## 198   PAN_HAL    BOS          4
## 199   PAN_HAL    BOS          3
## 200   PAN_HAL    BOS          5
## 201   PAN_HAL    BOS          8
## 202   PAN_HAL    BOS         12
## 203   PAN_HAL    BOS          2
## 204   PAN_HAL    BOS          7
## 205   PAN_HAL    BOS         11
## 206   PAN_HAL    BOS          1
## 207   PAN_HAL    BOS         11
## 208   PAN_HAL    BOS         11
## 209   PAN_HAL    BOS         12
## 210   PAN_HAL    BOS          8
## 211   BUT_JAM    BOS         10
## 212   BUT_JAM    BOS          9
## 213   BUT_JAM    BOS          1
## 214   BUT_JAM    BOS          2
## 215   BUT_JAM    BOS          5
## 216   BUT_JAM    BOS          8
## 217   BUT_JAM    BOS          1
## 218   BUT_JAM    BOS          5
## 219   BUT_JAM    BOS          3
## 220   BUT_JAM    BOS          5
## 221   BUT_JAM    BOS          1
## 222   BUT_JAM    BOS         10
## 223   BUT_JAM    BOS          1
## 224   BUT_JAM    BOS          1
## 225   BUT_JAM    BOS          8
## 226   BUT_JAM    BOS          8
## 227   BUT_JAM    BOS          8
## 228   BUT_JAM    BOS         10
## 229   BUT_JAM    BOS         10
## 230   BUT_JAM    BOS          8
## 231   BUT_JAM    BOS          9
## 232   BUT_JAM    BOS          6
## 233   BUT_JAM    BOS          2
## 234   BUT_JAM    BOS          5
## 235   BUT_JAM    BOS         10
## 236   BUT_JAM    BOS          4
## 237   BUT_JAM    BOS          9
## 238   BUT_JAM    BOS          6
## 239   BUT_JAM    BOS          8
## 240   BUT_JAM    BOS          1
## 241   FAL_SPA    BOS         33
## 242   FAL_SPA    BOS         20
## 243   FAL_SPA    BOS         16
## 244   FAL_SPA    BOS         28
## 245   FAL_SPA    BOS         16
## 246   FAL_SPA    BOS          4
## 247   FAL_SPA    BOS         10
## 248   FAL_SPA    BOS         17
## 249   FAL_SPA    BOS         27
## 250   FAL_SPA    BOS         24
## 251   FAL_SPA    BOS          3
## 252   FAL_SPA    BOS         23
## 253   FAL_SPA    BOS         26
## 254   FAL_SPA    BOS         33
## 255   FAL_SPA    BOS         28
## 256   FAL_SPA    BOS         15
## 257   FAL_SPA    BOS          4
## 258   FAL_SPA    BOS          3
## 259   FAL_SPA    BOS          6
## 260   FAL_SPA    BOS         24
## 261   FAL_SPA    BOS         20
## 262   FAL_SPA    BOS         14
## 263   FAL_SPA    BOS         14
## 264   FAL_SPA    BOS          8
## 265   FAL_SPA    BOS          9
## 266   FAL_SPA    BOS         18
## 267   FAL_SPA    BOS         28
## 268   FAL_SPA    BOS         30
## 269   FAL_SPA    BOS          2
## 270   FAL_SPA    BOS         20
## 271   FAL_PER    BOS         15
## 272   FAL_PER    BOS         13
## 273   FAL_PER    BOS         12
## 274   FAL_PER    BOS          4
## 275   FAL_PER    BOS          7
## 276   FAL_PER    BOS          7
## 277   FAL_PER    BOS          7
## 278   FAL_PER    BOS         16
## 279   FAL_PER    BOS          7
## 280   FAL_PER    BOS          5
## 281   FAL_PER    BOS          4
## 282   FAL_PER    BOS          4
## 283   FAL_PER    BOS          6
## 284   FAL_PER    BOS          3
## 285   FAL_PER    BOS          3
## 286   FAL_PER    BOS          6
## 287   FAL_PER    BOS          4
## 288   FAL_PER    BOS          2
## 289   FAL_PER    BOS         14
## 290   FAL_PER    BOS         11
## 291   FAL_PER    BOS         14
## 292   FAL_PER    BOS         11
## 293   FAL_PER    BOS         15
## 294   FAL_PER    BOS         10
## 295   FAL_PER    BOS          8
## 296   FAL_PER    BOS         13
## 297   FAL_PER    BOS         14
## 298   FAL_PER    BOS          5
## 299   FAL_PER    BOS          6
## 300   FAL_PER    BOS         12
## 301   PAN_HAL    BOS         12
## 302   PAN_HAL    BOS          9
## 303   PAN_HAL    BOS          8
## 304   PAN_HAL    BOS          8
## 305   PAN_HAL    BOS          7
## 306   PAN_HAL    BOS          3
## 307   PAN_HAL    BOS          1
## 308   PAN_HAL    BOS          3
## 309   PAN_HAL    BOS          2
## 310   PAN_HAL    BOS          8
## 311   PAN_HAL    BOS          7
## 312   PAN_HAL    BOS          6
## 313   PAN_HAL    BOS          7
## 314   PAN_HAL    BOS          7
## 315   PAN_HAL    BOS          3
## 316   PAN_HAL    BOS          8
## 317   PAN_HAL    BOS         10
## 318   PAN_HAL    BOS          3
## 319   PAN_HAL    BOS         12
## 320   PAN_HAL    BOS          1
## 321   PAN_HAL    BOS          5
## 322   PAN_HAL    BOS         12
## 323   PAN_HAL    BOS          6
## 324   PAN_HAL    BOS          8
## 325   PAN_HAL    BOS          1
## 326   PAN_HAL    BOS         11
## 327   PAN_HAL    BOS         12
## 328   PAN_HAL    BOS          1
## 329   PAN_HAL    BOS          6
## 330   PAN_HAL    BOS          2
## 331   BUT_JAM    BOS          4
## 332   BUT_JAM    BOS          8
## 333   BUT_JAM    BOS         10
## 334   BUT_JAM    BOS         10
## 335   BUT_JAM    BOS          2
## 336   BUT_JAM    BOS         10
## 337   BUT_JAM    BOS          6
## 338   BUT_JAM    BOS          1
## 339   BUT_JAM    BOS          7
## 340   BUT_JAM    BOS          8
## 341   BUT_JAM    BOS          1
## 342   BUT_JAM    BOS          4
## 343   BUT_JAM    BOS          3
## 344   BUT_JAM    BOS          3
## 345   BUT_JAM    BOS          1
## 346   BUT_JAM    BOS          6
## 347   BUT_JAM    BOS          2
## 348   BUT_JAM    BOS          2
## 349   BUT_JAM    BOS          7
## 350   BUT_JAM    BOS          6
## 351   BUT_JAM    BOS          4
## 352   BUT_JAM    BOS          2
## 353   BUT_JAM    BOS          2
## 354   BUT_JAM    BOS          9
## 355   BUT_JAM    BOS          4
## 356   BUT_JAM    BOS          5
## 357   BUT_JAM    BOS          5
## 358   BUT_JAM    BOS          7
## 359   BUT_JAM    BOS          3
## 360   BUT_JAM    BOS          2
## 361   FAL_SPA    BOS          8
## 362   FAL_SPA    BOS         18
## 363   FAL_SPA    BOS         11
## 364   FAL_SPA    BOS         16
## 365   FAL_SPA    BOS         30
## 366   FAL_SPA    BOS         32
## 367   FAL_SPA    BOS          4
## 368   FAL_SPA    BOS         22
## 369   FAL_SPA    BOS         24
## 370   FAL_SPA    BOS          1
## 371   FAL_SPA    BOS         12
## 372   FAL_SPA    BOS         20
## 373   FAL_SPA    BOS          6
## 374   FAL_SPA    BOS         23
## 375   FAL_SPA    BOS         27
## 376   FAL_SPA    BOS         20
## 377   FAL_SPA    BOS          4
## 378   FAL_SPA    BOS         26
## 379   FAL_SPA    BOS         20
## 380   FAL_SPA    BOS         16
## 381   FAL_SPA    BOS         28
## 382   FAL_SPA    BOS         17
## 383   FAL_SPA    BOS         24
## 384   FAL_SPA    BOS         23
## 385   FAL_SPA    BOS         13
## 386   FAL_SPA    BOS         19
## 387   FAL_SPA    BOS          8
## 388   FAL_SPA    BOS          2
## 389   FAL_SPA    BOS          2
## 390   FAL_SPA    BOS         32
## 391   FAL_PER    BOS          7
## 392   FAL_PER    BOS         16
## 393   FAL_PER    BOS          3
## 394   FAL_PER    BOS          9
## 395   FAL_PER    BOS         12
## 396   FAL_PER    BOS         13
## 397   FAL_PER    BOS          4
## 398   FAL_PER    BOS          4
## 399   FAL_PER    BOS          3
## 400   FAL_PER    BOS          7
## 401   FAL_PER    BOS         16
## 402   FAL_PER    BOS         13
## 403   FAL_PER    BOS          4
## 404   FAL_PER    BOS         12
## 405   FAL_PER    BOS         14
## 406   FAL_PER    BOS          4
## 407   FAL_PER    BOS          3
## 408   FAL_PER    BOS         15
## 409   FAL_PER    BOS         14
## 410   FAL_PER    BOS         14
## 411   FAL_PER    BOS          6
## 412   FAL_PER    BOS         13
## 413   FAL_PER    BOS          5
## 414   FAL_PER    BOS          5
## 415   FAL_PER    BOS          5
## 416   FAL_PER    BOS         12
## 417   FAL_PER    BOS         16
## 418   FAL_PER    BOS          7
## 419   FAL_PER    BOS         10
## 420   FAL_PER    BOS         14
## 421   PAN_HAL    BOS          6
## 422   PAN_HAL    BOS          8
## 423   PAN_HAL    BOS          9
## 424   PAN_HAL    BOS         12
## 425   PAN_HAL    BOS          4
## 426   PAN_HAL    BOS          6
## 427   PAN_HAL    BOS          5
## 428   PAN_HAL    BOS         11
## 429   PAN_HAL    BOS          9
## 430   PAN_HAL    BOS          9
## 431   PAN_HAL    BOS          2
## 432   PAN_HAL    BOS         12
## 433   PAN_HAL    BOS          3
## 434   PAN_HAL    BOS          5
## 435   PAN_HAL    BOS          2
## 436   PAN_HAL    BOS          8
## 437   PAN_HAL    BOS          5
## 438   PAN_HAL    BOS          8
## 439   PAN_HAL    BOS          2
## 440   PAN_HAL    BOS          7
## 441   PAN_HAL    BOS          4
## 442   PAN_HAL    BOS          1
## 443   PAN_HAL    BOS         10
## 444   PAN_HAL    BOS          2
## 445   PAN_HAL    BOS          5
## 446   PAN_HAL    BOS          9
## 447   PAN_HAL    BOS          7
## 448   PAN_HAL    BOS         11
## 449   PAN_HAL    BOS          4
## 450   PAN_HAL    BOS          3
## 451   BUT_JAM    BOS          4
## 452   BUT_JAM    BOS          7
## 453   BUT_JAM    BOS          7
## 454   BUT_JAM    BOS          3
## 455   BUT_JAM    BOS          2
## 456   BUT_JAM    BOS          8
## 457   BUT_JAM    BOS          2
## 458   BUT_JAM    BOS          5
## 459   BUT_JAM    BOS          1
## 460   BUT_JAM    BOS          4
## 461   BUT_JAM    BOS          5
## 462   BUT_JAM    BOS          5
## 463   BUT_JAM    BOS          2
## 464   BUT_JAM    BOS         10
## 465   BUT_JAM    BOS          7
## 466   BUT_JAM    BOS          5
## 467   BUT_JAM    BOS          7
## 468   BUT_JAM    BOS          7
## 469   BUT_JAM    BOS          8
## 470   BUT_JAM    BOS          7
## 471   BUT_JAM    BOS         10
## 472   BUT_JAM    BOS          8
## 473   BUT_JAM    BOS          1
## 474   BUT_JAM    BOS          5
## 475   BUT_JAM    BOS          8
## 476   BUT_JAM    BOS          1
## 477   BUT_JAM    BOS          8
## 478   BUT_JAM    BOS          9
## 479   BUT_JAM    BOS          4
## 480   BUT_JAM    BOS          4
## 481   FAL_SPA    PAS         24
## 482   FAL_SPA    PAS         32
## 483   FAL_SPA    PAS         18
## 484   FAL_SPA    PAS         26
## 485   FAL_SPA    PAS          8
## 486   FAL_SPA    PAS         28
## 487   FAL_SPA    PAS         21
## 488   FAL_SPA    PAS         16
## 489   FAL_SPA    PAS         19
## 490   FAL_SPA    PAS         30
## 491   FAL_SPA    PAS          2
## 492   FAL_SPA    PAS         31
## 493   FAL_SPA    PAS          8
## 494   FAL_SPA    PAS         11
## 495   FAL_SPA    PAS         20
## 496   FAL_SPA    PAS         23
## 497   FAL_SPA    PAS          8
## 498   FAL_SPA    PAS         10
## 499   FAL_SPA    PAS         12
## 500   FAL_SPA    PAS         31
## 501   FAL_SPA    PAS         14
## 502   FAL_SPA    PAS          1
## 503   FAL_SPA    PAS         17
## 504   FAL_SPA    PAS         10
## 505   FAL_SPA    PAS         15
## 506   FAL_SPA    PAS         19
## 507   FAL_SPA    PAS         10
## 508   FAL_SPA    PAS          6
## 509   FAL_SPA    PAS         26
## 510   FAL_SPA    PAS         18
## 511   FAL_PER    PAS          5
## 512   FAL_PER    PAS          3
## 513   FAL_PER    PAS          5
## 514   FAL_PER    PAS          4
## 515   FAL_PER    PAS         11
## 516   FAL_PER    PAS         14
## 517   FAL_PER    PAS          9
## 518   FAL_PER    PAS          9
## 519   FAL_PER    PAS          9
## 520   FAL_PER    PAS          9
## 521   FAL_PER    PAS          4
## 522   FAL_PER    PAS          3
## 523   FAL_PER    PAS         11
## 524   FAL_PER    PAS          5
## 525   FAL_PER    PAS          5
## 526   FAL_PER    PAS         16
## 527   FAL_PER    PAS         10
## 528   FAL_PER    PAS          2
## 529   FAL_PER    PAS          1
## 530   FAL_PER    PAS          3
## 531   FAL_PER    PAS         10
## 532   FAL_PER    PAS         10
## 533   FAL_PER    PAS         10
## 534   FAL_PER    PAS          1
## 535   FAL_PER    PAS          7
## 536   FAL_PER    PAS         14
## 537   FAL_PER    PAS          5
## 538   FAL_PER    PAS         13
## 539   FAL_PER    PAS          3
## 540   FAL_PER    PAS         10
## 541   PAN_HAL    PAS          2
## 542   PAN_HAL    PAS          8
## 543   PAN_HAL    PAS          4
## 544   PAN_HAL    PAS          9
## 545   PAN_HAL    PAS         12
## 546   PAN_HAL    PAS          4
## 547   PAN_HAL    PAS          9
## 548   PAN_HAL    PAS          6
## 549   PAN_HAL    PAS          7
## 550   PAN_HAL    PAS         12
## 551   PAN_HAL    PAS          5
## 552   PAN_HAL    PAS          6
## 553   PAN_HAL    PAS         10
## 554   PAN_HAL    PAS         12
## 555   PAN_HAL    PAS          3
## 556   PAN_HAL    PAS          6
## 557   PAN_HAL    PAS          8
## 558   PAN_HAL    PAS          5
## 559   PAN_HAL    PAS          7
## 560   PAN_HAL    PAS          2
## 561   PAN_HAL    PAS          5
## 562   PAN_HAL    PAS         11
## 563   PAN_HAL    PAS          1
## 564   PAN_HAL    PAS          9
## 565   PAN_HAL    PAS         10
## 566   PAN_HAL    PAS          7
## 567   PAN_HAL    PAS          5
## 568   PAN_HAL    PAS          6
## 569   PAN_HAL    PAS          4
## 570   PAN_HAL    PAS          8
## 571   BUT_JAM    PAS          9
## 572   BUT_JAM    PAS          7
## 573   BUT_JAM    PAS          7
## 574   BUT_JAM    PAS          1
## 575   BUT_JAM    PAS          5
## 576   BUT_JAM    PAS          8
## 577   BUT_JAM    PAS          4
## 578   BUT_JAM    PAS          3
## 579   BUT_JAM    PAS          8
## 580   BUT_JAM    PAS          3
## 581   BUT_JAM    PAS          5
## 582   BUT_JAM    PAS         10
## 583   BUT_JAM    PAS          1
## 584   BUT_JAM    PAS          6
## 585   BUT_JAM    PAS          8
## 586   BUT_JAM    PAS          8
## 587   BUT_JAM    PAS          4
## 588   BUT_JAM    PAS          8
## 589   BUT_JAM    PAS          5
## 590   BUT_JAM    PAS          5
## 591   BUT_JAM    PAS          3
## 592   BUT_JAM    PAS         10
## 593   BUT_JAM    PAS          1
## 594   BUT_JAM    PAS         10
## 595   BUT_JAM    PAS         10
## 596   BUT_JAM    PAS          3
## 597   BUT_JAM    PAS          1
## 598   BUT_JAM    PAS          3
## 599   BUT_JAM    PAS          3
## 600   BUT_JAM    PAS          2
## 601   FAL_SPA    PAS         13
## 602   FAL_SPA    PAS         27
## 603   FAL_SPA    PAS          2
## 604   FAL_SPA    PAS         28
## 605   FAL_SPA    PAS         30
## 606   FAL_SPA    PAS         19
## 607   FAL_SPA    PAS         10
## 608   FAL_SPA    PAS         31
## 609   FAL_SPA    PAS          8
## 610   FAL_SPA    PAS          9
## 611   FAL_SPA    PAS         10
## 612   FAL_SPA    PAS         19
## 613   FAL_SPA    PAS         14
## 614   FAL_SPA    PAS          3
## 615   FAL_SPA    PAS          9
## 616   FAL_SPA    PAS         26
## 617   FAL_SPA    PAS         21
## 618   FAL_SPA    PAS          8
## 619   FAL_SPA    PAS         10
## 620   FAL_SPA    PAS         23
## 621   FAL_SPA    PAS         12
## 622   FAL_SPA    PAS         18
## 623   FAL_SPA    PAS         14
## 624   FAL_SPA    PAS         30
## 625   FAL_SPA    PAS         29
## 626   FAL_SPA    PAS         29
## 627   FAL_SPA    PAS         11
## 628   FAL_SPA    PAS         12
## 629   FAL_SPA    PAS         12
## 630   FAL_SPA    PAS          4
## 631   FAL_PER    PAS         11
## 632   FAL_PER    PAS          9
## 633   FAL_PER    PAS         14
## 634   FAL_PER    PAS         13
## 635   FAL_PER    PAS         14
## 636   FAL_PER    PAS         16
## 637   FAL_PER    PAS         12
## 638   FAL_PER    PAS          9
## 639   FAL_PER    PAS          7
## 640   FAL_PER    PAS          3
## 641   FAL_PER    PAS          6
## 642   FAL_PER    PAS          8
## 643   FAL_PER    PAS          2
## 644   FAL_PER    PAS          4
## 645   FAL_PER    PAS         12
## 646   FAL_PER    PAS         13
## 647   FAL_PER    PAS          3
## 648   FAL_PER    PAS         11
## 649   FAL_PER    PAS         13
## 650   FAL_PER    PAS         15
## 651   FAL_PER    PAS         10
## 652   FAL_PER    PAS          7
## 653   FAL_PER    PAS         13
## 654   FAL_PER    PAS         15
## 655   FAL_PER    PAS          8
## 656   FAL_PER    PAS          9
## 657   FAL_PER    PAS          4
## 658   FAL_PER    PAS         15
## 659   FAL_PER    PAS          4
## 660   FAL_PER    PAS         13
## 661   PAN_HAL    PAS         12
## 662   PAN_HAL    PAS          1
## 663   PAN_HAL    PAS          2
## 664   PAN_HAL    PAS          3
## 665   PAN_HAL    PAS          9
## 666   PAN_HAL    PAS         10
## 667   PAN_HAL    PAS          4
## 668   PAN_HAL    PAS          6
## 669   PAN_HAL    PAS          6
## 670   PAN_HAL    PAS         10
## 671   PAN_HAL    PAS          5
## 672   PAN_HAL    PAS          1
## 673   PAN_HAL    PAS          4
## 674   PAN_HAL    PAS          5
## 675   PAN_HAL    PAS          6
## 676   PAN_HAL    PAS         11
## 677   PAN_HAL    PAS          8
## 678   PAN_HAL    PAS          5
## 679   PAN_HAL    PAS         11
## 680   PAN_HAL    PAS          6
## 681   PAN_HAL    PAS          3
## 682   PAN_HAL    PAS         11
## 683   PAN_HAL    PAS          1
## 684   PAN_HAL    PAS          8
## 685   PAN_HAL    PAS          5
## 686   PAN_HAL    PAS          1
## 687   PAN_HAL    PAS          6
## 688   PAN_HAL    PAS          1
## 689   PAN_HAL    PAS          2
## 690   PAN_HAL    PAS          3
## 691   BUT_JAM    PAS          1
## 692   BUT_JAM    PAS          4
## 693   BUT_JAM    PAS          5
## 694   BUT_JAM    PAS          5
## 695   BUT_JAM    PAS          2
## 696   BUT_JAM    PAS         10
## 697   BUT_JAM    PAS          6
## 698   BUT_JAM    PAS          8
## 699   BUT_JAM    PAS          4
## 700   BUT_JAM    PAS          1
## 701   BUT_JAM    PAS          6
## 702   BUT_JAM    PAS          6
## 703   BUT_JAM    PAS          3
## 704   BUT_JAM    PAS          7
## 705   BUT_JAM    PAS          4
## 706   BUT_JAM    PAS          3
## 707   BUT_JAM    PAS          5
## 708   BUT_JAM    PAS          5
## 709   BUT_JAM    PAS          2
## 710   BUT_JAM    PAS          4
## 711   BUT_JAM    PAS          5
## 712   BUT_JAM    PAS          9
## 713   BUT_JAM    PAS          8
## 714   BUT_JAM    PAS          1
## 715   BUT_JAM    PAS          8
## 716   BUT_JAM    PAS         10
## 717   BUT_JAM    PAS          9
## 718   BUT_JAM    PAS          9
## 719   BUT_JAM    PAS          3
## 720   BUT_JAM    PAS         10
## 721   FAL_SPA    PAS         24
## 722   FAL_SPA    PAS         21
## 723   FAL_SPA    PAS          8
## 724   FAL_SPA    PAS         17
## 725   FAL_SPA    PAS         17
## 726   FAL_SPA    PAS          5
## 727   FAL_SPA    PAS         18
## 728   FAL_SPA    PAS         26
## 729   FAL_SPA    PAS         27
## 730   FAL_SPA    PAS          7
## 731   FAL_SPA    PAS          7
## 732   FAL_SPA    PAS         12
## 733   FAL_SPA    PAS         17
## 734   FAL_SPA    PAS         13
## 735   FAL_SPA    PAS          9
## 736   FAL_SPA    PAS         21
## 737   FAL_SPA    PAS          6
## 738   FAL_SPA    PAS         14
## 739   FAL_SPA    PAS         30
## 740   FAL_SPA    PAS         31
## 741   FAL_SPA    PAS          2
## 742   FAL_SPA    PAS         14
## 743   FAL_SPA    PAS          3
## 744   FAL_SPA    PAS         27
## 745   FAL_SPA    PAS         32
## 746   FAL_SPA    PAS         28
## 747   FAL_SPA    PAS         23
## 748   FAL_SPA    PAS         32
## 749   FAL_SPA    PAS         19
## 750   FAL_SPA    PAS         33
## 751   FAL_PER    PAS          2
## 752   FAL_PER    PAS          5
## 753   FAL_PER    PAS         16
## 754   FAL_PER    PAS          7
## 755   FAL_PER    PAS          9
## 756   FAL_PER    PAS          3
## 757   FAL_PER    PAS         14
## 758   FAL_PER    PAS         11
## 759   FAL_PER    PAS         14
## 760   FAL_PER    PAS         10
## 761   FAL_PER    PAS          9
## 762   FAL_PER    PAS          4
## 763   FAL_PER    PAS          3
## 764   FAL_PER    PAS          5
## 765   FAL_PER    PAS          5
## 766   FAL_PER    PAS          5
## 767   FAL_PER    PAS         16
## 768   FAL_PER    PAS          3
## 769   FAL_PER    PAS         10
## 770   FAL_PER    PAS         14
## 771   FAL_PER    PAS          3
## 772   FAL_PER    PAS          7
## 773   FAL_PER    PAS          6
## 774   FAL_PER    PAS          6
## 775   FAL_PER    PAS         13
## 776   FAL_PER    PAS          6
## 777   FAL_PER    PAS         15
## 778   FAL_PER    PAS          6
## 779   FAL_PER    PAS          5
## 780   FAL_PER    PAS          9
## 781   PAN_HAL    PAS          8
## 782   PAN_HAL    PAS          2
## 783   PAN_HAL    PAS          1
## 784   PAN_HAL    PAS          5
## 785   PAN_HAL    PAS          3
## 786   PAN_HAL    PAS          9
## 787   PAN_HAL    PAS          6
## 788   PAN_HAL    PAS          8
## 789   PAN_HAL    PAS          5
## 790   PAN_HAL    PAS          4
## 791   PAN_HAL    PAS          6
## 792   PAN_HAL    PAS          1
## 793   PAN_HAL    PAS          2
## 794   PAN_HAL    PAS          6
## 795   PAN_HAL    PAS         11
## 796   PAN_HAL    PAS          7
## 797   PAN_HAL    PAS          9
## 798   PAN_HAL    PAS          1
## 799   PAN_HAL    PAS         11
## 800   PAN_HAL    PAS          9
## 801   PAN_HAL    PAS         12
## 802   PAN_HAL    PAS          5
## 803   PAN_HAL    PAS          8
## 804   PAN_HAL    PAS          3
## 805   PAN_HAL    PAS          5
## 806   PAN_HAL    PAS         12
## 807   PAN_HAL    PAS          9
## 808   PAN_HAL    PAS          2
## 809   PAN_HAL    PAS          8
## 810   PAN_HAL    PAS         11
## 811   BUT_JAM    PAS          1
## 812   BUT_JAM    PAS          9
## 813   BUT_JAM    PAS          5
## 814   BUT_JAM    PAS          4
## 815   BUT_JAM    PAS          8
## 816   BUT_JAM    PAS          8
## 817   BUT_JAM    PAS         10
## 818   BUT_JAM    PAS          3
## 819   BUT_JAM    PAS          2
## 820   BUT_JAM    PAS          3
## 821   BUT_JAM    PAS          5
## 822   BUT_JAM    PAS          6
## 823   BUT_JAM    PAS          4
## 824   BUT_JAM    PAS         10
## 825   BUT_JAM    PAS          5
## 826   BUT_JAM    PAS          4
## 827   BUT_JAM    PAS          6
## 828   BUT_JAM    PAS          1
## 829   BUT_JAM    PAS          9
## 830   BUT_JAM    PAS          2
## 831   BUT_JAM    PAS          1
## 832   BUT_JAM    PAS          4
## 833   BUT_JAM    PAS          9
## 834   BUT_JAM    PAS          1
## 835   BUT_JAM    PAS          7
## 836   BUT_JAM    PAS          3
## 837   BUT_JAM    PAS          8
## 838   BUT_JAM    PAS          3
## 839   BUT_JAM    PAS          1
## 840   BUT_JAM    PAS          1
## 841   FAL_SPA    PAS          8
## 842   FAL_SPA    PAS         21
## 843   FAL_SPA    PAS         11
## 844   FAL_SPA    PAS         12
## 845   FAL_SPA    PAS         11
## 846   FAL_SPA    PAS          3
## 847   FAL_SPA    PAS         22
## 848   FAL_SPA    PAS         12
## 849   FAL_SPA    PAS         33
## 850   FAL_SPA    PAS          8
## 851   FAL_SPA    PAS         26
## 852   FAL_SPA    PAS          4
## 853   FAL_SPA    PAS         19
## 854   FAL_SPA    PAS         23
## 855   FAL_SPA    PAS         30
## 856   FAL_SPA    PAS         32
## 857   FAL_SPA    PAS          7
## 858   FAL_SPA    PAS         19
## 859   FAL_SPA    PAS         33
## 860   FAL_SPA    PAS         26
## 861   FAL_SPA    PAS         14
## 862   FAL_SPA    PAS         24
## 863   FAL_SPA    PAS          4
## 864   FAL_SPA    PAS         12
## 865   FAL_SPA    PAS         18
## 866   FAL_SPA    PAS         25
## 867   FAL_SPA    PAS         15
## 868   FAL_SPA    PAS          7
## 869   FAL_SPA    PAS         28
## 870   FAL_SPA    PAS          9
## 871   FAL_PER    PAS          9
## 872   FAL_PER    PAS          2
## 873   FAL_PER    PAS         10
## 874   FAL_PER    PAS          9
## 875   FAL_PER    PAS          9
## 876   FAL_PER    PAS         10
## 877   FAL_PER    PAS         15
## 878   FAL_PER    PAS          4
## 879   FAL_PER    PAS          2
## 880   FAL_PER    PAS         15
## 881   FAL_PER    PAS         14
## 882   FAL_PER    PAS         10
## 883   FAL_PER    PAS          9
## 884   FAL_PER    PAS         14
## 885   FAL_PER    PAS          2
## 886   FAL_PER    PAS          7
## 887   FAL_PER    PAS          6
## 888   FAL_PER    PAS         15
## 889   FAL_PER    PAS         10
## 890   FAL_PER    PAS          1
## 891   FAL_PER    PAS         12
## 892   FAL_PER    PAS          1
## 893   FAL_PER    PAS         12
## 894   FAL_PER    PAS          3
## 895   FAL_PER    PAS         15
## 896   FAL_PER    PAS          9
## 897   FAL_PER    PAS          7
## 898   FAL_PER    PAS          4
## 899   FAL_PER    PAS          5
## 900   FAL_PER    PAS         11
## 901   PAN_HAL    PAS         10
## 902   PAN_HAL    PAS          1
## 903   PAN_HAL    PAS          2
## 904   PAN_HAL    PAS          2
## 905   PAN_HAL    PAS          7
## 906   PAN_HAL    PAS         12
## 907   PAN_HAL    PAS          6
## 908   PAN_HAL    PAS         11
## 909   PAN_HAL    PAS          6
## 910   PAN_HAL    PAS          5
## 911   PAN_HAL    PAS          7
## 912   PAN_HAL    PAS         12
## 913   PAN_HAL    PAS         11
## 914   PAN_HAL    PAS          7
## 915   PAN_HAL    PAS          6
## 916   PAN_HAL    PAS          4
## 917   PAN_HAL    PAS          3
## 918   PAN_HAL    PAS          7
## 919   PAN_HAL    PAS          1
## 920   PAN_HAL    PAS          2
## 921   PAN_HAL    PAS          2
## 922   PAN_HAL    PAS         11
## 923   PAN_HAL    PAS         11
## 924   PAN_HAL    PAS         12
## 925   PAN_HAL    PAS          4
## 926   PAN_HAL    PAS         11
## 927   PAN_HAL    PAS          6
## 928   PAN_HAL    PAS          8
## 929   PAN_HAL    PAS          1
## 930   PAN_HAL    PAS          9
## 931   BUT_JAM    PAS         10
## 932   BUT_JAM    PAS          8
## 933   BUT_JAM    PAS          1
## 934   BUT_JAM    PAS          3
## 935   BUT_JAM    PAS          6
## 936   BUT_JAM    PAS          8
## 937   BUT_JAM    PAS         10
## 938   BUT_JAM    PAS          9
## 939   BUT_JAM    PAS          4
## 940   BUT_JAM    PAS          2
## 941   BUT_JAM    PAS          9
## 942   BUT_JAM    PAS          6
## 943   BUT_JAM    PAS          3
## 944   BUT_JAM    PAS          5
## 945   BUT_JAM    PAS         10
## 946   BUT_JAM    PAS          6
## 947   BUT_JAM    PAS          5
## 948   BUT_JAM    PAS          7
## 949   BUT_JAM    PAS          9
## 950   BUT_JAM    PAS         10
## 951   BUT_JAM    PAS          1
## 952   BUT_JAM    PAS          7
## 953   BUT_JAM    PAS          1
## 954   BUT_JAM    PAS          8
## 955   BUT_JAM    PAS          4
## 956   BUT_JAM    PAS          1
## 957   BUT_JAM    PAS          1
## 958   BUT_JAM    PAS          7
## 959   BUT_JAM    PAS          6
## 960   BUT_JAM    PAS          1
## 961   FAL_SPA    HUM         28
## 962   FAL_SPA    HUM         14
## 963   FAL_SPA    HUM          8
## 964   FAL_SPA    HUM         19
## 965   FAL_SPA    HUM         12
## 966   FAL_SPA    HUM         18
## 967   FAL_SPA    HUM         17
## 968   FAL_SPA    HUM         29
## 969   FAL_SPA    HUM          4
## 970   FAL_SPA    HUM         10
## 971   FAL_SPA    HUM          9
## 972   FAL_SPA    HUM         17
## 973   FAL_SPA    HUM          7
## 974   FAL_SPA    HUM          1
## 975   FAL_SPA    HUM          2
## 976   FAL_SPA    HUM         12
## 977   FAL_SPA    HUM         16
## 978   FAL_SPA    HUM         11
## 979   FAL_SPA    HUM         32
## 980   FAL_SPA    HUM          8
## 981   FAL_SPA    HUM          7
## 982   FAL_SPA    HUM         23
## 983   FAL_SPA    HUM         18
## 984   FAL_SPA    HUM         15
## 985   FAL_SPA    HUM         16
## 986   FAL_SPA    HUM          3
## 987   FAL_SPA    HUM         30
## 988   FAL_SPA    HUM         26
## 989   FAL_SPA    HUM         30
## 990   FAL_SPA    HUM          8
## 991   FAL_PER    HUM          1
## 992   FAL_PER    HUM          4
## 993   FAL_PER    HUM          8
## 994   FAL_PER    HUM         11
## 995   FAL_PER    HUM          1
## 996   FAL_PER    HUM         11
## 997   FAL_PER    HUM          8
## 998   FAL_PER    HUM         11
## 999   FAL_PER    HUM         14
## 1000  FAL_PER    HUM          7
## 1001  FAL_PER    HUM          6
## 1002  FAL_PER    HUM         12
## 1003  FAL_PER    HUM          2
## 1004  FAL_PER    HUM          3
## 1005  FAL_PER    HUM         13
## 1006  FAL_PER    HUM         10
## 1007  FAL_PER    HUM          8
## 1008  FAL_PER    HUM         16
## 1009  FAL_PER    HUM          1
## 1010  FAL_PER    HUM          5
## 1011  FAL_PER    HUM          5
## 1012  FAL_PER    HUM          5
## 1013  FAL_PER    HUM         16
## 1014  FAL_PER    HUM          3
## 1015  FAL_PER    HUM          1
## 1016  FAL_PER    HUM         14
## 1017  FAL_PER    HUM          9
## 1018  FAL_PER    HUM          4
## 1019  FAL_PER    HUM          4
## 1020  FAL_PER    HUM          2
## 1021  PAN_HAL    HUM          5
## 1022  PAN_HAL    HUM          7
## 1023  PAN_HAL    HUM          8
## 1024  PAN_HAL    HUM          9
## 1025  PAN_HAL    HUM          5
## 1026  PAN_HAL    HUM          5
## 1027  PAN_HAL    HUM          6
## 1028  PAN_HAL    HUM          8
## 1029  PAN_HAL    HUM          3
## 1030  PAN_HAL    HUM          2
## 1031  PAN_HAL    HUM          4
## 1032  PAN_HAL    HUM         11
## 1033  PAN_HAL    HUM         11
## 1034  PAN_HAL    HUM         12
## 1035  PAN_HAL    HUM          6
## 1036  PAN_HAL    HUM          2
## 1037  PAN_HAL    HUM          2
## 1038  PAN_HAL    HUM          1
## 1039  PAN_HAL    HUM          5
## 1040  PAN_HAL    HUM          2
## 1041  PAN_HAL    HUM          2
## 1042  PAN_HAL    HUM         12
## 1043  PAN_HAL    HUM          8
## 1044  PAN_HAL    HUM         12
## 1045  PAN_HAL    HUM          7
## 1046  PAN_HAL    HUM         10
## 1047  PAN_HAL    HUM          8
## 1048  PAN_HAL    HUM          2
## 1049  PAN_HAL    HUM          9
## 1050  PAN_HAL    HUM          4
## 1051  BUT_JAM    HUM          9
## 1052  BUT_JAM    HUM          4
## 1053  BUT_JAM    HUM          8
## 1054  BUT_JAM    HUM          1
## 1055  BUT_JAM    HUM          1
## 1056  BUT_JAM    HUM          9
## 1057  BUT_JAM    HUM          8
## 1058  BUT_JAM    HUM          7
## 1059  BUT_JAM    HUM          2
## 1060  BUT_JAM    HUM          2
## 1061  BUT_JAM    HUM          4
## 1062  BUT_JAM    HUM          5
## 1063  BUT_JAM    HUM          8
## 1064  BUT_JAM    HUM          1
## 1065  BUT_JAM    HUM          7
## 1066  BUT_JAM    HUM          7
## 1067  BUT_JAM    HUM          7
## 1068  BUT_JAM    HUM          4
## 1069  BUT_JAM    HUM          6
## 1070  BUT_JAM    HUM          4
## 1071  BUT_JAM    HUM          1
## 1072  BUT_JAM    HUM          8
## 1073  BUT_JAM    HUM          4
## 1074  BUT_JAM    HUM          3
## 1075  BUT_JAM    HUM          4
## 1076  BUT_JAM    HUM          5
## 1077  BUT_JAM    HUM          2
## 1078  BUT_JAM    HUM          2
## 1079  BUT_JAM    HUM          8
## 1080  BUT_JAM    HUM          6
## 1081  FAL_SPA    HUM         15
## 1082  FAL_SPA    HUM         11
## 1083  FAL_SPA    HUM         16
## 1084  FAL_SPA    HUM         16
## 1085  FAL_SPA    HUM         18
## 1086  FAL_SPA    HUM         24
## 1087  FAL_SPA    HUM         18
## 1088  FAL_SPA    HUM         26
## 1089  FAL_SPA    HUM          8
## 1090  FAL_SPA    HUM          8
## 1091  FAL_SPA    HUM         22
## 1092  FAL_SPA    HUM          9
## 1093  FAL_SPA    HUM          3
## 1094  FAL_SPA    HUM         21
## 1095  FAL_SPA    HUM         31
## 1096  FAL_SPA    HUM          9
## 1097  FAL_SPA    HUM         32
## 1098  FAL_SPA    HUM          1
## 1099  FAL_SPA    HUM         33
## 1100  FAL_SPA    HUM         29
## 1101  FAL_SPA    HUM         12
## 1102  FAL_SPA    HUM         18
## 1103  FAL_SPA    HUM         29
## 1104  FAL_SPA    HUM         31
## 1105  FAL_SPA    HUM          8
## 1106  FAL_SPA    HUM         17
## 1107  FAL_SPA    HUM         27
## 1108  FAL_SPA    HUM         12
## 1109  FAL_SPA    HUM         12
## 1110  FAL_SPA    HUM         20
## 1111  FAL_PER    HUM          4
## 1112  FAL_PER    HUM          5
## 1113  FAL_PER    HUM          7
## 1114  FAL_PER    HUM          6
## 1115  FAL_PER    HUM          5
## 1116  FAL_PER    HUM         12
## 1117  FAL_PER    HUM         13
## 1118  FAL_PER    HUM         10
## 1119  FAL_PER    HUM         14
## 1120  FAL_PER    HUM         11
## 1121  FAL_PER    HUM         14
## 1122  FAL_PER    HUM          9
## 1123  FAL_PER    HUM         14
## 1124  FAL_PER    HUM         16
## 1125  FAL_PER    HUM         14
## 1126  FAL_PER    HUM         15
## 1127  FAL_PER    HUM          9
## 1128  FAL_PER    HUM          1
## 1129  FAL_PER    HUM         13
## 1130  FAL_PER    HUM         12
## 1131  FAL_PER    HUM          5
## 1132  FAL_PER    HUM          7
## 1133  FAL_PER    HUM          5
## 1134  FAL_PER    HUM          1
## 1135  FAL_PER    HUM          8
## 1136  FAL_PER    HUM          4
## 1137  FAL_PER    HUM         14
## 1138  FAL_PER    HUM         15
## 1139  FAL_PER    HUM          1
## 1140  FAL_PER    HUM          6
## 1141  PAN_HAL    HUM          7
## 1142  PAN_HAL    HUM         11
## 1143  PAN_HAL    HUM         10
## 1144  PAN_HAL    HUM          7
## 1145  PAN_HAL    HUM          9
## 1146  PAN_HAL    HUM         11
## 1147  PAN_HAL    HUM          4
## 1148  PAN_HAL    HUM         12
## 1149  PAN_HAL    HUM          1
## 1150  PAN_HAL    HUM          9
## 1151  PAN_HAL    HUM         11
## 1152  PAN_HAL    HUM          8
## 1153  PAN_HAL    HUM         12
## 1154  PAN_HAL    HUM          3
## 1155  PAN_HAL    HUM          3
## 1156  PAN_HAL    HUM          8
## 1157  PAN_HAL    HUM         10
## 1158  PAN_HAL    HUM          5
## 1159  PAN_HAL    HUM          9
## 1160  PAN_HAL    HUM          8
## 1161  PAN_HAL    HUM          2
## 1162  PAN_HAL    HUM          2
## 1163  PAN_HAL    HUM          7
## 1164  PAN_HAL    HUM          5
## 1165  PAN_HAL    HUM          6
## 1166  PAN_HAL    HUM         11
## 1167  PAN_HAL    HUM         10
## 1168  PAN_HAL    HUM          9
## 1169  PAN_HAL    HUM         10
## 1170  PAN_HAL    HUM          1
## 1171  BUT_JAM    HUM          2
## 1172  BUT_JAM    HUM          5
## 1173  BUT_JAM    HUM         10
## 1174  BUT_JAM    HUM          8
## 1175  BUT_JAM    HUM          3
## 1176  BUT_JAM    HUM          4
## 1177  BUT_JAM    HUM          2
## 1178  BUT_JAM    HUM          3
## 1179  BUT_JAM    HUM          8
## 1180  BUT_JAM    HUM          1
## 1181  BUT_JAM    HUM          3
## 1182  BUT_JAM    HUM          4
## 1183  BUT_JAM    HUM          9
## 1184  BUT_JAM    HUM          5
## 1185  BUT_JAM    HUM          5
## 1186  BUT_JAM    HUM          9
## 1187  BUT_JAM    HUM          9
## 1188  BUT_JAM    HUM          4
## 1189  BUT_JAM    HUM          3
## 1190  BUT_JAM    HUM          2
## 1191  BUT_JAM    HUM         10
## 1192  BUT_JAM    HUM          1
## 1193  BUT_JAM    HUM          5
## 1194  BUT_JAM    HUM          4
## 1195  BUT_JAM    HUM          9
## 1196  BUT_JAM    HUM          2
## 1197  BUT_JAM    HUM          4
## 1198  BUT_JAM    HUM          6
## 1199  BUT_JAM    HUM          1
## 1200  BUT_JAM    HUM          9
## 1201  FAL_SPA    HUM         32
## 1202  FAL_SPA    HUM          8
## 1203  FAL_SPA    HUM         14
## 1204  FAL_SPA    HUM         23
## 1205  FAL_SPA    HUM         20
## 1206  FAL_SPA    HUM         23
## 1207  FAL_SPA    HUM         12
## 1208  FAL_SPA    HUM         29
## 1209  FAL_SPA    HUM         25
## 1210  FAL_SPA    HUM         26
## 1211  FAL_SPA    HUM         24
## 1212  FAL_SPA    HUM          1
## 1213  FAL_SPA    HUM         27
## 1214  FAL_SPA    HUM         13
## 1215  FAL_SPA    HUM         12
## 1216  FAL_SPA    HUM         26
## 1217  FAL_SPA    HUM         29
## 1218  FAL_SPA    HUM         24
## 1219  FAL_SPA    HUM         18
## 1220  FAL_SPA    HUM         17
## 1221  FAL_SPA    HUM         19
## 1222  FAL_SPA    HUM         28
## 1223  FAL_SPA    HUM          9
## 1224  FAL_SPA    HUM         20
## 1225  FAL_SPA    HUM          1
## 1226  FAL_SPA    HUM         29
## 1227  FAL_SPA    HUM         32
## 1228  FAL_SPA    HUM         16
## 1229  FAL_SPA    HUM         32
## 1230  FAL_SPA    HUM         29
## 1231  FAL_PER    HUM         16
## 1232  FAL_PER    HUM          6
## 1233  FAL_PER    HUM         14
## 1234  FAL_PER    HUM          9
## 1235  FAL_PER    HUM          1
## 1236  FAL_PER    HUM          4
## 1237  FAL_PER    HUM         16
## 1238  FAL_PER    HUM         13
## 1239  FAL_PER    HUM         15
## 1240  FAL_PER    HUM         10
## 1241  FAL_PER    HUM          4
## 1242  FAL_PER    HUM          9
## 1243  FAL_PER    HUM          4
## 1244  FAL_PER    HUM          6
## 1245  FAL_PER    HUM         12
## 1246  FAL_PER    HUM         16
## 1247  FAL_PER    HUM          2
## 1248  FAL_PER    HUM          7
## 1249  FAL_PER    HUM         15
## 1250  FAL_PER    HUM          4
## 1251  FAL_PER    HUM          2
## 1252  FAL_PER    HUM         16
## 1253  FAL_PER    HUM          5
## 1254  FAL_PER    HUM         13
## 1255  FAL_PER    HUM          6
## 1256  FAL_PER    HUM         12
## 1257  FAL_PER    HUM         16
## 1258  FAL_PER    HUM         16
## 1259  FAL_PER    HUM         13
## 1260  FAL_PER    HUM          8
## 1261  PAN_HAL    HUM          5
## 1262  PAN_HAL    HUM         12
## 1263  PAN_HAL    HUM         10
## 1264  PAN_HAL    HUM          9
## 1265  PAN_HAL    HUM          7
## 1266  PAN_HAL    HUM          4
## 1267  PAN_HAL    HUM         10
## 1268  PAN_HAL    HUM          5
## 1269  PAN_HAL    HUM         10
## 1270  PAN_HAL    HUM          7
## 1271  PAN_HAL    HUM         11
## 1272  PAN_HAL    HUM          8
## 1273  PAN_HAL    HUM          4
## 1274  PAN_HAL    HUM          8
## 1275  PAN_HAL    HUM         12
## 1276  PAN_HAL    HUM          8
## 1277  PAN_HAL    HUM          3
## 1278  PAN_HAL    HUM          9
## 1279  PAN_HAL    HUM          6
## 1280  PAN_HAL    HUM          6
## 1281  PAN_HAL    HUM          5
## 1282  PAN_HAL    HUM          2
## 1283  PAN_HAL    HUM         10
## 1284  PAN_HAL    HUM          6
## 1285  PAN_HAL    HUM          8
## 1286  PAN_HAL    HUM         12
## 1287  PAN_HAL    HUM          6
## 1288  PAN_HAL    HUM          6
## 1289  PAN_HAL    HUM         11
## 1290  PAN_HAL    HUM          2
## 1291  BUT_JAM    HUM          4
## 1292  BUT_JAM    HUM          2
## 1293  BUT_JAM    HUM          7
## 1294  BUT_JAM    HUM          2
## 1295  BUT_JAM    HUM          6
## 1296  BUT_JAM    HUM          1
## 1297  BUT_JAM    HUM          9
## 1298  BUT_JAM    HUM          6
## 1299  BUT_JAM    HUM          5
## 1300  BUT_JAM    HUM          6
## 1301  BUT_JAM    HUM          6
## 1302  BUT_JAM    HUM          5
## 1303  BUT_JAM    HUM         10
## 1304  BUT_JAM    HUM         10
## 1305  BUT_JAM    HUM         10
## 1306  BUT_JAM    HUM          6
## 1307  BUT_JAM    HUM          4
## 1308  BUT_JAM    HUM          5
## 1309  BUT_JAM    HUM          1
## 1310  BUT_JAM    HUM         10
## 1311  BUT_JAM    HUM          5
## 1312  BUT_JAM    HUM          7
## 1313  BUT_JAM    HUM          9
## 1314  BUT_JAM    HUM          4
## 1315  BUT_JAM    HUM          9
## 1316  BUT_JAM    HUM          5
## 1317  BUT_JAM    HUM          1
## 1318  BUT_JAM    HUM          2
## 1319  BUT_JAM    HUM          6
## 1320  BUT_JAM    HUM         10
## 1321  FAL_SPA    HUM         29
## 1322  FAL_SPA    HUM         22
## 1323  FAL_SPA    HUM         14
## 1324  FAL_SPA    HUM          6
## 1325  FAL_SPA    HUM          8
## 1326  FAL_SPA    HUM          8
## 1327  FAL_SPA    HUM         27
## 1328  FAL_SPA    HUM         19
## 1329  FAL_SPA    HUM         19
## 1330  FAL_SPA    HUM         23
## 1331  FAL_SPA    HUM         33
## 1332  FAL_SPA    HUM         21
## 1333  FAL_SPA    HUM         12
## 1334  FAL_SPA    HUM         30
## 1335  FAL_SPA    HUM         24
## 1336  FAL_SPA    HUM         11
## 1337  FAL_SPA    HUM         17
## 1338  FAL_SPA    HUM         14
## 1339  FAL_SPA    HUM         15
## 1340  FAL_SPA    HUM         17
## 1341  FAL_SPA    HUM          4
## 1342  FAL_SPA    HUM         29
## 1343  FAL_SPA    HUM          7
## 1344  FAL_SPA    HUM          6
## 1345  FAL_SPA    HUM         33
## 1346  FAL_SPA    HUM         26
## 1347  FAL_SPA    HUM         13
## 1348  FAL_SPA    HUM         23
## 1349  FAL_SPA    HUM         14
## 1350  FAL_SPA    HUM         32
## 1351  FAL_PER    HUM         14
## 1352  FAL_PER    HUM          7
## 1353  FAL_PER    HUM         14
## 1354  FAL_PER    HUM         16
## 1355  FAL_PER    HUM          6
## 1356  FAL_PER    HUM         13
## 1357  FAL_PER    HUM         14
## 1358  FAL_PER    HUM          3
## 1359  FAL_PER    HUM         11
## 1360  FAL_PER    HUM          8
## 1361  FAL_PER    HUM         10
## 1362  FAL_PER    HUM          4
## 1363  FAL_PER    HUM          1
## 1364  FAL_PER    HUM          3
## 1365  FAL_PER    HUM          8
## 1366  FAL_PER    HUM         14
## 1367  FAL_PER    HUM         12
## 1368  FAL_PER    HUM          3
## 1369  FAL_PER    HUM          6
## 1370  FAL_PER    HUM          6
## 1371  FAL_PER    HUM          2
## 1372  FAL_PER    HUM          6
## 1373  FAL_PER    HUM          7
## 1374  FAL_PER    HUM          5
## 1375  FAL_PER    HUM          8
## 1376  FAL_PER    HUM          7
## 1377  FAL_PER    HUM         14
## 1378  FAL_PER    HUM         12
## 1379  FAL_PER    HUM         10
## 1380  FAL_PER    HUM         16
## 1381  PAN_HAL    HUM         11
## 1382  PAN_HAL    HUM         11
## 1383  PAN_HAL    HUM          8
## 1384  PAN_HAL    HUM          8
## 1385  PAN_HAL    HUM          7
## 1386  PAN_HAL    HUM          9
## 1387  PAN_HAL    HUM          1
## 1388  PAN_HAL    HUM          5
## 1389  PAN_HAL    HUM          6
## 1390  PAN_HAL    HUM          1
## 1391  PAN_HAL    HUM          3
## 1392  PAN_HAL    HUM          5
## 1393  PAN_HAL    HUM          7
## 1394  PAN_HAL    HUM          5
## 1395  PAN_HAL    HUM          6
## 1396  PAN_HAL    HUM          9
## 1397  PAN_HAL    HUM         11
## 1398  PAN_HAL    HUM         10
## 1399  PAN_HAL    HUM         11
## 1400  PAN_HAL    HUM          8
## 1401  PAN_HAL    HUM          9
## 1402  PAN_HAL    HUM          4
## 1403  PAN_HAL    HUM         10
## 1404  PAN_HAL    HUM          4
## 1405  PAN_HAL    HUM         12
## 1406  PAN_HAL    HUM          4
## 1407  PAN_HAL    HUM         12
## 1408  PAN_HAL    HUM          5
## 1409  PAN_HAL    HUM          4
## 1410  PAN_HAL    HUM          7
## 1411  BUT_JAM    HUM          5
## 1412  BUT_JAM    HUM          1
## 1413  BUT_JAM    HUM          1
## 1414  BUT_JAM    HUM          1
## 1415  BUT_JAM    HUM          4
## 1416  BUT_JAM    HUM          8
## 1417  BUT_JAM    HUM         10
## 1418  BUT_JAM    HUM         10
## 1419  BUT_JAM    HUM          8
## 1420  BUT_JAM    HUM          9
## 1421  BUT_JAM    HUM          4
## 1422  BUT_JAM    HUM          1
## 1423  BUT_JAM    HUM          8
## 1424  BUT_JAM    HUM          6
## 1425  BUT_JAM    HUM          2
## 1426  BUT_JAM    HUM          8
## 1427  BUT_JAM    HUM          4
## 1428  BUT_JAM    HUM          4
## 1429  BUT_JAM    HUM          3
## 1430  BUT_JAM    HUM          5
## 1431  BUT_JAM    HUM         10
## 1432  BUT_JAM    HUM          1
## 1433  BUT_JAM    HUM          6
## 1434  BUT_JAM    HUM          4
## 1435  BUT_JAM    HUM         10
## 1436  BUT_JAM    HUM          2
## 1437  BUT_JAM    HUM          7
## 1438  BUT_JAM    HUM          5
## 1439  BUT_JAM    HUM          1
## 1440  BUT_JAM    HUM         10
## 1441  FAL_SPA     ZU         19
## 1442  FAL_SPA     ZU          3
## 1443  FAL_SPA     ZU          1
## 1444  FAL_SPA     ZU          4
## 1445  FAL_SPA     ZU         28
## 1446  FAL_SPA     ZU          2
## 1447  FAL_SPA     ZU         16
## 1448  FAL_SPA     ZU         17
## 1449  FAL_SPA     ZU          2
## 1450  FAL_SPA     ZU         32
## 1451  FAL_SPA     ZU         17
## 1452  FAL_SPA     ZU          5
## 1453  FAL_SPA     ZU         33
## 1454  FAL_SPA     ZU         29
## 1455  FAL_SPA     ZU         15
## 1456  FAL_SPA     ZU         17
## 1457  FAL_SPA     ZU         21
## 1458  FAL_SPA     ZU         27
## 1459  FAL_SPA     ZU         33
## 1460  FAL_SPA     ZU          2
## 1461  FAL_SPA     ZU          8
## 1462  FAL_SPA     ZU         14
## 1463  FAL_SPA     ZU         33
## 1464  FAL_SPA     ZU         28
## 1465  FAL_SPA     ZU          1
## 1466  FAL_SPA     ZU          6
## 1467  FAL_SPA     ZU         16
## 1468  FAL_SPA     ZU         17
## 1469  FAL_SPA     ZU          6
## 1470  FAL_SPA     ZU         33
## 1471  FAL_PER     ZU         14
## 1472  FAL_PER     ZU          2
## 1473  FAL_PER     ZU          5
## 1474  FAL_PER     ZU          9
## 1475  FAL_PER     ZU         15
## 1476  FAL_PER     ZU          3
## 1477  FAL_PER     ZU         13
## 1478  FAL_PER     ZU         14
## 1479  FAL_PER     ZU          1
## 1480  FAL_PER     ZU         10
## 1481  FAL_PER     ZU          4
## 1482  FAL_PER     ZU          6
## 1483  FAL_PER     ZU          8
## 1484  FAL_PER     ZU          3
## 1485  FAL_PER     ZU          5
## 1486  FAL_PER     ZU         16
## 1487  FAL_PER     ZU          1
## 1488  FAL_PER     ZU         16
## 1489  FAL_PER     ZU         15
## 1490  FAL_PER     ZU          9
## 1491  FAL_PER     ZU          4
## 1492  FAL_PER     ZU         14
## 1493  FAL_PER     ZU          2
## 1494  FAL_PER     ZU         13
## 1495  FAL_PER     ZU          4
## 1496  FAL_PER     ZU         14
## 1497  FAL_PER     ZU         16
## 1498  FAL_PER     ZU          8
## 1499  FAL_PER     ZU          9
## 1500  FAL_PER     ZU          2
## 1501  PAN_HAL     ZU         11
## 1502  PAN_HAL     ZU          7
## 1503  PAN_HAL     ZU         10
## 1504  PAN_HAL     ZU          1
## 1505  PAN_HAL     ZU         12
## 1506  PAN_HAL     ZU          4
## 1507  PAN_HAL     ZU          6
## 1508  PAN_HAL     ZU         10
## 1509  PAN_HAL     ZU          6
## 1510  PAN_HAL     ZU          9
## 1511  PAN_HAL     ZU         11
## 1512  PAN_HAL     ZU          3
## 1513  PAN_HAL     ZU          5
## 1514  PAN_HAL     ZU         12
## 1515  PAN_HAL     ZU          2
## 1516  PAN_HAL     ZU          6
## 1517  PAN_HAL     ZU          9
## 1518  PAN_HAL     ZU          3
## 1519  PAN_HAL     ZU          5
## 1520  PAN_HAL     ZU          3
## 1521  PAN_HAL     ZU          1
## 1522  PAN_HAL     ZU         12
## 1523  PAN_HAL     ZU          6
## 1524  PAN_HAL     ZU         10
## 1525  PAN_HAL     ZU          3
## 1526  PAN_HAL     ZU         10
## 1527  PAN_HAL     ZU         12
## 1528  PAN_HAL     ZU          7
## 1529  PAN_HAL     ZU         12
## 1530  PAN_HAL     ZU          3
## 1531  BUT_JAM     ZU          2
## 1532  BUT_JAM     ZU          4
## 1533  BUT_JAM     ZU          3
## 1534  BUT_JAM     ZU          9
## 1535  BUT_JAM     ZU          2
## 1536  BUT_JAM     ZU          9
## 1537  BUT_JAM     ZU          7
## 1538  BUT_JAM     ZU          8
## 1539  BUT_JAM     ZU          8
## 1540  BUT_JAM     ZU          4
## 1541  BUT_JAM     ZU          6
## 1542  BUT_JAM     ZU          3
## 1543  BUT_JAM     ZU          3
## 1544  BUT_JAM     ZU          5
## 1545  BUT_JAM     ZU          6
## 1546  BUT_JAM     ZU          8
## 1547  BUT_JAM     ZU          2
## 1548  BUT_JAM     ZU          6
## 1549  BUT_JAM     ZU          6
## 1550  BUT_JAM     ZU          1
## 1551  BUT_JAM     ZU          4
## 1552  BUT_JAM     ZU         10
## 1553  BUT_JAM     ZU          5
## 1554  BUT_JAM     ZU         10
## 1555  BUT_JAM     ZU          8
## 1556  BUT_JAM     ZU          2
## 1557  BUT_JAM     ZU          6
## 1558  BUT_JAM     ZU          4
## 1559  BUT_JAM     ZU          3
## 1560  BUT_JAM     ZU          6
## 1561  FAL_SPA     ZU         23
## 1562  FAL_SPA     ZU         21
## 1563  FAL_SPA     ZU         14
## 1564  FAL_SPA     ZU         10
## 1565  FAL_SPA     ZU         32
## 1566  FAL_SPA     ZU         13
## 1567  FAL_SPA     ZU         30
## 1568  FAL_SPA     ZU         29
## 1569  FAL_SPA     ZU          6
## 1570  FAL_SPA     ZU         22
## 1571  FAL_SPA     ZU         23
## 1572  FAL_SPA     ZU          2
## 1573  FAL_SPA     ZU         13
## 1574  FAL_SPA     ZU          8
## 1575  FAL_SPA     ZU         24
## 1576  FAL_SPA     ZU         29
## 1577  FAL_SPA     ZU          2
## 1578  FAL_SPA     ZU         29
## 1579  FAL_SPA     ZU         29
## 1580  FAL_SPA     ZU         17
## 1581  FAL_SPA     ZU          5
## 1582  FAL_SPA     ZU         23
## 1583  FAL_SPA     ZU         24
## 1584  FAL_SPA     ZU         33
## 1585  FAL_SPA     ZU         25
## 1586  FAL_SPA     ZU         30
## 1587  FAL_SPA     ZU         12
## 1588  FAL_SPA     ZU         12
## 1589  FAL_SPA     ZU         23
## 1590  FAL_SPA     ZU          9
## 1591  FAL_PER     ZU         16
## 1592  FAL_PER     ZU          9
## 1593  FAL_PER     ZU          5
## 1594  FAL_PER     ZU         16
## 1595  FAL_PER     ZU          1
## 1596  FAL_PER     ZU         12
## 1597  FAL_PER     ZU         15
## 1598  FAL_PER     ZU          9
## 1599  FAL_PER     ZU          1
## 1600  FAL_PER     ZU          6
## 1601  FAL_PER     ZU         14
## 1602  FAL_PER     ZU          9
## 1603  FAL_PER     ZU         12
## 1604  FAL_PER     ZU         16
## 1605  FAL_PER     ZU          8
## 1606  FAL_PER     ZU         16
## 1607  FAL_PER     ZU          4
## 1608  FAL_PER     ZU         15
## 1609  FAL_PER     ZU          8
## 1610  FAL_PER     ZU          8
## 1611  FAL_PER     ZU          9
## 1612  FAL_PER     ZU         16
## 1613  FAL_PER     ZU         15
## 1614  FAL_PER     ZU         16
## 1615  FAL_PER     ZU         11
## 1616  FAL_PER     ZU          7
## 1617  FAL_PER     ZU          3
## 1618  FAL_PER     ZU         16
## 1619  FAL_PER     ZU          8
## 1620  FAL_PER     ZU          9
## 1621  PAN_HAL     ZU          9
## 1622  PAN_HAL     ZU          8
## 1623  PAN_HAL     ZU          8
## 1624  PAN_HAL     ZU          6
## 1625  PAN_HAL     ZU          4
## 1626  PAN_HAL     ZU          9
## 1627  PAN_HAL     ZU          5
## 1628  PAN_HAL     ZU         10
## 1629  PAN_HAL     ZU          8
## 1630  PAN_HAL     ZU          9
## 1631  PAN_HAL     ZU          1
## 1632  PAN_HAL     ZU         10
## 1633  PAN_HAL     ZU          1
## 1634  PAN_HAL     ZU         11
## 1635  PAN_HAL     ZU          8
## 1636  PAN_HAL     ZU          3
## 1637  PAN_HAL     ZU          3
## 1638  PAN_HAL     ZU          9
## 1639  PAN_HAL     ZU          6
## 1640  PAN_HAL     ZU          5
## 1641  PAN_HAL     ZU          2
## 1642  PAN_HAL     ZU          6
## 1643  PAN_HAL     ZU          6
## 1644  PAN_HAL     ZU          5
## 1645  PAN_HAL     ZU          4
## 1646  PAN_HAL     ZU          4
## 1647  PAN_HAL     ZU          7
## 1648  PAN_HAL     ZU          3
## 1649  PAN_HAL     ZU          8
## 1650  PAN_HAL     ZU          3
## 1651  BUT_JAM     ZU          8
## 1652  BUT_JAM     ZU          5
## 1653  BUT_JAM     ZU          7
## 1654  BUT_JAM     ZU          1
## 1655  BUT_JAM     ZU          2
## 1656  BUT_JAM     ZU          1
## 1657  BUT_JAM     ZU         10
## 1658  BUT_JAM     ZU          5
## 1659  BUT_JAM     ZU          1
## 1660  BUT_JAM     ZU          4
## 1661  BUT_JAM     ZU          8
## 1662  BUT_JAM     ZU          6
## 1663  BUT_JAM     ZU          9
## 1664  BUT_JAM     ZU          2
## 1665  BUT_JAM     ZU          2
## 1666  BUT_JAM     ZU         10
## 1667  BUT_JAM     ZU          2
## 1668  BUT_JAM     ZU          9
## 1669  BUT_JAM     ZU          5
## 1670  BUT_JAM     ZU          3
## 1671  BUT_JAM     ZU          4
## 1672  BUT_JAM     ZU          5
## 1673  BUT_JAM     ZU          6
## 1674  BUT_JAM     ZU          6
## 1675  BUT_JAM     ZU          9
## 1676  BUT_JAM     ZU          5
## 1677  BUT_JAM     ZU         10
## 1678  BUT_JAM     ZU          6
## 1679  BUT_JAM     ZU          1
## 1680  BUT_JAM     ZU          4
## 1681  FAL_SPA     ZU         29
## 1682  FAL_SPA     ZU          5
## 1683  FAL_SPA     ZU         27
## 1684  FAL_SPA     ZU         23
## 1685  FAL_SPA     ZU          4
## 1686  FAL_SPA     ZU          5
## 1687  FAL_SPA     ZU          2
## 1688  FAL_SPA     ZU         13
## 1689  FAL_SPA     ZU         14
## 1690  FAL_SPA     ZU          2
## 1691  FAL_SPA     ZU         28
## 1692  FAL_SPA     ZU         13
## 1693  FAL_SPA     ZU         24
## 1694  FAL_SPA     ZU         32
## 1695  FAL_SPA     ZU         31
## 1696  FAL_SPA     ZU         17
## 1697  FAL_SPA     ZU         18
## 1698  FAL_SPA     ZU          8
## 1699  FAL_SPA     ZU         22
## 1700  FAL_SPA     ZU          7
## 1701  FAL_SPA     ZU         33
## 1702  FAL_SPA     ZU          1
## 1703  FAL_SPA     ZU         10
## 1704  FAL_SPA     ZU         24
## 1705  FAL_SPA     ZU         27
## 1706  FAL_SPA     ZU         19
## 1707  FAL_SPA     ZU         18
## 1708  FAL_SPA     ZU          5
## 1709  FAL_SPA     ZU         13
## 1710  FAL_SPA     ZU         30
## 1711  FAL_PER     ZU          9
## 1712  FAL_PER     ZU          2
## 1713  FAL_PER     ZU         11
## 1714  FAL_PER     ZU          6
## 1715  FAL_PER     ZU          2
## 1716  FAL_PER     ZU         12
## 1717  FAL_PER     ZU         14
## 1718  FAL_PER     ZU          4
## 1719  FAL_PER     ZU         15
## 1720  FAL_PER     ZU         13
## 1721  FAL_PER     ZU          5
## 1722  FAL_PER     ZU          8
## 1723  FAL_PER     ZU         15
## 1724  FAL_PER     ZU         10
## 1725  FAL_PER     ZU          8
## 1726  FAL_PER     ZU         11
## 1727  FAL_PER     ZU         13
## 1728  FAL_PER     ZU          1
## 1729  FAL_PER     ZU          5
## 1730  FAL_PER     ZU          9
## 1731  FAL_PER     ZU          1
## 1732  FAL_PER     ZU          2
## 1733  FAL_PER     ZU          6
## 1734  FAL_PER     ZU          8
## 1735  FAL_PER     ZU         11
## 1736  FAL_PER     ZU         13
## 1737  FAL_PER     ZU          3
## 1738  FAL_PER     ZU          4
## 1739  FAL_PER     ZU          4
## 1740  FAL_PER     ZU          8
## 1741  PAN_HAL     ZU          6
## 1742  PAN_HAL     ZU         12
## 1743  PAN_HAL     ZU         11
## 1744  PAN_HAL     ZU         12
## 1745  PAN_HAL     ZU          2
## 1746  PAN_HAL     ZU         12
## 1747  PAN_HAL     ZU          1
## 1748  PAN_HAL     ZU          3
## 1749  PAN_HAL     ZU         11
## 1750  PAN_HAL     ZU         12
## 1751  PAN_HAL     ZU          7
## 1752  PAN_HAL     ZU         11
## 1753  PAN_HAL     ZU         10
## 1754  PAN_HAL     ZU          1
## 1755  PAN_HAL     ZU          6
## 1756  PAN_HAL     ZU          6
## 1757  PAN_HAL     ZU          6
## 1758  PAN_HAL     ZU         12
## 1759  PAN_HAL     ZU         12
## 1760  PAN_HAL     ZU          8
## 1761  PAN_HAL     ZU         11
## 1762  PAN_HAL     ZU          6
## 1763  PAN_HAL     ZU          6
## 1764  PAN_HAL     ZU          5
## 1765  PAN_HAL     ZU         10
## 1766  PAN_HAL     ZU          7
## 1767  PAN_HAL     ZU         11
## 1768  PAN_HAL     ZU          7
## 1769  PAN_HAL     ZU         11
## 1770  PAN_HAL     ZU         12
## 1771  BUT_JAM     ZU          2
## 1772  BUT_JAM     ZU         10
## 1773  BUT_JAM     ZU          2
## 1774  BUT_JAM     ZU          3
## 1775  BUT_JAM     ZU          6
## 1776  BUT_JAM     ZU          2
## 1777  BUT_JAM     ZU          7
## 1778  BUT_JAM     ZU          2
## 1779  BUT_JAM     ZU          1
## 1780  BUT_JAM     ZU          9
## 1781  BUT_JAM     ZU          6
## 1782  BUT_JAM     ZU          5
## 1783  BUT_JAM     ZU          7
## 1784  BUT_JAM     ZU          3
## 1785  BUT_JAM     ZU         10
## 1786  BUT_JAM     ZU          4
## 1787  BUT_JAM     ZU         10
## 1788  BUT_JAM     ZU          2
## 1789  BUT_JAM     ZU          3
## 1790  BUT_JAM     ZU          4
## 1791  BUT_JAM     ZU          2
## 1792  BUT_JAM     ZU          9
## 1793  BUT_JAM     ZU         10
## 1794  BUT_JAM     ZU          2
## 1795  BUT_JAM     ZU          5
## 1796  BUT_JAM     ZU          5
## 1797  BUT_JAM     ZU          9
## 1798  BUT_JAM     ZU          7
## 1799  BUT_JAM     ZU          1
## 1800  BUT_JAM     ZU          5
## 1801  FAL_SPA     ZU          2
## 1802  FAL_SPA     ZU         29
## 1803  FAL_SPA     ZU          1
## 1804  FAL_SPA     ZU         21
## 1805  FAL_SPA     ZU          3
## 1806  FAL_SPA     ZU         14
## 1807  FAL_SPA     ZU          4
## 1808  FAL_SPA     ZU         16
## 1809  FAL_SPA     ZU         23
## 1810  FAL_SPA     ZU         20
## 1811  FAL_SPA     ZU         28
## 1812  FAL_SPA     ZU          1
## 1813  FAL_SPA     ZU         29
## 1814  FAL_SPA     ZU          7
## 1815  FAL_SPA     ZU         30
## 1816  FAL_SPA     ZU         33
## 1817  FAL_SPA     ZU          9
## 1818  FAL_SPA     ZU         14
## 1819  FAL_SPA     ZU         21
## 1820  FAL_SPA     ZU         14
## 1821  FAL_SPA     ZU          9
## 1822  FAL_SPA     ZU         17
## 1823  FAL_SPA     ZU         32
## 1824  FAL_SPA     ZU         11
## 1825  FAL_SPA     ZU          2
## 1826  FAL_SPA     ZU         25
## 1827  FAL_SPA     ZU         16
## 1828  FAL_SPA     ZU         20
## 1829  FAL_SPA     ZU         32
## 1830  FAL_SPA     ZU         30
## 1831  FAL_PER     ZU          4
## 1832  FAL_PER     ZU          6
## 1833  FAL_PER     ZU          7
## 1834  FAL_PER     ZU          1
## 1835  FAL_PER     ZU          5
## 1836  FAL_PER     ZU         16
## 1837  FAL_PER     ZU          3
## 1838  FAL_PER     ZU          8
## 1839  FAL_PER     ZU         14
## 1840  FAL_PER     ZU          8
## 1841  FAL_PER     ZU         10
## 1842  FAL_PER     ZU         11
## 1843  FAL_PER     ZU          4
## 1844  FAL_PER     ZU          7
## 1845  FAL_PER     ZU         14
## 1846  FAL_PER     ZU         10
## 1847  FAL_PER     ZU          3
## 1848  FAL_PER     ZU         14
## 1849  FAL_PER     ZU         15
## 1850  FAL_PER     ZU          4
## 1851  FAL_PER     ZU          1
## 1852  FAL_PER     ZU          2
## 1853  FAL_PER     ZU          2
## 1854  FAL_PER     ZU         15
## 1855  FAL_PER     ZU         16
## 1856  FAL_PER     ZU          1
## 1857  FAL_PER     ZU          1
## 1858  FAL_PER     ZU          9
## 1859  FAL_PER     ZU          3
## 1860  FAL_PER     ZU         10
## 1861  PAN_HAL     ZU          4
## 1862  PAN_HAL     ZU          9
## 1863  PAN_HAL     ZU          1
## 1864  PAN_HAL     ZU          8
## 1865  PAN_HAL     ZU          6
## 1866  PAN_HAL     ZU          6
## 1867  PAN_HAL     ZU          6
## 1868  PAN_HAL     ZU          1
## 1869  PAN_HAL     ZU         11
## 1870  PAN_HAL     ZU          6
## 1871  PAN_HAL     ZU          2
## 1872  PAN_HAL     ZU          9
## 1873  PAN_HAL     ZU          9
## 1874  PAN_HAL     ZU         10
## 1875  PAN_HAL     ZU         11
## 1876  PAN_HAL     ZU          3
## 1877  PAN_HAL     ZU          8
## 1878  PAN_HAL     ZU          7
## 1879  PAN_HAL     ZU         10
## 1880  PAN_HAL     ZU          3
## 1881  PAN_HAL     ZU          5
## 1882  PAN_HAL     ZU          8
## 1883  PAN_HAL     ZU         10
## 1884  PAN_HAL     ZU          1
## 1885  PAN_HAL     ZU          9
## 1886  PAN_HAL     ZU          7
## 1887  PAN_HAL     ZU         11
## 1888  PAN_HAL     ZU          5
## 1889  PAN_HAL     ZU          3
## 1890  PAN_HAL     ZU          2
## 1891  BUT_JAM     ZU          3
## 1892  BUT_JAM     ZU         10
## 1893  BUT_JAM     ZU          4
## 1894  BUT_JAM     ZU          4
## 1895  BUT_JAM     ZU          4
## 1896  BUT_JAM     ZU          3
## 1897  BUT_JAM     ZU          1
## 1898  BUT_JAM     ZU          1
## 1899  BUT_JAM     ZU          6
## 1900  BUT_JAM     ZU          1
## 1901  BUT_JAM     ZU          3
## 1902  BUT_JAM     ZU          7
## 1903  BUT_JAM     ZU          3
## 1904  BUT_JAM     ZU         10
## 1905  BUT_JAM     ZU          8
## 1906  BUT_JAM     ZU          1
## 1907  BUT_JAM     ZU          2
## 1908  BUT_JAM     ZU          8
## 1909  BUT_JAM     ZU          5
## 1910  BUT_JAM     ZU          1
## 1911  BUT_JAM     ZU          3
## 1912  BUT_JAM     ZU          4
## 1913  BUT_JAM     ZU          8
## 1914  BUT_JAM     ZU          8
## 1915  BUT_JAM     ZU          1
## 1916  BUT_JAM     ZU          3
## 1917  BUT_JAM     ZU          5
## 1918  BUT_JAM     ZU          1
## 1919  BUT_JAM     ZU          7
## 1920  BUT_JAM     ZU          9
ggplot(data=dataf, aes(x=reorder(Especies,Abundancia), y=Abundancia, fill=Sitios)) + geom_bar(stat="identity", position="dodge") + labs(x="Especies",y="Abundancia") + scale_fill_manual("Sitios", values = c("BOS" = "#27AE60", "HUM" = "#3498DB", "PAS" = "#F1C40F", "ZU" = "#7F8C8D"))

# SE TOMA LA DECISION DE USAR OTRO GRAF XQ ESCALA NO COINCIDE

# PARA REPORTAR RESULTADOS DE EXPERIMENTO:
T10G<-c(0.446698161,0.179728779,0.627716076,1.010205394,0.130444902,1.413243637,0.302696537,0.835317793,0,0.799502599)
mean(T10G)  # 0.5745554
## [1] 0.5745554
sd(T10G)
## [1] 0.4457334
T15G<-c(0.751361251,0.344687488,0.101090408,0.728284237,1.406586705,0.883504893,1.241286268,0,0.19890516,0.998550348)
mean(T15G) # 0.6654257
## [1] 0.6654257
sd(T15G)
## [1] 0.4867787
T25G<-c(2.448914392,2.916482249,2.442398114,1.555197822,2.713667911,3.225865253,1.901868629,1.819783841,2.207417681,0.864375646)
mean(T25G) # 2.209597
## [1] 2.209597
sd(T25G)
## [1] 0.6978188
T35G<-c(0.784284516,0.86464155,0.589636935,1.482825868,0,0.595385042,0.799215124,1.049051366,1.518858427,0.715094445)
mean(T35G)  # 0.8398993
## [1] 0.8398993
sd(T35G)
## [1] 0.4430789
T10P<-c(3.139374257,1.304773372,0.086640707,3.458546482,3.614733365,2.245183599,0.010503563,0,4.570463289,0.307572441)
mean(T10P)  # 1.873779
## [1] 1.873779
sd(T10P)
## [1] 1.748258
T15P<-c(2.626811205,2.966476409,2.961043711,4.068104976,2.278065886,4.695521743,4.151473974,3.865028751,3.186247818,2.066820612)
mean(T15P)  # 3.28656
## [1] 3.28656
sd(T15P)
## [1] 0.871605
T25P<-c(5.916183182,5.911558375,3.870983715,6.075054829,3.958550918,3.975439264,3.769403852,5.007383662,5.386924372,4.900727216)
mean(T25P)  # 4.877221
## [1] 4.877221
sd(T25P)
## [1] 0.9284594
T35P<-c(4.2974188,4.296696882,3.638162131,4.854464373,4.112298984,4.82375161,3.513894431,4.033807044,3.962681611,3.919699049)
mean(T35P)  # 4.145287
## [1] 4.145287
sd(T35P)
## [1] 0.4425806
# PARA REPORTAR RESULTADOS DE MUESTREO:
AbSpaBOS<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="FAL_SPA" & BD_AvesRapaces$Sitio == "BOS"] 
mean(AbSpaBOS)   # 17.48333
## [1] 17.48333
sd(AbSpaBOS)  # 9.462872
## [1] 9.462872
AbSpaHUM<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="FAL_SPA" & BD_AvesRapaces$Sitio == "HUM"] 
mean(AbSpaHUM)  # 18
## [1] 18
sd(AbSpaHUM)  # 9.013529
## [1] 9.013529
AbSpaPAS<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="FAL_SPA" & BD_AvesRapaces$Sitio == "PAS"] 
mean(AbSpaPAS) # 17.2
## [1] 17.2
sd(AbSpaPAS) # 9.075408
## [1] 9.075408
AbSpaZU<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="FAL_SPA" & BD_AvesRapaces$Sitio == "ZU"] 
mean(AbSpaZU)  # 17.28333
## [1] 17.28333
sd(AbSpaZU)  # 10.37708
## [1] 10.37708
AbPerBOS<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="FAL_PER" & BD_AvesRapaces$Sitio == "BOS"] 
mean(AbPerBOS)  # 8.575
## [1] 8.575
sd(AbPerBOS) # 4.412858
## [1] 4.412858
AbPerHUM<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="FAL_PER" & BD_AvesRapaces$Sitio == "HUM"] 
mean(AbPerHUM) # 8.625
## [1] 8.625
sd(AbPerHUM)  # 4.773989
## [1] 4.773989
AbPerPAS<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="FAL_PER" & BD_AvesRapaces$Sitio == "PAS"] 
mean(AbPerPAS)  # 8.391667
## [1] 8.391667
sd(AbPerPAS)  # 4.352717
## [1] 4.352717
AbPerZU<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="FAL_PER" & BD_AvesRapaces$Sitio == "ZU"] 
mean(AbPerZU)  #  8.516667
## [1] 8.516667
sd(AbPerZU)  # 4.970473
## [1] 4.970473
AbJamBOS<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="BUT_JAM" & BD_AvesRapaces$Sitio == "BOS"] 
mean(AbJamBOS) # 5.55
## [1] 5.55
sd(AbJamBOS)  # 2.981314
## [1] 2.981314
AbJamHUM<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="BUT_JAM" & BD_AvesRapaces$Sitio == "HUM"] 
mean(AbJamHUM) # 5.233333
## [1] 5.233333
sd(AbJamHUM)  # 2.932671
## [1] 2.932671
AbJamPAS<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="BUT_JAM" & BD_AvesRapaces$Sitio == "PAS"] 
mean(AbJamPAS) # 5.291667
## [1] 5.291667
sd(AbJamPAS)  # 2.996906
## [1] 2.996906
AbJamZU<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="BUT_JAM" & BD_AvesRapaces$Sitio == "ZU"] 
mean(AbJamZU)  # 5.025
## [1] 5.025
sd(AbJamZU)  # 2.879842
## [1] 2.879842
AbHalBOS<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="PAN_HAL" & BD_AvesRapaces$Sitio == "BOS"] 
mean(AbHalBOS)  # 6.3
## [1] 6.3
sd(AbHalBOS)  # 3.647493
## [1] 3.647493
AbHalHUM<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="PAN_HAL" & BD_AvesRapaces$Sitio == "HUM"] 
mean(AbHalHUM)  # 7.033333
## [1] 7.033333
sd(AbHalHUM)  # 3.243438
## [1] 3.243438
AbHalPAS<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="PAN_HAL" & BD_AvesRapaces$Sitio == "PAS"] 
mean(AbHalPAS)  # 6.291667
## [1] 6.291667
sd(AbHalPAS)  # 3.438281
## [1] 3.438281
AbHalZU<-BD_AvesRapaces$Abundancia[BD_AvesRapaces$Especies=="PAN_HAL" & BD_AvesRapaces$Sitio == "ZU"] 
mean(AbHalZU)  # 6.916667
## [1] 6.916667
sd(AbHalZU)  # 3.419134
## [1] 3.419134
# -------- #
THal<-BD_AvesRapaces$Talla[BD_AvesRapaces$Especies=="PAN_HAL"] 
mean(THal) # 61.11042
## [1] 61.11042
sd(THal)
## [1] 3.066809
TJam<-BD_AvesRapaces$Talla[BD_AvesRapaces$Especies=="BUT_JAM"] 
mean(TJam)  # 58.4875
## [1] 58.4875
sd(TJam)
## [1] 4.034798
TSpa<-BD_AvesRapaces$Talla[BD_AvesRapaces$Especies=="FAL_SPA"] 
mean(TSpa)  #  18.63958
## [1] 18.63958
sd(TSpa)
## [1] 6.459549
TPer<-BD_AvesRapaces$Talla[BD_AvesRapaces$Especies=="FAL_PER"] 
mean(TPer)   # 42.36042
## [1] 42.36042