title: Use mapply to fit list of lm models to list of data frames [R]
body: <p>Is this possible? I can get <code>mapply</code> to work with the <code>help</code> examples, but I can't get a trivial example with <code>lm</code> to work. Here's my attempt which returns a matrix, instead of a list of <code>lm</code> objects.</p>
  
  <pre><code>temp.df &lt;- list(
                  data.frame(a = rep(1:10, each = 10), b = 1:100, c = rnorm(100), d = rnorm(100, 2))
                             )
  temp.df[[2]] &lt;- subset(temp.df[[1]], a &gt; 2)
  temp.mod &lt;- list(a ~ b,
                a ~ b + c,
               a ~ b + c + d)
  temp.lm &lt;- mapply(lm, formula = temp.mod, data = temp.df[c(1,1,2)])
  temp.sum &lt;- lapply(temp.lm, summary)
  </code></pre>
  
  <p>Should I just stick with <code>lapply</code> and specifying <code>data =</code> each time? Thanks!</p>