当前，采用空间模型比较普遍，但现有的模型有很多缺陷，主要有二：（1）矩阵式0-1变量矩阵；（2）个体间的影响是线性的、对称的。
这两点导致现有文献的结论相互冲突，备受争议。
直观地说，第一，0-1变量会导致，比如，上海和广州，或上海和北京，或北京和广州之间是没有影响的，但事实上，他们之间是有相互影响的，这是现有模型无法做到的，不论是以是否相邻来去0或1值，或者以个体（城市）间的距离来去0或1值；
第二，线性的、对称的相互影响导致的问题更加严重。比如，上海和嘉兴（杭州的一个地级城市，靠近上海，位于上海和杭州之间），现有的模型假定上海对嘉兴的影响，等于嘉兴对上海的影响，这显然是不可能的；而且，现有模型假定上海对嘉兴的影响在样本期内是保持不变的，这也是不符合客观事实的。
这两点是严重与Lucas批评相违背的，因此必须对模型和计量方法进行调整。
特别声明：以上的问题主要发生在国内文献中，国外英文文献（近期）对基础模型进行了修改，比如AER（应该是2015第8期）允许是连续型的的矩阵。
但是现在国内很多学者之间套用当前软件提供的最基本的模型和计量方法，直接放入数据，run......

不过，要将threshold effect和TV effect放入spatial model中一起考虑，难度是有的，比如用国内地级市数据（至少260个左右城市），那矩阵就是260×260×T（T表示样本期时间跨度，允许非对称效应存在（aij不等于aji），故矩阵元素（aij）不除2），这个数据大家可以想象一下其庞大性。

因此，stata恐怕有点勉为其难了，需要用R。当然，用R来写也更方便一点。
R比stata快，除了两者数据库结构不一致外，stata只能用ado来完成，这与官方内置的command是有很大差距的，ado的command肯定要慢于官方内置command，当数据少时，感觉不到，但大数据时，立马显现了。

grid with bs方法需要用
再次之前linearity test方法与PSTR，PTR类似（先做Taylor展开，取1阶，最多取3阶（理由：参见现有文献））
大概思路就是这样。

自己动手，丰衣足食！