Hello,

I'm running a few POPAN models; the models use age as a time-varying covariate (I have multiple initial ages). When I get to using popan.derived to calculate pop abundance by group/occasion, I get the error "Computations for mlogit parameters with covariate values cannot be specified with index column in data; use separate indices argument".

It looks like I can get the group / occasion N-hat values if I load the RMark results into R. However, I can't figure out how to properly sum N-hat estimates across groups within occasion (since I'm not sure where the vcv matrix is hiding). I also thought of exporting the models into MARK, re-running them as mcmc models, then follow this post: http://www.phidot.org/forum/viewtopic.php?f=1&t=3744 to re-export the results into R and sum them up that way. But that still seems unreasonably complicated.

Is there anywhere a worked example of how to calculate N-hat values with individual covariates?

Thank you!