add_M2: Attach the margin-based fit indices to a fitted model
Description
Computes M2 for the model and for the independence baseline, and
returns the fitted object with a TestFitIndicesM2 component added. The
print method then shows the response-pattern indices and the margin-based ones
side by side.
This is a separate step rather than part of the fit because it is expensive:
the cost is the Cholesky factorisation of a dense matrix whose size grows with
the square of the item count (see M2).