Just looking for some advice on running a dead recovery model...
Some background info, I have 7991 tagged and released individuals with 727 captures over three years.
I binned the data into even weekly time intervals (nocc=149) but am now interested in getting an annual estimate of survival.
- Code: Select all
time.intervalsweek<-c(rep(0.019230769, 149))
wall17.processed=process.data(data,model="Brownie", time.intervals = time.intervalsweek)
wall17.ddl=make.design.data(wall17.processed)
wall17.ddl$S$Year=factor(c(rep("spring74Spring75", 52), rep ("spring75spring76", 52), rep("spring76march3177", 45)))
S.year=list(formula=~Year)
f.dot=list(formula=~1)
model.list1=create.model.list("Brownie")
I altered the time interval to be equal to one year (1/52) and binned the data into the yearly groupings and ran the brownie dead recovery model
Im getting annual survival of roughly 7% (one survival interval is at boundary) which seems unrealistically small.
Could this possibly due to the small initial tags released and <10% recovery, or am I interpreting something wrong with the time.intervals code set up??
Thanks in advance!