Perhaps this one: Modeling spatial point patterns with a marked log-Gaussian Cox process — PyMC example gallery