pyphotomol.utils package#

Submodules#

pyphotomol.utils.data_import module#

Data import utilities for mass photometry analysis.

This module provides functions to import data from various file formats commonly used in mass photometry experiments, including HDF5 files from Refeyn instruments and CSV files with custom data formats.

Functions#

import_file_h5 : Import data from Refeyn HDF5 files import_csv : Import data from CSV files import_movie_h5 : Import movie data from HDF5 files

The import functions automatically handle different data structures and perform basic preprocessing like NaN removal and calibration conversions.

pyphotomol.utils.data_import.import_file_h5(filename)[source]#

Import mass photometry data from HDF5 files generated by Refeyn instruments.

This function reads contrast and mass data from standard Refeyn HDF5 file formats. It automatically handles different data structures and performs calibration conversions when necessary. NaN values are filtered out automatically.

Parameters:

filename (str) – Path to the HDF5 file to import

Returns:

  • contrasts (np.ndarray) – Array of contrast values with NaN values removed

  • masses_kDa (np.ndarray or None) – Array of mass values in kDa units, or None if not available. If masses are not directly available but calibration data exists, they will be computed from contrasts using the calibration parameters.

Notes

The function searches for data in the following order: 1. Direct ‘contrasts’ and ‘masses_kDa’ datasets 2. ‘calibrated_values’ with calibration parameters 3. ‘per_movie_events’ for movie-based data

Raises:

  • FileNotFoundError – If the specified file does not exist

  • KeyError – If required datasets are not found in the HDF5 file

  • ValueError – If the file format is not recognized or data is corrupted

Examples

Import standard Refeyn data:

>>> contrasts, masses = import_file_h5('experiment.h5')
>>> print(f"Loaded {len(contrasts)} events")
>>> if masses is not None:
...     print(f"Mass range: {masses.min():.1f} - {masses.max():.1f} kDa")
pyphotomol.utils.data_import.import_csv(filename)[source]#

Import data from a CSV file generated by a Refeyn instrument.

The CSV file should contain columns ‘contrasts’ and optionally ‘masses_kDa’. NaN values are automatically filtered out from the imported data.

Parameters:

filename (str) – Path to the CSV file to import

Returns:

  • contrasts (np.ndarray) – Array of contrast values with NaN values removed

  • masses_kDa (np.ndarray or None) – Array of mass values in kDa, or None if not available in the CSV file

Raises:

  • FileNotFoundError – If the specified file does not exist

  • KeyError – If the required ‘contrasts’ column is not found in the CSV file

  • ValueError – If the file format is not recognized or data is corrupted

Examples

Import CSV data with contrasts only:

>>> contrasts, masses = import_csv('contrasts_only.csv')
>>> print(f"Loaded {len(contrasts)} contrast measurements")
>>> print(f"Masses available: {masses is not None}")

Import CSV data with both contrasts and masses:

>>> contrasts, masses = import_csv('full_data.csv')
>>> if masses is not None:
...     print(f"Mass range: {masses.min():.1f} - {masses.max():.1f} kDa")

pyphotomol.utils.helpers module#

pyphotomol.utils.helpers.contrasts_to_mass(contrasts, slope, intercept)[source]#

Function to convert masses from contrasts using known calibration parameters.

Caution! slope and intercept are based on f(mass) = contrast !!!! In other words, contrast = slope*mass + intercept

Parameters:

  • contrasts (np.ndarray) – Contrasts to convert.

  • slope (float) – Slope of the calibration line.

  • intercept (float) – Intercept of the calibration line.

Returns:

Converted masses in kDa.

Return type:

np.ndarray

pyphotomol.utils.helpers.create_histogram(vector, window=[0, 2000], bin_width=10)[source]#

Creates an histogram of the provided vector within a specified window and bin width.

Parameters:

  • vector (np.ndarray) – The data to create the histogram from.

  • window (list, default [0, 2000]) – The range of values to include in the histogram [min, max].

  • bin_width (float, default 10) – The width of each bin in the histogram.

Returns:

  • histogram_centers (np.ndarray) – The x-coordinates of the histogram bins.

  • hist_counts (np.ndarray) – The counts of values in each bin (Y-axis of the histogram).

  • hist_nbins (int) – The number of bins in the histogram.

Examples

For contrast data:

>>> centers, counts, nbins = create_histogram(contrasts, window=[-1, 0], bin_width=0.0004)

For mass data:

