Note that ``fit()`` can be used to calculate mean and marginal effects at arbitrary points, while r_squared is calculated at ``exog``. So there may not be a previously calculated result.

That said, the regression interface is not very well thought out compared to that of density estimation. Needs more testing as well, especially for mixed data and of the marginal effects.