OSCR

Population coding under the scale invariance of high-dimensional noise.

Code ↔ Paper

6 matches between paragraphs of the paper and lines of its authors' code, computed by the harvester (lexical-v1). Click a colored paragraph or line to see its counterpart.

The 6 matches
  1. [1] § MATERIALS AND METHODS › Simulation of neural activity with bounded and unbounded LFI ↔ data_simulated_gamma2_manyrep/gen_simdata.py, lines 12 to 99 · score 0.85 · inverse Wishart distribution, lognormal distribution, firing rates, mouse TX42, Simulation, matrix · Report a mismatch
  2. [2] § MATERIALS AND METHODS › Bias correcting the eigenvalues and angular occupations › Eigenvalue bias correction ↔ code/code_janelia_trialeffect_new/bias_correct_morevals_allrec.ipynb, lines 12 to 106 · score 0.58 · initial guess, Bias corrected eigenvalues, iteratively, naive · Report a mismatch
  3. [3] § MATERIALS AND METHODS › Bias correcting the eigenvalues and angular occupations › Eigenvalue bias correction ↔ code/code_simdata_alphapos_gamma2/bias_correct.ipynb, lines 32 to 82 · score 0.58 · initial guess, Bias corrected eigenvalues, iteratively, naive · Report a mismatch
  4. [4] § MATERIALS AND METHODS › Bias correcting the eigenvalues and angular occupations › Angular occupation bias correction › Eigenvector bias correction. ↔ code/code_janelia_trialeffect_new/bias_correct_morevals_allrec.ipynb, lines 12 to 106 · score 0.55 · initial guess, bias corrected, algorithm, naive, cosines, squared · Report a mismatch
  5. [5] § MATERIALS AND METHODS › Bias correcting the eigenvalues and angular occupations › Angular occupation bias correction › Eigenvector bias correction. ↔ code/code_simdata_alphapos_gamma2/bias_correct.ipynb, lines 32 to 82 · score 0.55 · initial guess, bias corrected, algorithm, naive, cosines, squared · Report a mismatch
  6. [6] § MATERIALS AND METHODS › Analysis of mouse V1 data ↔ code/functions/data_extract.py, lines 26 to 75 · score 0.51 · 43–47 deg, static, orientations, 43 deg, stimuli · Report a mismatch

Paper

Loaded from Europe PMC by your browser; OSCR does not store it: doi.org · Europe PMC

The paper is loaded when this pane is shown.

The authors' code

