Marcar y recuperar

Marcar y recuperar es un método comúnmente usado en ecología para estimar el tamaño de una población animal donde no es práctico contar a cada individuo.[1]​ Una parte de la población es capturada, marcada y liberada. Posteriormente, se capturará otra porción y se contará el número de individuos marcados dentro de la muestra. Dado que el número de individuos marcados dentro de la segunda muestra debe ser proporcional al número de individuos marcados en toda la población, se puede obtener una estimación del tamaño total de la población dividiendo el número de individuos marcados por la proporción de individuos marcados en la segunda muestra. El método es más útil cuando no es práctico contar a todos los individuos de la población. Otros nombres para este método, o métodos estrechamente relacionados, incluyen captura-recaptura, captura-marca-recaptura, marca-recaptura, avistamiento de nuevo, marca-liberación-recaptura, estimación de sistemas múltiples, recuperación de banda, el método Petersen,[2]​ y el método Lincoln.

Collar etiquetado rock hyrax
Grajilla con un anillo de aluminio numerado en su tarso izquierdo
El biólogo está marcando un caracol ámbar ovado de Chittenango para monitorear la población.
right side view of a snail with a number 87 on its shell
Caracol de color ámbar ovado de Chittenango marcado.

Otra aplicación importante de estos métodos es la epidemiología,[3]​ donde se utilizan para estimar la integridad de la verificación de registros de enfermedades. Las aplicaciones típicas incluyen estimar el número de personas que necesitan servicios particulares (es decir, servicios para niños con discapacidades de aprendizaje, servicios para ancianos médicamente frágiles que viven en la comunidad) o con condiciones particulares (es decir, drogadictos ilegales, personas infectadas con el VIH, etc.).[4]

Trabajo de campo relacionado con marca-recaptura

editar

Por lo general, un investigador visita un área de estudio y usa trampas para capturar a un grupo de individuos vivos. Cada uno de estos individuos se marca con un identificador único (por ejemplo, una etiqueta o banda numerada) y luego se libera ileso de regreso al medio ambiente. En 1896, C. G. Johannes Petersen utilizó por primera vez un método de marcado-recaptura para un estudio ecológico para estimar las poblaciones de solla, Pleuronectes platesa.[5]

Se deja pasar suficiente tiempo para que los individuos marcados se redistribuyan entre la población no marcada.[5]

A continuación, el investigador regresa y captura otra muestra de individuos. Algunas personas de esta segunda muestra habrán sido marcadas durante la visita inicial y ahora se conocen como recapturas.[6]​ Otros animales capturados durante la segunda visita, no habrán sido capturados durante la primera visita al área de estudio. Estos animales sin marcar generalmente reciben una etiqueta o banda durante la segunda visita y luego son liberados.[5]

El tamaño de la población se puede estimar a partir de tan solo dos visitas al área de estudio. Por lo general, se realizan más de dos visitas, especialmente si se desean estimaciones de supervivencia o movimiento. Independientemente del número total de visitas, el investigador simplemente registra la fecha de cada captura de cada individuo. Las "historias de captura" generadas se analizan matemáticamente para estimar el tamaño, la supervivencia o el movimiento de la población.[5]

Al capturar y marcar organismos, los ecólogos deben considerar el bienestar de los organismos. Si el identificador elegido daña al organismo, entonces su comportamiento puede volverse irregular.

Condiciones

editar

Independientemente del método escogido, este debería asegurar que[7]​:

  • Todos los animales marcados deber poder ser reconocidos como marcados durante la totalidad del estudio.
  • Todos los animaless de la población deben tener la misma probabilidad de ser capturados a lo largo del estudio
  • Los animales no deben de presentar un aprendizaje hacia las trampas utilizadas para su captura.
  • El proceso de marcado no debe afectar su supervivencia natural del animal.
  • Se libere a los animales marcados en el mismo lugar en que fueron capturados.

En algunos métodos se requiere la existencia de un periodo de tiempo apreciable entre el marcado y la recaptura para garantizar que los animlaes marcados se distribuyan entre la población.[7]

Notación

editar

Dejando

N = Número de animales en la población
n = Número de animales marcados en la primera visita
K = Número de animales capturados en la segunda visita
k = Número de animales recapturados que fueron marcados

Un biólogo quiere estimar el tamaño de una población de tortugas en un lago. Captura 10 tortugas en su primera visita al lago y les marca el lomo con pintura. Una semana después regresa al lago y captura 15 tortugas. Cinco de estas 15 tortugas tienen pintura en la espalda, lo que indica que son animales recapturados. Este ejemplo es (n, K, k) = (10, 15, 5). El problema es estimar N.

 