>>> centers, counts, nbins = create_histogram(masses, window=[0, 2000], bin_width=10)
pyphotomol.utils.helpers.count_binding_events_from_masses(masses_kDa)[source]#

Count binding events based on the provided mass data.

Parameters:

masses_kDa (np.ndarray) – Array of masses in kDa.

Returns:

  • n_binding (int) – Number of binding events.

  • n_unbinding (int) – Number of unbinding events.

pyphotomol.utils.helpers.count_binding_events_from_contrasts(contrasts)[source]#

Count binding events based on the provided contrast data.

Parameters:

contrasts (np.ndarray) – Array of contrasts.

Returns:

  • n_binding (int) – Number of binding events (negative contrasts).

  • n_unbinding (int) – Number of unbinding events (positive contrasts).

pyphotomol.utils.helpers.guess_peaks(x, histogram_centers, height=14, distance=4, prominence=8, masses=True)[source]#

Try to find peaks in the histogram data.

Automatically finds peaks in histogram data with adaptive parameters based on the data range. For mass data, different distance thresholds are used for different mass ranges to account for peak spacing variations.

Parameters:

  • x (np.ndarray) – The histogram counts data to find peaks in.

  • histogram_centers (np.ndarray) – The centers of the histogram bins corresponding to x values.

  • height (int, default 14) – Minimum height of peaks.

  • distance (int, default 4) – Minimum distance between peaks (will be scaled for different mass ranges).

  • prominence (int, default 8) – Minimum prominence of peaks.

  • masses (bool, default True) – If True, find peaks in mass data; if False, find peaks in contrast data.

Returns:

Histogram centers of the found peaks.

Return type:

np.ndarray

Examples

For contrast data:

>>> peaks = guess_peaks(hist_counts, hist_centers, height=10, distance=4, prominence=4, masses=False)

For mass data:

>>> peaks = guess_peaks(hist_counts, hist_centers, height=14, distance=4, prominence=8, masses=True)
pyphotomol.utils.helpers.multi_gauss(x, *params)[source]#

Multiple gaussian function Inputs x values and gaussian parameters Outputs y values

Parameters:

  • x (np.ndarray) – The x-coordinates where the Gaussian function is evaluated.

  • *params (tuple) – Parameters for the Gaussian functions, where each set of parameters consists of: - Center (mean) of the Gaussian - Amplitude (height) of the Gaussian - Standard deviation (width) of the Gaussian

Returns:

The y-coordinates of the multiple-Gaussian function evaluated at x.

Return type:

np.ndarray

Examples

>>> x = np.linspace(0, 10, 100)
>>> params = [5, 1, 0.5, 7, 0.5, 0.2]  # Two Gaussians
>>> y = multi_gauss(x, *params)
pyphotomol.utils.helpers.truncated_multi_gauss_with_baseline(x, lower_limit_of_resolution, baseline, *params)[source]#

Multiple gaussian function with a baseline and a lower limit of resolution. Inputs x values, gaussian parameters, lower limit of resolution, and baseline. Outputs y values.

Parameters:

  • x (np.ndarray) – The x-coordinates where the Gaussian function is evaluated.

  • lower_limit_of_resolution (float) – The lower limit of resolution below which the function is zero

  • baseline (float) – The baseline value to be added to the Gaussian function.

  • *params (tuple) – Parameters for the Gaussian functions, where each set of parameters consists of: - Center (mean) of the Gaussian - Amplitude (height) of the Gaussian - Standard deviation (width) of the Gaussian

Returns:

The y-coordinates of the multiple-Gaussian function evaluated at x, truncated below the lower limit of resolution and with a baseline added.

Return type:

np.ndarray

pyphotomol.utils.helpers.fit_histogram(hist_counts, hist_centers, guess_positions=[66, 148, 480], mean_tolerance=None, std_tolerance=None, threshold=40, baseline=0, masses=True, fit_baseline=False)[source]#

Fit a histogram with multiple truncated gaussians.

Parameters:

  • hist_counts (np.ndarray) – The counts of values in each bin of the histogram.

  • hist_centers (np.ndarray) – The centers of the histogram bins.

  • guess_positions (list, default [66,148,480]) – Initial guesses for the positions of the peaks.

  • mean_tolerance (int, default 100) – Tolerance for the peak positions. If None, it will be copied from guess_positions.

  • std_tolerance (int, default 200) – Maximum standard deviation for the peaks. If None, it will be copied from guess_positions.

  • threshold (int, default 40 for masses in kDa units) – For masses, minimum value that can be observed. For contrasts, it is be the max value that can be observed. It should be a negative value.

  • baseline (float, default 0) – Baseline value to be added to the fit.

  • masses (bool, default True) – If True, the fit is for mass data; if False, it is for contrast data.

  • fit_baseline (bool, default False) – If True, the fit will include a baseline parameter. The ‘baseline’ argument will be ignored.

