只有两个ID的plm中的重复对(id-time)错误

问题描述 投票:1回答:1

我正在尝试使用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) 

我是对的。据我了解,这表明没有国家/年份组合的重复项。

r duplicates dataset plm
1个回答
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)

类似地,如果我们有idtime和某些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
© www.soinside.com 2019 - 2024. All rights reserved.