用R从多元线性回归中提取β
我有34个回归,对于不同的城市都是一样的,使用函数lm()完成,我想制作一个表,列出城市的名称和每个回归的第三个beta 我用这个循环创建了回归用R从多元线性回归中提取β,r,lm,R,Lm,我有34个回归,对于不同的城市都是一样的,使用函数lm()完成,我想制作一个表,列出城市的名称和每个回归的第三个beta 我用这个循环创建了回归 for ( i in 1:34 ){ #create a subset data data_sub <- subset(dif_dif, Comuna_subida == comunas[i]) assign(paste("lm", comunas[i] , sep = "_"), lm
for ( i in 1:34 ){
#create a subset data
data_sub <- subset(dif_dif, Comuna_subida == comunas[i])
assign(paste("lm", comunas[i] , sep = "_"), lm(log(viajes) ~
# Explanatory variables
cierre_colegio + Cuarentena.x +
dias_en_cuarentena +
# Control variables
n_comunas_cuarentena +
personas_en_cuarentena +
# Time Variables
dia_semana_Tuesday + dia_semana_Wednesday +
dia_semana_Thursday + dia_semana_Friday +
mes_4 + mes_5 + mes_6 + mes_7 + semana +
semana_corta + ano,
data = data_sub))
}
非常感谢 在for循环外部用两列初始化data.frame。我们称之为df。然后在for循环中,添加一行代码,为每个lm模型创建新行。它应该类似于:
df[i,]而不是将结果分配给具有不同名称的变量,您可以创建一个列表。
例如:
库(purrr)
图书馆(dplyr)
#计算回归
lmcalc市政A 2.023770
#>B市1.981352
#>市政C 2.015938
由(v0.3.0)于2020年9月22日创建
imap
函数允许使用表示列表中当前对象的.x
和.y
名称在列表上循环。请使用dput
添加数据,并显示相同的预期输出。阅读并了解如何给出建议。
Municipality beta_Cuarentena.x
X1 0.08
X2 0.15
X3 0.42
...etc