一、引言
考虑到单区制空间杜宾模型主要用于估计平均空间溢出效应，而无法捕捉到不同区制之间的非对称空间溢出效应。进一步，又考虑到传统的分区估计仅仅是单区制空间杜宾模型的重复，没有考虑到不同区制之间的非对称空间溢出效应，不具有严格的可对比性，难以准确衡量不同区制之间的空间互动关系以及溢出效应差异。因此，在传统单区制空间杜宾模型进出上扩展而来的两区制空间杜宾模型（two-regime spatial Durbin model）就应运而生了。
二、演示
我们会把《两区制空间杜宾模型的设定和估计演示案例》中的讲义、代码以及原始数据在“第八届空间计量研讨会”上与广大学员分享。需要说明的是，本代码参考了Elhorst and Fréret（2009）的做法，展示如下：
%清空命令窗口clcclear% Read data（读取数据）A=xlsread('C:\Users\user\Desktop\两区制空间杜宾模型案例\两区制空间杜宾模型回归演示\teenchee_data.xlsx');W=xlsread('C:\Users\user\Desktop\两区制空间杜宾模型案例\两区制空间杜宾模型回归演示\teenchee_swmatrix.xlsx');% Row-normalize WW=normw(W);% Model parameters and y and x variablesN=30;T=16;nobs=N*T;K=14;y=A(:,4);dum=A(:,11);xh=A(:,[5,6,7,8,9,10]);% Create wx variablesfor t=1:Tt1=1+(t-1)*N;t2=t*N;wx(t1:t2,:)= W*xh(t1:t2,:);endx=[dum xh wx];info.model=1;results = sarregime_panel(y,x,dum,W,T,info);vnames=strvcat('rcei','dum','ti','ts','gdpi','popi','eii','fdii','W*ti','W*ts','W*gdpi','W*popi','W*eii','W*fdii');prt_spreg(results,vnames,1)
图1
图2-1
图2-2
