mdl.c2.thaw() report("mdl") from sherpa.fit import Fit f = Fit(d, mdl) report("f") report("f.data") report("f.model") dump("f.model.c2.val") mdl.c2 = 1 dump("f.model.c2.val") print("Starting statistic: {:.3f}".format(f.calc_stat())) sinfo1 = f.calc_stat_info() report("sinfo1") d.ignore(0, 5) sinfo2 = f.calc_stat_info() d.notice() dump("sinfo1.numpoints") dump("sinfo2.numpoints") res = f.fit() if res.succeeded: print("Fit succeeded") if not res.succeeded: print("**** ERRRR, the fit failed folks") report("res.format()") report("res")
""" residplot.prepare(d, mdl, f.stat) /home/djburke/miniconda2/envs/sherpa410-py35/lib/python3.5/site-packages/sherpa-4.10.0-py3.5-linux-x86_64.egg/sherpa/plot/__init__.py:1128: RuntimeWarning: divide by zero encountered in true_divide return (ylist[0] - ylist[1]) / staterr Traceback (most recent call last): File "example.py", line 125, in <module> residplot.prepare(d, mdl, f.stat) File "/home/djburke/miniconda2/envs/sherpa410-py35/lib/python3.5/site-packages/sherpa-4.10.0-py3.5-linux-x86_64.egg/sherpa/plot/__init__.py", line 1140, in prepare self.yerr = staterr / staterr TypeError: unsupported operand type(s) for /: 'list' and 'list' """ d.ignore(None, 1) statinfo = f.calc_stat_info() report("statinfo") dump("statinfo.rstat == fitres3.rstat") dump("f.estmethod.name") coverrs = f.est_errors() report("coverrs.format()") dump("f.estmethod.sigma") f.estmethod.sigma = 1.6 coverrs90 = f.est_errors() report("coverrs90.format()") from sherpa.estmethods import Confidence