dataprob 0.9.0

Simple example

The following code generates noisy linear data and uses dataprob to find the maximum likelihood estimate of its slope and intercept.

import dataprob
import numpy as np

# Generate "experimental" linear data (slope = 5, intercept = 5.7) that has
# random noise on each point.
x_array = np.linspace(0,10,25)
noise = np.random.normal(loc=0,scale=0.5,size=x_array.shape)
y_obs = 5*x_array + 5.7 + noise

# 1. Define a linear model
def linear_model(m=1,b=1,x=[]):
    return m*x + b

# 2. Set up the analysis. 'method' can be "ml", "mcmc", or "bootstrap"
f = dataprob.setup(linear_model,
                   method="ml",
                   non_fit_kwargs={"x":x_array})

# 3. Fit the parameters of linear_model model to y_obs, assuming uncertainty
#    of 0.5 on each observed point.
f.fit(y_obs=y_obs,
      y_std=0.5)

# 4. Access results
print(f.fit_df)
fig = dataprob.plot_summary(f)
fig = dataprob.plot_corner(f)

The f.fit_df dataframe will look something like:

The plots will be:

Installation

We recommend installing dataprob with pip:

pip install dataprob

To install from source and run tests:

git clone https://github.com/harmslab/dataprob.git
cd dataprob
pip install .

# to run test-suite
pytest --runslow

Examples

A good way to learn how to use the library is by working through examples. The following notebooks are included in the dataprob/examples/ directory. They are self-contained demonstrations in which dataprob is used to analyze various classes of experimental data. The links below launch each notebook in Google colab:

• linear.ipynb: fit a linear model to noisy data (2 parameter, linear)

• binding.ipynb: a single-site binding interaction (2 parameter, sigmoidal curve)

• michaelis-menten.ipynb: Michaelis-Menten model of enzyme kinetics (2 parameter, sigmoidal curve)

• lagged-exponential.ipynb: bacterial growth curve with initial lag phase (3 parameter, exponential)

• multi-gaussian.ipynb: two overlapping normal distributions (6 parameter, Gaussian)

• periodic.ipynb: periodic data (3 parameter, sine)

• polynomial.ipynb: nonlinear data with no obvious form (5 parameter, polynomial)

• linear-extrapolation-folding.ipynb: protein equilibrium unfolding data (6 parameter, linear embedded in sigmoidal)

Documentation

Full documentation is on readthedocs.