MATLAB-空間計量模型詳細步驟

1、I.excel與MATLA利接: Excel: 選項一一加載項一一 COM 加載項一一轉(zhuǎn)到一一沒有勾選項 2 .MATLAB?裝目錄中尋找toolboxexlink點擊，啟 用宏 E:MATLABtoolboxexlink 芒砂 .薩比口現(xiàn) 9UM *Comterits 2013/12/282:14 MATLABCede 1KB 醒exdlink 2013/10/933:48 MicrosoftExcelB. 108KE QeMcllink2QO3 2013/10/923:40 MicrosoftExcel 305KB 網(wǎng)ExliSamp 2009/11/519:58 MicrosoftExc

2、el 72KB cetfunctionhelp 2010/12/42:17 MATLABCode 2KB getfunctionlkt 2010/12417 MATLABCede 3KB cetfunrti3n5igna：ure 2OL2/5/42:59 MATLABCede 4KB 國inf口 2007/7/240:22 XML登 1KB 然后，Excel中就出現(xiàn) MATLAB工具 （注意 Excel中的數(shù)據(jù)：）帶樣式 靜事 刪賒格式 W 目即求知，國岷.之清除. 單元格 打了4V推用蹄選且持周kJ皤TLA肛其 COM 加蓑項. 、rmrafni一】i山Hn1i 上4 4UU/JJJ U1比

3、W4 UIUU 4353曲3 IB141 IdI4.b 2 29 9 4 41 1 48402424840242 01801810180181 0.1535190.153519 47757624775762 579.579.2 2 45214521 29 42 4683981 0.118179 -0.0326 率 8M899 6354.1 4965.1 2 29 9 4 43 3 4.3966834.396683 0.3505260.350526 0.17S4080.17S408 4.8564634.856463 5152251522 3950.3950.2 2 Z9 4 44 4 45Mg6

4、7 0120077 -00326 4568586 10621 1465.4 2 29 9 4 45 5 456642g456642g 02295330229533 02408060240806 4 4758261758261 5017350173 3853338533 2 29 9 4 46 6 4 4707727707727 01002670100267 0.1313910.131391 47116924711692 43334333 36053605 曷量數(shù)陶權(quán)重期居 WdataBWdataBwlwlK KSheet3Sheet3 3 .啟動matlab （1） 點擊 startMATLA

5、B （2） senddatatomatlab，并對變量矩陣變量進行命名（注意：選取變量為數(shù)值，不包括 各變量） （data 表中數(shù)據(jù)進行命名） （空間權(quán)重進行命名） 2 2j j 1 1 4512i4512itotnjtotnj ChCh閣外 - -iii 皿1 1 3 3零處珈 Ul5Ul5 W W ) 4 4 q q 1 1 1 1心椀H H 川陋H H iiTiiTffi!ffi! 泗姒 1?M71?M7 Q j j .妣雅 T T例 4W”4W” 4E4E iWiW 例點 G G q q 5 5 5MI1P5MI1P 1132111321 4 4i i跑* iKJf.iKJf. a i

6、 i 口1 1；猊 4 4，如 4 4的 Lift2MLift2M 31,31,P P q ? ? 專 如 師1 1 1714217142 Illllll WiWi q q I I 0堆璉 - -fl161fl1614444 / / inoiHinoiH 的， Q P P 1 1附期 4 4的 4lrlr淑 刑*1*1 A1 1! ! I 4 4閑 川ilftllilftll 4.u4.u躋五 4414411111 0 0 i ii i 1 1購M M， 卓璃2b2b 4I4I加 m m測 niltnilt- - I I 1 12 2 - -022H2022H2 - -aZ3BnaZ3Bn 4

7、.27Tm4.27Tm snirsnir q q U U 爾1 1 4 4皿中 - -flfl 詢r r a a 1 14 4 - -tmumtmum - -0 01 1硒 - -a aITW，25ir25ir MlMl q q 1 13 3 A&TOIA&TOI 用 fli13311331 3123121 1 MMMMt Q 1 11 1 4rni4rni 41fliH41fliHr r - -a aIBH *E3l2f*E3l2f MWMW Z1MZ1M a 1 1, , 用E3WiE3Wi - -fl3fl3fl1fl1 j j- -mimi MSFMSF g g 1 1

