Question

我是r的新手，正在考虑使用plm（）的固定效果进行回归分析。我选择双向方法来考虑时间和个体影响。但是，在运行下面的代码后，我仍然收到以下消息：

Error in pdata.frame(data, index) : 
  variable id does not exist (individual index)

代码在这里：

pdata<-DATABASE[,c(2:4,13:21)]
pdata$id<-group_indices(pdata,ISO3.p,Productcode)
coutnin<-dcast.data.table(pdata,ISO3.p+Productcode~.,value.var = "id")
setcolorder(pdata,neworder=c("id","Year"))
pdata<-pdata.frame(pdata,index=c("id","Year"))
reg<-plm(pdata,diff(TV,1)~diff(RERcp,1)+diff(GDPR.p,1)-diff(GDPR.r,1),effect="twoways",model="within",index = c("id","Year"))

请注意，pdata结构表明id变量中存在数字形式的多个级别，我最初尝试使用字符串类型的变量，但始终收到相同的结果：

Classes ‘data.table’ and 'data.frame':  1211800 obs. of  13 variables:
 $ id         : int  4835 6050 13158 15247 17164 18401 19564 23553 24895 27541 ...
 $ Year       : int  1996 1996 1996 1996 1996 1996 1996 1996 1996 1996 ...
 $ Productcode: chr  "101" "101" "101" "101" ...
 $ ISO3.p     : Factor w/ 171 levels "ABW","AFG","AGO",..: 8 9 20 22 27 28 29 34 37 40 ...
 $ e          : num  0.245 -0.238 1.624 0.693 0.31 ...
 $ RERcp      : num  -0.14073 -0.16277 1.01262 0.03908 -0.00243 ...
 $ RERpp      : num  -0.1712 NA NA NA -0.0952 ...
 $ RER_GVC    : num  -3.44 NaN NA NA NaN ...
 $ GDPR.p     : num  27.5 26.6 23.5 20.3 27.8 ...
 $ GDPR.r     : num  30.4 30.4 30.4 30.4 30.4 ...
 $ GVCPos     : num  0.141 0.141 0.141 0.141 0.141 ...
 $ GVCPar     : num  0.436 0.436 0.436 0.436 0.436 ...
 $ TV         : num  17.1 17.1 17.1 17.1 17.1 ...
 - attr(*, ".internal.selfref")=<externalptr>

当我将data.table转换为pdata.frame时，我没有收到任何警告，它仅在运行plm函数之后发生。从运行View（table（index（index（pdata），useNA =“ ifany”））开始，它显示的值不大于1，因此我假设我的数据中没有重复的obs。

非常感谢您的帮助！

Answer 1

尝试将数据参数放在plm语句的第二位。如果已经将pdata转换为pdata.frame，请在index语句中省略plm参数，即尝试以下操作：

reg <- plm(diff(TV,1) ~ diff(RERcp,1)+diff(GDPR.p,1)-diff(GDPR.r,1), data = pdata, effect = "twoways", model = "within")

PLM无法识别我的ID变量名称

1 个答案: