Journal article icon

Journal article

Profile likelihood analysis for a stochastic model of diffusion in heterogeneous media

Abstract:
We compute profile likelihoods for a stochastic model of diffusive transport motivated by experimental observations of heat conduction in layered skin tissues. This process is modelled as a random walk in a layered one-dimensional material, where each layer has a distinct particle hopping rate. Particles are released at some location, and the duration of time taken for each particle to reach an absorbing boundary is recorded. To explore whether this data can be used to identify the hopping rates in each layer, we compute various profile likelihoods using two methods: first, an exact likelihood is evaluated using a relatively expensive Markov chain approach; and, second we form an approximate likelihood by assuming the distribution of exit times is given by a Gamma distribution whose first two moments match the expected moments from the continuum limit description of the stochastic model. Using the exact and approximate likelihoods we construct various profile likelihoods for a range of problems. In cases where parameter values are not identifiable, we make progress by re-interpreting those data with a reduced model with a smaller number of layers.
Publication status:
Published
Peer review status:
Peer reviewed

Actions

Access Document

Files:
Publisher copy:
10.1098/rspa.2021.0214

Authors


Publisher:
The Royal Society
Journal:
Proceedings of the Royal Society A: Mathematical, Physical and Engineering Sciences More from this journal
Volume:
447
Issue:
2250
Article number:
20210214
Publication date:
2021-06-09
Acceptance date:
2021-05-10
DOI:
EISSN:
1471-2946
ISSN:
1364-5021


Language:
English
Keywords:
Pubs id:
1147859
Local pid:
pubs:1147859
Deposit date:
2021-05-27
ARK identifier:

Terms of use


Views and Downloads






If you are the owner of this record, you can report an update to it here: Report update to this record

TO TOP