amici.petab
PEtab import related code.
- amici.petab.import_petab_problem(petab_problem, model_output_dir=None, model_name=None, compile_=None, non_estimated_parameters_as_constants=True, **kwargs)[source]
Create an AMICI model for a PEtab problem.
- Parameters:
petab_problem (
petab.v1.problem.Problem
) – A petab problem containing all relevant information on the model.model_output_dir (
str
|pathlib.Path
|None
) – Directory to write the model code to. It will be created if it doesn’t exist. Defaults to current directory.model_name (
str
) – Name of the generated model module. Defaults to the ID of the model or the model file name without the extension.compile – If
True
, the model will be compiled. IfFalse
, the model will not be compiled. IfNone
, the model will be compiled if it cannot be imported.non_estimated_parameters_as_constants – Whether parameters marked as non-estimated in PEtab should be considered constant in AMICI. Setting this to
True
will reduce model size and simulation times. If sensitivities with respect to those parameters are required, this should be set toFalse
.kwargs – Additional keyword arguments to be passed to
amici.sbml_import.SbmlImporter.sbml2amici()
oramici.pysb_import.pysb2amici()
, depending on the model type.
- Return type:
- Returns:
The imported model.
- amici.petab.rdatas_to_measurement_df(rdatas, model, measurement_df)[source]
Create a measurement dataframe in the PEtab format from the passed
rdatas
and own information.- Parameters:
rdatas (
collections.abc.Sequence
[amici.amici.ReturnData
]) – A sequence of rdatas with the ordering ofpetab.get_simulation_conditions()
.model (
typing.Union
[amici.amici.Model
,amici.amici.ModelPtr
]) – AMICI model used to generaterdatas
.measurement_df (
pandas.core.frame.DataFrame
) – PEtab measurement table used to generaterdatas
.
- Return type:
- Returns:
A dataframe built from the rdatas in the format of
measurement_df
.
- amici.petab.rdatas_to_simulation_df(rdatas, model, measurement_df)[source]
Create a PEtab simulation dataframe from
amici.amici.ReturnData
s.See
rdatas_to_measurement_df()
for details, only that model outputs will appear in columnsimulation
instead ofmeasurement
.- Return type:
- amici.petab.simulate_petab(petab_problem, amici_model, solver=None, problem_parameters=None, simulation_conditions=None, edatas=None, parameter_mapping=None, scaled_parameters=False, log_level=30, num_threads=1, failfast=True, scaled_gradients=False)[source]
Simulate PEtab model.
Note
Regardless of scaled_parameters, unscaled sensitivities are returned, unless scaled_gradients=True.
- Parameters:
petab_problem (
petab.v1.problem.Problem
) – PEtab problem to work on.amici_model (
typing.Union
[amici.amici.Model
,amici.amici.ModelPtr
]) – AMICI Model assumed to be compatible withpetab_problem
.solver (
amici.amici.Solver
|None
) – An AMICI solver. Will use default options if None.problem_parameters (
dict
[str
,float
] |None
) – Run simulation with these parameters. IfNone
, PEtabnominalValues
will be used. To be provided as dict, mapping PEtab problem parameters to SBML IDs.simulation_conditions (
pandas.core.frame.DataFrame
|dict
) – Result ofpetab.get_simulation_conditions()
. Can be provided to save time if this has be obtained before. Not required ifedatas
andparameter_mapping
are provided.edatas (
list
[typing.Union
[amici.amici.ExpData
,amici.amici.ExpDataPtr
]]) – Experimental data. Parameters are inserted in-place for simulation.parameter_mapping (
amici.petab.parameter_mapping.ParameterMapping
) – Optional precomputed PEtab parameter mapping for efficiency, as generated bycreate_parameter_mapping()
withscaled_parameters=True
.scaled_parameters (
bool
|None
) – IfTrue
,problem_parameters
are assumed to be on the scale provided in the PEtab parameter table and will be unscaled. IfFalse
, they are assumed to be in linear scale. If parameter_mapping is provided, this must match the value of scaled_parameters used to generate the mapping.log_level (
int
) – Log level, seeamici.logging
module.num_threads (
int
) – Number of threads to use for simulating multiple conditions (only used if compiled with OpenMP).failfast (
bool
) – Returns as soon as an integration failure is encountered, skipping any remaining simulations.scaled_gradients (
bool
) – Whether to compute gradients on parameter scale (True
) or not (False
).
- Return type:
- Returns:
Dictionary of
cost function value (
LLH
),list of
amici.amici.ReturnData
(RDATAS
),list of
amici.amici.ExpData
(EDATAS
),
corresponding to the different simulation conditions. For ordering of simulation conditions, see
petab.Problem.get_simulation_conditions_from_measurement_df()
.