Hi all,
I'm working on fitting a multi-species, multi-season occupancy model. One of the things I'd like to be able to incorporate is if average occurrence probabilities vary through seasons by some trait intrinsic to species e.g. species habitat preference, to evaluate whether more broad species groups share a common temporal trend.
I can see in tMsPGOcc() the occurrence covariates can vary by site or site/primary time period level, and the detection covariates can vary the same or additionally at the observation level. Does that mean it's not currently possible to include species-level covariates (traits etc.) in the occurrence model? Or have I missed something?
If that is the case, would an alternative be to extract this from the posterior post-fitting, e.g. if I have a trait with two levels, I could feasibly compare the average species-level trends for species of each level or the proportion with x-directional temporal change per trait level or the like.
Thank you in advance for any guidance or support!
All the best,
Oscar