“…Each parameter had a uniform prior distribution as follows: frequency of introduction ϕ ~Unif(0,60), β d ~Unif(0.0001,1), β c ~Unif(0.0001,0.99), ρ h ~Unif(0.0005,0.1), ρ c ~Unif(0.0005,0.8), π ~Unif(0.1,1.5), θ ~Unif(0.5,6), λ ~Unif(0.1, 2.5), β wd ~Unif(0.01,1), and β wc ~Unif(0.001,1). Prior distribution ranges were informed by movement and contact data (Kay et al, ; Pepin et al, ; Podgórski et al, , ). As part of the parameter generation process we implemented the following constraints for each parameter set: β d > β c , β wd > β d , β wc > β c ; to further inform prior distributions with biologically realistic knowledge.…”