`__
for more about the data and model used to generate this figure).
.. image:: ./bootstrap_fig.png
**Figure 1:** Scatterplot matrix of parameter values from 500 rounds of
bootstrap resampling fits to an SDSS *r*-band image of the galaxy IC
3478 (single-Sérsic model, no PSF convolution). Note the clear
correlations between the Sérsic model parameters (*n*, *r_e*, *I_e*).
See `here <./pyimfit_bootstrap_BtoT.html>`__ for an example of using
bootstrap output to estimate uncertainties on derived quantities, such
as bulge/total values.
.. raw:: html
Using MCMC
----------
Estimates of parameter uncertainties and correlations can also be
obtained via Markov-Chain Monte Carlo (MCMC) approaches. Although the
MCMC option of **Imfit** (``imfit-mcmc``) is not part of PyImfit, you
*can* use instances of the Imfit class with Python-based MCMC codes,
such as `emcee `__; see
`here <./pyimfit_emcee.html>`__ for an example.