How do I create a “macro” for regressors in R?

by

Here are some alternatives. No packages are used in the first 3.

1) reformulate

fo <- reformulate(regressors, response = "income")
lm(fo, Duncan)

or you may wish to write the last line as this so that the formula that is shown in the output looks nicer:

do.call("lm", list(fo, quote(Duncan)))

in which case the Call: line of the output appears as expected, namely:

Call:
lm(formula = income ~ education + prestige, data = Duncan)

2) lm(dataframe)

lm( Duncan[c("income", regressors)] )

The Call: line of the output look like this:

Call:
lm(formula = Duncan[c("income", regressors)])

but we can make it look exactly as in the do.call solution in (1) with this code:

fo <- formula(model.frame(income ~., Duncan[c("income", regressors)]))
do.call("lm", list(fo, quote(Duncan)))

3) dot

An alternative similar to that suggested by @jenesaisquoi in the comments is:

lm(income ~., Duncan[c("income", regressors)])

The approach discussed in (2) to the Call: output also works here.

4) fn$ Prefacing a function with fn$ enables string interpolation in its arguments. This solution is nearly identical to the desired syntax shown in the question using $ in place of @ to perform substitution and the flexible substitution could readily extend to more complex scenarios. The quote(Duncan) in the code could be written as just Duncan and it will still run but the Call: shown in the lm output will look better if you use quote(Duncan).

library(gsubfn)

rhs <- paste(regressors, collapse = "+")
fn$lm("income ~ $rhs", quote(Duncan))

The Call: line looks almost identical to the do.call solutions above — only spacing and quotes differ:

Call:
lm(formula = "income ~ education+prestige", data = Duncan)

If you wanted it absolutely the same then:

fo <- fn$formula("income ~ $rhs")
do.call("lm", list(fo, quote(Duncan)))

More Related Contents:

  1. How to create a loop that will make regression models in R? [closed]
  2. How to debug “contrasts can be applied only to factors with 2 or more levels” error?
  3. Predict() – Maybe I’m not understanding it
  4. How to succinctly write a formula with many variables from a data frame?
  5. Fitting a linear model with multiple LHS
  6. Getting Warning: ” ‘newdata’ had 1 row but variables found have 32 rows” on predict.lm
  7. `lm` summary not display all factor levels
  8. How does predict.lm() compute confidence interval and prediction interval?
  9. Extract regression coefficient values
  10. Pass a vector of variables into lm() formula
  11. Examples of the perils of globals in R and Stata
  12. How to get name of variable in R (substitute)?
  13. Why is using update on a lm inside a grouped data.table losing its model data?
  14. Plot polynomial regression curve in R
  15. How do I extract just the number from a named number (without the name)?
  16. How to identify/delete non-UTF-8 characters in R
  17. Mutate multiple columns in a dataframe
  18. predict.lm() in a loop. warning: prediction from a rank-deficient fit may be misleading
  19. Read Stata 13 file in R
  20. How `poly()` generates orthogonal polynomials? How to understand the “coefs” returned?
  21. lm function in R does not give coefficients for all factor levels in categorical data
  22. Is there a faster lm function
  23. Reading multiple files into multiple data frames
  24. Showing string in formula and not as variable in lm fit
  25. predict.lm() with an unknown factor level in test data
  26. Set one or more of coefficients to a specific integer
  27. R error which says “Models were not all fitted to the same size of dataset”
  28. Messy plot when plotting predictions of a polynomial regression using lm() in R
  29. Linear model function lm() error: NA/NaN/Inf in foreign function call (arg 1)
  30. Different Robust Standard Errors of Logit Regression in Stata and R

Leave a Comment