用MATLAB做空间计量时，panel_effects_sar(results,vnames,W);显示出错

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

clear

A=xlsread('cigarette.xls');

W1=xlsread('Spat-Sym-US.xls');

T=30; % numberof time periods

N=46; % numberof regions

% row-normalizeW

W=normw(W1); %function of LeSage

y=A(:,[3]); % columnnumber in the data matrix that corresponds to the dependent variable

x=A(:,[4,6]); %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);

[nobsK]=size(x);

%----------------------------------------------------------------------------------------

% No fixedeffects + spatially lagged dependent variable

info.lflag=0; %required for exact results

info.model=0;

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

% New routinesto calculate effects estimates

results=sar_panel_FE(y,[xconstantx],W,T,info);

vnames=strvcat('logcit','intercept','logp','logy');

% Print outcoefficient estimates

prt_sp(results,vnames,1);

% Print outeffects estimates

spat_model=0;

direct_indirect_effects_estimates(results,W,spat_model);

panel_effects_sar(results,vnames,W);（一直都是这个代码有问题）

引用了不存在的字段 'hessi'。

出错 panel_effects_sar (line 49)

sigmat = results.hessi - diag(diag(results.hessi)) + diag(diag(abs(results.hessi)));

但是文件夹里有一个名为hessian.m的函数了，我一开始以为函数不一致，我把函数名和变量改成hessi还是有问题！真心不知道问题出在哪里？

另外：panel_effects_sar(results,vnames,W);函数代码出错的地方：

% cheat here to fix the numerical hessian

% Use MCMC to get good results results.xwith

sigmat = results.hessi - diag(diag(results.hessi)) + diag(diag(abs(results.hessi)));

sigmatt = sigmat(1:end-1,1:end-1);

[R,posdef] = chol(sigmatt);

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：results effects Effect MATLAB Result

你好，请问你解决了吗

黄騰騰 发表于 2018-9-13 10:18
你好，请问你解决了吗

解决了，文件包的问题，用LeSage提供的就好

theory1993 发表于 2018-9-13 14:16
解决了，文件包的问题，用LeSage提供的就好

可以告诉我怎么解决的吗，谢谢了

黄騰騰 发表于 2018-9-13 10:18
你好，请问你解决了吗

您好，您解决了吗

陈崇云 发表于 2018-10-25 14:12
您好，您解决了吗

没有，我用R做了

theory1993 发表于 2018-9-13 14:16
解决了，文件包的问题，用LeSage提供的就好

你好，请问可以把你包给我一份吗？