|
非常感谢楼主的无私奉献,前几个模型都跑出来了,只是在跑下面的两个程序时:
% SAC / SARAR / Clifford-Ord/ Kelejian-Prucha model
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
info.bc=1;
results=sac_panel_FE(y,x,W,T,info);
vnames=strvcat('logcit','logp','logy');
prt_spnew(results,vnames,1);
% Print out effects estimates
spat_model=0;
direct_indirect_effects_estimates(results,W,spat_model);
% ----------------------------------------------------------------------------------------
% Full model
info.model=3;
info.fe=0; % Do not print intercept and fixed effects; use info.fe=1 to turn on
info.bc=1;
results=sac_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);
遇到到问题:??? Error using ==> prt_spnew
results structure not known by prt_sp function
请问大神,这个问题怎么解决?
|