Nuttallia wrote:Hi everyone,
I'm using mark recapture data collected at 5 and 7 month time intervals. I am interested in calcuating annual survival so I adjusted the time intervals from 5 and 7 months to yearly fractions (i.e. 0.417 and 0.583).
Its important to understand what tweaking the time intervals does in MARK. Suppose you had put in the number of months: 5, 7, 5, 7,...
What MARK does is calculate the Lth root of the estimate value over that parameter. So, for example, if for a 5 month interval the probability of survivng from start to end is 0.7738, then MARK would report the 5th root of this, which is 0.95. And, this value would be the estimated monthly survival rate. What MARK is doing is trying to transform the estimates suchthat they all have the same interval. The only way it can do this is to the the Lth root of the estimate.
If you had used 5, 7, 5, 7, you could have taken your estimates of monthly survival for each part of each year, and then used the products of them to determine the annual estimate.
Alternatively, if you use 0.41667 and 0.5833 as you did (MARK handles that just fine), instead of monthly estimates, you'd get estimates of annual survival - but (and this is where you need to be careful), you'd get annual estimates based on the survival over the two different parts of the season, and they mght not be the same (i.e., MARK would estimate annual survival based on first 5 months, then another estimate of annual survival for next 7 months, and so on...).
If you want annual survival, then all you really need to do is take that data from one particular point of the year, and build your encounter history that way. Estimation of survival is not dependent on when you start the calendar (e.g., Jan -> Jan estimates will be the same as April -> April estimates). But, you have data at two times of the year. Having data at two discrete times of the year lets you partition annual survival among seasons, which is potentially much more interesting than grinding out annual survival. I'd suggest using 5, 7, 5, 7..., and then look at whether or not there is evidence of seasonal variation in estimated monthly survival. This is straightforward. Imagine 3 years of data (real years, sampled twice each year). A time dependent PIM would look like
- Code: Select all
1 2 3 4 5 6
2 3 4 5 6
3 4 5 6
4 5 6
5 6
6
The key is remembering that 1, 3, 5 form a set (say, first 5 months of each of the three years), and that 2, 4, 6 for another set (say, second 7 months of each of the three years). If you didn't think there was annual vairation in monthly survival, but that monthly survival varied between seasons (within years) you would compare
- Code: Select all
1 2 1 2 1 2
2 1 2 1 2
1 2 1 2
2 1 2
1 2
2
with
- Code: Select all
1 1 1 1 1 1
1 1 1 1 1
1 1 1 1
1 1 1
1
1
The former PIM allows for seasonal variation in monthly survival, whereas the latter PIM constrains them to be the same.
This sort of approach, with some care, will let you partition annual survival into seasonal components, which is invariably a good thing.