步学会用MATLAB做空间计量回归详细步骤.docx

与MATLAB链接：

Excel ：

选项——加载项——COM 加载项——转到——没有勾选项

2.MATLAB安装目录中寻找 toolbox —— exlink ——点击 , 启用宏

E:\MATLAB\toolbox\exlink

然后， Excel 中就出现MATLAB 工具

（注意 Excel 中的数据：）

3. 启动 matlab

（1）点击 start MATLAB

（2 ）senddata to matlab，并对变量矩阵变量进行命名（注意：选取变量为数值，不包括各变量）

(data 表中数据进行命名)

（空间权重进行命名）

（3 ）导入MATLAB中的两个矩阵变量就可以看见

4.将 elhorst 和 jplv7 两个程序文件夹复制到MATLAB 安装目录的 toolbox文件夹

5.设置路径：

6.输入程序，得出结果

T=30;

N=46;

W=normw(W1);

y=A(:,3);

x=A(:,[4,6]);

xconstant=ones(N*T,1);

[nobs K]=size(x);

results=ols(y,[xconstant x]);

vnames=strvcat('logcit','intercept','logp','logy');

prt_reg(results,vnames,1);

sige=*((nobs-K)/nobs);

loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*'*

% The (robust)LM tests developed by Elhorst

LMsarsem_panel(results,W,y,[xconstant x]); % (Robust) LM tests 解释

每一行分别表示：

该面板数据的时期数为30（ T=30 ），

该面板数据有30 个地区（ N=30 ），

将空间权重矩阵标准化（W=normw(w1)），

将名为 A（以矩阵形式出现在MATLABA 中）的变量的第 3 列数据定义为被解释变量y,

将名为 A 的变量的第4、 5、 6 列数据定义为解释变量矩阵x，

定义一个有N*T 行， 1 列的全 1 矩阵，该矩阵名为：xconstant ，（ ones 即为全

1 矩阵）

说明解释变量矩阵x 的大小：有nobs 行， K 列。（ size 为描述矩阵的大小）。

附录：

静态面板空间计量经济学

一、 OLS 静态面板编程

1、普通面板编程

T=30;

N=46;

W=normw(W1);

y=A(:,3);

x=A(:,[4,6]);

xconstant=ones(N*T,1);

[nobs K]=size(x);

results=ols(y,[xconstant x]);

vnames=strvcat( 'logcit' ,'intercept' ,'logp' ,'logy' ); prt_reg(results,vnames,1);

sige=*((nobs-K)/nobs);

loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*'*

% The (robust)LM tests developed by Elhorst

LMsarsem_panel(results,W,y,[xconstant x]); % (Robust) LM tests

2、空间固定 OLS （spatial-fixed effects ）

T=30;

N=46;

W=normw(W1);

y=A(:,3);

x=A(:,[4,6]);

xconstant=ones(N*T,1);

[nobs K]=size(x);

model=1;

[ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,x,N,T,model );

results=ols(ywith,xwith);

vnames=strvcat('logcit','logp','logy'); % should be changed if x is

prt_reg(results,vnames);

sfe=meanny-meannx*; % including the constant term

yme = y - mean(y);

et=ones(T,1);

error=y-kron(et,sfe)-x*;

rsqr1 = error'*error;

rsqr2 = yme'*yme;

FE_rsqr2 =- rsqr1/rsqr2 % r-squared including fixed effects sige=*((nobs-K)/nobs);

logliksfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*'*

LMsarsem_panel(results,W,ywith,xwith); % (Robust) LM tests

3、时期固定 OLS（time-period fixed effects ）

T=30;

N=46;

W=normw(W1);