关于运用matlab运算产业集聚测度方法中的DO指数

DO指数应该只需要用到企业经纬度坐标即可，想请问大家帮忙梳理DO指数matlab复现具体步骤，代码思路，尤其是全局及局部置信区间的蒙特卡洛模拟等，感谢大家！