Papers › Modeling Latent Neural Dynamics with Gaussian Process Switching Linear Dynamical Systems

Modeling Latent Neural Dynamics with Gaussian Process Switching Linear Dynamical Systems

19 Jul 2024arXiv:2408.03330archive 2025-07-28

Amber Hu, David Zoltowski, Aditya Nair, David Anderson, Lea Duncker, Scott Linderman

Understanding how the collective activity of neural populations relates to computation and ultimately behavior is a key goal in neuroscience. To this end, statistical methods which describe high-dimensional neural time series in terms of low-dimensional latent dynamics have played a fundamental role in characterizing neural systems. Yet, what constitutes a successful method involves two opposing criteria: (1) methods should be expressive enough to capture complex nonlinear dynamics, and (2) they should maintain a notion of interpretability often only warranted by simpler linear models. In this paper, we develop an approach that balances these two objectives: the Gaussian Process Switching Linear Dynamical System (gpSLDS). Our method builds on previous work modeling the latent state evolution via a stochastic differential equation whose nonlinear dynamics are described by a Gaussian process (GP-SDEs). We propose a novel kernel function which enforces smoothly interpolated locally linear dynamics, and therefore expresses flexible -- yet interpretable -- dynamics akin to those of recurrent switching linear dynamical systems (rSLDS). Our approach resolves key limitations of the rSLDS such as artifactual oscillations in dynamics near discrete state boundaries, while also providing posterior uncertainty estimates of the dynamics. To fit our models, we leverage a modified learning objective which improves the estimation accuracy of kernel hyperparameters compared to previous GP-SDE fitting approaches. We apply our method to synthetic data and data recorded in two neuroscience experiments and demonstrate favorable performance in comparison to the rSLDS.

PaperPDFCodeCode Syntology ran

In Syntology Open this paper in Syntology's Atlas, the map of the papers in Syntology's graph and their citations.

For agents, Syntology's MCP tool lists every function and class Syntology harvested from this paper and whether it ran (how to connect): get_harvested_code_for_paper(arxiv_id="2408.03330")

Code

Syntology Ran 7 of 17 code samples harvested from 1 repository linked to this paper; 10 have no recorded run. Of those that ran: 3 ran · honoured contract; 3 ran · our draft was wrong; 1 ran with no contract checked.

By repository: official repository: 17 samples from 1 repository, 7 ran. The run record, sample by sample. “Ran” means executed on a synthesized input, not that the code is correct or reproduces the paper.

lindermanlab/gpslds officialmentioned in papermentioned on GitHubjax report

Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.

Code Syntology ran Syntology

17 samples harvested; 7 ran; 3 honoured the contract we drafted; 10 have no recorded run. Read from Syntology's graph 2026-09-24; that is when this build read the record, not when the samples ran.

3ran · honoured contract
3ran · our draft was wrong
1ran
10unverified

Licence: 17 of the 17 samples are pointer only, meaning Syntology does not serve that copy's text. This page shows no code text for any sample; each one links to its file in the repository.

Harvested from lindermanlab/gpslds. “Ran” means the sample executed on a synthesized input. It does not mean the output is correct, and nothing here reproduces the paper's results. “Honoured” and “violated” refer to a contract Syntology drafted from the code itself; “our draft was wrong” and “fixture could not drive it” are failures of Syntology's instrument, not of the code.

Each sample ends with its code_sha256, Syntology's identity for that exact code. An agent fetches the stored sample with Syntology's MCP tool get_code(code_sha256="…") (how to connect); click an identity to copy that call.

Repository labels, per sample. official repository: The archive marks this repository official for the paper. named in the paper: The archive records that the paper mentions this repository; it is not marked official. community (archive-listed): In the archive's code links for this paper, not marked official and not recorded as mentioned in the paper. found in paper text by Syntology: Syntology found this repository in the paper's own text; whether it is the authors' implementation is not asserted. community: Not in the archive's code links for this paper; a community repository Syntology harvested. Samples from a repository marked official are listed first. Licence labels name the repository's licence as recorded at harvest. “Pointer only” means Syntology does not serve that copy's text, for one of four reasons: no licence file was found; the licence was not identified; the licence is recorded as permissive but that copy's record is not marked cleared; or the licence is outside the permissive list Syntology serves text under (MIT, Apache-2.0, BSD and similar). Some licences outside that list permit redistribution, such as WTFPL, and GPL-3.0 under its conditions; they are simply not on the list. Hover a licence label for the reason. File links open the file on GitHub at the default branch, which may have changed since the harvest.

backward_pass lindermanlab/gpslds/gpslds/em.py official repository ran no licence file found · pointer only · 01539fa5868b696a · report
forward_pass lindermanlab/gpslds/gpslds/em.py official repository ran · honoured contract fingerprinted no licence file found · pointer only · eccd8179d36b008e · report
initialize_vem lindermanlab/gpslds/gpslds/em.py official repository ran · our draft was wrong fingerprinted no licence file found · pointer only · a7c5f160da2cfdda · report
kl lindermanlab/gpslds/gpslds/em.py official repository ran · our draft was wrong no licence file found · pointer only · f9158e7f06fbf8cd · report
kl_over_time lindermanlab/gpslds/gpslds/em.py official repository ran · our draft was wrong no licence file found · pointer only · 84e9a7c03f3d8afc · report
update_init_params lindermanlab/gpslds/gpslds/em.py official repository ran · honoured contract fingerprinted no licence file found · pointer only · 076890d8c1ad8081 · report
update_init_variational_params lindermanlab/gpslds/gpslds/em.py official repository ran · honoured contract fingerprinted no licence file found · pointer only · 516131e165058691 · report
compute_elbo lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · 5f0281db933e4f38 · report
compute_elbo_all_trials lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · 1e6fb8e7a60902a6 · report
compute_elbo_per_trial lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · ddb0c7fb04a66e53 · report
e_step lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · 862793aea2eb9de2 · report
fit_variational_em lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · 3833cbf98f80abc0 · report
make_gram lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · 1889949c23a245d1 · report
sgd lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · f7e343250fd26034 · report
update_B lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · 4668d8d91c4fa7a2 · report
update_q_u lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · b877229bff463484 · report
variational_step lindermanlab/gpslds/gpslds/em.py official repository unverified no licence file found · pointer only · c83c2cc902919ac8 · report

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

Methods

Gaussian Process

Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections