Vpop calibration
Description
A set of Python tools to allow for virtual population calibration, using a non-linear mixed effects (NLME) model approach, combined with surrogate models in order to speed up the simulation of QSP models.
The approach was mainly inspired from 1.
Here is an overview of the proposed workflow:
Currently available features
- Surrogate modeling using gaussian processes, implemented using GPyTorch
- Synthetic data generation using ODE models. The current implementation uses scipy.integrate.solve_ivp, parallelized with multiprocessing
- Non-linear mixed effect models, see the dedicated doc:
- Log-distributed parameters
- Additive or multiplicative error model
- Covariates handling
- Known individual patient descriptors (i.e. covariates with no effect on other descriptors outside of the structural model)
- SAEM: see the dedicated doc
- Optimization of random and fixed effects using repeated longitudinal data
- Individual parameter estimation
- Conditional distribution sampling: the posterior conditional distribution may directly be sampled using Metropolis-Hastings algorithm
- Empirical Bayesian Estimators (EBEs) or Maximum A Posteriori (MAP) estimates: these values represent the mode of the conditional distribution for each parameter and each patient. The current implementation is sub-optimal for
torchsurrogate models as it does not leverage the gradient predictions.
- Diagnostics:
- Post-inference diagnostics are available, and different type sof weighted residuals may be plotted (IWRES, PWRES, NPDE)
Getting started
- Tutorial: this notebook demonstrates step-by-step how to create and train a surrogate model, using a reference ODE model and a GP surrogate model. It then showcases how to optimize the surrogate model on synthetic data using SAEM
In-depth examples
- Surrogate model training
- Running SAEM without surrogate model
- Running SAEM on a reference ODE model. Note: the current implementation is notably under-optimized for running SAEM directly on an ODE structural model. This is implemented for testing purposes mostly
- Benchmarking examples
Support
For any issue or comments, please reach out to paul.lemarre@novainsilico.ai, or feel free to open an issue in the repo directly.
Authors
- Paul Lemarre
- Eléonore Dravet
- Hugo Alves
Acknowledgements
- Adeline Leclercq-Samson
- Eliott Tixier
- Louis Philippe
QSPC26 poster
This work was presented at QSPC2026 in Leiden, and all the corresponding material is introduced in this document. The benchmark notebooks for standard data sets are available directly for orange trees and theophylline.
The poster itself is available in this document. The full list of references is available in this document.
Roadmap
Here are a few planned features, in an unordered and non-exhaustive list:
- NLME:
- Support additional error models (additive-multiplicative, power, etc...)
- Support additional covariate models (categorical covariates)
- Compute likelihood via importance sampling following population parameters optimization
- Structural models:
- Integrate with SBML models (e.g. Roadrunner)
- Surrogate models:
- Support additional surrogate models in PyTorch
- Optimizer:
- Add preconditioned Stochastic-Gradient-Descent (SGD) method for surrogate model optimization
References
-
Grenier et al. 2018: Grenier, E., Helbert, C., Louvet, V. et al. Population parametrization of costly black box models using iterations between SAEM algorithm and kriging. Comp. Appl. Math. 37, 161–173 (2018). https://doi.org/10.1007/s40314-016-0337-5 ↩
3 maintainers
Metadata
3 maintainers
Release files for vpop-calibration 4.7.1
For a detailed explanation of source distributions (sdists) and built distributions (wheels), please see the package formats documentation.
Source distribution (sdist)
| File | Size | Uploaded | |
|---|---|---|---|
| vpop_calibration-4.7.1.tar.gz | 278.8 kB | Details |
Built distribution (wheel)
| File | Interpreter | ABI | Platform | Reset |
|---|---|---|---|---|
| vpop_calibration-4.7.1-py3-none-any.whl | Python 3 | none | any | Details |
Total release size: 602.2 kB
Release files / vpop_calibration-4.7.1.tar.gz
| Download URL | vpop_calibration-4.7.1.tar.gz |
|---|---|
| Size | 278.8 kB |
| Tags | Source |
|
SHA-256 checksum How to use checksums |
783c21cb81117f3c90640d7b0ab382dbc6fc6331b299a9a5e2eb722fb888118b
|
|
BLAKE2b-256 checksum How to use checksums |
8d0b6ad2a509b291d5490fda64d517ce2acdaa8a6c45a3bd23bdb6c72ccdca36
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
poetry/2.4.1 CPython/3.13.13 Linux/6.18.33
|
Release files / vpop_calibration-4.7.1-py3-none-any.whl
| Download URL | vpop_calibration-4.7.1-py3-none-any.whl |
|---|---|
| Size | 323.4 kB |
| Tags | Python 3 |
|
SHA-256 checksum How to use checksums |
885b953583865c65e3897b50d83871fea7ee5d5db17bca59f8ef7c42b6089122
|
|
BLAKE2b-256 checksum How to use checksums |
39c9c1c82580bf59782bc87bdd596a0d36ef28336155aa1d07ab1c1fbc51d930
|
| Upload date | |
|
Uploaded using Trusted Publishing? What is trusted publishing? |
No |
| Uploaded via |
poetry/2.4.1 CPython/3.13.13 Linux/6.18.33
|
