including demographic stochasticity?

[Aariq]

I'm using IPMs to do population viability analysis for the first time (I've done PVA with individual-based models before) and struggling a bit to understand how different sources of stochasticity can be included. I'm using a stochastic IPM with parameter selection and random environmental states for each year. The pop_state vector will be different for every run, but this is entirely due to environmental stochasticity. Is there a way to include demographic stochasticity in an IPM with ipmr?

I saw the suggestion in Ellner, Childs, and Rees section 10.3 to do something like rpois(length(nt), K %*% nt), but I'm not sure if this is appropriate or something that can be implemented in ipmr.