Step 1: extract the data into a data frame. new_DF<-foo_df

Step 2: Fit initial model initial_model <- lm(WHAT_YOU'RE_TRYING_TO_GET ~ columns + by + name + from + dataframe, data = new_DF)

Step 3: Review model summary, remove non-significant predictors, and refit as needed summary(initial_model) #Replace 'significant_predictors' accordingly reduced_model <- lm(WHAT_YOU'RE_TRYING_TO_GET ~ less + columns + by + than + before, data = new_DF)