8、1 1 4M4Mr r4242 - -4HI4I4HI4I - -07H3107H31 4I?4I?r r q q 1 1J J IPHIIPHI - -azMiazMi 1dOEZd.1dOEZd. 網(wǎng))M)M q q 1 1 i叫 司4G4G明 中 4 4網(wǎng)】 例口 a a 2 21 1 ?2TIM2?2TIM2 11 - -a a皿0 0 1 1 5UZ5UZ q q l 3FWTH3FWTH 1 1HMfHMf a a Z Z .traiiM.traiiMZ Z - -aeraer】 - -aZJT7aZJT7T T iraiaiiraiai 3Mif3Mif 八 w w jriiHM

9、i 用MwrMwr flfl 1 1 4A.4A.If.If. -uo （3） 導(dǎo)入 MATLAB中的兩個矩陣變量就可以看見詳細恬息 4 .將elhorst和jplv7兩個程序文件夾復(fù)制到MATLAB安裝目 錄的toolbox文件夾 2017/1/1415；08 文件央 5 .設(shè)置路徑: 工作區(qū) |mJ 值 IJSOJCSd口血46x46double elhorstelhorst 2017/1/15&46 文件夾 6 .輸入程序，得出結(jié)果 1=30;1=30; N=40;N=40; W=norJILVW1); y=A(y=A(：j3)j3)； K=A(：J4,6); xconsxcon

10、s- -tant=arietant=arieE1)：nobsK=sizenobsK=size; ; T=30; N=46; W=norm(W1); y=A(:,3); x=A(:,4,6); xconstant=ones(N*T,1); nobsK=size(x); results=ols(y,xconstantx); 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

11、/(2*sige)*results.resid*results.resid %The(robust)LMtestsdevelopedbyElhorst LMsarsem_panel(results,W,y,xconstantx);%(Robust)LMtests 解釋 每一行分別表示： 該面板數(shù)據(jù)的時期數(shù)為 30(T=30), 該面板數(shù)據(jù)有 30 個地區(qū)(N=30), 將空間權(quán)重矩陣標準化(W=normw(w1), 將名為 A(以矩陣形式出現(xiàn)在 MATLABA 中)的變量的第 3 列數(shù)據(jù)定義為被解釋變量 y, 將名為 A 的變量的第 4、5、6 列數(shù)據(jù)定義為解釋變量矩陣 x, 定義一個有 N*

12、T行，1 列的全 1 矩陣，該矩陣名為：xconstant,(ones 即為全 1 矩陣)說明解釋變量矩陣 x的大?。河?nobs 行，K列。(size為描述矩陣的大小)。 附錄： 靜態(tài)面板空間計量經(jīng)濟學(xué) 一、OLS靜態(tài)面板編程 1、普通面板編程 T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); xconstant=ones(N*T,1); nobsK=size(x); results=ols(y,xconstantx); vnames=strvcat(logcit,intercept,logp,logy); prt_reg(results,vnam

13、es,1); sige=results.sige*(nobs-K)/nobs); loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid*results.resid %The(robust)LMtestsdevelopedbyElhorst LMsarsem_panel(results,W,y,xconstantx);%(Robust)LMtests 2、空間固定OLS(spatial-fixedeffects) T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); xconstant=ones

14、(N*T,1); nobsK=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);%shouldbechangedifxis changed prt_reg(results,vnames); sfe=meanny-meannx*results.beta;%includingtheconstantterm yme=y-mean(y); et=ones(T,1); error

