eval（expr、envir、enclos）のエラー：オブジェクト 'mpg'が見つかりません

Question

data.tableパッケージを使用して2変数のグループワウンド回帰係数を計算しようとしています。

ここでは、ダミーデータでサンプルコードを投稿しました。

#Model Dependent varaible 
reg_dep_vars<-"mpg" 
#Model independent variable 
reg_ind_vars<-c("cyl","drat") 

reg_data<-as.data.table(mtcars) 

#creating a formula with depedent and independent variables which going to be used in the model. 
reg_formula<-as.formula(paste(paste("reg_data$",reg_dep_vars,sep=""),"~",paste(paste("reg_data$",reg_ind_vars,sep=""),collapse="+"))) 

OUT<-reg_data[,.(intercept=coef(lm(reg_formula))[1],cyl=coef(lm(reg_formula))[2],drat=coef(lm(reg_formula))[3],P=glance(lm(reg_formula))$p.value,F=summary(lm(reg_formula))$fstatistic[1]),by=.(am,gear)] 

上記のコードでは、cyl変数とdrat変数の見積もりを調べようとしています。グループ別の変数はamとgearです。

上記のコードを使用すると、次のエラーが発生します。 "eval（expr、envir、enclos）のエラー：オブジェクト 'mpg'が見つかりません"

誰でもこのことを教えてもらえますか？

Answer 1

使用

library(tidyverse) 

mtcars_model <- function(df) { 
    lm.fit(y = df[[reg_dep_vars]], x = as.matrix(df[reg_ind_vars])) 
}  

test <- mtcars %>% 
    group_by(am, gear) %>% 
    nest() %>% 
    mutate(model = map(data, mtcars_model)) 

はhttp://r4ds.had.co.nz/many-models.htmlで詳細をご覧ください。

PS：申し訳ありませんが、私は数式に立つことはできません。