plm封装的权值

Question

我的数据框架如下所示：

unique.groups<- letters[1:5]
unique_timez<- 1:20
groups<- rep(unique.groups, each=20)
my.times<-rep(unique_timez, 5)

play.data<- data.frame(groups, my.times, y= rnorm(100), x=rnorm(100), POP= 1:100)

我想进行以下加权回归：

plm(y~x + factor(my.times) , 
data=play.data, 
index=c('groups','my.times'), model='within', weights= POP)

但我不相信plm包允许重量。我在下面的模型中寻找系数的答案：

fit.regular<- lm(y~x + factor(my.times) + factor(my.groups),
weights= POP, data= play.data)
desired.answer<- coefficients(fit.regular)

然而，我正在寻找一个与plm软件包的答案，因为它是更快地得到系数的内部估计与plm与更大的数据集和许多组。

Answer 1

编辑：这个问题已经不存在了，因为plm现在使用了一个权重函数(见上面的@Helix123 123注释)。

尽管我不知道plm包有什么解决方案，但lfe包中的felm函数在固定效果上下文中正确地处理权重(从示例代码的语法来看，这似乎是您所需要的)。特别是在有许多观察和小组在场的情况下，它的重点是速度。

lfe包只关注固定的效果，因此如果您需要随机效果，lme4包可能会根据您的需要进行be more suited。

Answer 2

我要找的正是这个信息。我找到了一个软件包的作者http://r.789695.n4.nabble.com/Longitudinal-Weights-in-PLM-package-td3298823.html的答案，这似乎表明在plm包中没有直接使用权重的方法。