Returns:

  • popt (np.ndarray) – Optimized parameters for the fit.

  • fit (np.ndarray) – Fitted values for the histogram. The first column is the x-coordinates, followed by the individual Gaussian fits and the total fit.

  • fit_error (np.ndarray) – Errors of the fitted parameters.

Examples

For contrast data:

>>> popt, fit, errors = fit_histogram(counts, centers, mean_tolerance=0.05, std_tolerance=0.1)
pyphotomol.utils.helpers.create_fit_table(popt, popt_error, fit, n_binding, n_unbinding, hist_centers, masses=True, include_errors=True)[source]#

Generate a pandas DataFrame that summarizes fit results

Parameters:

  • popt (np.ndarray) – Optimized parameters from the fit.

  • popt_error (np.ndarray) – Errors of the fitted parameters.

  • fit (np.ndarray) – Fitted values for the histogram.

  • n_binding (int) – Number of binding events.

  • n_unbinding (int) – Number of unbinding events.

  • hist_centers (np.ndarray) – The centers of the histogram bins.

  • masses (bool, default True) – If True, the fit is for mass data; if False, it is for contrast data.

  • include_errors (bool, default True) – If True, include errors in the fit table.

Returns:

fit_table – DataFrame containing the fit results.

Return type:

pd.DataFrame

pyphotomol.utils.helpers.r_squared(data, fit)[source]#

Compute the R-squared value for the fit.

Parameters:

  • data (np.ndarray) – The original data.

  • fit (np.ndarray) – The fitted values.

Returns:

The R-squared value.

Return type:

float

pyphotomol.utils.helpers.calibrate(calib_floats, fit_table)[source]#

Calibration based on contrasts histogram

Parameters:

  • calib_floats (list) – List of calibration standards in kDa (e.g. [66, 146, 480]).

  • fit_table (pd.DataFrame) – DataFrame containing the fit results. Created by create_fit_table.

Returns:

calibration_dic – Dictionary containing the calibration results: - ‘standards’: Calibration standards used. - ‘exp_points’: Expected points from the fit. - ‘fit_params’: Parameters of the fit. - ‘fit_r2’: R-squared value of the fit.

Return type:

dict

pyphotomol.utils.helpers.evaluate_fitting(fitted_values, low_bounds, high_bounds, tolerance=0.02)[source]#

Evaluate the quality of a constrained fitting by checking if the fitted values are within the specified bounds.

Parameters:

  • fitted_values (np.ndarray) – The values obtained from the fitting process.

  • low_bounds (np.ndarray) – The lower bounds for the fitted values.

  • high_bounds (np.ndarray) – The upper bounds for the fitted values.

Returns:

True if all fitted values are within the specified bounds False otherwise.

Return type:

bool

pyphotomol.utils.palette module#

Color palettes for PyPhotoMol plotting.

pyphotomol.utils.plotting module#

Plotting utilities for mass photometry data visualization.

class pyphotomol.utils.plotting.PlotConfig(plot_width: int = 1000, plot_height: int = 400, plot_type: str = 'png', font_size: int = 14, normalize: bool = False, contrasts: bool = False, cst_factor_for_contrast: float = 1, x_range: List[float] | None = None)[source]#

Bases: object

General plot configuration

plot_width: int = 1000#
plot_height: int = 400#
plot_type: str = 'png'#
font_size: int = 14#
normalize: bool = False#
contrasts: bool = False#
cst_factor_for_contrast: float = 1#
x_range: List[float] | None = None#
__init__(plot_width: int = 1000, plot_height: int = 400, plot_type: str = 'png', font_size: int = 14, normalize: bool = False, contrasts: bool = False, cst_factor_for_contrast: float = 1, x_range: List[float] | None = None) → None#
class pyphotomol.utils.plotting.AxisConfig(showgrid_x: bool = True, showgrid_y: bool = True, n_y_axis_ticks: int = 3, axis_linewidth: int = 1, axis_tickwidth: int = 1, axis_gridwidth: int = 1)[source]#

Bases: object

Axis styling configuration

