def chi_squared(self, p=None):
chi2s = self.chi_squareds(p)<EOL>if chi2s == None: return None<EOL>return sum(self.chi_squareds(p))<EOL>
Returns the total chi squared (summed over all massaged data sets).

p=None means use the fit results.