We introduce a flexible and scalable class of Bayesian geostatistical models for discrete data, based on nearest-neighbor mixture processes (NNMP), referred to as NNMP. To define the joint probability mass function (pmf) over set spatial locations, we build from local mixtures conditional pmfs using directed graphical model, with acyclic graph that summarizes nearest neighbor structure. The app...