固定效应or随机效应？如何用matlab进行豪斯曼检验

不论是空间计量还是传统意义上的计量经济学，涉及到面板数据都可能需要进行豪斯曼检验，本文将给大家提供一份快捷的豪斯曼检验代码，在直接使用该代码之前，还需要导入一些自己使用的数据源：主要如下：

导入代码为：

A为文件读取，读取名字为‘guiyihua’的CSV文件

第二条是提取y，第三条是提取x，W为空间权重矩阵，不知道权重矩阵如何制作的好兄弟可以参考我上一个文章（该矩阵也可以不需要，如果做的并非空间计量，同时该矩阵最好已经标准化过）。

T为时间，N为数据行，xconstant为全部由1组成的T*N的列向量。在输入完需要的数据之后，得到结果如下：

可以看到需要的数据已经导入。

写入脚本：

xconstant=ones(N*T,1);
[nobs, K]=size(x);
%接下来计算空间固定效应模型的对数似然值；
fprintf('空间固定效应对数似然值');
model=1;
[ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,x,N,T,model);
results=ols(ywith,xwith);
sige=results.sige*((nobs-K)/nobs);
logliksfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid
%计算时间固定效应模型的对数似然值
fprintf('时间固定效应对数似然值');
model=2;
[ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,x,N,T,model);
results=ols(ywith,xwith);
sige=results.sige*((nobs-K)/nobs);
logliktfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid
%计算空间与时间双固定效应模型的对数似然值
fprintf('空间与时间双固定效应对数似然值');
model=3;
[ywith,xwith,meanny,meannx,meanty,meantx]=demean(y,x,N,T,model);
results=ols(ywith,xwith);
sige=results.sige*((nobs-K)/nobs);
loglikstfe=-nobs/2*log(2*pi*sige)-1/(2*sige)*results.resid'*results.resid
%空间联合性显著检验
LR=-2*(logliktfe-loglikstfe);
dof=N;
pvalue=1-chis_prb(LR,dof);
fprintf('空间固定效应联合不显著')
fprintf(1,'=%9.4f,%6d,%9.4f\n',LR,dof,pvalue)
%时间联合性显著检验
LR=-2*(logliksfe-loglikstfe);
dof=T;
pvalue=1-chis_prb(LR,dof);
fprintf('时间固定效应联合不显著')
fprintf(1,'=%9.4f,%6d,%9.4f\n',LR,dof,pvalue)

得到如下结果：

可以看到，所有检验已经完成，模型支持双固定效应模型。

该代码的具体教程，我在blibli上传了制作的视频讲解，如果因我在CSDN中讲解的不到位导致各位看官没有理解的可以去blibli搜索科研虾观看详细讲解。同时如果各位想更多的交流与探讨的，可以加扣扣群：936269384，同时我也将代码上传至群文件。

固定效应or随机效应？如何用matlab进行豪斯曼检验相关推荐

豪斯曼检验matlab,面板工具变量与内生性操作及应用（固定随机、豪斯曼检验、过度识别检验、预测等）...
原标题:面板工具变量与内生性操作及应用(固定随机.豪斯曼检验.过度识别检验.预测等)
豪斯曼检验matlab,豪斯曼检验、空间面板模型选择等问题
各位大神,本人空间计量小白,最近看了姜磊老师的书,自己也进行了实证分析.有一些问题不太懂: 1.空间面板数据模型分为混合空间面板数据模型.空间固定效应空间面板数据模型.时间固定效应空间面板模型.空间和 ...
时间固定效应和个体固定效应的选择_互助问答第31期：固定效应与随机效应选择和面板数据处理...
样本描述:各位老师好,我的论文采用的是微观面板非平衡数据(合并了3波数据,总观测值6万左右),每波观测之间约有20%的样本不同(约10%的样本流失,10%的新样本补入),因变量是连续变量,核心自变量是 ...
随机效应估算与固定效应估算_一般混合线性模型固定效应、随机效应与另一随机向量的联合估计...
一般混合线性模型固定效应.随机效应与另一随机向量的联合估计周永正 [期刊名称] <数学的实践与认识> [年 ( 卷 ), 期] 2011(041)019 [摘要] 讨论一般混合线性模型 ...
面板数据回归模型（固定效应、随机效应、混合回归、变系数）、面板数据AR、VAR模型
Stata——固定效应模型、随机效应模型、混合效应模型（区别、实例）
目录 1 固定效应模型概念(Fixed Effects Model) 1.1 stata命令 1.1.1 LSDV法(Least squares dummy variable) 1.1.2 固定效应 ...
R语言中固定与随机效应Meta分析 - 效率和置信区间覆盖
Meta分析是综合现有证据的重要工具. 我们围绕Meta技术进行一些咨询,帮助客户解决独特的业务问题.它通常用于医学和临床环境中,以评估关于治疗或暴露对感兴趣结果的影响的现有证据.基本思想是将先前研究 ...
R语言线性混合效应模型（固定效应随机效应）和交互可视化3案例
最近我们被客户要求撰写关于线性混合效应模型的研究报告,包括一些图形和统计输出. 视频:线性混合效应模型(LMM,Linear Mixed Models)和R语言实现案例线性混合效应模型(LMM,Li ...
R语言非线性混合效应 NLME模型(固定效应随机效应)对抗哮喘药物茶碱动力学研究
最近我们被客户要求撰写关于非线性混合效应 NLME模型的研究报告,包括一些图形和统计输出. 相关视频:线性混合效应模型(LMM,Linear Mixed Models)和R语言实现线性混合效应模型( ...
随机效应估算与固定效应估算_固定效应模型的估计原理说明.pdf
固定效应模型的估计原理说明在面板数据线性回归模型中,如果对于不同的截面或不同的时间序列,只是模型的截距项是不同的,而模型的斜率系数是相同的,则称此模型为固定效应模型.固定效应模型分为三类: 1. ...

固定效应or随机效应？如何用matlab进行豪斯曼检验

固定效应or随机效应？如何用matlab进行豪斯曼检验相关推荐

最新文章

热门文章