lumicks.pylake.parameter_trace

parameter_trace(model, params, inverted_parameter, independent, dependent, **kwargs)

Invert a model with respect to one parameter. This function fits a unique parameter for every data point in this data-set while keeping all other parameters fixed. This can be used to for example invert the model with respect to the contour length or some other parameter.

Parameters:

model : Model

Fitting model.

params : Params

Model parameters.

inverted_parameter : str

Parameter to invert.

independent : array_like

Vector of values for the independent variable

dependent : array_like

Vector of values for the dependent variable

**kwargs :

Forwarded to scipy.optimize.least_squares

Examples

# Define the model to be fitted
model = pylake.inverted_odijk("model") + pylake.force_offset("model")

# Fit the overall model first
model.add_data("dataset1", f=force_data, d=distance_data)
current_fit = pylake.FdFit(model)
data_handle = current_fit.add_data("my data", force, distance)
current_fit.fit()

# Calculate a per data point contour length
lcs = parameter_trace(model, current_fit[data_handle], "model/Lc", distance, force)