步學(xué)會(huì)用MATLAB做空間計(jì)量回歸詳細(xì)步驟.docx

1、與MATLAB鏈接:Excel :選項(xiàng)一一加載項(xiàng)一一 COM加載項(xiàng)一一轉(zhuǎn)到一一沒(méi)有勾選項(xiàng)2. MATLAB安裝目錄中尋找toolboxexlink點(diǎn)擊，啟用宏E:MATLABtoolboxexlink務(wù)次口朗Contents2013/12/28 2d4MATLAB Code1 KBQL excllink2013/10/9 23:48WtcfOMil Excel 108 KBC| exdlink2OO32013/10/9 29：朝Microsoft Excel 屮305 KBB' ExliSamp2009/11/5 19;58MJcrosoft Excel .72 KBgetfuncti

2、onheip2010/12/4 2:17MAT LA冬 Code2 K8' getfunctionlist2010/12/4 2tl7MATLAB Code3 KB松getfLinctionsiriaiure2012/5/4 2:59N?ATLAB Code4 Kfi邕 info2007/7/24 0:22XML1 KB然后，Excel屮就出現(xiàn) MATLAB工具FCQ循入CXL壬蘭動(dòng)環(huán)印塞靜筆元格zTATLABIB(注意Excel中的數(shù)據(jù)：)zy4U4 uurjjjU 詫bU4U1UUZ0/MT4 11Z14 b29414 8402420 1801810 1585194 7757825

3、792452 129424.6839810.118179-0.03264.8548996354.14966.129434.3956830.3505260.1754084 86645351S223950.229444 5649670 120077-0 03264 5685861862 11465429454 56642902295380 2408064 7582615017 33853.829464 7077270.10028701313914 711692483.3L 360 5I權(quán):重?cái)?shù)據(jù)1 5heet3 j3啟動(dòng) matlab(1 ) 點(diǎn)擊 start MATLAB(2 ) senddat

4、a to matlab ，并對(duì)變量矩陣變量進(jìn)行命名(注意：選取變量為數(shù)值，不包括各變量)(HMk1woe>t0£nuonMP畑巧141 MH41叫紳SB?2W fQ2< ttfcuW跡MMl034 «7»)1fl 17M?)sw:W'T?H30d42i»U1WTWI叫懶席Q!曲如H.S617：F*iMHaf1 g就4 15214 15211 Jtt2lV血surV?«翊4 Zt4.4%如Q1皿軸剛5：出如: /jmth0t如T創(chuàng)町詢1* 40it4 6UM14的越血竹XMfeJS*對(duì)讓；Qu11小蝴4 ai)M<KM巧

5、0124T3JJ"-ozjm:-azwn丄2?跖廣4T5sw巧U01)4 WMK*4砂IfMimr/Q14i e*Bw4J17W?4JY7W：A?51M：'叫QuIJWM記 JJin陽(yáng)F1W r,i?a114 WM-qffiM*AW3IZ1MW0it1 WMI可l»U砒則4. H Z*Cfl* f 'giiimvaz4KNB1M4 V5W»WM砂E0»4 95 mi4 WJI砒3W!yzzrM：EH !y*0JI1如82FXM：磅理啓d »MUMM!M2：2?Ji4 JZlMrHFcrur4J u-nuA71ZW1如Z3Ufr0

6、22QWF4 UAH1WW皆323丄 fl9 «M2-flFT處5EiS«!M?wrAj4 r?MflTrtMiV.41 tv.udmiun e亠NUMJUK* i <（data表屮數(shù)據(jù)進(jìn)行命名）（空間權(quán)重進(jìn)行命名）（3 ） 導(dǎo)入MATLAB中的兩個(gè)矩陣變量就可以看見(jiàn)八 MATLAB R2*.亠. 主頁(yè)*應(yīng)馬稈序由資逹轄運(yùn)mott的二囚打開(kāi)麗C 'Simulink齊信我逹新逹導(dǎo)入4將elhorst和jplv7兩個(gè)程序文件夾復(fù)制到MATLAB安裝目錄的toolbox文件夾5- elhorst2017/1/W 1908文樣與心 jplv72017/1/15 &am