Jupyter notebook · 109 lines · 4.3 KB · no license · 2 matches

  1. # %%
  2. import sys
  3. sys.path
  4. sys.path.append("..")
  5. sys.path.append("/mnt/disk1/home/sumedh/.local/lib/python3.8/site-packages")
  6. import numpy as np, matplotlib.pyplot as plt, time, math, sklearn
  7. import random
  8. %matplotlib inline
  9. from tqdm.notebook import tqdm, trange
  10. import matplotlib
  11. # %%
  12. from functions.pls_analysis import get_cos2_eig
  13. from functions.bias_correction_functions import bias_correct_eigs, bias_correct_cos2
  14. mouse_ids = ['TX42', 'TX41', 'TX40', 'TX39', 'TX38'] #names of the mice
  15. max_trials = [2654, 2246, 2231, 1999, 1855]
  16. nvals_plot = [150*i for i in range(1, 10)]
  17. # eigenvalue bias correction params
  18. r=100 # Number of max iterations
  19. beta=0.05 # learning rate in log scale
  20. Nsub=200 # Number of subpopulations
  21. #cos2 bias correction params
  22. M=1000 # Number of iterations for cosine bias correction
  23. eta=0.01 # learning rate in log scale
  24. q=10 # number of repeats of the algorithm with different random seed
  25. ################################ Loop over all recordings (mice) ################################
  26. biggesttic = time.perf_counter()
  27. for rec_id in [1, 2, 3, 4]:
  28. numtrials_all = [100*i for i in range(5, 26) if 100*i<max_trials[rec_id]] #multiples of 100
  29. T = numtrials_all[-1]
  30. # get values only for specific N indices to perform bias correction
  31. #load results with the largest number of trials
  32. results_i = np.load(f'./results_janelia_trialeffect/rec_{rec_id}/janelia_rec{rec_id}_{T}tr.npz')
  33. cos2, eig = get_cos2_eig(results_i) #averaged over repeats (subsets)
  34. N_vals = results_i['Neuron_subsize']
  35. norm2_delmu = (results_i['norm_delmu']**2).mean(axis=0)
  36. inds_plot = [int(np.argwhere(N_vals==nvals_plot[i])) for i in range(len(nvals_plot))]
  37. eigs_x = eig[:, inds_plot]
  38. cos2_x = cos2[:, inds_plot]
  39. norm2dmu_x = norm2_delmu[inds_plot]
  40. eigs_bc = np.empty((N_vals[inds_plot[-1]], len(inds_plot)))
  41. eigs_bc[:] = np.nan
  42. cos2_bc = np.empty((N_vals[inds_plot[-1]], len(inds_plot)))
  43. cos2_bc[:] = np.nan
  44. bigtic = time.perf_counter()
  45. for i in range(len(inds_plot)):
  46. smalltic = time.perf_counter()
  47. N_curr = N_vals[inds_plot[i]]
  48. w=np.ones(N_curr)*0.1 # Used for smoothing the naive estimates of cosine squares to avoid negative values
  49. #eigenvalue bias correction
  50. eh = eigs_x[:N_curr, i]
  51. ei=eh # initial guess equal to the naive estimate
  52. ebs , cost = bias_correct_eigs(eh,ei,T,r,beta,Nsub)
  53. idx=np.argmin(cost)
  54. eb=ebs[idx,:] # This is the bias corrected eigenvalues
  55. eigs_bc[:N_curr,i] = eb
  56. # print(f"Eig bc for N={N_curr} done in {(time.perf_counter()-smalltic):.2f}s")
  57. with open('status_bc_morevals_allrec.txt', 'a') as f:
  58. f.write(f"\nRec={rec_id}: Eig bc for N={N_curr} done in {(time.perf_counter()-smalltic):.2f}s")
  59. smalltic2 = time.perf_counter()
  60. #cos2 bias correction
  61. cos2h = cos2_x[:N_curr,i]
  62. cosbcc, cosbc , error = bias_correct_cos2(cos2h,cos2h,eb,norm2dmu_x[i],T,Nsub,M,eta,q,w)
  63. cos2_bc[:N_curr,i] = cosbcc
  64. # print(f"Cos2 bc for N={N_curr} done in {(time.perf_counter()-smalltic2):.2f}s")
  65. with open('status_bc_morevals_allrec.txt', 'a') as f:
  66. f.write(f"\nRec={rec_id}: Cos2 bc for N={N_curr} done in {(time.perf_counter()-smalltic2):.2f}s")
  67. f.write(f"\nTotal time elapsed: {(time.perf_counter()-biggesttic):.2f}s")
  68. np.savez(f'./results_janelia_trialeffect/bc_morevals_eig_cos2_rec{rec_id}.npz', eigx=eigs_x, \
  69. eigbc=eigs_bc, cos2x=cos2_x, cos2bc=cos2_bc, N_vals=N_vals)
  70. # for i in range(len(inds_plot)):
  71. # ex = eigs_x[:N_vals[inds_plot[i]],i]
  72. # cx = cos2_x[:N_vals[inds_plot[i]],i]
  73. # eb = eigs_bc[:N_vals[inds_plot[i]],i]
  74. # cb = cos2_bc[:N_vals[inds_plot[i]],i]
  75. # coefs_x = np.polyfit(np.log(ex), \
  76. # np.log(cx), deg=1)
  77. # coefs_bc = np.polyfit(np.log(eb), \
  78. # np.log(cb), deg=1)
  79. # gamma_x[i] = coefs_x[0]
  80. # gamma_bc[i] = coefs_bc[0]
  81. # alpha_nv = alpha_from_eigs2(eigs_x, np.array(nvals_plot))
  82. # alpha_bc = alpha_from_eigs2(eigs_bc, np.array(nvals_plot))
  83. # np.savez(f'./results_janelia_trialeffect/bc_morevals_eig_cos2_rec{rec_id}.npz', eigx=eigs_x, \
  84. # eigbc=eigs_bc, cos2x=cos2_x, cos2bc=cos2_bc, N_vals=N_vals, gammax=gamma_x, gammabc=gamma_bc, alphax=alpha_nv, alphabc=alpha_bc)
  85. # %%

