第一步：
areg y x1 x2....., absorb(city)
city量很大，所以不适合用reg ....i.city

第二步：
画图：各个city 的b_[x1]（纵坐标）对 x3（横坐标）画图

非常感谢！