We fit our models with Stan 2.19.1 [30, 31] and R 3.6.2 [32]; Stan implements the No-U-Turn Hamiltonian Markov chain Monte Carlo algorithm [33] for Bayesian statistical inference.