Driscoll-Kraay 标准误

Driscoll-Kraay 标准误是什么？

Driscoll-Kraay 标准误处理的是面板误差中的横截面相关、异方差和序列相关。它不改变回归系数，只改变推断方式，让标准误对更复杂的误差结构稳健。

先用一句话理解

把 Driscoll-Kraay 想成承认不同地区或企业会一起受冲击。普通标准误假设误差相互独立，但金融危机、宏观政策或行业冲击会让很多单位同时波动。

什么时候用

把它和相近方法分清楚

核心直觉和模型公式

横截面平均上的 HAC

Driscoll-Kraay 可以理解为对横截面平均得分做 HAC 修正：

它修正的是协方差矩阵，因此系数点估计通常不变，变化的是标准误和显著性。

公式里的符号怎么读

专业使用口径

专业 DK 报告要明确它只是标准误替代口径。

一步一步做：从问题到结果的五步

结果怎么读：先读口径，再读数字

本页案例读表时，先记住这个口径：kernel = Driscoll-Kraay (Bartlett)。

初学者最容易错在哪里

先看这个案例的结论

• kernel = Driscoll-Kraay (Bartlett)。
• lag = 2。
• 状态 = xtscc 失败 (需 ssc install xtscc) rc=199。
• 这些数字来自页面里的结果表；写论文时先解释数值含义，再讨论理论含义。

变量要求和案例口径

实际代码

下面是本页案例里的最小 Stata 代码。它的作用是帮助你看清变量怎么进入模型、结果文件怎么生成，而不是要求你在网页里手动运行。

log using "/root/workspace/empirical-wizard/workspace/d2b9f5c9/analysis.log", replace text
global JOB_DIR "/root/workspace/empirical-wizard/workspace/d2b9f5c9"
set more off
adopath + "/root/ado/plus"
global DATA_PATH "/root/workspace/empirical-wizard/workspace/test_e2e/csmar_innovation.csv"
import delimited "/root/workspace/empirical-wizard/workspace/test_e2e/csmar_innovation.csv", clear case(preserve)
capture confirm global JOB_DIR
if _rc global JOB_DIR "."
* 自动去除完全重复行（同列同值），避免 N 虚增与 xtset 失败
quietly duplicates drop
local idvar ""
local timevar ""
capture confirm variable stkcd
if !_rc {
    capture confirm numeric variable stkcd
    if _rc {
        tempvar __ewiz_id
        capture encode stkcd, gen(`__ewiz_id')
        if !_rc local idvar "`__ewiz_id'"
    }
    else {
        local idvar "stkcd"
    }
}
else {
    di as text "面板ID变量不存在，跳过 xtset ID：stkcd"
}
capture confirm variable year
if !_rc {
    capture confirm numeric variable year
    if _rc {
        tempvar __ewiz_time
        capture encode year, gen(`__ewiz_time')
        if !_rc local timevar "`__ewiz_time'"
    }
    else {
        local timevar "year"
    }
}
else {
    di as text "时间变量不存在，跳过 xtset time：year"
}
if "`idvar'" != "" & "`timevar'" != "" {
    capture xtset `idvar' `timevar'
}

* ── Driscoll-Kraay 标准误 (xtscc) ──
capture noisily xtscc patent_count dfi_index roa lev size growth cashflow tobinq top1 dual board indep soe age, lag(2) fe
local _rc = _rc
tempname fh
capture file close `fh'
file open `fh' using "$JOB_DIR/driscoll_kraay_results.csv", write replace
file write `fh' "项目,值" _n
file write `fh' "因变量,patent_count" _n
file write `fh' "核心解释变量,dfi_index" _n
file write `fh' "kernel,Driscoll-Kraay (Bartlett)" _n
file write `fh' "lag,2" _n
if `_rc' == 0 {
    capture local _b : display %9.4f _b[dfi_index]
    capture local _se : display %9.4f _se[dfi_index]
    capture local _t = _b[dfi_index]/_se[dfi_index]
    capture local _t_s : display %9.3f `_t'
    capture local _p : display %9.4f 2*ttail(e(df_r), abs(`_t'))
    capture local _n = e(N)
    file write `fh' "系数,`_b'" _n
    file write `fh' "DK 标准误,`_se'" _n
    file write `fh' "t 值,`_t_s'" _n
    file write `fh' "p 值,`_p'" _n
    file write `fh' "N,`_n'" _n
}
else {
    file write `fh' "状态,xtscc 失败 (需 ssc install xtscc) rc=`_rc'" _n
}
file close `fh'
di "Driscoll-Kraay 完成"
log close

实际输出表

这张表是本方法页配套案例的输出结果，用来把前面的公式、变量口径和代码对应到实际数字。

案例图

这是一张由同一份案例数据生成的页面内诊断图。

论文里怎么写

本文在共用企业面板样本上报告Driscoll-Kraay 标准误，核心输出见 driscoll_kraay_results.csv。结果解释时同时关注样本口径、变量构造、系数方向、标准误和适用前提，避免只凭单个 p 值完成方法选择。

检查清单

• 先确认面板存在横截面相关或序列相关风险，DK 是标准误修正。
• 说明时间维度是否足够支撑 Driscoll-Kraay 推断。
• DK 改变的是标准误和显著性，不改变回归系数本身。

返回教程库 · 查看共用案例