  We describe how to consistently incorporate solar model uncertainties, along
with experimental errors and correlations, when analyzing solar neutrino data
to derive confidence limits on parameter space for proposed solutions of the
solar neutrino problem. Our work resolves ambiguities and inconsistencies in
the previous literature. As an application of our methods we calculate the
masses and mixing angles allowed by the current data for the proposed MSW
solution using both Bayesian and frequentist methods, allowing purely for solar
model flux variations, to compare with previous work. We consider the effects
of including metal diffusion in the solar models and also discuss implications
for future experiments.
