question-mark
Stuck on an issue?

Lightrun Answers was designed to reduce the constant googling that comes with debugging 3rd party libraries. It collects links to all the places you might be looking at while hunting down a tough bug.

And, if you’re still stuck at the end, we’re happy to hop on a call to see how we can help out.

Pypesto Optimization ignores bounds

See original GitHub issue

I recently found that the parameter estimation of the ICG PEtab model resulted in parameters that violate the bounds given in the PEtab parameter table.

The PEtab problem I used is: icg_sd.omex.zip

And I used the following code that produced the error:

import pypesto import pypesto.petab import pypesto.optimize as optimize import pypesto.store as store

pypesto_importer = pypesto.petab.PetabImporter.from_yaml(yaml_file) pypesto_problem = pypesto_importer.create_problem()

optimizer = optimize.ScipyDifferentialEvolutionOptimizer(options = {'strategy':'best1bin', 'maxiter':8000 ,'mutation':(0.2, 1), 'recombination':0.7, 'popsize':70, 'init':'latinhypercube', 'polish':True, 'workers' : 8})

n_starts = 10 result = optimize.minimize(pypesto_problem, optimizer=optimizer, n_starts=n_starts)

store.save_to_hdf5.write_result(result, os.path.join(Path , M + '_Pypesto_Parameters.hdf5'), overwrite= True)

Issue Analytics

  • State:closed
  • Created a year ago
  • Comments:6 (4 by maintainers)

github_iconTop GitHub Comments

1reaction
shoepflcommented, Jun 24, 2022

Okay thanks for the help - maybe that needs more clairification in the PEtab documentation as the description for the bounds state there: “Must be provided in linear space, independent of parameterScale.”

from this at least it is confusing when I have to give the bounds in linear scale explicitly but indeed the log of them are taken as boundaries. But thats probably a topic for the PEtab implementation so I will shift the discussion to there.

Thanks for your help!

0reactions
yannikschaeltecommented, Aug 1, 2022

Closing here as this concerns the petab documentation

Read more comments on GitHub >

github_iconTop Results From Across the Web

Scipy.optimize.minimize ignores the bounds I specify while ...
I am trying to minimize a function of 76 parameters. I want every parameter value to be between 0 and 5. The default...
Read more >
API reference — pyPESTO 0.2.14 documentation
lower_bound – array of potential lower bounds for the parameters ... result – A pyPESTO result that contains an optimization result.
Read more >
pyPESTO/optimizer.py at main · ICB-DCM/pyPESTO - GitHub
Contribute to ICB-DCM/pyPESTO development by creating an account on GitHub. ... bounds = scipy.optimize. ... minimize will ignore hessp otherwise.
Read more >
pypesto Changelog - pyup.io
Fix waterfall plot limits for non-offsetted log-plots (891) ... Optimize `pypesto/objective/amici_util.py::par_index_slices` (845) * Remove Solver.
Read more >
Fides: Reliable trust-region optimization for parameter ... - PLOS
Optimization -based calibration of ODE models on is often challenging, ... an interface to fides in the parameter estimation toolbox pyPESTO, ...
Read more >

github_iconTop Related Medium Post

No results found

github_iconTop Related StackOverflow Question

No results found

github_iconTroubleshoot Live Code

Lightrun enables developers to add logs, metrics and snapshots to live code - no restarts or redeploys required.
Start Free

github_iconTop Related Reddit Thread

No results found

github_iconTop Related Hackernoon Post

No results found

github_iconTop Related Tweet

No results found

github_iconTop Related Dev.to Post

No results found

github_iconTop Related Hashnode Post

No results found