bias_correct_morevals_allrec.ipynb at commit ec17015, no license · at the source

Overview

  1. Department of Neurobiology, David Geffen School of Medicine, University of California, Los Angeles, Los Angeles, CA 90095, USA
  2. School of Engineering and Applied Sciences, Harvard University, Cambridge, MA 02138, USA
  3. Graduate School of Informatics, Kyoto University, Kyoto 606-8501, Japan
  4. Center for Human Nature, Artificial Intelligence, and Neuroscience (CHAIN), Hokkaido University, Sapporo 060-0812, Japan
Institutions: University of California, Los Angeles (United States); Harvard University (United States); Hokkaido University (Japan); Kyoto University (Japan)
Journal: Science advances, volume 12, issue 39, article eadz9632
Dates: received 20 June 2025; accepted 19 August 2026; published online 25 September 2026; in print September 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.1126/sciadv.adz9632 · PMID 42789721 · PMCID PMC13614404 · OpenAlex W7214290456
Open access: gold, a free copy (OpenAlex)
Status: code verified (how the code is found and verified)
Categories: mouse (organism), systems (subfield)
Methods: Connectivity, Smoothing, state filtering, decompositions, Single-unit activity, calcium imaging, Spectral & time-frequency
MeSH: Models, Neurological*, Neurons*, Visual Cortex*, Animals, Mice (* major topic)
Topic: Neural dynamics and brain function (Cognitive Neuroscience, Neuroscience), according to OpenAlex
Citations: not cited yet (Europe PMC); 44 references in the paper

Abstract

High-dimensional scale-invariant neural activity is ubiquitous across brain regions and species, but its implications for information coding remain unclear. Here, we ask how stimulus information in the high-dimensional activity of mouse V1 scales with neuron number: Does it saturate due to noise correlations or increase without bound as subpopulations grow? Contrary to previous reports, we find that leading noise components that scale linearly with population size, and thus can limit information, are not sufficiently aligned with the signal to impose a bound. This conclusion follows from two scale-invariant power-law properties of neuronal responses in mouse V1: the noise eigenspectrum and alignment of noise components with the signal. We show that population subsampling links the observed power-law exponents to information boundedness and that information scaling depends on the full eigenspectrum rather than its leading modes. Last, we prove that, under subsampling, information-limiting correlations, if present, are differential correlations. Our findings clarify how information scales in high-dimensional neuronal activity under scale-invariant noise.

Reproduced under the paper's license (CC BY), from the paper cited above.

Repositories

Its files are read in the Code ↔ Paper reader above, with 6 matches between paragraphs and lines of code.

Zenodo 19504768

License: CC-BY-4.0
State: the link answers, verified on 28 September 2026
Evidence: files inventoried
Size: 1 file
Software Heritage: not checked
Found in: “Data, code, and materials availability:”
Not found: README, license file, CITATION.cff, environment file, tests, continuous integration, documentation
Tools: NumPy (36 files), Matplotlib (29 files), scikit-learn (23 files), SciPy (15 files), Numba (1 file)
Availability: 1 check, the latest on 28 September 2026: the link answers (HTTP 200)
  • 28 September 2026: the link answers (HTTP 200)
38 files

Sai-Sumedh/scaling-population-coding

License: none: the authors keep all their rights
State: the link answers, verified on 28 September 2026
Evidence: files inventoried
Commit: ec170158e146aba4c94296cf2a412d4b4dae8b2f, 24 July 2026
Languages: Jupyter (25), Python (13)
Size: 45 files, 38 scripts
Software Heritage: not archived
Found in: “Data, code, and materials availability:”
Holds: README, 25 notebooks
Not found: license file, CITATION.cff, environment file, tests, continuous integration, documentation
Tools: NumPy (37 files), Matplotlib (29 files), scikit-learn (23 files), SciPy (15 files), Numba (1 file)
Availability: 1 check, the latest on 28 September 2026: the link answers
  • 28 September 2026: the link answers