7、p;46女忖5 設(shè)置路徑:唱K甌址1_密密縝再文件裘渤國(guó)澤徑, v'卞 程序(£:) > MATLAB toolbox 勺；攬箕P1 el4A詠*修8田期F"云.:上訃/2017/1/16 12:19-jmsyil local2017/1/14 9 24/ Ire2017/1/14 9 22文舛滅臚ap2G17/1/14 Q 21matlab2017/1/1 9 24文枠rMthb prr>gram200-?G147017/1/15 259丈<4JB洌 matlabxl2017/1/1423如喪J J®b t 事再 If、mbc2017/

8、V14 9 25丈<4仏 iBS#【G)e怕pc2017/1/14 24a專業(yè)(氏nnet2017/1/14 925穩(wěn)序(&)A op2O17/V1 925d資科Vh)注 oo:irr?017/14 925文除J i-«1 CD敗祈谿M< 1現(xiàn) 文件舟jpl"6 輸入程序，得出結(jié)果4 C5一 . T 丁 c徐宀 r . D 1 d.tac - 門 - Q C .f "cwsT=30;N=46;W=normw(W1);y=A(:,3);x=A(:,4,6);xconstant=ones(N*T,1);nobs K=size(x);results=

9、ols(y,xconstant x);vnames=strvcat(logcitVi nterceptTlogp','logy)；prt_reg(results,vnames,1);sige=*(nobs-K )/no bs);loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*1*% The (robust)LM tests developed by ElhorstLMsarsem_panel(results,W5y,xconstant x); % (Robust) LM tests解釋每一行分別表示：該面板數(shù)據(jù)的時(shí)期數(shù)為30 ( T=30

10、 ),該面板數(shù)據(jù)有30個(gè)地區(qū)(N=30 ),將空間權(quán)重矩陣標(biāo)準(zhǔn)化(W=normw(w1),將名為A (以矩陣形式出現(xiàn)在 MATLABA中)的變量的第 3列數(shù)據(jù)定義為被解釋變量y，將名為A的變量的第4、5、6列數(shù)據(jù)定義為解釋變量矩陣 X,定義一個(gè)有N*T行，1列的全1矩陣，該矩陣名為：xconstant , ( ones即為全1矩陣)說(shuō)明解釋變量矩陣x的大小：有nobs行，K列。（size為描述矩陣的大?。?。附錄:靜態(tài)面板空間計(jì)量經(jīng)濟(jì)學(xué)一、OLS靜態(tài)面板編程1、普通面板編程T=30;N=46;W=normw(W1);y=A(：,3)；x=A(:,4,6);xconstant=ones(N*T,1

11、);nobs K=size(x);results=ols(y,xconstant x);vnames=strvcat( logcit* /intercept' ,'logp* ,'logy1);prt_reg(results,vnames,1);sige=*(nobs-K )/no bs);Ioglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*% The (robust)LM tests developed by ElhorstLMsarsem_panel(results,W,y,xconstant x); % (Robust) LM

12、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,mea nn y,mea nn x,meanty,meantx=demean(y,x,N,TJmodel);results=ols(ywith,xwith);vnames=strvcat(,logcit,/logp,logy,); % should be changed if x ischangedprt_reg(res

13、ults,vnames);sfe=meanny-meannx*; % including the constant term yme = y - mean(y);et=ones(T,1);error=y-kron(et,sfe)-x*;rsqrl = error" error;rsqr2 = yme'*yme;FE_rsqr2 = - rsqrl/rsqr2 % r-squared including fixed effectssige=*(nobs-K)/nobs);logliksfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*1*LMsarsem

14、_panel(results,W,ywith,xwith); % (Robust) LM tests3、時(shí)期固定 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,mea nn y,mea nn x,meanty,meantx=demean(y,x,N,T,model);results=ols(ywith,xwith);vnames=strvcat(logcit7logp,/logy

15、,); % should be changed if x is changedprt_reg(results,vnames);tfe=meanty-meantx*; % including the constant termyme = y - mean(y);en=ones(N,1);error=y-kron(tfe,en)-x*;rsqrl = errorf*error;rsqr2 = yme'*yme;FE_rsqr2 = - rsqr1/rsqr2 % r-squared including fixed effects sige=*(nobs-K )/no bs);loglikt

16、fe=-nobs/2*log(2*pi*sige)-1/(2*sige)*1*LMsarsem_panel(results,W,ywithJxwith); % (Robust) LM tests4、空間與時(shí)間雙固定模型T=30;N=46;W=normw(W1);y=A(：,3)；x=A(:,4,6);xconstant=ones(N*T,1);nobs K=size(x);model=3;ywith,xwith,mea nn y,meannx,meanty,meantx=demean(y,x,N5T,model)；results=ols(ywith,xwith);vnames=strvcat(

17、,logcit7logp,logy,); % should be changed if x is changedprt_reg(results,vnames)en=ones(N,1);et=ones(T,1);intercept=mean(y)-mean(x)*;sfe=meanny-mea nn x*(en,intercept);tfe=meanty-meantx*(et,intercept);yme = y - mean(y);ent=ones(N*T,1);error=y-kron(tfe,en)-kron(et,sfe)-x*(ent,intercept);rsqrl = errorf

18、*error;rsqr2 = yme'*yme;FE_rsqr2 = - rsqrl/rsqr2 % r-squared including fixed effectssige=*(nobs-K )/no bs);loglikstfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*1*LMsarsem_panel(results,W,ywith,xwith); % (Robust) LM tests二、靜態(tài)面板SAR模型1 > 無(wú)固定效應(yīng)(No fixed effects)T=30;N=46;W=normw(W1);y=A(:,3);x=A(:,4,6);f

19、or t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0;=0;=0;results=sar_panel_FE(y,xc onstant x,W,T,i nfo);vnames=strvcat( 'logcit', 'intercept', 'logp' , 'logy' );prt_sp new(results,v names, 1)% Print out effectsestimates spa

20、t_model=0;direct_indirect_effects_estimates(results,W,spat_model);pan el_effects_sar(results,vnames5W);2、空間固定效應(yīng)(Spatial fixed effects)T=30;N=46;W=normw(W1);y=A(:,3)；x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0；=1;=0；results=sar_pa nel_FE

21、(y,x,W,T,info);vnames=strvcat('logcit1 , logp1 , 'logy1 );prt_sp new(results，v names J)% Print out effectsestimates spat_model=0;pan el_effects_sar(results,vnames,W);direct_indirect_effects_estimates(results,W,spat_model);3、時(shí)點(diǎn)固定效應(yīng)(Time period fixed effectsT=30;N=46;W=normw(W1);y=A(:,3)；x=A(:

22、,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0;% required for exact results=2;=0;% Do not print intercept and fixed effects; use =1 to turn onresults=sar_pa nel_FE(y,x,W,T,info);vnames=strvcat('logcif , logp' , logy* );prt_sp new(result

23、s，v names J)% Print out effectsestimates spat_model=0;direct_indirect_effects_estimates(results,W,spat_model);4、雙固定效應(yīng) (Spatial and time period fixed effectspan el_effects_sar(results,vnames3W);T=30;N=46;W=normw(W1);y=A(:,3)；x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstan

24、t=ones(N*T,1);nobs K=size(x);=0;% required for exact results=3;=0;% Do not print intercept and fixed effects; use =1 to turn onresults=sar_pa nel_FE(y,x,W,T,info);vnames=strvcat('logcit1 , logp1 , 'logy1 );prt_sp new(results,v names J)% Print out effectsestimates spat_model=0;direct_indirect

25、_effects_estimates(results,W5spat_model);pan el_effects_sar(results,vnames3W);三、靜態(tài)面板SDM模型1、無(wú)固定效應(yīng)(No fixed effectsT=30;N=46;W=normw(W1);y=A(:,3)；x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0；=0；=0；results=sar_pa nel_FE(y,xc on stant x wx,W

26、,T,info);vnames=strvcat( 'logcit1 , Intercept1, 'logp1 , logy1 , 'WJogp' , 'WJogy' );prt_sp new(results,v names J)% Print out effectsestimates spat_model=1;direct_indirect_effects_estimates(results,W5spat_model);pan el_effects_sdm(results,vnames,W);2> 空間固定效應(yīng)(Spatial fixed

27、effects)T=30;N=46;W=normw(W1);y=A(:,3)；x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0;% required for exact results=1;=0;% Do not print intercept and fixed effects; use =1 to turn on,WJogy');results=sar_pa nel_FE(y,x wx,W,T,i nfo);vname

28、s=strvcat( 'logcit1 , logp1 , 'logy1 JWJogp' prt_sp new(results，v names)% Print out effectsestimates spat_model=1;direct_indirect_effects_estimates(results,W5spat_model);pan e l_ef f ects_sd m (res u I ts, v names,W);3、時(shí)點(diǎn)固定效應(yīng)(Time period fixed effectsT=30;N=46;x=A(:,4,6);for t=1:Tt1=(t-1

29、)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0;% required for exact results=2；=0;% Do not print intercept and fixed effects; use =1 to turn on% New routines to calculate effects estimatesresults=sar_panel_FE(y,x wx,W,T,info);,W*logy );vnames=strvcat( 'logcit1 ,

30、logp' , 'logy* JWJogp'% Print out coefficient estimatesprt_spnew(results,vnames, 1)% Print out effectsestimates spat_model=1;direct_indirect_effects_estimates(results,W5spat_model);pan e l_ef f ects_sd m (res u I ts, v names,W)4、雙固定效應(yīng) (Spatial and time period fixed effectsT=30;N=46;W=nor

31、mw(W1);y=A(:,3)；x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0；=0;% required for exact results=3;=0;% Do not print intercept and fixed effects; use =1 to turn onresults=sar_panel_FE(y,x wx,W,T,info);vnames=strvcat( 'logcit' , logp&

32、#39; , 'logy* , 'W*logp* , 'W*logy' );prt_sp new(results，v names, 1)% Print out effectsestimates spat_model=1;direct_indirect_effects_estimates(results,W5spat_model);pan el_effects_sdm(results,v names,W)wald test spatial lag% Wald test for spatial Durbin model against spatial lagmode

33、l btemp=;varcov=;Rafg=zeros(K,2*K+2);for k=1:KRafg(k,K+k)=1;end% R(1,3)=0 and R4)=0;Wald_spatial_lag=(Rafg*btemp),*inv(Rafg*varcov*Rafgl)*Rafg*btempprob_spatial_lag=1 -chis_cdf (Wald_spatial_lag, K)wald test spatial error% Wald test spatial Durbin model against spatial errormodel R=zeros(K,1);for k=

34、1:KR(k)=btemp(2*K+1 )*btemp(k)+btemp(K+k);% k changed in 1,7/12/2010% R(1 )=btemp(5)*btemp(1 )+btemp(3);% R(2)=btemp(5)*btemp(2)+btemp(4);endRafg=zeros(K,2*K+2);for k=1:KRafg(k,k) =btemp(2*K+1);% k changed in 1,7/12/2010Rafg(k,K+k)=1;Rafg(k,2*K+1 )=btemp(k);% Rafg(1,1 )=btemp(5);Rafg(1,3)=1 ;Rafg(1,

35、5)=btemp(1);% Rafg(2,2)=btemp(5);Rafg(2,4)=1 ;Rafg(2,5)=btemp(2);endWald_spatial_error=R'*inv(Rafg*varcov*Rafg')*Rprob_spatial_error=1 -chis_cdf (Wald_spatial_error,K)LR test spatial lagresultssar=sar_pa nel_FE(y,x,W,T,info);LR_spatial_lag=-2* (LR_spatial_lag,K)LR test spatial errorresultsse

36、m=sem_pa nel_FE(y,x，W,T,info);LR_spatial_error=-2* (LR_spatial_error,K)5、空間隨機(jī)效應(yīng)與時(shí)點(diǎn)固定效應(yīng)模型T=30;N=46;W=normw(W1);y=A(:,3);x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);ywith,xwith,meanny,meannx,meanty5meantx=demean(y,x wx,N,T,2);2=time dummies=

37、1；results=sar_pa nel_RE(ywith,xwith,W,T,i nfo);prt_sp new(results，v names J) spat_model=1;direct_indirect_effects_estimates(results,W5spat_model);pan e l_ef f ects_sd m (res u I ts, v names,W)wald test spatial lagbtemp=(1:2*K+2);varcov=(1:2*K+2,1:2*K+2);Rafg=zeros(K,2*K+2);for k=1:KRafg(k,K+k)=1;% R

38、(1,3)=0 and R(2,4)=0;endWald_spatial_lag=(Rafg*btemp)'*inv(Rafg*varcov*Rafg')*Rafg*btempprob_spatial_lag= chis_cdf (Wald_spatial_lag, K)wald test spatial errorR=zeros(K,1);for k=1:K% k changed in 1,R(k)=btemp(2*K+1)*btemp(k)+btemp(K+k);7/12/2010% R(1 )=btemp(5)*btemp(1 )+btemp(3);% R(2)=btem

39、p(5)*btemp(2)+btemp(4);endRafg=zeros(K,2*K+2);for k=1:KRafg(k,k) =btemp(2*K+1);% k changed in 1,7/12/2010Rafg(k,K+k)=1;Rafg(k,2*K+1 )=btemp(k);% Rafg(1,1)=btemp(5);Rafg(1,3)=1 ;Rafg(1,5)=btemp(1);% Rafg(2,2)=btemp(5);Rafg(2,4)=1 ;Rafg(2,5)=btemp(2);endWald_spatial_error=R'*inv(Rafg*varcov*Rafg&#

40、39;)*Rprob_spatial_error= 1-chis_cdf (Wald_spatial_error,K)LR test spatial lagresultssar=sar_panel_RE(ywith,xwith(:,1:K),W,T,info);LR_spatial_lag=-2* (LR_spatial_lag,K)LR test spatial errorresultssem=sem_pa nel_RE(ywith,xwith(:,1:K),W5T,info);LR_spatial_error=-2* (LR_spatial_error,K)四、靜態(tài)面板SEM模型1 >

41、; 無(wú)固定效應(yīng)(No fixed effects)T=30;N=46;W=normw(W1);y=A(:,3)；x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0；=0；=0；results=sem_panel_FE(y，xconstant x,W,T,info);vnames=strvcat('logcif, Sntercept1, logp* , logy1 );prt_sp new(results，v names, 1

42、)% Print out effectsestimates spat_model=0;direct_indirect_effects_estimates(results,W5spat_model);pan el_effects_sar(results,vnames,W);2、空間固定效應(yīng)(Spatial fixed effectsT=30;N=46;W=normw(W1);y=A(:,3)；x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x

43、);=0；=1;=0；results=sem_pa nel_FE(y，x,W,T,i nfo);vnames=strvcat('logcif , logp' , logy* );prt_sp new(results，v names J)% Print out effectsestimates spat_model=0;direct_in direct_effects_estimates(results,W,spat_model);pan el_effects_sar(results,v names,W);3、時(shí)點(diǎn)固定效應(yīng)(Time period fixed effects)T=

44、30；N=46;W=normw(W1);y=A(:,3)；x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxconstant=ones(N*T,1);nobs K=size(x);=0;% required for exact results=2;=0;% Do not print intercept and fixed effects; use =1 to turn onresults=sem_pa nel_FE(y,x,W,T,i nfo);vnames=strvcat('logcif , lo

45、gp' , logy* );prt_sp new(results，v names, 1)% Print out effectsestimates spat_model=0;direct_indirect_effects_estimates(results,W5spat_model);pan el_effects_sar(results,vnames,W);4、雙固定效應(yīng) (Spatial and time period fixed effects)T=30;N=46;W=normw(W1);y=A(:,3)； x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*

46、N;wx(t1:t2,:)=W*x(t1:t2,:);endxcon stant=ones(N*T,1);nobs K=size(x);=0;% required for exact results=3;=0;% Do not print intercept and fixed effects; use =1 to turn onresults=sem_pa nel_FE(y，x,W,T，i nfo);vnames=strvcat('logcit1 , logp' , 'logy* );prt_sp new(results，v names, 1)% Print out

47、effectsestimates spat_model=0;direct_indirect_effects_estimates(results,W5spat_model);pan el_effects_sar(results,vnames3W);五、靜態(tài)面板SDEM模型1 > 無(wú)固定效應(yīng)(No fixed effects)T=30;N=46;y=A(:,3)； x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxcon stant=ones(N*T,1);nobs K=size(x);=0；=0；=0；

48、results=sem_panel_FE(y，xconstant x wx,W,T,info);vnames=strvcat( 'logcit1 , Intercept1, 'logp* , logy1 , 'WJogp' , 'WJogy' );prt_sp new(results，v names J)% Print out effectsestimates spat_model=1;direct_indirect_effects_estimates(results,W5spat_model);pan el_effects_sdm(results,v names,W);2> 空間固定效應(yīng)(Spatial fixed effects)T=30;N=46;W=normw(W1);y=A(:,3)； x=A(:,4,6);for t=1:Tt1=(t-1)*N+1;t2=t*N;wx(t1:t2,:)=W*x(t1:t2,:);endxcon stant