showgrid_x: bool = True#
showgrid_y: bool = True#
n_y_axis_ticks: int = 3#
axis_linewidth: int = 1#
axis_tickwidth: int = 1#
axis_gridwidth: int = 1#
__init__(showgrid_x: bool = True, showgrid_y: bool = True, n_y_axis_ticks: int = 3, axis_linewidth: int = 1, axis_tickwidth: int = 1, axis_gridwidth: int = 1) → None#
class pyphotomol.utils.plotting.LayoutConfig(stacked: bool = False, show_subplot_titles: bool = False, vertical_spacing: float = 0.1, shared_yaxes: bool = True, extra_padding_y_label: float = 0)[source]#

Bases: object

Layout and spacing configuration

stacked: bool = False#
show_subplot_titles: bool = False#
vertical_spacing: float = 0.1#
shared_yaxes: bool = True#
extra_padding_y_label: float = 0#
__init__(stacked: bool = False, show_subplot_titles: bool = False, vertical_spacing: float = 0.1, shared_yaxes: bool = True, extra_padding_y_label: float = 0) → None#
class pyphotomol.utils.plotting.LegendConfig(add_masses_to_legend: bool = True, add_percentage_to_legend: bool = False, add_labels: bool = True, add_percentages: bool = True, line_width: int = 3)[source]#

Bases: object

Legend and labeling configuration

add_masses_to_legend: bool = True#
add_percentage_to_legend: bool = False#
add_labels: bool = True#
add_percentages: bool = True#
line_width: int = 3#
__init__(add_masses_to_legend: bool = True, add_percentage_to_legend: bool = False, add_labels: bool = True, add_percentages: bool = True, line_width: int = 3) → None#
pyphotomol.utils.plotting.config_fig(fig, plot_width=800, plot_height=600, plot_type='png', plot_title_for_download='plot')[source]#

Configure plotly figure with download options and toolbar settings.

Parameters:

  • fig (go.Figure) – Plotly figure object

  • plot_width (int, default 800) – Width of the plot in pixels

  • plot_height (int, default 600) – Height of the plot in pixels

  • plot_type (str, default "png") – Format for downloading the plot (e.g., “png”, “jpeg”)

  • plot_title_for_download (str, default "plot") – Title for the downloaded plot file

Returns:

Configured plotly figure

Return type:

go.Figure

pyphotomol.utils.plotting.create_axis_config(title=None, font_size=14, show_grid=True, axis_linewidth=1, axis_tickwidth=1, axis_gridwidth=1, nticks=None)[source]#

Create a standardized axis configuration dictionary for Plotly figures.

Parameters:

  • title (str or None) – Title for the axis

  • font_size (int, default 14) – Font size for axis title and tick labels

  • show_grid (bool, default True) – Whether to show grid lines

  • axis_linewidth (int, default 1) – Width of the axis line

  • axis_tickwidth (int, default 1) – Width of the tick marks

  • axis_gridwidth (int, default 1) – Width of the grid lines

  • nticks (int or None) – Number of ticks on the axis (optional guideline)

Returns:

Axis configuration dictionary for Plotly

Return type:

dict

pyphotomol.utils.plotting.add_histogram(fig, values, bin_info, normalize=False, hex_color='#377EB8')[source]#

Add mass histogram trace to plotly figure.

Parameters:

  • fig (go.Figure) – Plotly figure object

  • values (np.ndarray) – Array of mass or contrast values

  • bin_info (dict) – Dictionary containing histogram bin information

  • normalize (bool) – Whether to normalize the histogram

  • hex_color (str) – Hex color code for the histogram bars

Returns:

Updated plotly figure with histogram trace

Return type:

go.Figure

pyphotomol.utils.plotting.add_labels_to_fig(fig, fit_table, contrasts=False, scaling_factor=1.0, font_size=14, sels=None, stacked=False, max_y=None, show_values=True, show_counts=True, cst_factor_for_contrast=1000, subplot_row=None)[source]#

Add peak labels to plotly figure.

Parameters:

  • fig (go.Figure) – Plotly figure object

  • fit_table (pd.DataFrame) – DataFrame containing fit results

  • contrasts (bool) – Whether data represents contrasts

  • scaling_factor (float) – Factor to scale y values

  • font_size (int) – Font size for axes

  • sels (list or None) – Boolean list for selecting which peaks to label

  • stacked (bool) – Whether this is a stacked plot

  • max_y (float or None) – Maximum y value for label positioning, required if the y-axis range is shared

  • show_values (bool) – Whether to include mass / contrast values in labels

  • show_counts (bool) – Whether to include count percentages in labels

  • cst_factor_for_contrast (float) – Factor to convert contrasts for display

  • subplot_row (int or None) – Row number for stacked subplots (1-indexed), None for single plot

