Hasta ahora hemos utilizado unos ejemplos que vienen pre-formateada en una trama de datos de las historias de captura. Por lo general, sin embargo, nuestros datos de campo van a estar en un formato menos estructurado, por lo general una serie de días u otras ocasiones de encuentro y una lista de los animales (por número de identificación u otra etiqueta) que se capturan cada ocasión. Así que esta estructura típica es algo así como
Ocasión Id
1 7
1 18
1 22
....
2 6
2 18
2 3
2 22
...
etc
Lo que queremos en cambio, es una matriz de las historias de captura, donde las filas son el animal indiviual (pueden o no ser necesarias las etiquetas de identificación), y las columnas son las entradas de 1 indica que el animal fue capturado en esa ocasión, o 0, no . Por ejemplo, para 5 ocasiones algo así como
10100
00100
10000
etc
Por definición, por supuesto, nunca vamos a observar
00000
que son animales que nunca fueron capturados.
Hay un número de maneras de lograr este incluyendo
La construcción de la matriz con la mano (sólo se recomienda para problemas muy pequeños)
La construcción de la matriz usando Excel o Access tablas dinámicas (no es mi departamento)
R u otro programa para hacer el trabajo (enfoque aquí)
He escrito algo de código en R que lee datos de un archivo csv en ocasiones / formato de id, la convierte en una matriz de historias de captura usando las funciones en la biblioteca R reshape y vuelva a ejecutar una serie de modelos de captura cerrados con los datos formateados como entrada . Los datos se encuentran aquí y el programa ya está aquí. El formato de los datos se compone de 82 historias de captura de animales capturados al menos una vez en 5 ocasiones. Los datos se leen por la función de mark en la biblioteca de RMark y se ejecutan a través de una serie de modelos cerrados. La estimación de la estimación del modelo promediada de unos 106 y p ~ 0.25 indican que hubo más de 20 animales de la población que nunca fueron capturados, por lo que la estimación ingenua de 82 sería muy baja.
Siguiente: Introducción a MARK y RMark