I have a combination of live resightings and dead recoveries. I have 4 groups and defined 2 states (live and dead).
Individuals in all groups can be resighted alive, but only members of groups 1 & 2 can be recovered dead.
There is only 1 age class. I'm using a CAS model with separate survival/transition parameters in M-Surge.
I coded survival as 'f' and interpret transition as survival (I fixed the transition from(dead)to(live) to 0). I changed the transitions pattern matrix to
|1 0|
|1 0|
I want to run a model with additive effects of group and time (no interaction).
I coded p in Gemaco as
[g+t].to(1)+[g(1,2)+t].to(2)
I get separate estimates of p for live resightings for the 4 groups that vary in parallel over time, as expected.
I expected to only get estimates of p for dead recoveries for groups 1 & 2. However, I also get estimates for group 3 in the reduced set of parameters. Occasions 1-11 are 0, and I fixed the last occasion (13) in both p and psi, but occasion 12 is non-zero.
How should I code p to get dead recovery probabilities just for groups 1 & 2?
Thanks and merci beaucoup in advance for any advice.
Ginger