Returns:

Updated plotly figure with peak labels

Return type:

go.Figure

pyphotomol.utils.plotting.add_values_in_legend(legends, fit_table, contrasts=False, cst_factor_for_contrast=1000)[source]#

Add mass / contrast values to legend labels.

Parameters:

  • legends (list) – List of legend labels

  • fit_table (pd.DataFrame) – DataFrame with fit table data

  • contrasts (bool) – Whether data represents contrasts

  • cst_factor_for_contrast (float) – Factor to convert contrasts for display

Returns:

Updated legend labels with mass values

Return type:

list

pyphotomol.utils.plotting.add_gaussian_traces(fig, arr_fitted_values, fitted_params_table, legends, color_palette, sels, baseline=0, scaling_factor=1.0, show_values_in_legend=True, show_perc_in_legend=True, contrasts=False, cst_factor_for_contrast=1000, line_width=3, show_legend_vec=[])[source]#

Add individual Gaussian traces to plotly figure.

Parameters:

  • fig (go.Figure) – Plotly figure object

  • arr_fitted_values (np.ndarray) – Array of fitted values (x, individual Gaussians, sum)

  • fitted_params_table (pd.DataFrame) – DataFrame with fit table results

  • legends (list) – List of legend labels - Same length as columns in arr_fitted_values

  • color_palette (list) – List of hex color codes

  • sels (list) – Boolean list for selecting which traces to show

  • baseline (float) – Baseline value for filtering

  • scaling_factor (float) – Factor to scale y values

  • show_values_in_legend (bool) – Whether to add masses / contrast values to legend labels

  • show_perc_in_legend (bool) – Whether to add percentages to legend labels

  • contrasts (bool) – Whether data represents contrasts

  • add_labels (bool) – Whether to add peak labels

  • add_percentages (bool) – Whether to add percentage labels

  • stacked (bool) – Whether this is a stacked plot

  • cst_factor_for_contrast (float) – Factor to convert contrasts for display

  • subplot_row (int or None) – Row number for stacked subplots (1-indexed), None for single plot

  • line_width (int, default 3) – Width of the Gaussian trace lines

  • show_legend_vec (list, default []) – Boolean list for showing Gaussian traces in legend

Returns:

Updated plotly figure with Gaussian traces

Return type:

go.Figure

pyphotomol.utils.plotting.add_percentages_to_legend(legends, fit_table)[source]#

Add percentage values to legend labels.

Parameters:

  • legends (list) – List of legend labels

  • fit_table (pd.DataFrame) – DataFrame with fit table data

Returns:

Updated legend labels with percentage values

Return type:

list

pyphotomol.utils.plotting.add_gaussian_simulation_to_plot(fig, mean, std, amplitude, left_bound)[source]#

Add simulated truncated Gaussian trace to plotly figure.

Parameters:

  • fig (go.Figure) – Plotly figure object

  • mean (float) – Mean of the Gaussian distribution

  • std (float) – Standard deviation of the Gaussian distribution

  • amplitude (float) – Amplitude of the Gaussian distribution

  • left_bound (float) – Left boundary for truncation

Returns:

Updated plotly figure with simulation trace

Return type:

go.Figure

pyphotomol.utils.plotting.plot_histograms_and_fits(analyzer, legends_df=None, colors_hist=None, plot_config: PlotConfig = None, axis_config: AxisConfig = None, layout_config: LayoutConfig = None, legend_config: LegendConfig = None)[source]#

Create a comprehensive plot of PhotoMol fit data with histograms and Gaussian traces.

Parameters:

  • analyzer (pyphotomol.MPAnalyzer or pyphotomol.PyPhotoMol) – MPAnalyzer instance containing multiple PyPhotoMol models - or a single PyPhotoMol instance

  • legends_df (pd.DataFrame, optional) – DataFrame containing legends, colors, and selections with columns [‘legends’, ‘color’, ‘select’, ‘show_legend’] This dataframe affects the fitted curves only, not the histograms.

  • colors_hist (list, str, or pd.DataFrame, optional) – List of colors for histograms (one per model) If a string, it will be used for all histograms. If a DataFrame, it should have a column ‘color’ with hex color codes.

  • plot_config (PlotConfig, optional) – General plot configuration (dimensions, format, contrasts, etc.)

  • axis_config (AxisConfig, optional) – Axis styling configuration (grid, line widths, etc.)

  • layout_config (LayoutConfig, optional) – Layout configuration (stacked, spacing, etc.)

  • legend_config (LegendConfig, optional) – Legend and labeling configuration

