面板数据的Pesaran CD检验

lnxxn5zx  于 2023-05-26  发布在  其他
关注(0)|答案(1)|浏览(736)

我正在尝试对意大利各省的一些面板数据进行Pesaran CD测试。对于每个单位,我有22年,每个单位分为四分之一(因此每个单位有92个观察结果)。
变量“lag_lingua”在一个省的每个季度和年份具有相同的值,但它在省之间变化。
这是我使用的代码:

  1. library(plm)
  2. A <- data[-2]
  3. AA <- pdata.frame(A, index = c("id","time"))
  4. pcdtest(AA$lag_lingua)

我得到这个结果:

  1. Pesaran CD test for cross-sectional dependence in panels
  2. data:
  3. AA$lag_lingua
  4. z = NA, p-value = NA
  5. alternative hypothesis: cross-sectional dependence
  6. Warning message:
  7. In cor(wideres, use = "pairwise.complete.obs") :
  8. la deviazione standard è zero.

我可以想象,这是由于变量只在单元之间变化,而不是在单元“内部”变化,但我如何解决这个问题并测试CD?

r6vfmomb

r6vfmomb1#

您关于变量是时间常数(仅随单位变化)的假设是正确的。测试使用单位之间的相关性,如果这些单位是常数,则不定义相关性,因此测试返回NA
为了说明,让我们创建一个名为var的非时变变量并运行测试:

  1. library(plm)
  2. data("Grunfeld", package = "plm")
  3. # create time-constant variable 'var' which varies per individual
  4. Grunfeld$var <- Grunfeld$firm
  5. pGrun <- pdata.frame(Grunfeld)
  6. pcdtest(pGrun$var) # same as pcdtest(var ~ 1, data = pGrun)
  7. #> Warning in cor(wideres, use = "pairwise.complete.obs"): the standard deviation
  8. #> is zero
  9. #>
  10. #> Pesaran CD test for cross-sectional dependence in panels
  11. #>
  12. #> data: pGrun$var
  13. #> z = NA, p-value = NA
  14. #> alternative hypothesis: cross-sectional dependence
  15. # can check for variation across unit and time with pvar():
  16. pvar(pGrun$var)
  17. #> no time variation: var

pcdtest中的所有测试方法均基于不同单位(rho)的观测值之间的相关性。所以,基本上,调用cor(),输入几对常量向量,所以类似于这样:

  1. cor(c(1, 1), c(2, 2))
  2. #> Warning in cor(c(1, 1), c(2, 2)): the standard deviation is zero
  3. #> [1] NA

我们看到NA和警告,如pcdtest中的cor()函数,用于这种数据输入。pcdtest中最简单的相关性实现是通过设置test = "rho"给出的,pcdtest中的所有其他统计方法都基于rho值。
当变量在单位内不显示随时间的变化时,您根本无法测试横截面依赖性。

展开查看全部

相关问题