Here is the code I've use for my capture history (foxCH), as well as a summary of foxCH:
- Code: Select all
> foxCH = read.capthist("2021gridcapt.txt","2021gridtrap.txt", detector = "single")
No errors found :-)
Warning messages:
1: In make.capthist(captures = capturelist[[i]], traps = trps, fmt = fmt, :
dropping repeat detections within occasions at exclusive detectors (traps)
2: In make.capthist(captures = capturelist[[i]], traps = trps, fmt = fmt, :
dropping repeat detections within occasions at exclusive detectors (traps)
3: In make.capthist(captures = capturelist[[i]], traps = trps, fmt = fmt, :
dropping repeat detections within occasions at exclusive detectors (traps)
There should be 18 detectors associated with each session (G01, G02, ... G12), but they are each showing all 198 detectors used throughout the season.
- Code: Select all
> summary(foxCH, terse = TRUE)
G01 G02 G03 G04 G05 G06 G08 G09 G10 G11 G12
Occasions 5 5 5 5 5 5 5 5 5 5 5
Detections 59 26 22 26 25 33 25 22 23 45 23
Animals 25 14 16 14 9 13 8 10 12 19 11
[b]Detectors 198 198 198 198 198 198 198 198 198 198 198[/b]
Here is a brief summary of my session capture data and associated detector data for reference:
- Code: Select all
> foxcapt = read.table("C:/Users/Hospital/Desktop/2021_Grid_Densities/2021gridcapt.txt")
> head(foxcapt)
V1 V2 V3 V4
1 G01 619591 5 G01T15
2 G01 450500 5 G01T36
3 G01 619947 4 G01T33
4 G01 605234 2 G01T13
5 G01 622471 2 G01T18
6 G01 647185 2 G01T17
> foxtrap = read.table("C:/Users/Hospital/Desktop/2021_Grid_Densities/2021gridtrap.txt")
> head(foxtrap)
V1 V2 V3
1 G01T13 352594 3653294
2 G01T14 352444 3653094
3 G01T15 352294 3652894
4 G01T16 352494 3652744
5 G01T17 352644 3652944
6 G01T18 352794 3653144
...
193 G12T31 363713 3639826
194 G12T32 363515 3639673
195 G12T33 363316 3639521
196 G12T38 363469 3639323
197 G12T39 363667 3639475
198 G12T40 363865 3639627
Is there a way I can fix my capture history to assign the 18 detector names (i.e. G01T13, or G12T31) to their associated session (i.e., G01, or G12)? The fitted secr.fit model with the current foxCH and a buffer of 800m lists all 11 unique sessions, but provides the same incorrect estimate of density for each one. If I use the same code but input capture and trap text files for one grid (i.e., Grid 1), I get a reasonable density estimate. I know I can upload 11 capture and trap text files and run separate models as they've done in previous years on this project, but I'm hoping I can fix my capture history to allow for multi-sessions.
Thank you all so much for your help!