我正在尝试引导面板数据,然后运行池化面板回归,最后收集每个估计的系数。这是我到目前为止的代码:
df <- data.frame(id=c(rep('1',57),rep('2',57),rep('3',57),rep('4',57),rep('5',57)),YEAR=c(1961:2017), data=Project2)
N=length(Project2) # Counts the number of observations
B=10 # Number of times to recompute the estimate
stor.r2=matrix(c(rep(0,B)), nrow=B, ncol=20) #Storage matrix for the coefficients
for(i in 1:B) {
newdata.df = df[sample(nrow(df), 10, replace=TRUE), ]
wols.pop = plm(pop ~ Il+In+Io+Oh+Mis+TDum
+lag(cmt.)+lag(pop.)+lag(lv.)+lag(cr.)
+lag(cmt..1)+lag(pop..1)+lag(lv..1)+lag(cr..1)+lag(cpiplus)+lag(vlyeg)+0, data=newdata.df, weights=wls.pop/1000, model="pooling")
stor.r2[i] = summary(wols.pop)$r.squared
}
我不断收到以下错误消息:
row.names<-.data.frame
中的错误(*tmp*
,值= orig_rownames [as.numeric(row.names(data))]):
不允许重复'row.names'
另外:警告信息:
设置'row.names'时的非唯一值:
我不确定如何修复我的代码。