def _do_plot(x, y, dist=None, line=False, ax=None, fmt='bo', **kwargs):
    fig, ax = utils.create_mpl_ax(ax)
    ax.set_xmargin(0.02)
    ax.plot(x, y, fmt, **kwargs)
    if line:
        if line not in ['r','q','45','s']:
            msg = "%s option for line not understood" % line
            raise ValueError(msg)

        qqline(ax, line, x=x, y=y, dist=dist)

    return fig, ax
Boiler plate plotting function for the `ppplot`, `qqplot`, and
dataset/ETHPy150Open statsmodels/statsmodels/statsmodels/graphics/gofplots.py/_do_plot