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

MATLAB链接： Excel ： 选项——加载项—— COM 加载项——转到——没有勾选项 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 即为全 说明解释变量矩阵 x 的大小：有 nobs 行， K 列。（ size 为描述矩阵的大小） 。  1 矩阵） 附录： 静态面板空间计量经济学 一、 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 changed 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)