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]))