15、=y-kron(et,sfe)-x*results.beta; rsqrl=error*error; rsqr2=yme*yme; FE_rsqr2=1.0-rsqr1/rsqr2%r-squaredincludingfixedeffects sige=results.sige*(nobs-K)/nobs); logliksfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid*results.residLMsarsem_panel(results,W,ywith,xwith);%(Robust)LMtests 3、時期固定OLS(time-per

16、iodfixedeffects) T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); xconstant=ones(N*T,1); nobsK=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);%shouldbechangedifxis changed prt_reg(results,vnames); tfe=meanty-

17、meantx*results.beta;%includingtheconstantterm yme=y-mean(y); en=ones(N,1); error=y-kron(tfe,en)-x*results.beta; rsqrl=error*error; rsqr2=yme*yme; FE_rsqr2=1.0-rsqr1/rsqr2%r-squaredincludingfixedeffects sige=results.sige*(nobs-K)/nobs); logliktfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid*result

18、s.resid LMsarsem_panel(results,W,ywith,xwith);%(Robust)LMtests 4、空間與時間雙固定模型 T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); xconstant=ones(N*T,1); nobsK=size(x); model=3; ywith,xwith,meanny,meannx,meanty,meantx=demean(y,x,N,T,model ); results=ols(ywith,xwith); vnames=strvcat(logcit,logp,logy);%should

19、bechangedifxis changed prt_reg(results,vnames) en=ones(N,1); et=ones(T,1); intercept=mean(y)-mean(x)*results.beta; sfe=meanny-meannx*results.beta-kron(en,intercept); tfe=meanty-meantx*results.beta-kron(et,intercept); yme=y-mean(y); ent=ones(N*T,1); error=y-kron(tfe,en)-kron(et,sfe)-x*results.beta-kr

