I wrote a case study and would be happy to get feedback before finalizing it and putting it up on our web site.
Estimating Lotka-Volterra Predator-Prey Population Dynamics with Stan
The Lotka-Volterra equations define parametric differential equations for the fluctuation of predator and prey populations. To estimate the parameters of such a model, a model for measurement error and unexplained variation is layered on top of the deterministic dynamics. The model is coded in Stan and fit to data on Canadian lynxes and snowshoe hares based on numbers of pelts collected in the early 20th century by the Hudson Bay Company.
I had to post it here in .pdf form because the site won’t let me attach .html files.
- lotka-volterra-case-study-draft.pdf (1.7 MB)
The source is on GitHub in the
example-models repo if you feel inclined to interact via pull request:
- Source repo [GitHub]