我正在尝试使用plm软件包运行固定效果回归。回归代码如下:
fixed = plm(hp~crime,index=c('year','country'),data=data,model='within')
返回以下错误代码:
pdim.default(index [[1]],index [[2]]中的错误:重复对(id-time)
我已经搜索了网络,包括stackoverflow。我了解的是,plm只能使用两个ID来运行,因此,如果您有多个ID,您将在索引编制之前合并这些ID,以如何“欺骗” plm。但是,我的数据仅包含以下列:国家,年份,马力和犯罪,因此我不知道这是怎么可能的。
本质上是我要问的,我做错了什么吗?我是否仍需要合并这两个ID,或者这是我的行重复项中的错误,如果是这样的话,是否有可能通过编码找到重复项? (我手动尝试查看我的面板数据以查找ID的重复项,即国家1的第一年房价的多个值。
如果我跑步
any(table(data$country,data$year)!=1)
我是对的。据我了解,这表明没有国家/年份组合的重复项。
请考虑以下适当的数据。
set.seed(42)
(d1 <- transform(expand.grid(id=1:2, time=1:2), X=rnorm(4), y=rnorm(4)))
# id time X y
# 1 1 1 1.3709584 0.40426832
# 2 2 1 -0.5646982 -0.10612452
# 3 1 2 0.3631284 1.51152200
# 4 2 2 0.6328626 -0.09465904
library(plm)
plm(y ~ X, index=c("id", "time"), d1)
# works
现在让我们复制最后一行以模拟数据中的缺陷,
(d1 <- rbind(d1, d1[nrow(d1), ]))
# id time X y
# 1 1 1 1.3709584 0.40426832
# 2 2 1 -0.5646982 -0.10612452
# 3 1 2 0.3631284 1.51152200
# 4 2 2 0.6328626 -0.09465904
# 41 2 2 0.6328626 -0.09465904 ## duplicated (X and y may be different though)
我们收到错误的地方:
plm(y ~ X, index=c("id", "time"), d1)
# Error in pdim.default(index[[1]], index[[2]]) :
# duplicate couples (id-time)
类似地,如果我们有id
,time
和某些cond
ition的数据,则会发生错误:
(d2 <- transform(expand.grid(id=1:2, time=1:2, cond=0:1), X=rnorm(4), y=rnorm(4)))
# id time cond X y
# 1 1 1 0 2.0184237 -1.3888607
# 2 2 1 0 -0.0627141 -0.2787888
# 3 1 2 0 1.3048697 -0.1333213
# 4 2 2 0 2.2866454 0.6359504
# 5 1 1 1 2.0184237 -1.3888607
# 6 2 1 1 -0.0627141 -0.2787888
# 7 1 2 1 1.3048697 -0.1333213
# 8 2 2 1 2.2866454 0.6359504
plm(y ~ X, index=c("id", "time"), d2)
# Error in pdim.default(index[[1]], index[[2]]) :
# duplicate couples (id-time)
为了克服这个问题,我们可以从技术上合并两个索引,无论从统计意义上讲是什么:
(d2 <- transform(d2, id2=apply(d2[c("id", "cond")], 1, paste, collapse=".")))
# id time cond X y id2
# 1 1 1 0 2.0184237 -1.3888607 1.0
# 2 2 1 0 -0.0627141 -0.2787888 2.0
# 3 1 2 0 1.3048697 -0.1333213 1.0
# 4 2 2 0 2.2866454 0.6359504 2.0
# 5 1 1 1 2.0184237 -1.3888607 1.1
# 6 2 1 1 -0.0627141 -0.2787888 2.1
# 7 1 2 1 1.3048697 -0.1333213 1.1
# 8 2 2 1 2.2866454 0.6359504 2.1
plm(y ~ X, index=c("id2", "time"), d2)
# works
[最后,此stopifnot
应该不是产生错误,其中c("id", "time")
对应于您在plm(..., index=c("id", "time"))
中定义的内容:
stopifnot(!any(duplicated(d1[c("id", "time")])))
# Error: !any(duplicated(d1[c("id", "time")])) is not TRUE