species-level covariates

33 views
Skip to first unread message

Oscar Morton

unread,
Apr 25, 2025, 6:46:03 AMApr 25
to spOccupancy and spAbundance users
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

Jeffrey Doser

unread,
Apr 25, 2025, 5:02:05 PMApr 25
to Oscar Morton, spOccupancy and spAbundance users
Hi Oscar,

You are correct that spOccupancy does not directly support including species-level covariates in the multi-species model functions (including tMsPGOcc()). However, the alternative you propose there is a completely valid option, and the function postHocLM() in spOccupancy is designed exactly for that purpose. You could take the species-specific posterior samples from tMsPGOcc() and use those as the response in a subsequent linear model, testing the impact of species-traits on those species-specific effects. postHocLM() allows you to take the full posterior distribution of those effects and propagate that uncertainty into the subsequent regression. This is perhaps a bit less elegant than directly allowing for species-level effects when fitting the model, but it gets the job done.

Kind regards,

Jeff

--
You received this message because you are subscribed to the Google Groups "spOccupancy and spAbundance users" group.
To unsubscribe from this group and stop receiving emails from it, send an email to spocc-spabund-u...@googlegroups.com.
To view this discussion visit https://groups.google.com/d/msgid/spocc-spabund-users/d8cd9923-6859-42a9-ac05-4da5c424a282n%40googlegroups.com.


--
Jeffrey W. Doser, Ph.D.
Assistant Professor
Department of Forestry and Environmental Resources
North Carolina State University
Pronouns: he/him/his

Oscar Morton

unread,
Apr 28, 2025, 3:49:23 AMApr 28
to Jeffrey Doser, spOccupancy and spAbundance users
Hi Jeff, 

Thanks for the fast reply, and sorry I missed postHocLM(), that looks like it will work great.

Many thanks,
Oscar

Reply all
Reply to author
Forward
0 new messages