小编Zsl*_*ce 的帖子

采用plm包装的重量

我的数据框如下所示:

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)
Run Code Online (Sandbox Code Playgroud)

我想运行以下加权回归:

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

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

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

但是,我正在寻找plm包的答案,因为使用更大的数据集和许多组获得带有plm的估计器的系数要快得多.

r linear-regression panel-data plm

8
推荐指数
1
解决办法
1427
查看次数

将for循环中的结果保存为r中的向量

我有一个相当简单的问题.

我有以下代码:

for (i in 1:10) {print (i) } 
Run Code Online (Sandbox Code Playgroud)

如何将结果保存为矢量?

如果我试试

my_vector <- for (i in 1:10) {print (i) } 
Run Code Online (Sandbox Code Playgroud)

然后调用my_vector,我得到Null.

for-loop r vector

5
推荐指数
1
解决办法
3万
查看次数

标签 统计

r ×2

for-loop ×1

linear-regression ×1

panel-data ×1

plm ×1

vector ×1