You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Technically this would imply hiding such parameters from the minimizer and doing the optimization within the NLL function evaluation. There are probably a number of ways to try this out relatively easily (even without changes to pyhf itself). One could replace model.logpdf by manually evaluated Poisson + constraint terms, where the model prediction model.main_model.expected_data going into the Poisson term is performed with the optimal staterror parameter values. These parameters could be set to constant, but would just need to automatically get updated at each function evaluation. This implies no need for changes in the minimizer themselves.
I gave this a try, here's an implementation: https://gist.github.com/alexander-held/b020445d089c3ef79abe93d63543c5e7. Seems to be working, but I have not tested it extensively. There are a few to-do items and some cleanup needed. In addition, this will not be working correctly with (partially) pruned staterror configurations, so that needs some thought.
Summary
It would be great to add analytic staterror parameter optimization, as available in Combine: https://cms-analysis.github.io/HiggsAnalysis-CombinedLimit/part2/bin-wise-stats/#analytic-minimisation.
Technically this would imply hiding such parameters from the minimizer and doing the optimization within the NLL function evaluation. There are probably a number of ways to try this out relatively easily (even without changes to pyhf itself). One could replace
model.logpdf
by manually evaluated Poisson + constraint terms, where the model predictionmodel.main_model.expected_data
going into the Poisson term is performed with the optimal staterror parameter values. These parameters could be set to constant, but would just need to automatically get updated at each function evaluation. This implies no need for changes in the minimizer themselves.cc @nsmith-
Additional Information
n/a
Code of Conduct
The text was updated successfully, but these errors were encountered: