我正在尝试对意大利各省的一些面板数据进行Pesaran CD测试。对于每个单位,我有22年,每个单位分为四分之一(因此每个单位有92个观察结果)。
变量“lag_lingua”在一个省的每个季度和年份具有相同的值,但它在省之间变化。
这是我使用的代码:
library(plm)
A <- data[-2]
AA <- pdata.frame(A, index = c("id","time"))
pcdtest(AA$lag_lingua)
我得到这个结果:
Pesaran CD test for cross-sectional dependence in panels
data:
AA$lag_lingua
z = NA, p-value = NA
alternative hypothesis: cross-sectional dependence
Warning message:
In cor(wideres, use = "pairwise.complete.obs") :
la deviazione standard è zero.
我可以想象,这是由于变量只在单元之间变化,而不是在单元“内部”变化,但我如何解决这个问题并测试CD?
1条答案
按热度按时间r6vfmomb1#
您关于变量是时间常数(仅随单位变化)的假设是正确的。测试使用单位之间的相关性,如果这些单位是常数,则不定义相关性,因此测试返回
NA
。为了说明,让我们创建一个名为
var
的非时变变量并运行测试:pcdtest
中的所有测试方法均基于不同单位(rho)的观测值之间的相关性。所以,基本上,调用cor()
,输入几对常量向量,所以类似于这样:我们看到
NA
和警告,如pcdtest
中的cor()
函数,用于这种数据输入。pcdtest
中最简单的相关性实现是通过设置test = "rho"
给出的,pcdtest
中的所有其他统计方法都基于rho值。当变量在单位内不显示随时间的变化时,您根本无法测试横截面依赖性。