Returns:

Configured plotly figure object

Return type:

go.Figure

Examples

Simple plot with default settings:

>>> fig = plot_histograms_and_fits(analyzer, colors_hist=['blue', 'red'])
>>> fig.show()

Customized plot with configuration objects:

>>> plot_config = PlotConfig(plot_width=800, contrasts=True, x_range=[0, 500])
>>> layout_config = LayoutConfig(stacked=True, vertical_spacing=0.05)
>>> fig = plot_histograms_and_fits(analyzer, plot_config=plot_config,
...                               layout_config=layout_config)
>>> fig.show()

Plot with custom x-axis limits:

>>> plot_config = PlotConfig(x_range=[100, 800])  # Zoom to 100-800 kDa range
>>> fig = plot_histograms_and_fits(analyzer, plot_config=plot_config)
>>> fig.show()
pyphotomol.utils.plotting.plot_histogram(analyzer, colors_hist=None, plot_config: PlotConfig = None, axis_config: AxisConfig = None, layout_config: LayoutConfig = None)[source]#

Create a plot with only histograms from PhotoMol data (wrapper around plot_histograms_and_fits).

This function is a simplified wrapper that creates histogram-only plots without requiring fitted data or legend configuration.

Parameters:

  • analyzer (pyphotomol.MPAnalyzer or pyphotomol.PyPhotoMol) – MPAnalyzer instance containing multiple PyPhotoMol models or a single PyPhotoMol instance

  • colors_hist (list, optional) – List of colors for histograms (one per model)

  • plot_config (PlotConfig, optional) – General plot configuration (dimensions, format, contrasts, etc.)

  • axis_config (AxisConfig, optional) – Axis styling configuration (grid, line widths, etc.)

  • layout_config (LayoutConfig, optional) – Layout configuration (stacked, spacing, etc.)

Returns:

Configured plotly figure object with histograms only

Return type:

go.Figure

Examples

Create a simple histogram plot:

>>> fig = plot_histogram(analyzer, ['#FF5733', '#33C3FF'])
>>> fig.show()

Create stacked normalized histograms:

>>> plot_config = PlotConfig(normalize=True)
>>> layout_config = LayoutConfig(stacked=True)
>>> fig = plot_histogram(analyzer, ['blue', 'red'],
...                      plot_config=plot_config, layout_config=layout_config)
>>> fig.show()
pyphotomol.utils.plotting.plot_calibration(mass, contrast, slope, intercept, plot_config: PlotConfig = None, axis_config: AxisConfig = None)[source]#

Create a scatter plot of mass vs contrast with calibration line.

This function creates a visualization showing the relationship between mass and ratiometric contrast, with a fitted calibration line overlaid. This is useful for visualizing calibration quality and outliers.

Parameters:

  • mass (array-like) – Array of mass values in kDa

  • contrast (array-like) – Array of corresponding ratiometric contrast values

  • slope (float) – Slope of the calibration line (contrast = slope * mass + intercept)

  • intercept (float) – Intercept of the calibration line

  • plot_config (PlotConfig, optional) – General plot configuration (dimensions, format, axis size, etc.)

  • axis_config (AxisConfig, optional) – Axis styling configuration (grid, line widths, etc.)

Returns:

Plotly figure object containing the mass vs contrast calibration plot

Return type:

go.Figure

Examples

Plot mass vs contrast calibration:

>>> import numpy as np
>>> from pyphotomol.utils.plotting import plot_calibration, PlotConfig, AxisConfig
>>>
>>> # Simulated calibration data
>>> mass = np.array([66, 146, 480])
>>> contrast = np.array([-0.1, -0.2, -0.5])
>>> slope = -0.001
>>> intercept = 0.02
>>>
>>> # Simple plot with defaults
>>> fig = plot_calibration(mass, contrast, slope, intercept)
>>> fig.show()
>>>
>>> # Customized plot
>>> plot_config = PlotConfig(plot_width=600, plot_height=400, font_size=12)
>>> axis_config = AxisConfig(showgrid_x=False, n_y_axis_ticks=6)
>>> fig = plot_calibration(mass, contrast, slope, intercept,
...                       plot_config=plot_config, axis_config=axis_config)
>>> fig.show()