report("coverrs") dump("f.estmethod.sigma") f.estmethod.sigma = 1.6 coverrs90 = f.est_errors() report("coverrs90.format()") dump("coverrs90.percent") report("coverrs.extra_output") print([p.split('.')[1] for p in coverrs.parnames]) from sherpa.estmethods import Confidence f.estmethod = Confidence() print("*** about to start confidence run") conferrs = f.est_errors() print("*** finished confidence run") report("conferrs.format()") dump("f.estmethod.name") dump("f.estmethod.parallel") print("*** about to start confidence run") c1errs = f.est_errors(parlist=(mdl.c1, )) print("*** finished confidence run") report("c1errs") from sherpa.plot import IntervalProjection
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 f.estmethod = Confidence() print("*** start confidence errors") conferrs = f.est_errors() print("*** end confidence errors") report("conferrs.format()") print("*** start confidence errors") offseterrs = f.est_errors(parlist=(mdl.pars[1], )) print("*** end confidence errors") report("offseterrs") fmt = "{:13s} covar=±{:4.2f} conf={:+5.2f} {:+5.2f}" for i in range(len(conferrs.parnames)): print(fmt.format(conferrs.parnames[i], coverrs.parmaxes[i], conferrs.parmins[i], conferrs.parmaxes[i]))