Ireland_quasi_symmetry function

Fit for quasi-symmetry model. Obtained by subtraction, so no model-based probabilities.