article
Path weight sampling: Exact Monte Carlo computation of the mutual information between stochastic trajectories
published
yes
Manuel
Reinhardt
author
Gašper
Tkačik
author 3D494DCA-F248-11E8-B48F-1D18A9856A870000-0002-6699-1455
Pieter Rein
Ten Wolde
author
GaTk
department
Most natural and engineered information-processing systems transmit information via signals that vary in time. Computing the information transmission rate or the information encoded in the temporal characteristics of these signals requires the mutual information between the input and output signals as a function of time, i.e., between the input and output trajectories. Yet, this is notoriously difficult because of the high-dimensional nature of the trajectory space, and all existing techniques require approximations. We present an exact Monte Carlo technique called path weight sampling (PWS) that, for the first time, makes it possible to compute the mutual information between input and output trajectories for any stochastic system that is described by a master equation. The principal idea is to use the master equation to evaluate the exact conditional probability of an individual output trajectory for a given input trajectory and average this via Monte Carlo sampling in trajectory space to obtain the mutual information. We present three variants of PWS, which all generate the trajectories using the standard stochastic simulation algorithm. While direct PWS is a brute-force method, Rosenbluth-Rosenbluth PWS exploits the analogy between signal trajectory sampling and polymer sampling, and thermodynamic integration PWS is based on a reversible work calculation in trajectory space. PWS also makes it possible to compute the mutual information between input and output trajectories for systems with hidden internal states as well as systems with feedback from output to input. Applying PWS to the bacterial chemotaxis system, consisting of 182 coupled chemical reactions, demonstrates not only that the scheme is highly efficient but also that the number of receptor clusters is much smaller than hitherto believed, while their size is much larger.
https://research-explorer.ista.ac.at/download/14515/14522/2023_PhysReviewX_Reinhardt.pdf
application/pdfno
https://creativecommons.org/licenses/by/4.0/
American Physical Society2023
eng
Physical Review X
2160-3308
2203.0346110.1103/PhysRevX.13.041017
134
Reinhardt, Manuel, Gašper Tkačik, and Pieter Rein Ten Wolde. “Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories.” <i>Physical Review X</i>. American Physical Society, 2023. <a href="https://doi.org/10.1103/PhysRevX.13.041017">https://doi.org/10.1103/PhysRevX.13.041017</a>.
Reinhardt M, Tkačik G, Ten Wolde PR. Path weight sampling: Exact Monte Carlo computation of the mutual information between stochastic trajectories. <i>Physical Review X</i>. 2023;13(4). doi:<a href="https://doi.org/10.1103/PhysRevX.13.041017">10.1103/PhysRevX.13.041017</a>
Reinhardt M, Tkačik G, Ten Wolde PR. 2023. Path weight sampling: Exact Monte Carlo computation of the mutual information between stochastic trajectories. Physical Review X. 13(4), 041017.
Reinhardt, M., Tkačik, G., & Ten Wolde, P. R. (2023). Path weight sampling: Exact Monte Carlo computation of the mutual information between stochastic trajectories. <i>Physical Review X</i>. American Physical Society. <a href="https://doi.org/10.1103/PhysRevX.13.041017">https://doi.org/10.1103/PhysRevX.13.041017</a>
M. Reinhardt, G. Tkačik, and P. R. Ten Wolde, “Path weight sampling: Exact Monte Carlo computation of the mutual information between stochastic trajectories,” <i>Physical Review X</i>, vol. 13, no. 4. American Physical Society, 2023.
M. Reinhardt, G. Tkačik, P.R. Ten Wolde, Physical Review X 13 (2023).
Reinhardt, Manuel, et al. “Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories.” <i>Physical Review X</i>, vol. 13, no. 4, 041017, American Physical Society, 2023, doi:<a href="https://doi.org/10.1103/PhysRevX.13.041017">10.1103/PhysRevX.13.041017</a>.
145152023-11-12T23:00:55Z2023-11-13T09:03:30Z