 <p>This function could probably be improved upon. Tip of the hat to @thelatemail.</p> <pre><code>lm.3Trans = function(y, x1, x2, transformations = c(log,sqrt,pwer)){ pwer = function(x, p = 2) poly(x,p) res = lapply(transformations, function(f) lm(y ~ f(x1) + x2)) res } </code></pre> <p>This will output your models into a list to which you could then do something like:</p> <pre><code>lapply(lm.3Trans(y,x,x1), summary) </code></pre> <p>To get more detailed summaries</p> <p>It might be worth expanding this function to take a formula, and perhaps another arguement to identify which predictor should be transformed. Please let me know if that is useful. </p>