Matthijs Hollanders

Matthijs Hollanders

@matthollanders.bsky.social

Quantitative ecology, statistics, wildlife, field herping. Post-doctoral Research Fellow @ U of Canberra Consultant @ Quantecol, https://quantecol.com.au Wildlife tour guide @ Australian Wildlife Encounters, https://www.australianwildlifeencounters.com

508 followers271 following253 posts

Matthijs Hollanders reposted

View full thread

The ecological parameters are then recovered as: logit(psi) = alpha[1] + lambda[1] * eta log(gamma) = alpha[2] + lambda[2] * eta log(epsilon) = alpha[3] + lambda[3] * eta Maybe not everyone's cup of tea, but drastically reduces model complexity; nice for multispecies models! 5/n

View full thread

To avoid linear predictors on initial occupancy and both transition rates (which are usually highly correlated), occARU uses a structural equation modeling (SEM) approach using a latent variable `eta` (by site, season, and species), representing "site suitability" for predictors + random effect. 4/n

Constructing the transitions in continuous time, P = matrix_exp(Q * tau), accounts for potentially unequal time intervals between seasons. We would not expect occupancy status to be likely to change between short intervals, whereas longer gaps increase the likelihood. 2/n

Plotting functions have been updated with `by_region` arguments to plot i.e. intercepts and realised by region (and species), and an additional reduce_sum() function has been added that does within-chain parallisation over regions instead of sites. Now GPs can be computed in parallel. 2/2

View full thread

occARU accommodates site-level predictors for occupancy and detection and additionally site-by-survey level predictors for detection rates. Continuous, categorical, and ordinal predictors are currently supported. Model complexity is handled via global-local shrinkage priors. 4/n

occARU uses hierarchical multispecies Gaussian processes to infer spatial and temporal variation in detection rates, and estimates interspecific correlations in these processes. These effects are orthogonally projected to retain inference on fixed effects. 2/n

occARU takes in standardised CamtrapDB data by default, optionally thins detections with a user-specific value, and aggregates observations to surveys of arbitrary lengths. This balances temporal resolution with computational demands. 3/n

New pre-print! I cover a range of open capture-recapture models (single survey/robust design, multistate/multievent, in (Cormack-)Jolly-Seber variants) in Stan, and provide efficient log likelihood functions. I also introduce a method to account for unequal survey intervals in the entry process.

MC Stan@mc-stan.org

"Efficient Bayesian implementations of capture-recapture models with Stan" by @matthollanders.bsky.social discourse.mc-stan.org/t/bayesian-c...

In (Gaussian) linear models, we usually put normal priors for regression coefficients. How do we feel about putting gamma(a, a) priors in something like Poisson regression? Using rates for rates seems nice; the geometric means are still 1, meaning the priors are centered on no multiplicative effects

Hazard rates FTW 🎉

Quantecol@quantecol.bsky.social

By default, mark-recapture models parameterise the survival process with survival probabilities. This is unfortunate, because logit-linear survival functions aren't time invariant, and you can't just exponentiate multistate transition probabilities like you can survival. Instead, we should... 1/2

brms is one of the most amazing R packages ever. Nevertheless I can't get myself to use it for my own work because nothing beats coding up your own models where you have ultimate control and insight into what's going on. There remains a "black box" element unless you are deeply familiar with it.