|
1 | 1 | [↩ back to index](../README.md) |
2 | 2 | # Estimating parameters of ion current models from whole-cell voltage-clamp data |
3 | 3 |
|
4 | | -In this tutorial, we look at the problem of estimating the parameters of an ion current model from whole-cell voltage-clamp data. |
| 4 | +In these notebooks, we look at the problem of estimating the parameters of an ion current model from whole-cell voltage-clamp data. |
5 | 5 |
|
6 | 6 | As we go along, we'll create some classes and utility functions that may be useful in general. |
7 | 7 | These are all stored in [library.py](./library.py). |
8 | 8 |
|
9 | 9 | The follow topics are covered: |
10 | 10 |
|
11 | | -## [Introduction](introduction.ipynb) |
| 11 | +## [Introduction](introduction.ipynb) - [nbviewer](https://nbviewer.jupyter.org/github/CardiacModelling/fitting-notebooks/blob/main/ion-currents/introduction.ipynb) |
12 | 12 |
|
13 | | -This tutorial provides some background on the model we'll use throughout the tutorial. |
| 13 | +This notebook provides some background on the model we'll use in all examples. |
14 | 14 | It also introduces a first voltage-protocol (a simplified variant of the "staircase protocol"). |
15 | 15 |
|
16 | | -## [Basic simulations](basic-simulations.ipynb) |
| 16 | +## [Basic simulations](basic-simulations.ipynb) - [nbviewer](https://nbviewer.jupyter.org/github/CardiacModelling/fitting-notebooks/blob/main/ion-currents/basic-fitting.ipynb) |
17 | 17 |
|
18 | | -This tutorial shows how Myokit can be used to simulate patch-clamp experiments. |
| 18 | +This notebook shows how Myokit can be used to simulate patch-clamp experiments. |
19 | 19 | It shows you how to create a simulation from a model and protocol stored on disk, and discusses how to change model parameters. |
20 | 20 | Finally, it shows how steady-states can be calculated and set as initial conditions. |
21 | 21 |
|
22 | 22 | ## [Basic fitting](basic-fitting.ipynb) |
23 | 23 |
|
24 | | -In this tutorial we link Myokit to PINTS. |
| 24 | +In this notebook we link Myokit to PINTS. |
25 | 25 | Noise models are discussed and synthetic data is generated, after which an error measure is defined and minimised. |
26 | 26 | Inspecting the results, we show how tight solver tolerances are needed for fitting, and how the finite size of our experimental time series can cause a slight bias in the results. |
27 | 27 |
|
28 | 28 | ## Fitting to different voltage protocols |
29 | 29 |
|
30 | 30 | The next four notebooks discuss different voltage protocols, and the simulation methods appropriate to each one. |
31 | 31 |
|
32 | | -- [Combining step protocols with sine waves or ramps](more-protocols-1-steps-and-ramps.ipynb) |
| 32 | +- [Combining step protocols with sine waves or ramps](more-protocols-1-steps-and-ramps.ipynb) - [nbviewer](https://nbviewer.jupyter.org/github/CardiacModelling/fitting-notebooks/blob/main/ion-currents/more-protocols-1-steps-and-ramps.ipynb) |
33 | 33 | - [Simulating an AP protocol with "data clamp"](more-protocols-2-data-clamp.ipynb) |
34 | 34 | - [Analytical solvers for simple step protocols](more-protocols-3-analytic-solvers.ipynb) |
35 | 35 | - [Fitting to multiple simple step protocols](more-protocols-4-multiple-protocols.ipynb) |
36 | 36 |
|
37 | 37 | ## [Setting boundaries on model parameters](boundaries.ipynb) |
38 | 38 |
|
39 | | -In this tutorial we show how some parameters can cause numerical issues during simulation, and how we can catch and report these errors. |
| 39 | +In this notebook we show how some parameters can cause numerical issues during simulation, and how we can catch and report these errors. |
40 | 40 | We then inspect the model equations and use previous estimates of our parameters (or quantities related to the parameters) to define some very wide boundaries, or "prior estimates". |
41 | 41 | Finally, we show how we can use this kind of reasoning to define univariate boundaries (one on each parameter), and multivariate boundaries (which restrict the maximum rate coefficients seen during a simulation). |
42 | 42 |
|
43 | 43 | ## [Selecting starting points for an optimisation](starting-points.ipynb) |
44 | 44 |
|
45 | | -Continuing from the previous tutorial, this tutorial shows how we can sample from within the (univariate and multivariate) boundaries to select starting points for an optisiation. |
46 | | -At the end of this tutorial we briefly discuss a "repeated-fits" strategy which allows you to test the reliability of results obtained on real data, where the "true" parameters are not known. |
| 45 | +Continuing from the previous chapters, this notebook shows how we can sample from within the (univariate and multivariate) boundaries to select starting points for an optisiation. |
| 46 | +At the end of this notebook we briefly discuss a "repeated-fits" strategy which allows you to test the reliability of results obtained on real data, where the "true" parameters are not known. |
47 | 47 |
|
48 | 48 | ## [Searching in a transformed space](transformations.ipynb) |
49 | 49 |
|
50 | | -This tutorial shows how you can create wrappers around models and boundaries to run optimisations on a transformed parameter space. |
| 50 | +This notebook shows how you can create wrappers around models and boundaries to run optimisations on a transformed parameter space. |
51 | 51 |
|
52 | 52 | ## [Running big fitting experiments](big-fitting.ipynb) |
53 | 53 |
|
54 | | -This tutorial focusses on the practical side of fitting. |
| 54 | +This notebook focusses on the practical side of fitting. |
55 | 55 | It introduces methods to store simulation results to disk, load and analyse them, and shows a way to "reserve" filenames when multiple processes are running at once. |
56 | 56 | It ends with a brief note on multiprocessing. |
57 | 57 |
|
|
0 commit comments