三层：
第一层： individual events  (negative binomial, not Gaussian)
第二层： urban facilities (e.g. grocery stores, bus stops),
第三层： neighborhood conditions, such as unemployment rate and poverty rate.
目的：主要是看哪些因素影响第一层。

注：
1、第一层的individual events是negative binomial分布的，大部分是0。
2、第一层和第二层没有包含关系，都是单个spatial points. 但是第一层的分布由第二层的分布决定。

建这个multilevel spatial model，该用哪个统计方法？   是否要上Bayesian hierarchical spatial model? 还是HLM就行？

有现成的package或code吗？

多谢。