Fit a spatial occupancy model via Gibbs sampling