Исправленные панелью стандартные ошибки для регрессии данных поперечного сечения временных рядов

Question

Я выполнил регрессию для моих данных поперечного сечения временных рядов. Теперь я хотел бы включить исправленные панелью стандартные ошибки. PLM регрессия работает без проблем. К сожалению, я не уверен, как это сделать, поскольку функция pcse не работает.

library(plm)
reg1 <- plm(Var2 ~ Var1 + Var3 + Var4, data=testdf, index=c("country", "year"), model="random")

library(pcse)
pcse(reg1, country, year, pairwise=FALSE)

Является ли pcse неправильной функцией, и есть ли возможность включить вычисление исправленные панелью стандартные ошибки в регрессии plm уже?

Answer 1

Для пакета pcse требуется объект класса lm. Поскольку вы оценили модель с plm, pcse выдает ошибку. Тем не менее, plm имеет встроенную функцию для оценки стандартных ошибок, исправленных панелями, vcovBK():

library(plm)
library(lmtest)

data("Grunfeld", package = "plm")
plm_fit <- plm(inv ~ value + capital, data = Grunfeld, model = "random", index = c("firm", "year"))
coeftest(plm_fit, vcov=vcovBK(plm_fit, cluster = "time"))
#> 
#> t test of coefficients:
#> 
#>               Estimate Std. Error t value  Pr(>|t|)    
#> (Intercept) -57.834415  30.396538 -1.9027   0.05854 .  
#> value         0.109781   0.016230  6.7643 1.491e-10 ***
#> capital       0.308113   0.024574 12.5380 < 2.2e-16 ***
#> ---
#> Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1

^{Создано в 2020-04-16 пакетом Представить (v0.3.0)}