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.​ 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,​ 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.

Caracol de color ámbar ovado de Chittenango marcado.

Otra aplicación importante de estos métodos es la epidemiología,​ 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.).​

Trabajo de campo relacionado con marca-recaptura

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.​

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

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.​ 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.​

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.​

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

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

• 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.​

Notación

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.

$N=n*K/k$

Estimador Lincoln-Petersen

El método Lincoln-Petersen​ (también conocido como índice Petersen-Lincoln​ 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:

${\hat {N}}={\frac {Kn}{k}},$

Derivación

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).​

Esto implica que, en la segunda muestra, la proporción de individuos marcados que son capturados ($k/K$) debe ser igual a la proporción de la población total marcada ($n/N$). 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,

${\frac {k}{K}}={\frac {n}{N}}.$

Una reordenación de esto da

${\hat {N}}={\frac {Kn}{k}},$

la fórmula utilizada para el método Lincoln-Petersen.​

Cálculo de muestra

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

${\hat {N}}={\frac {Kn}{k}}={\frac {10\times 15}{5}}=30$

Estimador de Chapman

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.​ El estimador de Chapman proporciona un estimador alternativo menos sesgado del tamaño de la población:

${\hat {N}}_{C}={\frac {(K+1)(n+1)}{k+1}}-1$

Cálculo de muestra

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

${\hat {N}}_{C}={\frac {(K+1)(n+1)}{k+1}}-1={\frac {11\times 16}{6}}-1=28.3$

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."​

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

Un aproximado $100(1-\alpha )\%$, el intervalo de confianza para el tamaño de la población N se puede obtener como:$K+n-k-0.5+{\frac {(K-k+0.5)(n-k+0.5)}{(k+0.5)}}\exp(\pm z_{\alpha /2}{\hat {\sigma }}_{0.5})$ ,

donde ${\textstyle z_{\alpha /2}}$ corresponde a la $1-\alpha /2$ cuantil de una variable aleatoria normal estándar, y

${\hat {\sigma }}_{0.5}={\sqrt {{\frac {1}{k+0.5}}+{\frac {1}{K-k+0.5}}+{\frac {1}{n-k+0.5}}+{\frac {k+0.5}{(n-k+0.5)(K-k+0.5)}}}}$ .

Se ha demostrado que este intervalo de confianza tiene probabilidades de cobertura reales cercanas a la nominal. $100(1-\alpha )\%$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.​

Estimación bayesiana

El valor medio ± desviación estándar es

$N\approx \mu \pm {\sqrt {\mu \epsilon }}$

donde

$\mu ={\frac {(K-1)(n-1)}{k-2}}$ para $k>2$

$\epsilon ={\frac {(K-k+1)(n-k+1)}{(k-2)(k-3)}}$ para $k>3$

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

Probabilidad de captura

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,​ y se ha utilizado tanto en ecología como en epidemiología para detectar enfermedades animales o humanas,​ 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.​ 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).​ 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:

$P={\frac {5}{10}}fq+{\frac {5}{90}}(1-f)(1-q),$

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 :

$P=\left({\frac {5}{10}}q-{\frac {5}{90}}(1-q)\right)f+{\frac {5}{90}}(1-q).$

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:

$\left({\frac {5}{10}}q-{\frac {5}{90}}(1-q)\right)>0,$

que se simplifica a:

$q>{\frac {1}{10}}.$

Este es un ejemplo de optimización lineal.​ 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

Existen modelos estadísticos muy elaborados disponibles para el análisis de estos experimentos.​ 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)",​ "Modelos Loglineales para Experimentos de Captura-Recaptura (Rcapture)",​ y "Muestreo de Distancia de Marca-Recaptura (mrds)".​ Estos modelos también pueden adaptarse a programas especializados como MARK​ o M-SURGE.​

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​ (descritos anteriormente como una expansión del método Lincoln-Petersen para poblaciones cerradas).​

Enfoques integrados

El modelado de datos de marcado-recaptura tiende hacia un enfoque más integrador,​ 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.​

Véase también

• Problema del tanque alemán, para la estimación del tamaño de la población cuando los elementos están numerados.
• Etiquetar y liberar
• Estimación de abundancia
• Seguimiento de vida silvestre por GPS

Bibliografía

• Besbeas, P; Freeman, S. N.; Morgan, B. J. T.; Catchpole, E. A. (2002). «Integrating mark-recapture-recovery and census data to estimate animal abundance and demographic parameters.». Biometrics 58 (3): 540-547. PMID 12229988. doi:10.1111/j.0006-341X.2002.00540.x. 
• Martin-Löf, P. (1961). «Mortality rate calculations on ringed birds with special reference to the Dunlin Calidris alpina». Arkiv för Zoologi (Zoology Files), Kungliga Svenska Vetenskapsakademien (The Royal Swedish Academy of Sciences) Serie 2. Band 13 (21). 
• Maunder, M. N. (2004). «Population viability analysis, based on combining integrated, Bayesian, and hierarchical analyses». Acta Oecologica 26 (2): 85-94. Bibcode:2004AcO....26...85M. doi:10.1016/j.actao.2003.11.008. 
• Phillips, C. A.; M. J. Dreslik; J. R. Johnson; J. E. Petzing (2001). «Application of population estimation to pond breeding salamanders». Transactions of the Illinois Academy of Science 94 (2): 111-118. 
• Royle, J. A.; R. M. Dorazio (2008). Hierarchical Modeling and Inference in Ecology. Elsevier. ISBN 978-1-930665-55-2. 
• Seber, G.A.F. (2002). The Estimation of Animal Abundance and Related Parameters. Caldwel, New Jersey: Blackburn Press. ISBN 1-930665-55-5. 
• Schaub, M; Gimenez, O.; Sierro, A.; Arlettaz, R (2007). «Use of Integrated Modeling to Enhance Estimates of Population Dynamics Obtained from Limited Data». Conservation Biology 21 (4): 945-955. PMID 17650245. doi:10.1111/j.1523-1739.2007.00743.x. 
• Williams, B. K.; J. D. Nichols; M. J. Conroy (2002). Analysis and Management of Animal Populations. San Diego, California: Academic Press. ISBN 0-12-754406-2. 
• Chao, A; Tsay, P. K.; Lin, S. H.; Shau, W. Y.; Chao, D. Y. (2001). «The applications of capture-recapture models to epidemiological data». Statistics in Medicine 20 (20): 3123-3157. PMID 11590637. doi:10.1002/sim.996. 

Otras lecturas

• 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

• Una introducción histórica a los métodos de captura-recaptura
• Análisis de datos de captura-recaptura

• Datos: Q796147