title: statsmodels: printing summary of more than one regression models together
body: <p>In the Python library <code>Statsmodels</code>, you can print out the regression results with <code>print(results.summary())</code>, how can I print out the summary of more than one regressions in one table, for better comparison? </p>
  
  <p>A linear regression, code taken from <code>statsmodels</code> documentation:</p>
  
  <pre><code>nsample = 100
  x = np.linspace(0, 10, 100)
  X = np.column_stack((x, x**2))
  beta = np.array([0.1, 10])
  e = np.random.normal(size=nsample)
  y = np.dot(X, beta) + e
  
  model = sm.OLS(y, X)
  results_noconstant = model.fit()
  </code></pre>
  
  <p>Then I add a constant to the model and run the regression again:</p>
  
  <pre><code>beta = np.array([1, 0.1, 10])
  X = sm.add_constant(X)
  y = np.dot(X, beta) + e 
  
  model = sm.OLS(y, X)
  results_withconstant = model.fit()
  </code></pre>
  
  <p>I'd like to see the summaries of <code>results_noconstant</code> and <code>results_withconstant</code> printed out in one table. This should be a very useful function, but I didn't find any instruction about this in the <code>statsmodels</code> documentation.</p>
  
  <p><strong>EDIT</strong>: The regression table I had in mind would be something like <a href="http://repec.org/bocode/e/estout/esttab012c.png" rel="nofollow">this</a>, I wonder whether there is ready-made functionality to do this.</p>