Elhorst编写的程序与Eviews跑出来的结果不一致

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

本人用elhorst的程序包demoLMsarsem_panel跑出来了教材《空间计量经济学-从横截面数据到空间面板》的结果，其中包括一个个体固定效应模型和一个时间固定效应模型。但是我用同样的数据用Eviews跑个体固定效应模型后发现给出了大相径庭的结果，为什么会这样呢？这是Eviews的结果


这是Matlab的结果


只有普通OLS两者的结果差不多，其他的大相径庭呢~

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：Elhorst EVIEWS Views Horst Eview

有兴趣的朋友可以自己试一试~所以我也不清楚哪出了问题……

一个有常数，一个没有常数项

magicsun 发表于 2017-10-5 10:11
一个有常数，一个没有常数项

那个常数只是matlab没有报告出来吧

应该都在里面

magicsun 发表于 2017-10-6 09:27
应该都在里面

都在里面是指？

就是所有的参数和t值都列出来吧？

magicsun 发表于 2017-10-6 09:57
就是所有的参数和t值都列出来吧？

都列出来了~在elhorst的SDM程序里是可以选择是否把固定效应的截距项显示出来，但是LM检验这里没法选择，但是看代码它是有做固定效应的，就是不知道问题出在哪了。（是的其实我就是看不懂固定效应的算法……）

magicsun 发表于 2017-10-6 09:57
就是所有的参数和t值都列出来吧？

1. fprintf(1,'****************************OLS****************************\n');
   
2. results=ols(y,[xconstant x]);
   
3. vnames=strvcat('ln(lp)','intercept','p','ln(e)','ln(m)','h');%备用为:'ln(d)','ln(d)''city','xx'
   
4. prt_reg(results,vnames,1);
   
5. sige=results.sige*((nobs-K)/nobs);
   
6. loglikols=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid
   
7. % The (robust)LM tests developed by Elhorst
   
8. LMsarsem_panel(results,W,y,[xconstant x]); % (Robust) LM tests
   
9. % The lm tests developed by Donald Lacombe
   
10. % see http://www.rri.wvu.edu/lacombe/~lacombe.htm
    
11. % ----------------------------------------------------------------------------------------
    
12. % spatial fixed effects + (robust) LM tests for spatial lag and spatial error model
    
13. % fixed effects, within estimator
    
14. % demeaning of the y and x variables
    
15. fprintf(1,'****************************个体固定效应OLS****************************\n');
    
16. model=1;
    
17. [ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,x,N,T,model);
    
18. results=ols(ywith,xwith);
    
19. vnames=strvcat('ln(lp)','p','ln(e)','ln(m)','h'); % should be changed if x is changed
    
20. prt_reg(results,vnames);
    
21. sfe=meanny-meannx*results.beta; % including the constant term
    
22. yme = y - mean(y);
    
23. et=ones(T,1);
    
24. error=y-kron(et,sfe)-x*results.beta;

复制代码

代码如上