由于样本量有点多，直接上附件了。代码如下附注：region == 1 的情况和region == 0 回归的观察值相加，不等于excel表里的全样本1612个，同时不region分组的观测值，logit回归以后也会丢失两个观测值。

因为我用excel筛选已经剔除各变量的NA了，不知道这两个（或者分组情况下更多的）样本是咋丢的