六步学会：用MATLAB做空间计量回归详细步骤（每一步都有截图）.docx

1.excel与MATLAB链接：Excel：选项——加载项——COM加载项——转到——没有勾选项2.MATLAB安装目录中寻找toolbox——exlink——点击,启用宏E:\MATLAB\toolbox\exlink然后，Excel中就出现MATLAB工具（注意Excel中的数据：）3.启动matlab点击start MATLABsenddata to matlab ，并对变量矩阵变量进行命名（注意：选取变量为数值，不包括各变量）(data表中数据进行命名)（空间权重进行命名）导入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=results.sige*((nobs-K)/nobs);loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid*results.resid% The (robust)LM tests developed by ElhorstLMsarsem_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=results.sige*((nobs-K)/nobs);loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid*results.resid% The (robust)LM tests developed by ElhorstLMsarsem_panel(results,W,y,[xconstant x]); % (Robust) LM tests2、空间固定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 changedprt_reg(results,vnames);sfe=meanny-meannx*results.beta; % including the constant termyme = y - mean(y);et=ones(T,1);error=y-kron(et,sfe)-x*results.beta;rsqr1 = error*error;rsqr2 = yme*yme;FE_rsqr2 = 1.0 - rsqr1/rsqr2 % r-squared including fixed effectssige=results.sige*((nobs-K)/nobs);logliksfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*resu