39 files

The paper's code and data availability statement is in the Data section.

Tracing map

Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.

What the map holds:

  • 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
  • 75 scripts, each with its path and the digest of its content;
  • 6 matches between paragraphs of the paper and lines of the code (method lexical-v1);
  • neither the text of the paper nor the code itself.

Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.

Data

No dataset and no data link were found in the paper.

Data, code, and materials availability

All data and code needed to evaluate and reproduce the results in the paper are present in the paper and/or the Supplementary Materials. The two-photon calcium imaging recordings analyzed in this study are publicly available from Pachitariu et al. (44), “Recordings of 20,000 neurons from V1 in response to oriented stimuli,” https://doi.org/10.25378/janelia.8279387. The code for the data analysis is archived on Zenodo (https://doi.org/10.5281/zenodo.19504768) and linked to the GitHub repository: https://github.com/Sai-Sumedh/scaling-population-coding. This study did not generate new materials.

Reproduced under the paper's license (CC BY), from the paper cited above.

Versions

The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.

Version 1, 28 September 2026: the first record

Recorded: type, language, journal, volume, issue, pages, dates, 3 authors, 5 MeSH terms, 42 references.

Cite

This paper

Moosavi, S. A., Hindupur, S. S. R., & Shimazaki, H. (2026). Population coding under the scale invariance of high-dimensional noise. Science advances, 12(39), eadz9632. https://doi.org/10.1126/sciadv.adz9632

BibTeX

