8

我正在尝试在plm具有固定效果的 R 包中运行回归model = 'within'，同时具有聚集的标准错误。使用Cigar来自 的数据集plm，我正在运行：

require(plm)
require(lmtest)
data(Cigar)
model <- plm(price ~ sales + factor(state), model = 'within', data = Cigar)
coeftest(model, vcovHC(model, type = 'HC0', cluster = 'group'))

  Estimate Std. Error t value Pr(>|t|)    
sales  -1.21956    0.21136 -5.7701 9.84e-09

这与我使用 Stata（将雪茄文件编写为 .dta）得到的结果（略有）不同：

use cigar

xtset state year

xtreg price sales, fe vce(cluster state)


price   Coef.   Std. Err.   t   P>t [95% Conf.  Interval]

sales   -1.219563   .2137726    -5.70   0.000   -1.650124   -.7890033

即，标准误和 T 统计量不同。我尝试使用不同的“类型”重新运行 R 代码，但没有一个给出与 Stata 相同的结果。我错过了什么吗？

4

2 回答 2

9

Stata 使用有限样本校正来减少由于集群数量有限而导致的错误向下偏差。它是方差-协方差矩阵的乘法因子，$c=\frac{G}{G-1} \cdot \frac{N-1}{NK}$，其中 G 是组数，N 是观察次数，K 是参数的数量。我认为coeftest只使用 $c'=\frac{N-1}{NK}$ 因为如果我用 c 中第一项的平方来缩放 R 的标准误差，我得到的东西非常接近 Stata 的标准误差：

display 0.21136*(46/(46-1))^(.5)
.21369554

以下是我将如何复制 Stata 在 R 中所做的事情：

require(plm)
require(lmtest)
data(Cigar)
model <- plm(price ~ sales, model = 'within', data = Cigar)
G <- length(unique(Cigar$state))
c <- G/(G - 1)
coeftest(model,c * vcovHC(model, type = "HC1", cluster = "group"))

这产生：

t test of coefficients:

       Estimate Std. Error  t value   Pr(>|t|)    
sales -1.219563   0.213773 -5.70496 1.4319e-08 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

这与 Stata 的误差 0.2137726 和 -5.70 的 t-stat 一致。

这段代码可能并不理想，因为数据中的状态数可能与回归中的状态数不同，但我懒得弄清楚如何获得正确的面板数。

于 2015-10-15T23:41:08.413 回答
8

plmStata 使用在1.5中实现的特定小样本校正。

试试这个：

require(plm)
require(lmtest)
data(Cigar)
model <- plm(price ~ sales + factor(state), model = 'within', data = Cigar)
coeftest(model, function(x) vcovHC(x, type = 'sss'))

这将产生：

t test of coefficients:

      Estimate Std. Error t value  Pr(>|t|)    
sales  -1.2196     0.2137  -5.707 1.415e-08 ***
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

这给出了最多 3 位数的相同 SE 估计值：

x <- coeftest(model, function(x) vcovHC(x, type = 'sss'))
x[ , "Std. Error"]
## [1] 0.2136951
于 2016-03-05T09:19:11.540 回答