《六步学会用MATLAB做空间计量回归详细步骤》由会员分享，可在线阅读，更多相关《六步学会用MATLAB做空间计量回归详细步骤（24页珍藏版）》请在金锄头文库上搜索。
1、1. excel 与 MATLAB 链接:Excel：选项一一加载项一一OM加载项转到一一没有勾选项2. MATLAB安装目录中寻找toolboxexlink点击，启用宏E:MATLABtoolboxexlinkT悸取口朗固 Contents2013/12/2 & 2:14MATLAB Code1 KBglT eKclIink2013/10/9 23:48Microsoft Excel .108 KB阳 acdlink2DD32013/10/9 23:48IMiizosoft Excel 305 KBS ExliSamp2009/11/5 19:58Microsoft Excel 72 KB餐
2、| getfunctiDnhelp2010/12/4 2:17MATLAB Code2 KB擋 getfunctionlist2010/12/4 2:17MATLAE Code3 KB餐| getfunctionisign呂tuire*2012/5/4 2:59MATLAB Code4 KE圏 info2007/7/24 0:22XML注1 KB然后，Excel中就出现MATLAB工具注意 Excel 中的数据：）变量数扌dataJ vj1 权董数据,Sheet33.启动 matlab（1） 点击 start MATLAB（ 2） senddata to matlab ，并对变量矩阵变量进行命
3、名（注意：选取变量为数值，不包括各变量）5 W5 七咖* Hl晰f 1W2W2 ii i$?tni iurns11 imi?i -fli&as -4)238?a述is 羽风审 刃2zhiJ-2rrir?1 I iJUMZ -O unn -gjPMI J.MM1 u i raiMs -a a 乜iFH laiisa13刃2瞬H 叩JZH2 4恥1II 1 r-5!?l -OQIW -OJZMI 10CI3IZIV IMPIJI 乜 1 萍M -AJM3 i 11TMII 1 M3r42 -Ql!5!fll -C!2M2I 丄 XZ&ffl i 1! IMTZaS Hllftll -dJMi id
4、OZZd 21 1M1M2 -OflHKZ -0JZH3 4JM1321 1T2W2 也网勺叩记吋a 4-21HTS ZZ 142IK2E -0 01 -OJEfldi 鼻声如23 1W1M3 -0 1TM -0JJP3 lK5iaBV#laWe wme旨 汕哪1严=-0jjM Ma皿寸ESvii49 曲型E翼山:J砂*Ihb Ihpi ImwOD贷暑OISMf汕加辿 i 14183( 临畑1 a 542232i54)21712七 121i rv2 -1.132(-4JW13.17W3i556i?2ti & 11(12 Hi.im HBim 丄2币: ? b 诧il -42611 -I2M2
5、丄网22M I iHSai -fl if.iai 七 i卿 丄训喇（data 表中数据进行命名）MB4U4.UU/JJJU.iykbU4U.lUUkbf4.555dbJlk!4.529414.8402420.1801810 1585194.775782579.2452.129424.6839810.118179-0.03264.8648996364.14965.129434.3956830.3505260 1754084.8564535162.23960.229444.5849570.120077-0.03254.5685861862.11465.429464.6664290.2295380.
6、2408064.7582615017.33853.829464.7077270.1002870 1313914.711692483.3360.5空间权重进行命名)(3)导入MATLAB中的两个矩阵变量就可以看见S工柞区lSOxS double46x4& doubfeg O O 3打开腔本*. 凸打-工作区因;清除工忙区-./清除命令*立件A ESimulink库SIMULIMK帀局-*iii ri.TLAB R20148 _疼錘稈序4.将 elhorst和jplv7两个程序文件夹复制到MATLAB 安装 目录的toolbox文件夹=窝6寸I二E2器嘶寸H/V卜 M占 350 寸鹽6 寸T/TA
7、norvl 窝5寸5AU2 雹帯寸T3E0E BN mUHtE 寸輕6寸T/TM12 gm mrZKFE &6 寸T7TE2 垃心- odE -r uqE -xqE-4BE对nJmoorEEbmdqe-lew T q4EEdIBE -圣一召卫J d-u睡艸匸匸噩桁肾6 - Si3 巴噩j ca泉 3 s 肩 lim+I 駅超脚空 S 诚KleflEhl遇 瑕0 比画L10亠 xoqo4 SVILXE 4也型翼哩？吐蛍押-ss.g0 GmEE t.s冴 ss 岂&I寸T-I/VZJ 乂 一巴。上帀 电N=46;W=normw(W1);y=A(:,3);x=A(:,4,6);xconstant=o
8、nes(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,xco
9、nstant 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;
10、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*resu lts.resid% The (robust)LM tests deve
11、loped by Elhorst LMsarsem_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
12、(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
13、effects sige=results.sige*(nobs-K)/nobs);logliksfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid*resu lts.residLMsarsem_panel(results,W,ywith,xwith); % (Robust) LM tests3、时期固定OLS (time-period 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=2;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);tfe=meanty-meantx*results.beta; % including the constant termyme = y - mean(y);en=ones(N,1);
