Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories
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...
| Main Authors: | , , |
|---|---|
| Format: | Article |
| Language: | English |
| Published: |
American Physical Society
2023-10-01
|
| Series: | Physical Review X |
| Online Access: | http://doi.org/10.1103/PhysRevX.13.041017 |
| _version_ | 1797648770463170560 |
|---|---|
| author | Manuel Reinhardt Gašper Tkačik Pieter Rein ten Wolde |
| author_facet | Manuel Reinhardt Gašper Tkačik Pieter Rein ten Wolde |
| author_sort | Manuel Reinhardt |
| collection | DOAJ |
| description | 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. |
| first_indexed | 2024-03-11T15:36:49Z |
| format | Article |
| id | doaj.art-052d49c4982841afb3c0a33436e2993a |
| institution | Directory Open Access Journal |
| issn | 2160-3308 |
| language | English |
| last_indexed | 2024-03-11T15:36:49Z |
| publishDate | 2023-10-01 |
| publisher | American Physical Society |
| record_format | Article |
| series | Physical Review X |
| spelling | doaj.art-052d49c4982841afb3c0a33436e2993a2023-10-26T15:06:44ZengAmerican Physical SocietyPhysical Review X2160-33082023-10-0113404101710.1103/PhysRevX.13.041017Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic TrajectoriesManuel ReinhardtGašper TkačikPieter Rein ten WoldeMost 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.http://doi.org/10.1103/PhysRevX.13.041017 |
| spellingShingle | Manuel Reinhardt Gašper Tkačik Pieter Rein ten Wolde Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories Physical Review X |
| title | Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories |
| title_full | Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories |
| title_fullStr | Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories |
| title_full_unstemmed | Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories |
| title_short | Path Weight Sampling: Exact Monte Carlo Computation of the Mutual Information between Stochastic Trajectories |
| title_sort | path weight sampling exact monte carlo computation of the mutual information between stochastic trajectories |
| url | http://doi.org/10.1103/PhysRevX.13.041017 |
| work_keys_str_mv | AT manuelreinhardt pathweightsamplingexactmontecarlocomputationofthemutualinformationbetweenstochastictrajectories AT gaspertkacik pathweightsamplingexactmontecarlocomputationofthemutualinformationbetweenstochastictrajectories AT pieterreintenwolde pathweightsamplingexactmontecarlocomputationofthemutualinformationbetweenstochastictrajectories |