20、on(ent,intercept); rsqr1=error*error; rsqr2=yme*yme; FE_rsqr2=1.0-rsqr1/rsqr2%r-squaredincludingfixedeffects sige=results.sige*(nobs-K)/nobs); loglikstfe=-nobs/2*log(2*pi*sige)- 1/(2*sige)*results.resid*results.resid LMsarsem_panel(results,W,ywith,xwith);%(Robust)LMtests 二、靜態(tài)面板SAR模型 1、無固定效應(yīng)(Nofixede

21、ffects) T=30; N=46; W=normw(W1); y=A(：,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0; info.model=0; info.fe=0; results=sar_panel_FE(y,xconstantx,W,T,info); vnames=strvcat(logcit,intercept,logp,logy); prt_spnew(results,v

22、names,1) %Printouteffectsestimates spat_model=0; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sar(results,vnames,W); 2、空間固定效應(yīng)(Spatialfixedeffects) T=30; N=46; W=normw(W1); y=A(：,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1

23、); nobsK=size(x); info.lflag=0; info.model=1; info.fe=0; results=sar_panel_FE(y,x,W,T,info); vnames=strvcat(logcit,logp,logy); prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=0; direct_indirect_effects_estimates(results,W,spat_model);panel_effects_sar(results,vnames,W); 3、日寸點固定效應(yīng)(Ti

24、meperiodfixedeffects) T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0;%requiredforexactresults info.model=2; info.fe=0;%Donotprintinterceptandfixedeffects;useinfo.fe=1toturnon results=sar_pa

25、nel_FE(y,x,W,T,info); vnames=strvcat(logcit,logp,logy); prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=0; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sar(results,vnames,W); 4、雙固定效應(yīng)(Spatialandtimeperiodfixedeffects T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,

26、4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0;%requiredforexactresults info.model=3; info.fe=0;%Donotprintinterceptandfixedeffects;useinfo.fe=1toturnon results=sar_panel_FE(y,x,W,T,info); vnames=strvcat(logcit,logp,logy); prt_spn

27、ew(results,vnames,1) %Printouteffectsestimates spat_model=0; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sar(results,vnames,W); 三、靜態(tài)面板SDM模型 1、無固定效應(yīng)(Nofixedeffects) T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xc

28、onstant=ones(N*T,1); nobsK=size(x); info.lflag=0; info.model=0; info.fe=0; results=sar_panel_FE(y,xconstantxwx,W,T,info); vnames=strvcat(logcit,intercept,logp,logy,W*logp,W*logy);prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=1; direct_indirect_effects_estimates(results,W,spat_mode

29、l); panel_effects_sdm(results,vnames,W); 2、空間固定效應(yīng)(Spatialfixedeffects) T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0;%requiredforexactresults info.model=1; info.fe=0;%Donotprintinterceptan

30、dfixedeffects;useinfo.fe=1toturnon results=sar_panel_FE(y,xwx,W,T,info); vnames=strvcat(logcit,logp,logy,W*logp,W*logy);prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=1; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sdm(results,vnames,W); 3、時點固定效應(yīng)(Timeperio

31、dfixedeffects) T=30; N=46; W=norm(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0;%requiredforexactresults info.model=2; info.fe=0;%Donotprintinterceptandfixedeffects;useinfo.fe=1toturnon %Newroutinestocalculat

32、eeffectsestimates results=sar_panel_FE(y,xwx,W,T,info); vnames=strvcat(logcit,logp,logy,W*logp,W*logy); %Printoutcoefficientestimates prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=1; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sdm(results,vnames,W) 4、雙固定

33、效應(yīng)(Spatialandtimeperiodfixedeffects) T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.bc=0; info.lflag=0;%requiredforexactresults info.model=3; info.fe=0;%Donotprintinterceptandfixedeffects;useinfo.f

34、e=1toturnon results=sar_panel_FE(y,xwx,W,T,info); vnames=strvcat(logcit,logp,logy,W*logp,W*logy);prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=1; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sdm(results,vnames,W) waldtestspatiallag %WaldtestforspatialDurb

35、inmodelagainstspatiallagmodelbtemp=results.parm; varcov=results.cov; Rafg=zeros(K,2*K+2); fork=1:K Rafg(k,K+k)=1;%R(1,3)=0andR(2,4)=0; end Wald_spatial_lag=(Rafg*btemp)*inv(Rafg*varcov*Rafg)*Rafg*btempprob_spatial_lag=1-chis_cdf(Wald_spatial_lag,K) waldtestspatialerror %WaldtestspatialDurbinmodelaga

36、instspatialerrormodelR=zeros(K,1); fork=1:K R(k)=btemp(2*K+1)*btemp(k)+btemp(K+k);%kchangedin1, 7/12/2010%R(1)=btemp(5)*btemp(1)+btemp(3); %R(2)=btemp(5)*btemp(2)+btemp(4); end Rafg=zeros(K,2*K+2); fork=1:K Rafg(k,k)=btemp(2*K+1);%kchangedin1,7/12/2010 Rafg(k,K+k)=1; Rafg(k,2*K+1)=btemp(k); %Rafg(1,

37、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);end Wald_spatial_error=R*inv(Rafg*varcov*Rafg)*Rprob_spatial_error=1-chis_cdf(Wald_spatial_error,K) LRtestspatiallag resultssar=sar_panel_FE(y,x,W,T,info); LR_spatial_lag=-2*(resultssar.lik-results.lik) pr

38、ob_spatial_lag=1-chis_cdf(LR_spatial_lag,K) LRtestspatialerror resultssem=sem_panel_FE(y,x,W,T,info); LR_spatial_error=-2*(resultssem.lik-results.lik)prob_spatial_error=1-chis_cdf(LR_spatial_error,K) 5、空間隨機效應(yīng)與時點固定效應(yīng)模型 T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1

39、:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); ywith,xwith,meanny,meannx,meanty,meantx=demean(y,xwx,N,T,2); 2=timedummies info.model=1; results=sar_panel_RE(ywith,xwith,W,T,info); prt_spnew(results,vnames,1) spat_model=1; direct_indirect_effects_estimates(results,W,spat_model); panel

40、_effects_sdm(results,vnames,W) waldtestspatiallag btemp=results.parm(1:2*K+2); varcov=results.cov(1:2*K+2,1:2*K+2); Rafg=zeros(K,2*K+2); fork=1:K Rafg(k,K+k)=1;%R(1,3)=0andR(2,4)=0; end Wald_spatial_lag=(Rafg*btemp)*inv(Rafg*varcov*Rafg)*Rafg*btempprob_spatial_lag=1-chis_cdf(Wald_spatial_lag,K) wald

41、testspatialerror R=zeros(K,1); fork=1:K R(k)=btemp(2*K+1)*btemp(k)+btemp(K+k);%kchangedin1, 7/12/2010%R(1)=btemp(5)*btemp(1)+btemp(3); %R(2)=btemp(5)*btemp(2)+btemp(4); end Rafg=zeros(K,2*K+2);fork=1:KRafg(k,k)=btemp(2*K+1); %kchangedin1,7/12/2010Rafg(k,K+k)=1; Rafg(k,2*K+1)=btemp(k); %Rafg(1,1)=bte

42、mp(5);Rafg(1,3)=1;Rafg(1,5)=btemp(1); %Rafg(2,2)=btemp(5);Rafg(2,4)=1;Rafg(2,5)=btemp(2);end Wald_spatial_error=R*inv(Rafg*varcov*Rafg)*Rprob_spatial_error=1-chis_cdf(Wald_spatial_error,K) LRtestspatiallag resultssar=sar_panel_RE(ywith,xwith(:,1:K),W,T,info); LR_spatial_lag=-2*(resultssar.lik-result

43、s.lik) prob_spatial_lag=1-chis_cdf(LR_spatial_lag,K) LRtestspatialerror resultssem=sem_panel_RE(ywith,xwith(:,1:K),W,T,info); LR_spatial_error=-2*(resultssem.lik-results.lik) prob_spatial_error=1-chis_cdf(LR_spatial_error,K) 四、靜態(tài)面板SEM模型 1、無固定效應(yīng)（Nofixedeffects T=30; N=46; W=normw(W1); y=A(:,3); x=A(:

44、,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0; info.model=0; info.fe=0; results=sem_panel_FE(y,xconstantx,W,T,info); vnames=strvcat(logcitprt_spnew(results,vnames,1)%Printouteffectsestimatesspat_model=0; direct_indirect_effects_

45、estimates(results,W,spat_model); panel_effects_sar(results,vnames,W); 2、空間固定效應(yīng)(Spatialfixedeffects) T=30; N=46; ,intercept logp,logy); W=normw(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0; info.model=1; info

46、.fe=0; results=sem_panel_FE(y,x,W,T,info); vnames=strvcat(logcit,logp,logy); prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=0; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sar(results,vnames,W); 3、時點固定效應(yīng)(Timeperiodfixedeffects T=30; N=46; W=normw(W1); y=A(

47、:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0;%requiredforexactresults info.model=2;info.fe=0;%Donotprintinterceptandfixedeffects;useinfo.fe=1toturnon results=sem_panel_FE(y,x,W,T,info); vnames=strvcat(logcit,logp,log

48、y); prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=0; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sar(results,vnames,W); 4、雙固定效應(yīng)(Spatialandtimeperiodfixedeffects T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(

49、t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0;%requiredforexactresults info.model=3; info.fe=0;%Donotprintinterceptandfixedeffects;useinfo.fe=1toturnon results=sem_panel_FE(y,x,W,T,info); vnames=strvcat(logcit,logp,logy); prt_spnew(results,vnames,1) %Printouteffectsestimates spat_

50、model=0; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sar(results,vnames,W); 五、靜態(tài)面板SDEM模型 1、無固定效應(yīng)(Nofixedeffects) T=30; N=46; W=normw(W1); y=A(:,3); x=A(:,4,6); fort=1:T t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end xconstant=ones(N*T,1); nobsK=size(x); info.lflag=0; i

51、nfo.model=0; info.fe=0; results=sem_panel_FE(y,xconstantxwx,W,T,info); vnames=strvcat(logcit,intercept,logp,logy,W*logp,W*logy);prt_spnew(results,vnames,1) %Printouteffectsestimates spat_model=1; direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sdm(results,vnames,W); 2、空間固定效應(yīng)(Spatialfi