空间计量模型SAC模型matlab代码，sac_panel_FE

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

此代码是J.Paul. Elhorst空间计量经济学：从横截面数据到空间面板此书中sar_panel_FE和sem_panel_FE代码的延续，压缩文件内包含以下3段代码。代码运行同样需在matlab加载了经济学工具箱和Elhorst提供的jplv7工具箱之后运行。 SAC.zip (5.96 KB, 需要: 3 个论坛币) 本附件包括：

• f2_sacpanel.m
• f_sacpanel.m
• sac_panel_FE.m

2019-3-6 11:21:00 上传

需要: 3 个论坛币  [购买]

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：空间计量经济学 空间计量经济 计量经济学 横截面数据 经济学工具

请问知道怎么用吗？就是如何加载自己的数据

dengjhui 发表于 2019-3-6 13:02
请问知道怎么用吗？就是如何加载自己的数据

跟其他模型一样。

A=xlsread('D:\matlab\toolbox\SpatEcon_ren\Files_SLX_paper\cigarette.xls');

W1=xlsread('D:\matlab\toolbox\SpatEcon_ren\Files_SLX_paper\Spat-Sym-US.xls');

%

T=9; % number of time periods

N=31; % number of regions

% row-normalize W

W=normw(W1); % function of LeSage

y=A(:,1); % column number in the data matrix that corresponds to the dependent variable

x=A(:,[2,3,4,5,6,7,8,9]); % column numbers in the data matrix that correspond to the independent variables

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);           %size返回的是行和列

et=ones(T,1);

en=ones(N,1);

% All models include spatial and time period fixed effects

% 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);

可是我就是这样的代码 老是运行不了

结果是这样
>> a new day
未定义与 'char' 类型的输入参数相对应的函数 'a'。

楼主，您好。我想问一下，为啥我用这个无法格式化输出？显示结果：results structure not known by prt_sp function

挥霍你的迷恋 发表于 2019-12-20 21:22
楼主，您好。我想问一下，为啥我用这个无法格式化输出？显示结果：results structure not known by prt_sp  ...

同问，请问你解决了吗？

到这一步应该说明估计完成了，prt_spnew 函数是显示SAC模型估计结果的函数，可查看此函数中有没有关于SAC模型的结果显示的代码，参照其他模型的结果显示代码添加一段代码进去即可。

今古一相接8 发表于 2019-12-21 21:00
同问，请问你解决了吗？

用最新的工具包就可以了，老版本的工具包没有SAC格式化输出的函数