Stata：面板聚类标准误-自动确定最优聚类层级和数量-xtregcluster

发布时间：2021-11-06 阅读 546

Stata连享会主页 || 视频 || 推文 || 知乎 || Bilibili 站

温馨提示： 定期清理浏览器缓存，可以获得最佳浏览体验。

New！ lianxh 命令发布了：
随时搜索推文、Stata 资源。安装：
. ssc install lianxh
详情参见帮助文件 (有惊喜)：
. help lianxh
连享会新命令：cnssc, ihelp, rdbalance, gitee, installpkg

✌ 课程详情： https://gitee.com/lianxh/Course

⛳ 课程主页： https://gitee.com/lianxh/Course

⛳ Stata 系列推文：

全部 | Stata入门 | Stata教程 | Stata资源 | Stata命令
计量专题 | 论文写作 | 数据分享 | 专题课程
结果输出 | Stata绘图 | 数据处理 | Stata程序
回归分析 | 面板数据 | 交乘项-调节 | IV-GMM
内生性-因果推断 | 倍分法DID | 断点回归RDD | PSM-Matching | 合成控制法
Probit-Logit | 时间序列 | 空间计量 | 分位数回归 | 生存分析 | SFA-DEA
文本分析-爬虫 | Python-R-Matlab | 机器学习
Markdown | 工具软件 | 其它

☝ PDF下载 - 推文合集

作者：吕卓阳 (厦门大学)
邮箱：lvzy20@163.com

1. 研究背景

在针对上市公司的研究中，大家经常在如下两种估计命令之间摇摆不定：第一种是我们控制了时间固定效应、以及公司或行业固定效应，进行面板回归：

. xtreg y x1 x2 i.year, fe 
. reg   y x1 x2 i.year i.industry

第二种是在第一种的基础上，考虑标准误的聚类：

. xtreg y x1 x2 i.year, fe vce(cluster firmID)
. reg   y x1 x2 i.year i.industry, vce(cluster industry)

回归的方程如下：

y_{i t} = α_{i} + λ_{t} + x_{i t} β + ε_{i t}

y_{i t} = θ_{j} + λ_{t} + x_{i t} β + ε_{i t}

两种回归的方式不同在于第二种回归考虑了标准误的聚类问题，即假设标准误在公司或行业的层面聚类。这种假设相比于假定标准误差满足 $i . i . d$ 的基础上更加合理，因为实际问题中标准误在某一层面上允许其相关是合理的，譬如行业内的标准误存在相关性，而行业间标准误不相关。

我们平时常常是基于经济学的意义预先设定聚类的数量和分块大小，鲜少从统计意义上去估计最优聚类方法，本文将为大家介绍一种新的「聚类分析」方法及其 Stata 实现命令 xtregcluster。其优点在于不需要先验地假设聚类的数量和分块的大小，而是根据数据建立惩罚函数来降低总残差平方和，本着「让数据说话」的理念寻找最优聚类数量和分块大小。

2. 方法介绍

对于标准的面板数据分析，一般施加所有个体具有相同斜率系数的限制，个体不可观测的异质性取决于个体不随时间改变的效应，但从理论与实证上都很难证明该限制 (Burnside，1996；Baltagi 和 Griffin，1997；Pesaran 等，1999)。

由此，Sarafidis 和 Weber (2015) 开发了一种回归聚类方法，将个体分为不同的聚类。这种聚类不依赖于参数分布的先验假设。譬如，我们想探究学生成绩的影响因素，若在计算标准误时，选择在班级层面上进行聚类调整，则意味着我们假定同一个班级的同学的成绩更具有相关性 (可能源于班风、班主任个人特征等因素)。但这类假设往往是基于研究者的理论分析给出的假设，而非通过数据分析得到的。

Sarafidis 开发的回归聚类方法是基于数据确定的，假设聚类数量与分块的大小都是未知，基于最小化总残差平方和的最优化方法并设置严格增加聚类数量的惩罚函数，从而将个体分为不同的聚类。在每个聚类中，斜率系数是同质的，而在组间斜率的系数是异质的。在没有参数结构的先验信息的情况下，xtregcluster 可以用于探索斜率参数的异质性，适用于分析大样本个体和短时间序列样本。此外，从统计意义上也可以检验单个实体先验分类是否最优。

3. 命令介绍

3.1 理论部分

考虑如下面板回归模型：

y_{i_{ω} t} = β_{ω}^{'} x_{i_{ω} t} + ϵ_{i_{ω} t}

其中， $y_{i_{ω} t}$ 表示个体 $i$ 在 $w$ 标准误聚类情形下 $t$ 时刻被解释变量的观测值， $X_{i_{ω} t}$ 表示 $K \times 1$ 阶在标准误聚类情形下的解释变量， $β_{ω}$ 为 $K \times 1$ 阶的估计系数，并且由于在 $w$ 情形下聚类，系数在组内相同但在组间不同。残差项可以分解如下：

ϵ_{i_{ω} t} = η_{i_{ω}} + τ_{t} + e_{i_{ω} t}

其中， $e_{i_{ω} t}$ 为个体不可观测的残差项， $τ_{t}$ 表示时间固定效应， $e_{i_{ω} t}$ 可能与 $%热门资讯 ⛵ 连享会 - 课程主页 12/05 计量专题 35935 ⏩ 因果推断专题-RDD-DID-IV-合成控制 11/10 专题课程 950 连享会主页-推文列表-按类别 10/10 专题课程 33516 ☝ Stata Journal 来啦！随便看 (2001-2020) 02/02 Stata资源 11012 Stata：线性趋势检验-lintrend 12/07 Stata命令 24 LaTeX-如何使用Overleaf生成PPT 12/05 工具软件 63 经典计量经济学教材推荐 12/04 计量专题 331 ⏩ 连享会：2022 Stata寒假班-全新上线 12/03 计量专题 108 Stata：自选择偏误之双栏模型简介-(Double-hurdle-model) 12/03 内生性-因果推断 67 Stata结果输出：addest自己定制输出的统计量 12/02 结果输出 117关于我们友情链接联系我们版权所有：连享会技术支持：穗源网络科技备案号：陕ICP备19021402号-1$(document).ready(function(){
var _w = $("#article_box"). width();

$("#test-editormd-view2 img").each(function(i){
var img = $(this);
var realWidth;
var realHeight;
$("<img/>").attr("src", $(img).attr("src")).on('load',function() {
realWidth = this.width;
realHeight = this.height;
if(realWidth >= _w){
$(img).css("width","100%").css("height","auto");
}else{
$(img).css("width",realWidth+'px').css("height",realHeight+'px');
}

})
});

$.get("/index/api/collection_state",{id:"31be9ab09c5ff"},function(e){
if(e.code == 1){
$("#collection").text('取消收藏');
}else{
$("#collection").text('收藏文章');
}
});

});

function collection(){
$.get("/index/api/collection",{id:"31be9ab09c5ff"},function(e){
if(e.code == 1){
alert('收藏成功');
$("#collection").text('取消收藏');
return;
}else if (e.code==3) {
alert('没有登录！请登录后进行收藏');
return;
}else if (e.code==4) {
alert('取消收藏成功');
$("#collection").text('收藏文章');
return;
}else if (e.code==5) {
alert('取消收藏失败');
return;
}else{
alert('收藏失败');
}

});
}$