Papers › Inferring stochastic low-rank recurrent neural networks from neural data

Inferring stochastic low-rank recurrent neural networks from neural data

24 Jun 2024arXiv:2406.16749archive 2025-07-28

Matthijs Pals, A Erdem Sağtekin, Felix Pei, Manuel Gloeckler, Jakob H Macke

A central aim in computational neuroscience is to relate the activity of large populations of neurons to an underlying dynamical system. Models of these neural dynamics should ideally be both interpretable and fit the observed data well. Low-rank recurrent neural networks (RNNs) exhibit such interpretability by having tractable dynamics. However, it is unclear how to best fit low-rank RNNs to data consisting of noisy observations of an underlying stochastic system. Here, we propose to fit stochastic low-rank RNNs with variational sequential Monte Carlo methods. We validate our method on several datasets consisting of both continuous and spiking neural data, where we obtain lower dimensional latent dynamics than current state of the art methods. Additionally, for low-rank models with piecewise linear nonlinearities, we show how to efficiently identify all fixed points in polynomial rather than exponential cost in the number of units, making analysis of the inferred dynamics tractable for large RNNs. Our method both elucidates the dynamical systems underlying experimental recordings and provides a generative model whose trajectories match observed variability.

PaperPDFCodeCode Syntology ran

In Syntology View this paper on Syntology: its repositories, every harvested function with whether it ran, its licence and the call to fetch it.

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="2406.16749")

Code

Syntology Ran 23 of 25 code samples harvested from 1 repository linked to this paper; 2 have no recorded run. Of those that ran: 23 ran with no contract checked.

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

mackelab/smc_rnns officialmentioned in papermentioned on GitHubpytorchApache-2.0 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

25 samples harvested; 23 ran; 0 honoured the contract we drafted; 2 have no recorded run. Read from Syntology's graph 2026-09-25; that is when this build read the record, not when the samples ran.

23ran
2unverified

Licence: 0 of the 25 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 mackelab/smc_rnns. “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.

PL_Jacobian mackelab/smc_rnns/fixed_points/stability.py official repository ran Apache-2.0 (permissive) · e791be93f5d34005 · report
Relu_derivative mackelab/smc_rnns/fixed_points/stability.py official repository ran fingerprinted Apache-2.0 (permissive) · 580464f8c0ffe2eb · report
calc_isi_stats mackelab/smc_rnns/evaluation/calc_stats.py official repository ran Apache-2.0 (permissive) · bfed00e6f5d83eaa · report
calculate_correlation mackelab/smc_rnns/evaluation/calc_stats.py official repository ran fingerprinted Apache-2.0 (permissive) · 8e47cd6071cb66e5 · report
clean_from_outliers mackelab/smc_rnns/evaluation/kl_Gauss.py official repository ran Apache-2.0 (permissive) · 58c4df2377e1b0b8 · report
construct_relu_matrix mackelab/smc_rnns/fixed_points/constrained_scify.py official repository ran Apache-2.0 (permissive) · 0dd3d7e404b4a977 · report
construct_relu_matrix_list mackelab/smc_rnns/fixed_points/constrained_scify.py official repository ran Apache-2.0 (permissive) · 37f39f6d5b10b7a6 · report
ensure_length_is_even mackelab/smc_rnns/evaluation/pse.py official repository ran fingerprinted Apache-2.0 (permissive) · 28f6dfd84d1926e3 · report
estimate_cross_correlation mackelab/smc_rnns/evaluation/calc_stats.py official repository ran Apache-2.0 (permissive) · 1f4f2991b2c87ba3 · report
fft_smoothed mackelab/smc_rnns/evaluation/pse.py official repository ran fingerprinted Apache-2.0 (permissive) · 429eb4472c867f62 · report
find_fixed_points_analytic mackelab/smc_rnns/fixed_points/find_fixed_points_analytic.py official repository ran Apache-2.0 (permissive) · 7fe34e5b23ee89f6 · report
get_average_spectrum mackelab/smc_rnns/evaluation/pse.py official repository ran fingerprinted Apache-2.0 (permissive) · cb11c1c39000382f · report
get_default_params mackelab/smc_rnns/py_rnn/default_params.py official repository ran Apache-2.0 (permissive) · 93bac9f432ea3e20 · report
h5_to_dict mackelab/smc_rnns/evaluation/eval_random_search.py official repository ran Apache-2.0 (permissive) · ce13363aeb64a0f0 · report
initialize_loadings mackelab/smc_rnns/py_rnn/initializers.py official repository ran Apache-2.0 (permissive) · 0763d929a3c9130c · report
initialize_w_inp mackelab/smc_rnns/py_rnn/initializers.py official repository ran Apache-2.0 (permissive) · 7489a305dc0b3c32 · report
initialize_w_rec mackelab/smc_rnns/py_rnn/initializers.py official repository ran Apache-2.0 (permissive) · cec5b0ec7d0fab17 · report
kl_between_two_gaussians mackelab/smc_rnns/evaluation/kl_Gauss.py official repository ran fingerprinted Apache-2.0 (permissive) · 3ebb3b68437965d3 · report
mask_none mackelab/smc_rnns/py_rnn/model.py official repository ran fingerprinted Apache-2.0 (permissive) · 229da4a3e2f0ad09 · report
mean_rate mackelab/smc_rnns/evaluation/eval_kl_pse.py official repository ran fingerprinted Apache-2.0 (permissive) · 97e84c476d2b594b · report
powerset mackelab/smc_rnns/fixed_points/constrained_scify.py official repository ran Apache-2.0 (permissive) · 07159a7ad120cadd · report
powerset mackelab/smc_rnns/fixed_points/find_fixed_points_analytic.py official repository ran Apache-2.0 (permissive) · 5803a5532384606b · report
set_nonlinearity mackelab/smc_rnns/py_rnn/model.py official repository ran Apache-2.0 (permissive) · d9dd7f86ac752725 · report
eval_likelihood_gmm_for_diagonal_cov mackelab/smc_rnns/evaluation/kl_Gauss.py official repository unverified Apache-2.0 (permissive) · 654cda3e3112fd10 · report
predict mackelab/smc_rnns/py_rnn/model.py official repository unverified Apache-2.0 (permissive) · 69ff95878cf21a72 · report

Results from the paper archive 2025-07-28

No leaderboard rows for this paper in the archive.

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