7月2号的原创推文推送以来,我的公众号后台有不少留言。其中有关Stata的问题我已经给出专门回答,有读者留言,希望我推一期使用MATLAB做空间面板的推送。
本期以我国31个省级行政单位(除港、澳、台)2001-2010年的空间面板数据为例,所用数据及jplv7包文件下载方法请见文末留言区。
加载分析所须的jplv7包
找到MATLAB安装路径中的toolbox文件夹,例如我的电脑是C:\Program Files\MATLAB\R2012a\toolbox,把jplv7文件夹复制进去,然后在MATLAB中点击file-set path-add with subfolders,选择jplv7以设置路径:
执行如下程序代码
% 读取数据
[data,name]=xlsread('data',1);
% 读取权重矩阵
weight=xlsread('data',2);
% 提取变量
y=data(:,[4]);
x=data(:,[5,6,7,8,9,10,11,12,13]);
% 标准化权重矩阵
W=normw(weight);
% 设置面板
info.model = 1;
% SAR模型
results_sar=sar_panel_FE(y,x,W,T,info);
% SEM模型
results_sem=sem_panel_FE(y,x,W,T,info);
vnames=strvcat('house','urbanization','student','patent','foreigninvest','industry','research','fixedassets','lnpopulation','lngdp');
prt_spnew(results_sar,vnames,1)
prt_spnew(results_sem,vnames,1)
LMsarsem_panel(results_sar,W,y,x)
LMsarsem_panel(results_sem,W,y,x)
% SDM模型
t1=(t-1)*N+1;t2=t*N;
wx(t1:t2,:)=W*x(t1:t2,:);
results_sdm=sar_panel_FE(y,[x,wx],W,T,info);
vnames=strvcat('house','urbanization','student','patent','foreigninvest','industry','research','fixedassets','lnpopulation','lngdp',...
'W*urbanization','W*student','W*patent','W*foreigninvest','W*industry','W*research','W*fixedassets','W*lnpopulation','W*lngdp');
prt_spnew(results_sdm,vnames,1)
结果输出: