Supplementary data, and code, for
Preprint available on bioRxiv.
Raw data in the /Data
folder corresponds to experiment indices in the main document as follows:
Index | Folder |
---|---|
1 | 8000_E2 |
2 | 8000_G2 |
3 | 10000_H1 |
4 | 10000_H2 |
5 | 8000_F2 |
6 | 10000_A2 |
7 | 12000_H3 |
8 | 12000_F3 |
9 | 12000_D2 |
Experimental data is recorded 12 h after seeding. Hence, raw data indicated by filename ending corresponds to the following times:
Filename end | Time |
---|---|
"_12h" | 0 h |
"_30h" | 18 h |
"_48h" | 36 h |
For more details, see ExperimentalData.xlsx
- MATLAB R2018a or later
- MEX and C compiler
- (optional) Parallel computing toolbox
If you do not have this toolbox, you must replace
parfor
withfor
in the following files:ABC_SimulateModel.m
ABC_SMC.m
Before use, you must compile the MEX code contained in IBM.c
using the
Install.m
script or typing:
mex IBM.c -R2018a
into the command window.
See Example_IBM.m
for an example IBM realisation.
- Parallel computing toolbox highly recommend
- This function can take a large amount of time to run for
Nsamples > 100
- We have set
Nsamples = 100
andalpha = 0.1
for demonstrative purposes
- We have set
Run Example_ABCRejection.m
code (in the main document, we use Nsamples = 100000
)
- Parallel computing toolbox highly recommend
- This function can take a large amount of time to run for
N > 10
- We have set
N = 8
and not used the full sequence of thresholds for demonstrative purposes
- We have set
Run Example_ABCSMC.m
code (in the main document, we use Nsamples = 5000
)