经济学家

 找回密码
 快速注册

QQ登录

只需一步,快速开始

查看: 1143|回复: 4

空间计量SLM模型的问题

[复制链接]
发表于 2015-3-9 08:30 | 显示全部楼层 |阅读模式
10金币

求教大侠,以下是我用matlabSLM模型计量的结果,请问logL值是否就是那个“NaN"的值?如果是的那为什么会这样?

Pooled model with spatially lagged dependent variable, spatial and time period fixed effects
Dependent Variable =               y  
R-squared          =    0.8504   
corr-squared       =    0.4854   
sigma^2            =    0.1130
Nobs,Nvar,#FE      =    155,     7,    41  
log-likelihood     =              NaN
# of iterations    =      1   
min and max rho    =   -1.0000,   1.0000
total time in secs =    0.0250
time for optimiz   =    0.0060
time for lndet     =    0.0100
time for t-stats   =    0.0010
No lndet approximation used
***************************************************************
Variable        Coefficient  Asymptot t-stat    z-probability
x1                 0.428738         0.278092         0.780942
x2                -0.685447        -2.101475         0.035599
x3                 4.047355         1.065698         0.286560
x4                -1.387638        -1.207144         0.227377
x5                 6.405075         2.773594         0.005544
x6                 0.951532         9.305283         0.000000
W*dep.var.        -0.183058        -1.382462         0.166830


评分

参与人数 1银币 +5 收起 理由
cbw153053970 + 5 鼓励交流!

查看全部评分

 楼主| 发表于 2015-3-9 16:37 | 显示全部楼层
自己顶顶,求大神解疑啊!!
回复

使用道具 举报

发表于 2015-3-10 10:43 | 显示全部楼层
虽然不是很懂,但是是空间自回归模型么!求教大神是怎么在matlab中做的,贴点代码我看看
回复

使用道具 举报

 楼主| 发表于 2015-3-10 16:04 | 显示全部楼层
cbw153053970 发表于 2015-3-10 10:43
虽然不是很懂,但是是空间自回归模型么!求教大神是怎么在matlab中做的,贴点代码我看看

代码其实是公开的呀。这是固定效应中的无固定和空间固定的代码,都类似。。
T=5; % number of time periods
N=31; % number of regions
% row-normalize W
W=normw(ww); % function of LeSage
y=pc(:,[4]); % column number in the data matrix that corresponds to the dependent variable
x=pc(:,[12,14,15,18,20]); % 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);
% ----------------------------------------------------------------------------------------
% No fixed effects + spatially lagged dependent variable
info.lflag=0; % required for exact results
info.model=0;
info.fe=0; % Do not print intercept and fixed effects; use info.fe=1 to turn on
% New routines to calculate effects estimates
results=sar_panel_FE(y,[xconstant x],W,T,info);
% vnames=strvcat('lny','intercept','lnAgdp','lnAGr','lnPni','Pcl','lnEha','lnDlc');
vnames=strvcat('lny','intercept','lnAGr','lnPni','Pcl','lnEha','lnDlc');
% Print out coefficient estimates
prt_sp(results,vnames,1);
% ----------------------------------------------------------------------------------------
% Spatial fixed effects + spatially lagged dependent variable
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
% New routines to calculate effects estimates
results=sar_panel_FE(y,x,W,T,info);
% vnames=strvcat('lny','lnAgdp','lnAGr','lnPni','Pcl','lnEha','lnDlc');
vnames=strvcat('lny','lnAGr','lnPni','Pcl','lnEha','lnDlc');
% Print out coefficient estimates
prt_sp(results,vnames,1);
回复

使用道具 举报

 楼主| 发表于 2015-3-10 19:01 | 显示全部楼层
自己顶顶啊,求大神指点啊!!!
回复

使用道具 举报

关闭

广告上一条 /1 下一条

小黑屋|手机版|南康家具|家具批发|整体橱柜|经济学家 赣ICP备11002034号

GMT+8, 2017-6-26 04:46 , Processed in 0.536585 second(s), 23 queries .

Powered by Discuz! X3.3

© 2001-2017 Comsenz Inc.

快速回复 返回顶部 返回列表