@article{moosavi2026population,
author = {Moosavi, S Amin and Hindupur, Sai Sumedh R and Shimazaki, Hideaki},
title = {{Population coding under the scale invariance of high-dimensional noise}},
journal = {Science advances},
year = {2026},
month = sep,
volume = {12},
number = {39},
pages = {eadz9632},
publisher = {American Association for the Advancement of Science},
issn = {2375-2548},
doi = {10.1126/sciadv.adz9632},
url = {https://doi.org/10.1126/sciadv.adz9632},
pmid = {42789721},
pmcid = {PMC13614404}
}

RIS

TY - JOUR
AU - Moosavi, S Amin
AU - Hindupur, Sai Sumedh R
AU - Shimazaki, Hideaki
TI - Population coding under the scale invariance of high-dimensional noise
T2 - Science advances
J2 - Sci Adv
PY - 2026
DA - 2026/09/25
VL - 12
IS - 39
SP - eadz9632
SN - 2375-2548
PB - American Association for the Advancement of Science
DO - 10.1126/sciadv.adz9632
UR - https://doi.org/10.1126/sciadv.adz9632
LA - en
ER -

CSL-JSON

{
"id": "10.1126/sciadv.adz9632",
"type": "article-journal",
"title": "Population coding under the scale invariance of high-dimensional noise",
"container-title": "Science advances",
"author": [
{
"family": "Moosavi",
"given": "S Amin"
},
{
"family": "Hindupur",
"given": "Sai Sumedh R"
},
{
"family": "Shimazaki",
"given": "Hideaki"
}
],
"container-title-short": "Sci Adv",
"volume": "12",
"issue": "39",
"page": "eadz9632",
"DOI": "10.1126/sciadv.adz9632",
"PMID": "42789721",
"PMCID": "PMC13614404",
"ISSN": "2375-2548",
"publisher": "American Association for the Advancement of Science",
"URL": "https://doi.org/10.1126/sciadv.adz9632",
"language": "en",
"issued": {
"date-parts": [
[
2026,
9,
25
]
]
}
}

The tracing map gets a citation of its own once an author has validated it and it has a DOI.

Similar papers

The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.

[1] doi:10.1126/sciadv.abg5244 [code]
Internally generated population activity in cortical networks hinders information transmission.
Journal: Science advances
In common: systems, 16 references
[2] doi:10.1038/s41467-022-31254-y [code]
Noise correlations in neural ensemble activity limit the accuracy of hippocampal spatial representations.
Journal: Nature communications
In common: systems, mouse, 15 references
[3] doi:10.1038/s41467-024-54322-x [code]
Orthogonalization of spontaneous and stimulus-driven activity by hierarchical neocortical areal network in primates.
Journal: Nature communications
In common: systems, mouse, 13 references
[4] doi:10.1371/journal.pone.0271136 [code]
Targeted dimensionality reduction enables reliable estimation of neural population coding accuracy from trial-limited data.
Journal: PloS one
In common: SciPy, Matplotlib, NumPy, systems, 12 references
[5] doi:10.1371/journal.pcbi.1011667 [code]
Modeling the diverse effects of divisive normalization on noise correlations.
Journal: PLoS computational biology
In common: systems, mouse, 13 references
[6] doi:10.1038/s41467-025-62086-1 [code]
Relating natural image statistics to patterns of response covariability in macaque primary visual cortex.
Journal: Nature communications
In common: SciPy, Matplotlib, NumPy, systems, 11 references
[7] doi:10.7554/elife.97848 [code]
Visual information is broadcast among cortical areas in discrete channels.
Journal: eLife
In common: mouse, 11 references
[8] doi:10.1371/journal.pcbi.1012000 [code]
What does the mean mean? A simple test for neuroscience.
Journal: PLoS computational biology
In common: scikit-learn, SciPy, Matplotlib, 1 other tool, mouse, 9 references
[9] doi:10.7554/elife.89936 [code]
Task-specific invariant representation in auditory cortex.
Journal: eLife
In common: scikit-learn, SciPy, Matplotlib, 1 other tool, 9 references
[10] doi:10.1038/s41467-024-50563-y [code]
Flexible neural population dynamics govern the speed and stability of sensory encoding in mouse visual cortex.
Journal: Nature communications
In common: systems, mouse, 10 references

Contribute

The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.

Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.

Add the badge to the code's README

OSCR: the badge of this paper

If you are an author or a maintainer of this paper's code, you can add this badge. It says what OSCR knows of the paper (its code, its validated tracing map) and links the code to this page. Paste one of these snippets under the README's title. Keep the two marks around it, so that the badge can be replaced in place and is never added twice. OSCR changes nothing in your repository.

Markdown
<!-- oscr-badge:start -->
[![OSCR](https://openscicode.org/badge/doi/10.1126/sciadv.adz9632.svg)](https://openscicode.org/badge/doi/10.1126/sciadv.adz9632)
<!-- oscr-badge:end -->
reStructuredText
.. oscr-badge:start

.. image:: https://openscicode.org/badge/doi/10.1126/sciadv.adz9632.svg
   :target: https://openscicode.org/badge/doi/10.1126/sciadv.adz9632
   :alt: OSCR

.. oscr-badge:end
HTML
<!-- oscr-badge:start -->
<a href="https://openscicode.org/badge/doi/10.1126/sciadv.adz9632"><img src="https://openscicode.org/badge/doi/10.1126/sciadv.adz9632.svg" alt="OSCR"></a>
<!-- oscr-badge:end -->

On GitHub, open the README, choose Edit (the pencil), paste the snippet under the title, then "Commit changes".

If you choose to, you can have the badge kept up to date for you, with oscr badge --write from the command line or with the GitHub Action in your repository's workflows. Every other form and option (a card, AsciiDoc, LaTeX, shields.io) is on the badge page.

Request its removal

To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).

Discussion, reproductions, activity

Discussion: code errors and mismatches between the code and the paper about this paper and its code, from readers and its authors, beside the paper and its tracing map.

Report a mismatch between the code and the paper · Report a code error (signed in)

Reproductions: reports from readers who ran the authors' code and did not get the paper's result, what they ran, with which environment, commit and data.

Report a reproduction failure (signed in)

Versions of its code: the releases of its code tied to a version of this paper (the preprint, the accepted manuscript, the version of record…), each with the tracing map versioned with it.

Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.