Estimador Lincoln-Petersen

editar

El método Lincoln-Petersen[8]​ (también conocido como índice Petersen-Lincoln[5]​ o índice Lincoln) se puede utilizar para estimar el tamaño de la población si solo se realizan dos visitas al área de estudio. Este método asume que la población de estudio es "cerrada". En otras palabras, las dos visitas al área de estudio son lo suficientemente cercanas en el tiempo para que ninguna persona muera, nazca o se mude dentro o fuera del área de estudio entre visitas. El modelo también asume que no se caen marcas de los animales entre visitas al sitio de campo por parte del investigador, y que el investigador registra correctamente todas las marcas.

Dadas esas condiciones, el tamaño de población estimado es:

 

Derivación

editar

Se asume que todos los individuos tienen la misma probabilidad de ser capturados en la segunda muestra, independientemente de si fueron capturados previamente en la primera muestra (con solo dos muestras, esta suposición no se puede probar directamente).[9]

Esto implica que, en la segunda muestra, la proporción de individuos marcados que son capturados ( ) debe ser igual a la proporción de la población total marcada ( ). Por ejemplo, si la mitad de los individuos marcados fueran recapturados, se supondría que la mitad de la población total se incluyó en la segunda muestra.

En símbolos,

 

Una reordenación de esto da

 

la fórmula utilizada para el método Lincoln-Petersen.[9]

Cálculo de muestra

editar

En el ejemplo (n, K, k) = (10, 15, 5) el método Lincoln-Petersen estima que hay 30 tortugas en el lago.

 

Estimador de Chapman

editar

El estimador de Lincoln-Petersen es asintóticamente insesgado a medida que el tamaño de la muestra se acerca al infinito, pero está sesgado en tamaños de muestra pequeños.[10]​ El estimador de Chapman proporciona un estimador alternativo menos sesgado del tamaño de la población:

 

Cálculo de muestra

editar

El ejemplo (K, n, k) = (10, 15, 5) da

 

Tenga en cuenta que la respuesta proporcionada por esta ecuación debe ser truncada, no redondeada. Así, el método Chapman estima 28 tortugas en el lago.

Sorprendentemente, la estimación de Chapman fue una conjetura de un rango de posibles estimadores:

"En la práctica, el número entero inmediatamente menor que (K+1)(n+1)/(k+1) o incluso Kn/(k+1) será ser la estimación. La forma anterior es más conveniente para propósitos matemáticos."[10]

Chapman también encontró que el estimador podría tener un sesgo negativo considerable para Kn/N pequeño, pero no le preocupaba porque las desviaciones estándar estimadas eran grandes para estos casos.

Intervalo de confianza

editar

Un aproximado  , el intervalo de confianza para el tamaño de la población N se puede obtener como:  ,

donde   corresponde a la   cuantil de una variable aleatoria normal estándar, y

  .

Se ha demostrado que este intervalo de confianza tiene probabilidades de cobertura reales cercanas a la nominal.  nivel incluso para poblaciones pequeñas y probabilidades de captura extremas (cercanas a 0 o 1), en cuyo caso otros intervalos de confianza no alcanzan los niveles de cobertura nominal.[11]

Estimación bayesiana

editar

El valor medio ± desviación estándar es

 

donde

  para  
  para  

El ejemplo (K, n, k) = (10, 15, 5) da la estimación N ≈ 42 ± 21.5

Probabilidad de captura

editar
 
Campañol de banco, Myodes glareolus, en un estudio de captura y liberación de la población de pequeños mamíferos para London Wildlife Trust en la reserva natural local de Gunnersbury Triangle

La probabilidad de captura se refiere a la probabilidad de detectar un animal individual o una persona de interés,[12]​ y se ha utilizado tanto en ecología como en epidemiología para detectar enfermedades animales o humanas,[13]​ respectivamente.

La probabilidad de captura se define a menudo como un modelo de dos variables, en el que f se define como la fracción de un recurso finito dedicado a detectar el animal o la persona de interés de un sector de alto riesgo de una población animal o humana, y q es la frecuencia de tiempo en que el problema (por ejemplo, una enfermedad animal) ocurre en el sector de alto riesgo frente al de bajo riesgo.[14]​ Por ejemplo, una aplicación del modelo en la década de 1920 fue detectar portadores de fiebre tifoidea en Londres, que llegaban de zonas con altas tasas de tuberculosis (probabilidad q de que un pasajero con la enfermedad viniera de dicha zona, donde q>0,5), o tasas bajas (probabilidad 1- q).[15]​ Se postuló que solo se pudieron detectar 5 de cada 100 de los viajeros y que 10 de cada 100 procedían de la zona de alto riesgo. Entonces, la probabilidad de captura P se definió como:

 

