豪斯曼检验matlab,豪斯曼检验、空间面板模型选择等问题
各位大神,本人空间计量小白,最近看了姜磊老师的书,自己也进行了实证分析。有一些问题不太懂:
1、空间面板数据模型分为混合空间面板数据模型、空间固定效应空间面板数据模型、时间固定效应空间面板模型、空间和时间双固定效应的空间面板数据模型、随机效应的空间面板数据模型,而豪斯曼检验可以用来判断到底是选择固定效应模型还是随机效应模型。那豪斯曼统计量到底是如何来构造呢?是否需要按照空间时间双固定效应模型vs.随机效应模型、空间固定效应模型vs.空间随机效应模型、时间固定效应模型vs.时间随机效应模型分别来构造?如果是空间随机时间固定效应模型,那么豪斯曼检验统计量是不是按照空间和时间双固定效应模型vs.随机效应模型来构造?
2.我在应用自己的数据,运行MATLAB中的demopanelscompare代码时,到最后一个模型“空间随机时间固定效应模型”,就会出现错误:
部分代码:
% Spatial random effects and time period fixed effects + spatially lagged dependent variable + spatially
% independent variables
[ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,[x wx],N,T,2); % 2=time dummies
info.model=1;%不清楚此处为什么还是选择1,空间固定效应模型不才是1吗?
results=sar_panel_RE(ywith,xwith,W,T,info);
prt_sp(results,vnames,1);
% Print out effects estimates
spat_model=1;
direct_indirect_effects_estimates(results,W,spat_model);
panel_effects_sdm(results,vnames,W);
% Wald test for spatial lag model
btemp=results.parm(1:2*K+2);
varcov=results.cov(1:2*K+2,1:2*K+2);
Rafg=zeros(K,2*K+2);
for k=1:K
Rafg(k,K+k)=1; % R(1,3)=0 and R(2,4)=0;
end
Wald_spatial_lag=(Rafg*btemp)'*inv(Rafg*varcov*Rafg')*Rafg*btemp
prob_spatial_lag= 1-chis_cdf (Wald_spatial_lag, K) % probability greater than 0.05 points to insignificance
% LR test spatial Durbin model against spatial lag model (requires
% estimation results of the spatial lag model to be available)
resultssar=sar_panel_RE(ywith,xwith(:,1:K),W,T,info);
LR_spatial_lag=-2*(resultssar.lik-results.lik)
prob_spatial_lag=1-chis_cdf (LR_spatial_lag,K) % probability greater than 0.05 points to insignificance
% Wald test for spatial error model
R=zeros(K,1);
for k=1:K
R(k)=btemp(2*K+1)*btemp(k)+btemp(K+k); % k changed in 1, 7/12/2010
% R(1)=btemp(5)*btemp(1)+btemp(3);
% R(2)=btemp(5)*btemp(2)+btemp(4);
end
Rafg=zeros(K,2*K+2);
for k=1:K
Rafg(k,k) =btemp(2*K+1); % k changed in 1, 7/12/2010
Rafg(k,K+k) =1;
Rafg(k,2*K+1)=btemp(k);
% Rafg(1,1)=btemp(5);Rafg(1,3)=1;Rafg(1,5)=btemp(1);
% Rafg(2,2)=btemp(5);Rafg(2,4)=1;Rafg(2,5)=btemp(2);
end
Wald_spatial_error=R'*inv(Rafg*varcov*Rafg')*R
prob_spatial_error= 1-chis_cdf (Wald_spatial_error,K) % probability greater than 0.05 points to insignificance
% LR test spatial Durbin model against spatial error model (requires
% estimation results of the spatial error model to be available
resultssem=sem_panel_RE(ywith,xwith(:,1:K),W,T,info);
LR_spatial_error=-2*(resultssem.lik-results.lik)
prob_spatial_error=1-chis_cdf (LR_spatial_error,K) % probability greater than 0.05 points to insignificance
%运行到这一行时,就会提示错误。
错误使用 gammainc
输入必须为实数完全双精度值或单精度值。
出错 gamm_cdf (line 25)
cdf = gammainc(x,a);
出错 chis_cdf (line 28)
F = gamm_cdf(x/2,a*0.5);
出错 demopanelscompareorigin (line 310)
prob_spatial_error=1-chis_cdf (LR_spatial_error,K) % probability greater than 0.05 points to insignificance
我往前找错,发现出错的原因是resultssem.lik数值为:1031.40668903056 + 47.1238898038469i;
不知道为什么会出现这样的问题?该代码中前面还运行了无固定效应、空间固定效应、时间固定效应和时空双固定效应,都没有报错。用相同的代码,运行书中的数据,也没有问题。难道是我的数据存在什么问题吗?
3.我在运行SDM模型的几种固定效应、时间效应时,发现时间固定效应的相关系数R2(0.47)明显高于空间固定效应(0.10)、时空固定效应(0.12),而最大似然值时间固定效应模型(679.8)却显著低于空间固定效应(1121.4)和时空固定效应(1142.4),同时,空间固定效应和时空固定效应的变量显著性也好于时间固定效应。为什么不同指标差别那么大?可能是什么原因造成的呢?
豪斯曼检验matlab,豪斯曼检验、空间面板模型选择等问题相关推荐
- matlab独立t样本检验,matlab与单样本t检验
R语言,SPSS几EXCEL解决单样本T检验问题.母猪的怀孕期未114d,今抽测... spss单样本t检验_数学_自然科学_专业资料. 文档贡献者 920450863 贡献于2012-12-21 1 ...
- wald检验matlab,Wald检验和Wald-F检验一样么?高手帮帮忙
Wald检验只需要估计无约束模型 Wald-F检验需要估计有约束.跟无约束的两个模型 下面贴出matlab jplv7工具箱中的waldf函数,在\jplv7\regress目录下 function ...
- matlab likelihood,Matlab做空间面板模型log-likelihood为NAN怎么办
初学者,程序完全照抄,求大神纠错 ,SAR固定效应是程序包弄的不知道哪里出了问题... load 'pc.mat' pc; load 'wl.mat' wl; load 'pc1.mat' pc1; ...
- matlab2016a贝叶斯工具箱,matlab贝叶斯工具箱
matlab贝叶斯工具箱是由Kevin Murphy开发的编程学习工具,该工具采用MATLAB语言编制,可实现贝叶斯网络结构学习.参数学习.推理和构建贝叶斯分类器,此工具箱在贝叶斯学习编程方面非常灵活 ...
- 建立空间计量模型时对数据进行LM检验的Stata代码
通过LM检验可以确定需要建立的空间计量模型. use "(权重矩阵路路径)\W.dta" //注意保存成dta格式 //扩大矩阵,time后括号内的是倍数 spcs2xt (第一个 ...
- matlab贝叶斯优化工具箱_经济学人的神器——BEAR(贝叶斯估计、分析和回归工具包)...
武林至尊,宝刀屠龙, 号令天下,莫敢不从, 倚天不出,谁与争锋. --金庸·<倚天屠龙记> Bayesian Estimation, Analysis and Regression(简写为 ...
- 安装matlab贝叶斯网络工具箱
安装matlab贝叶斯网络工具箱 下载工具箱 安装工具箱 安装matlab贝叶斯网络工具箱 1. 下载工具箱 下载链接: 官方下载:http://www.cs.ubc.ca/~murphyk/Soft ...
- 变分推断 | MATLAB实现VBMC变分贝叶斯蒙特卡洛模拟的贝叶斯推断
变分推断 | MATLAB实现变分贝叶斯蒙特卡洛模拟的贝叶斯推断 目录 变分推断 | MATLAB实现变分贝叶斯蒙特卡洛模拟的贝叶斯推断 效果一览 基本介绍 研究内容 模型描述 模型设计 参考资料 效 ...
- 【图像重建】基于matlab布雷格曼迭代算法集合ART算法CT图像重建【含Matlab源码 1905期】
一.获取代码方式 获取代码方式1: 完整代码已上传我的资源:[图像重建]基于matlab布雷格曼迭代算法集合ART算法CT图像重建[含Matlab源码 1905期] 获取代码方式2: 通过订阅紫极神光 ...
最新文章
- linux vscode vim,VSCode Vim进阶操作
- 项目由于装运点不同交货拆分解决方案
- C#简易通讯录的开发试题
- 【C++11新特性】 C++11智能指针之weak_ptr
- scp和sftp常用操作
- 后端的日期类型赋值前端表单_Spring Boot实践--前端字符串日期自动转换成后台date类型。...
- 用c语言写的电话簿的程序,用C语言散列表实现电话薄
- Apache与Subversion的简单配置(Apache SVN版本控制)
- 61.Linux/Unix 系统编程手册(下) -- SOCKET: 高级主题
- Gson解析json数据
- 2021-07-02
- jpa原生query_SpringDataJpa使用原生sql的小坑
- TDMA解调接收机IP core
- ubuntu桌面美化mac_教程:为你的linux桌面带来Mac OS Mojave的体验
- 史上最全推广小程序实操方法
- 如何购买腾讯云服务器?腾讯云服务器购买教程
- python画笑脸表情_用Matplotlib,妈妈再也不担心我没有表情包斗图了
- JNI的替代者—使用JNA访问Java外部功能接口 DLL
- 如何开始学剪辑?零基础详细解析
- 用 22 张照片打开 23 年
热门文章
- 挪威访学2:SOLA UIS 学生宿舍
- 20层的试炼html5,Vue.js-02:新手村的试炼 - 新世界的武器(指令)
- xp系统outlook突然无法连接服务器,outlook邮箱打不开是什么原因
- JVM各内存区域存放内容
- self-supervised text erasing with controllable image synthesis
- gnuplot绘图程序中对线型(linetype)、点型(pointtype)、线条宽度(linewidth)、点大小(pointsize)、图样
- 解释一下label中的写法:plt.plot(t, sig, b-, linewidth=2, label=r$\sigma(t) = \frac{1}{1 + e^{-t}}$)...
- CentOS 7 安装 MariaDB 10.1
- Envoy架构理解--理解xDS/Listener/Cluster/Router/Filter
- matlab 去条带噪声,一种图像条带噪声及坏线消除方法