NBlat function

Fit Negative Binomial Spatial Model