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])
# HACK due to https://github.com/sherpa/sherpa/issues/342
f.thaw_indices = tuple(i for i, p in enumerate(f.model.pars) if not p.frozen)
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")
