µÚÒ»·¶ÎÄÍø - רҵÎÄÕ·¶ÀýÎĵµ×ÊÁÏ·ÖÏíÆ½Ì¨

Áù²½Ñ§»áÓÃMATLAB×ö¿Õ¼ä¼ÆÁ¿»Ø¹éÏêϸ²½Öè

À´Ô´£ºÓû§·ÖÏí ʱ¼ä£º2025/6/12 14:49:47 ±¾ÎÄÓÉloading ·ÖÏí ÏÂÔØÕâÆªÎĵµÊÖ»ú°æ
˵Ã÷£ºÎÄÕÂÄÚÈݽö¹©Ô¤ÀÀ£¬²¿·ÖÄÚÈÝ¿ÉÄܲ»È«£¬ÐèÒªÍêÕûÎĵµ»òÕßÐèÒª¸´ÖÆÄÚÈÝ£¬ÇëÏÂÔØwordºóʹÓá£ÏÂÔØwordÓÐÎÊÌâÇëÌí¼Ó΢ÐźÅ:xxxxxxx»òQQ£ºxxxxxx ´¦Àí£¨¾¡¿ÉÄܸøÄúÌṩÍêÕûÎĵµ£©£¬¸ÐлÄúµÄÖ§³ÖÓëÁ½⡣

wx(t1:t2,:)=W*x(t1:t2,:); end

xconstant=ones(N*T,1); [nobs K]=size(x); info.lflag=0; info.model=0; info.fe=0;

results=sem_panel_FE(y,[xconstant x wx],W,T,info);

vnames=strvcat('logcit','intercept','logp','logy','W*logp','W*logy'); prt_spnew(results,vnames,1) % Print out effects estimates spat_model=1;

direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sdm(results,vnames,W);

2¡¢¿Õ¼ä¹Ì¶¨Ð§Ó¦£¨Spatial fixed effects£©

T=30; N=46; W=normw(W1); y=A(:,[3]); x=A(:,[4,6]); for t=1:T

t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end

xconstant=ones(N*T,1); [nobs K]=size(x);

info.lflag=0; % required for exact results info.model=1;

info.fe=0; % Do not print intercept and fixed effects; use info.fe=1 to turn on

results=sem_panel_FE(y,[x wx],W,T,info);

vnames=strvcat('logcit','logp','logy','W*logp','W*logy'); prt_spnew(results,vnames,1) % Print out effects estimates spat_model=1;

direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sdm(results,vnames,W);

3¡¢Ê±µã¹Ì¶¨Ð§Ó¦£¨Time period fixed effects£©

T=30; N=46; W=normw(W1); y=A(:,[3]);

21

x=A(:,[4,6]); for t=1:T

t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end

xconstant=ones(N*T,1); [nobs K]=size(x);

info.lflag=0; % required for exact results info.model=2;

info.fe=0; % Do not print intercept and fixed effects; use info.fe=1 to turn on

% New routines to calculate effects estimates results=sem_panel_FE(y,[x wx],W,T,info);

vnames=strvcat('logcit','logp','logy','W*logp','W*logy'); % Print out coefficient estimates prt_spnew(results,vnames,1) % Print out effects estimates spat_model=1;

direct_indirect_effects_estimates(results,W,spat_model); panel_effects_sdm(results,vnames,W)

4¡¢Ë«¹Ì¶¨Ð§Ó¦£¨Spatial and time period fixed effects£©

T=30; N=46; W=normw(W1); y=A(:,[3]); x=A(:,[4,6]); for t=1:T

t1=(t-1)*N+1;t2=t*N; wx(t1:t2,:)=W*x(t1:t2,:); end

xconstant=ones(N*T,1); [nobs K]=size(x); info.bc=0;

info.lflag=0; % required for exact results info.model=3;

info.fe=0; % Do not print intercept and fixed effects; use info.fe=1 to turn on

results=sem_panel_FE(y,[x wx],W,T,info);

vnames=strvcat('logcit','logp','logy','W*logp','W*logy'); prt_spnew(results,vnames,1) % Print out effects estimates spat_model=1;

direct_indirect_effects_estimates(results,W,spat_model);

22

panel_effects_sdm(results,vnames,W)

23

Áù²½Ñ§»áÓÃMATLAB×ö¿Õ¼ä¼ÆÁ¿»Ø¹éÏêϸ²½Öè.doc ½«±¾ÎĵÄWordÎĵµÏÂÔØµ½µçÄÔ£¬·½±ã¸´ÖÆ¡¢±à¼­¡¢ÊղغʹòÓ¡
±¾ÎÄÁ´½Ó£ºhttps://www.diyifanwen.net/c3sems5t22j9vfqx3d4pq7px008twst015cz_6.html£¨×ªÔØÇë×¢Ã÷ÎÄÕÂÀ´Ô´£©
ÈÈÃÅÍÆ¼ö
Copyright © 2012-2023 µÚÒ»·¶ÎÄÍø °æÈ¨ËùÓÐ ÃâÔðÉùÃ÷ | ÁªÏµÎÒÃÇ
ÉùÃ÷ :±¾ÍøÕ¾×ðÖØ²¢±£»¤ÖªÊ¶²úȨ£¬¸ù¾Ý¡¶ÐÅÏ¢ÍøÂç´«²¥È¨±£»¤ÌõÀý¡·£¬Èç¹ûÎÒÃÇ×ªÔØµÄ×÷Æ·ÇÖ·¸ÁËÄúµÄȨÀû,ÇëÔÚÒ»¸öÔÂÄÚ֪ͨÎÒÃÇ£¬ÎÒÃǻἰʱɾ³ý¡£
¿Í·þQQ£ºxxxxxx ÓÊÏ䣺xxxxxx@qq.com
ÓåICP±¸2023013149ºÅ
Top