donde el primer término se refiere a la probabilidad de detección (probabilidad de captura) en una zona de alto riesgo, y el último término se refiere a la probabilidad de detección en una zona de bajo riesgo. Es importante destacar que la fórmula se puede reescribir como una ecuación lineal en términos de f :

 

Debido a que esta es una función lineal, se deduce que para ciertas versiones de q para las cuales la pendiente de esta línea (el primer término multiplicado por f) es positiva, todo el recurso de detección debe dedicarse a la población de alto riesgo (f debe establecerse en 1 para maximizar la probabilidad de captura), mientras que para otro valor de q, para el cual la pendiente de la línea es negativa, toda la detección debe dedicarse a la población de bajo riesgo (f debe establecerse en 0. Podemos resolver la ecuación anterior para los valores de q para los cuales la pendiente será positiva para determinar los valores para los cuales f debería establecerse en 1 para maximizar la probabilidad de captura:

 

que se simplifica a:

 

Este es un ejemplo de optimización lineal.[14]​ En casos más complejos, donde más de un recurso f se dedica a más de dos áreas, se suele utilizar la optimización multivariante, a través del algoritmo simplex o sus derivados.

Más de dos visitas

editar

Existen modelos estadísticos muy elaborados disponibles para el análisis de estos experimentos.[16]​ Un modelo simple que se adapta fácilmente a las tres fuentes, o al estudio de tres visitas, es ajustar un modelo de regresión de Poisson. Los modelos sofisticados de marcado y recaptura pueden adaptarse a varios paquetes para el lenguaje de programación Open Source R. Estos incluyen "Captura-Recaptura Espacialmente Explícita (secr)",[17]​ "Modelos Loglineales para Experimentos de Captura-Recaptura (Rcapture)",[18]​ y "Muestreo de Distancia de Marca-Recaptura (mrds)".[19]​ Estos modelos también pueden adaptarse a programas especializados como MARK[20]​ o M-SURGE.[21]

Otros métodos relacionados que se utilizan a menudo incluyen el modelo Jolly-Seber (utilizado en poblaciones abiertas y para estimaciones de censos múltiples) y los estimadores de Schnabel[22]​ (descritos anteriormente como una expansión del método Lincoln-Petersen para poblaciones cerradas).[23]

Enfoques integrados

editar

El modelado de datos de marcado-recaptura tiende hacia un enfoque más integrador,[24]​ que combina los datos de marcado-recaptura con modelos de dinámica de población y otros tipos de datos. El enfoque integrado es más exigente desde el punto de vista computacional, pero extrae más información de los datos mejorando los parámetros y las estimaciones de incertidumbre.[25]

Véase también

editar

Referencias

editar
  1. «Mark-Recapture». www2.nau.edu. 
  2. Krebs, Charles J. (2009). Ecology (6th edición). p. 119. ISBN 978-0-321-50743-3. 
  3. Chao, Anne; Tsay, P. K.; Lin, Sheng-Hsiang; Shau, Wen-Yi; Chao, Day-Yu (2001). «The applications of capture-recapture models to epidemiological data». Statistics in Medicine (en inglés) 20 (20): 3123-3157. ISSN 1097-0258. doi:10.1002/sim.996. 
  4. Allen (2019). «Estimating the Number of People Who Inject Drugs in A Rural County in Appalachia». American Journal of Public Health 109 (3): 445-450. PMC 6366498. PMID 30676803. doi:10.2105/AJPH.2018.304873. 
  5. a b c d e Southwood, T. R. E.; Henderson, P. (2000). Ecological Methods (3rd edición). Oxford: Blackwell Science. 
  6. «Definition of RECAPTURE». www.merriam-webster.com (en inglés). 
  7. a b Rojas González, Ramón Isaac; Zúñiga Vega, José Jaime (2005). Técnicas para el estudio de poblaciones de fauna silvestre. UNAM, Facultad de Estudios Superiores Iztacala, Unidad de Biología, Tecnología y Prototipos. ISBN 970-9000-32-2. OCLC 69656368. Consultado el 8 de febrero de 2022. 
  8. Seber, G. A. F. (1982). The Estimation of Animal Abundance and Related Parameters. Caldwel, New Jersey: Blackburn Press. ISBN 1-930665-55-5. 
  9. a b Charles J. Krebs (1999). Ecological Methodology (2nd edición). ISBN 9780321021731. 
  10. a b Chapman, D.G. (1951). Some properties of the hypergeometric distribution with applications to zoological sample censuses. 
  11. Sadinle, Mauricio (1 de octubre de 2009). «Transformed Logit Confidence Intervals for Small Populations in Single Capture–Recapture Estimation». Communications in Statistics - Simulation and Computation 38 (9): 1909-1924. ISSN 0361-0918. doi:10.1080/03610910903168595. 
  12. Drenner, Ray (1978). «Capture probability: the role of zooplankter escape in the selective feeding of planktivorous fish». Journal of the Fisheries Board of Canada 35 (10): 1370-1373. doi:10.1139/f78-215. 
  13. MacKenzie, Darryl (2002). «How should detection probability be incorporated into estimates of relative abundance?». Ecology 83 (9): 2387-2393. doi:10.1890/0012-9658(2002)083[2387:hsdpbi]2.0.co;2. 
  14. a b Bolker, Benjamin (2008). Ecological Models and Data in R. Princeton University Press. ISBN 9781400840908. 
  15. Unknown (1921). «The Health of London». Hosp Health Rev 1: 71-2. 
  16. McCrea, R.S. and Morgan, B.J.T. (2014) «Analysis of capture-recapture data.». Consultado el 19 de noviembre de 2014.  «Chapman and Hall/CRC Press». Consultado el 19 de noviembre de 2014. 
  17. Efford, Murray (2 de septiembre de 2016). «Spatially Explicit Capture-Recapture (secr)». Comprehensive R Archive Network (CRAN). Consultado el 2 de septiembre de 2016. 
  18. Rivest, Louis-Paul (1 de septiembre de 2014). «Loglinear Models for Capture-Recapture Experiments (Rcapture)». Comprehensive R Archive Network (CRAN). Consultado el 2 de septiembre de 2016. 
  19. Laake, Jeff (17 de agosto de 2015). «Mark-Recapture Distance Sampling (mrds)». Comprehensive R Archive Network (CRAN). 
  20. «Program MARK». Archivado desde el original el 21 de febrero de 2006. Consultado el 29 de mayo de 2013. 
  21. «Logiciels». Archivado desde el original el 24 de julio de 2009. 
  22. Schnabel, Z. E. (1938). «The Estimation of the Total Fish Population of a Lake». American Mathematical Monthly 45 (6): 348-352. doi:10.2307/2304025. 
  23. William J. Sutherland, ed. (1996). Ecological Census Techniques: A Handbook. Cambridge University Press. ISBN 0-521-47815-4. 
  24. Maunder M.N. (2003) Paradigm shifts in fisheries stock assessment: from integrated analysis to Bayesian analysis and back again. Natural Resource Modeling 16:465–475
  25. Maunder, M.N. (2001) Integrated Tagging and Catch-at-Age Analysis (ITCAAN). In Spatial Processes and Management of Fish Populations, edited by G.H. Kruse, N. Bez, A. Booth, M.W. Dorn, S. Hills, R.N. Lipcius, D. Pelletier, C. Roy, S.J. Smith, and D. Witherell, Alaska Sea Grant College Program Report No. AK-SG-01-02, University of Alaska Fairbanks, pp. 123–146.

Bibliografía

editar

Otras lecturas

editar
  • Bonett, D.G.; Woodward, J.A.; Bentler, P.M. (1986). «A Linear Model for Estimating the Size of a Closed Population». British Journal of Mathematical and Statistical Psychology 39: 28-40. PMID 3768264. doi:10.1111/j.2044-8317.1986.tb00843.x. 
  • Evans, M.A.; Bonett, D.G.; McDonald, L. (1994). «A General Theory for Analyzing Capture-recapture Data in Closed Populations». Biometrics 50 (2): 396-405. doi:10.2307/2533383. 
  • Lincoln, F. C. (1930). «Calculating Waterfowl Abundance on the Basis of Banding Returns». United States Department of Agriculture Circular 118: 1-4. 
  • Petersen, CGJ (1896). "La inmigración anual de solla joven al Limfjord desde el mar alemán", Informe de la estación biológica danesa (1895), 6, 5-84.
  • Schofield, JR (2007). "Más allá de la eliminación de defectos: estimación de defectos latentes con el método de captura-recaptura", Crosstalk, agosto de 2007; 27-29.

Enlaces externos

editar