  Context. Monte Carlo methods have enabled nuclear reaction rates from
uncertain inputs to be presented in a statistically meaningful manner. However,
these uncertainties are currently computed assuming no correlations between the
physical quantities that enter those calculations. This is not always an
appropriate assumption. Astrophysically important reactions are often dominated
by resonances, whose properties are usually normalized to a well-known
reference resonance. This insight provides a basis from which to develop a
flexible framework for including correlations in Monte Carlo reaction rate
calculations. Aims. The aim of this work is to develop and test a method for
including correlations in Monte Carlo reaction rate calculations when the input
has been normalized to a common reference. Methods. A mathematical framework is
developed for including correlations between input parameters in Monte Carlo
reaction rate calculations. The magnitude of those correlations is calculated
from the uncertainties typically reported in experimental papers, where full
correlation information is not available. The method is applied to four
illustrative examples: a fictional 3-resonance reaction,
$^{27}$Al(p,$\gamma$)$^{28}$Si, $^{23}$Na(p,$\alpha$)$^{20}$Ne, and
$^{23}$Na($\alpha$,p)$^{26}$Al. Results. Reaction rates at low temperatures
that are dominated by a few isolated resonances are found to minimally impacted
by correlation effects. However, reaction rates determined from many
overlapping resonances can be significantly affected. Uncertainties in the
$^{23}$Na($\alpha$,p)$^{26}$Al reaction, for example, increase by up to a
factor of 5. This highlights the need to take correlation effects into account
in reaction rate calculations, and provide insight into which cases are
expected to be most affected by them. The impact of correlation effects on
nucleosynthesis is also investigated.
