Datasets / DANDI:000774

DANDI:000774

Not yet verified

Spatial and temporal propagation of spiking activity in response to electrical stimulation in visual cortex

electrophysiological House mouse 1.6 GB 10 files
  • 11 reproduced
  • 3 partially reproduced
  • 7 not reproduced
  • 21 blocked

Tested on dataset version 0.260520.1753

View on DANDI Paper JSON How to read outcomes

Claims

From the paper Neural population dynamics of direct electrical stimulation of neocortex

Reproducing Hickman et al. (2026) from dandiset 000774 21

Notebook reproduce_hickman_2026.ipynb

Reproduced11

  • ReproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 3; Figure 4D; Figure S3Bquestion to the reader

    First-spike latencies after a pulse are bimodal, with an early peak, a pause at a median of 1.95 ms across parameters and a broader peak from about 2.5 ms, which sets the 2 ms direct window.

    the distribution of first-spike latencies, which exhibited a bimodal shape: an early peak immediately after stimulation, followed by a brief pause occurring at a median value of 1.95 ms across parameters

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 3; Figure 4D; Figure S3B

    The trough of the first-spike latency histogram (0.1 ms bins, centre in 0.5 to 5 ms) lies at 1.95, 1.75, 1.95 and 1.95 ms for 5 to 100 uA cathodal pulses, median 1.95 ms, with a higher bin before and after it in all 12 conditions. Over all 12 polarity and amplitude conditions the median is 1.45 ms (anodal and bipolar troughs 1.05 to 1.45 ms); the criterion takes the four cathodal amplitudes, as the code does. The later cathodal peaks are at 3.75 to 11.85 ms. With the first 75 trials and stable units only, the cathodal median is 1.85 ms. The verifier reads the paper's "median across parameters" as all 12 conditions (1.45 ms, outside tolerance); this is asked of a reader (questions/000774.yaml).

    How the test differed from the paper’s method

    • Latencies come from all trials of the condition's contact and all good units, as the code (figure4_figure5_figure6.ipynb cell 20), from the seven recordings that hold cathodal pulses.
  • ReproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 3; Figure S3C and S3D

    Gaussian mixtures of first-spike latencies give amplitude-specific direct windows (0.68, 1.14, 1.22 and 1.47 ms at 5, 25, 50 and 100 uA cathodal), and the proportion of direct responders does not differ from that of the 2 ms window.

    create a direct response window for each amplitude (5 μA: 0.68 ms; 25 μA: 1.14 ms; 50 μA: 1.22 ms, and 100 μA: 1.47 ms, Figure S3C ); there was no statistical difference in the overall proportion of direct responders compared to the 2-ms window used

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 3; Figure S3C and S3D

    Four-component mixtures give windows of 0.683, 1.139, 1.221 and 1.470 ms, and BIC chooses four components at every amplitude, so readings (a) and (b) agree. Over the seven recordings with cathodal pulses the responder proportion with the GMM windows does not differ from the 2 ms window at any amplitude (Wilcoxon, FDR-adjusted p 1.0, 0.667, 0.25 and 0.25) under both baseline readings. Pooled over amplitudes, the GMM windows give fewer responders (Wilcoxon p = 0.0022), a test the legend does not name.

    How the test differed from the paper’s method

    • Seven recordings, where the S3D legend implies nine; two recordings hold no cathodal pulses.
    • The crossover is taken as the first latency after the earliest component's mean where the second component's posterior reaches the first's; the Methods say only "crossover points where adjacent component responsibilities were equal" (STAR Methods, Gaussian mixture models).
  • ReproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 4; Figure 4H and 4I

    The spatial extent of direct responders stays relatively constant across amplitudes (1180 +- 559, 884 +- 652, 965 +- 451 and 522 +- 313 um), with no significant amplitude effect.

    the spatial extent of single-unit responders remained relatively constant across amplitudes (5 μ A: 1180 ± 559 μ m, 25 μ A: 884 ± 652 μ m, 50 μ A: 965 ± 451 μ m, and 100 μ A: 522 ± 313 μ m, One-way ANOVA, F = 1.265, p = 0.312)

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 4; Figure 4H and 4I

    What the paper reports

    • one-way ANOVA, F = 1.265, p = 0.312

    The per-recording extent of cathodal responders does not change with amplitude (one-way ANOVA F = 1.250, p = 0.317, 6, 6, 6 and 7 recordings). The means are 1121, 825, 813 and 452 um; those at 50 and 100 uA are outside 10% of 965 and 522 um. With the code's coordinates the means are 1179, 884, 868 and 522 um and F = 1.265, p = 0.312, the paper's F and p and its means at 5, 25 and 100 uA; 868 um at 50 uA (which the code also prints) is just outside 10% of the paper's 965 um.

    How the test differed from the paper’s method

    • The signed distance needs CCF coordinates of every channel, which the dandiset does not give; they are fitted per probe by a straight line through the units' coordinates on their channel number. The distances use the dataset's stimulation contact coordinates, which differ from the code's for jlh31 and jlh32.
  • ReproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 5; Figure 4I

    The proportion of direct responders among units within the evoked potential is small and increases with amplitude (0.014, 0.020, 0.031 and 0.030), so fewer than 5% of neurons are directly activated.

    The proportion of responders within the evoked potential increased with amplitude (5 μ A: 0.014 ± 0.12, 25 μ A: 0.020 ± 0.14, 50 μ A: 0.031 ± 0.172, and 100 μ A: 0.030 ± 0.170; linear ANOVA, F = 5.41, p = 0.020

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 5; Figure 4I

    What the paper reports

    • linear ANOVA, F = 5.41, p = 0.02 (n = 9 recordings)

    Under reading (c), the amplitude's own boundaries on the raw signed distance, the proportions of responders among units within the evoked potential are 0.0133, 0.0208, 0.0307 and 0.0301 (678 to 1427 units), all within 10% and below 0.05, rising with amplitude (OLS F = 4.99, p = 0.0256). With the code's coordinates reading (c) gives F = 5.41, p = 0.0201, the paper's values. Readings (a) and (b) give the trend (p = 0.0154 and 0.0015) with values below the paper's at 50 and 100 uA (0.0242 and 0.0254 under a). Reading (c) signs the raw distance as the code does (figure4_figure5_figure6.ipynb cell 7 makes it negative from the channel closest to the contact on); in the code's direct_df.csv distance <= 0 exactly where norm_distance <= 0 (all rows). The unsigned raw distance gives 0.0129, 0.0221, 0.0312 and 0.0293 (F = 4.55, p = 0.0330), with 0.0221 just outside 0.020 +- 0.002, so the outcome rests on the code's sign. The boundaries are the code's (outside input); the claim depends on fig2g-extent-increases, which is blocked.

    How the test differed from the paper’s method

    • The evoked-potential boundaries are the code's EP_distances.pkl (outside input).
  • ReproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 5; Figure 4J

    Among units responsive at any amplitude (99 units), direct response probability increases with amplitude (0.12, 0.14, 0.19 and 0.23), so higher amplitude both recruits more neurons and raises their spiking probability.

    among responders, direct-response probability increased with amplitude (5 μ A: 0.12 ± 0.20, 25 μ A: 0.14 ± 0.21, 50 μ A: 0.19 ± 0.25, and 100 μ A: 0.23 ± 0.31; F = 12.78, p = 0.0004, Figure 4J )

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 5; Figure 4J

    What the paper reports

    • one-way ANOVA, F = 12.78, p = 4 × 10⁻⁴ (n = 99 single units from 9 recordings)

    103 units respond (nDRP > 0.15) at some cathodal amplitude and 99 have a value at all four (from 7 recordings, 2 to 28 each). Mean nDRP is 0.116, 0.142, 0.194 and 0.233. Reading (a), the code's regression on numeric amplitude, gives F = 12.779, p = 0.00039 (paper F = 12.78, p = 0.0004); (b) a one-way ANOVA with amplitude categorical gives F = 4.44, p = 0.0044 and (c) a repeated-measures ANOVA F = 4.66, p = 0.0034. The number of good stable units per cathodal condition (2056, 2068, 2070, 2071) equals that of the code's intermediate table direct_df.csv (outside input, printed as a check). jlh31 and jlh34 ran most cathodal conditions twice at the one cathodal contact; the first run is used. With the second run, 91 units respond with values at all four amplitudes, mean nDRP 0.125, 0.092, 0.157 and 0.223, F = 12.41, p = 0.00048.

    How the test differed from the paper’s method

    • The stability filter is recomputed from the Methods (the paper's flags are not in the dataset); the time origin of the fit and the set of units the quartiles are taken over are not stated and were chosen as the first spike of the recording and all units (a standard choice). The unit counts per cathodal and anodal condition match the code's; bipolar ones differ by 1 to 5 units.
  • ReproducedResults, Stimulation amplitude and polarity tune unit-specific recruitment, paragraph 1; Figure 5A to 5C; Figure S5Aquestion to the reader

    Clustering the direct responders by response probability across amplitude gives six clusters, broadly dose-dependent (58.6% of responders) or amplitude-selective, responding almost only at one amplitude (19.2% at 5 uA, 15.2% at 25 uA, 7.1% at 100 uA).

    Hierarchical clustering produced six clusters ( Figure S5A ; cophenetic correlation coefficient (CCC) = 0.702, mean silhouette score = 0.352) that were broadly either dose-dependent or amplitude-selective. Surprisingly, while 58.6% of direct responders exhibited dose-dependent increases in response probability

    Results, Stimulation amplitude and polarity tune unit-specific recruitment, paragraph 1; Figure 5A to 5C; Figure S5A

    Under reading (a), the code's unscaled nDRP, Ward clustering of the 99 responders gives CCC = 0.702, silhouette = 0.352 and clusters of 15, 19, 40, 13, 7 and 5 units. Labelling them 25 uA, 5 uA, dose, dose, 100 uA and dose gives 15.2%, 19.2%, 58.6% and 7.1%, the paper's fractions, the assignment the criterion allows. By the claims file's half-peak rule the 7-unit cluster peaks at 50 uA (mean nDRP 0.651 against 0.343 at 100 uA) and the 13-unit cluster is the one selective for 100 uA, so no assignment that follows the clusters' peaks matches; the paper's "100 uA" cluster is better described as selective for 50 uA. Reading (b), vectors divided by their maximum, gives CCC 0.840, silhouette 0.501 and no matching assignment. Whether the criterion's "some assignment" clause is meant to pass clusters that do not respond almost only at one amplitude is asked of a reader.

  • ReproducedResults, Stimulation amplitude and polarity tune unit-specific recruitment, paragraph 2; Figure 5D to 5F; Figure S5Bverifier disagreed at first

    Clustering responders by their response to cathodal, anodal and bipolar pulses shows a polarity bias - 69.1% are polarity-specific, preferring bipolar (27.7%), anodal (23.6%) or cathodal (17.7%) pulses; smaller groups are dose-selective but polarity-agnostic (10.9%), cathodal-phase-led (12.7%) or monopolar-only (3.3%).

    revealed a striking polarity bias: 69.1% of directly activated units were polarity-specific, preferring bipolar (27.7%, Figure 5F , i ), anodal (23.6%, Figure 5F , ii ), or cathodal (17.7%, Figure 5F , iii ) stimulation

    Results, Stimulation amplitude and polarity tune unit-specific recruitment, paragraph 2; Figure 5D to 5F; Figure S5B

    Under reading (a), the code's population of units responsive to any polarity (220 units), CCC = 0.773 and silhouette = 0.538 are within 10%. Five assignments of the six clusters to the three polarities or none give all three fractions within 10%, and in two of them each labelled cluster's highest mean is at its label; the first of these gives bipolar 0.287, anodal 0.232 and cathodal 0.186 (polarity-specific 0.705), the criterion's "some assignment". It labels as bipolar both a cluster of 45 units (mean 1.0, 0.04, 0.03 for bipolar, cathodal, anodal) and one of 18 units (0.98, 0.40, 0.79), whose anodal response is 0.79 of its bipolar one. By the half-maximum rule only the first prefers bipolar, giving 0.205, 0.232 and 0.186 (total 0.623). With the bipolar contact nearest layer 5 units (235 units), the half-maximum rule gives 0.260, 0.213 and 0.196. Among the 99 cathodal responders (reading b) CCC 0.711 and silhouette 0.393 are outside tolerance. Whether the "some assignment" clause is meant this loosely is asked of a reader.

    How the test differed from the paper’s method

    • The contact nearest layer 5 is judged from the region of the nearest unit, as the dataset gives no region for the contacts.
  • ReproducedResults, Cell type correlates with direct activation distance distributions, paragraph 2; Figure 5I and 5Jrests on an assumption

    FS direct responders lie closer to the stimulation site than RS responders, and higher amplitude brings both closer.

    FS responders were clustered significantly closer to the stimulation site than RS responders (two-way ANOVA, F waveform = 8.04, p = 0.005)

    Results, Cell type correlates with direct activation distance distributions, paragraph 2; Figure 5I and 5J

    What the paper reports

    • waveform class: two-way ANOVA, F = 8.04, p = 0.005
    • amplitude: two-way ANOVA, F = 4.10, p = 0.008

    With the code's classes and the dataset's coordinates FS responders lie closer than RS (mean 111 against 267 um; waveform F = 7.02, p = 0.0090) and amplitude has an effect (F = 3.95, p = 0.0097); FS 205 um at 5 uA and RS 556 and 183 um at 5 and 100 uA are within 10%. With the code's coordinates F = 8.04, p = 0.0052 and F = 4.10, p = 0.0080, the paper's values. Under the Methods' k-means classes the waveform effect is not significant (F = 2.96, p = 0.0874). Over 20 k-means runs (that run and 19 from other seeds and initialisations; 296 to 470 FS units) the waveform effect is never significant (p 0.087 to 0.401); the run of lowest within-cluster sum of squares gives p = 0.127.

    How the test differed from the paper’s method

    • Assumed the code's FS and RS classes (a 0.4 ms duration threshold), where the Methods describe k-means. The code is what the paper ran, its classes give the paper's statistics exactly here, they are deterministic, and the k-means classes change with the initialisation (296 to 470 FS units over 20 runs); the k-means reading is shown. Asked of the authors (see findings).
    • Signed distances as in fig4h-extent-constant.
  • ReproducedResults, Cell type correlates with direct activation distance distributions, paragraph 3; Figure S6Brests on an assumption

    In the 10 ms after a pulse, the firing of FS and RS cortical units differs at every amplitude.

    we observed pronounced differences in spiking between waveform classes across amplitudes that evolved with stimulation strength (Kolmogorov-Smirnov tests at each amplitude

    Results, Cell type correlates with direct activation distance distributions, paragraph 3; Figure S6B

    What the paper reports

    • 5 uA: Kolmogorov-Smirnov, adjusted, p = 0.01
    • 25 uA: Kolmogorov-Smirnov, adjusted, p = 0.007
    • 50 uA: Kolmogorov-Smirnov, adjusted, p = 0.019
    • 100 uA: Kolmogorov-Smirnov, adjusted, p = 0.01 (n = 9 recordings)

    With the code's classes and 1 ms bins (the code's), FS and RS early firing differ at every amplitude (Holm-Sidak adjusted p 0.0056, 0.0038, 0.0127, 0.0127; 70 bins per class). With the legend's 2 ms bins no adjusted p is below 0.05 (0.0598 to 0.1151). Under the Methods' k-means classes only 100 uA differs (1 ms bins, adjusted p 0.0139). Over 20 k-means runs (that run and 19 from other seeds and initialisations) all four amplitudes differ in 9 and one in 11.

    How the test differed from the paper’s method

    • Assumed the code's FS and RS classes (see fig5j-fs-closer).
    • Seven recordings with cathodal pulses, where the legend gives nine.
  • ReproducedResults, Single pulse electrical stimulation leads to extensive circuit response, paragraph 3; Figure 6Frests on an assumptionquestion to the reader

    Response timing shows no propagating wave - time to peak inhibition and to peak rebound excitation do not correlate with distance from the stimulation site.

    This was confirmed statistically by the absence of a relationship between distance from stimulation and time to peak inhibition ( Figure 6F , ii ; Pearson r = 0.09) or time to peak rebound excitation ( Figure 6F , iii ; Pearson r = −0.03)

    Results, Single pulse electrical stimulation leads to extensive circuit response, paragraph 3; Figure 6F

    What the paper reports

    • time to peak inhibitory response, Figure 6F legend: Pearson correlation, r = -0.03, p = 0.51
    • time to peak excitatory response, Figure 6F legend: Pearson correlation, r = 0.09, p = 0.64

    With the code's peak rule, among 614 cortical circuit responders at 50 uA cathodal, time to peak excitation (median 205 ms) and to peak inhibition (median 70 ms) do not correlate with signed distance (Pearson r = -0.001, p = 0.987; r = 0.047, p = 0.244); with the cathodal contact for jlh32 too, r = 0.004 and 0.047 (p = 0.918 and 0.240); dropping the 19 units without a prominent peak, r = 0.002 and 0.050 (p = 0.969 and 0.219). With the criterion's 4 ms minimum width applied, both correlate (r = 0.131, p = 0.001 and r = -0.115, p = 0.004; with units without such a peak dropped, r = 0.153, p < 0.001 and r = -0.121, p = 0.004). jlh39 and jlh49v1 hold no cathodal pulses, so the readings with and without them are the same.

    How the test differed from the paper’s method

    • Assumed the code's peak rule. The criterion lists a minimum width of 4 ms and a centre of mass over +-6 ms, the values in the code's configuration (figure7.ipynb cell 11); the code converts them to 5 ms bins with round(), so the width (0.8 bins, rounded to 1) is not applied, as the code applies it only above 1 bin, and the centre of mass spans +-1 bin. The paper's r values come from this rule and the Methods give none. Applying the width filter changes the outcome; asked of a reader and of the authors.
  • ReproducedResults, Single pulse electrical stimulation leads to extensive circuit response, paragraph 5; Figure 6G

    After a pulse the population's dimensionality drops (participation ratio and effective dimensionality fall), even as firing rate fluctuates, with a brief rise in effective dimensionality within 20 ms at 5, 25 and 50 uA.

    the response window was associated with a pronounced reduction in neural dimensionality, as reflected by decreases in both PR and effective dimensionality

    Results, Single pulse electrical stimulation leads to extensive circuit response, paragraph 5; Figure 6G

    With 2 ms bins the normalised participation ratio in 12 to 150 ms is 0.982, 0.951, 0.765 and 0.459 and the effective dimensionality 0.986, 0.961, 0.778 and 0.476 at 5 to 100 uA, all below 1. Within 0 to 20 ms the largest bin of effective dimensionality is 1.024, 1.087 and 1.094 at 5, 25 and 50 uA (reading b; 1.043 at 100 uA as well), while its mean over the window is 0.929, 0.994 and 0.960 (reading a). With 5 ms bins the largest bin exceeds 1 only at 25 and 50 uA. The code's "layer" unit filter gives the same pattern. Against the pre-stimulus spread (-200 to -4 ms; largest normalised effective dimensionality 1.079, 1.056, 1.091 and 1.054 at 5 to 100 uA), the largest bin in 0 to 20 ms exceeds the pre-stimulus maximum at 25 uA (5 bins) and 50 uA (3 bins) but not at 5 uA (1.024) or 100 uA; the criterion does not compare with it.

    How the test differed from the paper’s method

    • Seven recordings with cathodal pulses, where the legend gives nine.

Partially reproduced3

  • Partially reproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 4; Figure 4H

    Higher amplitudes bring direct responders closer to the stimulation site (mean distance 426, 300, 137 and 146 um at 5, 25, 50 and 100 uA), driven by 5 against 50 and 5 against 100 uA.

    Stimulation amplitude significantly decreased the average distance from stimulation of responsive units ... one-way ANOVA, F = 4.15, p = 0.007), driven specifically by differences between 5 μ A and 50 μ A and between 5 μ A and 100 μ A (post hoc Tukey tests: 5 vs. 50 μ A, p = 0.015; 5 vs. 100 μ A, p = 0.019

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 4; Figure 4H

    What the paper reports

    • one-way ANOVA, F = 4.15, p = 0.007
    • 5 against 50 uA: Tukey HSD, p = 0.015
    • 5 against 100 uA: Tukey HSD, p = 0.019

    Mean signed distances of cathodal responders are 414, 264, 129 and 142 um (29, 34, 45 and 45 responders). The tests hold (ANOVA F = 3.97, p = 0.0094; Tukey 5 against 50 uA p = 0.0134, 5 against 100 uA p = 0.0199, other pairs above 0.42). The 25 uA mean is outside 10% of 300 um. With the code's coordinates the means are 426, 300, 137 and 146 um and F = 4.15, p = 0.0074, the paper's values; the difference is the dataset's contact coordinates for jlh31 and jlh32.

    How the test differed from the paper’s method

    • Signed distances as in fig4h-extent-constant (fitted channel coordinates, the dataset's contact coordinates).
  • Partially reproducedResults, Single pulse electrical stimulation leads to extensive circuit response, paragraph 2; Figure 6E; Figure S7C

    The proportion of cortical units in the circuit response grows with amplitude (0.14, 0.21, 0.52, 0.79), so more than half are recruited at the higher amplitudes; subcortical units respond less but also increase (0.08, 0.09, 0.17, 0.35).

    Further, the proportion of cortical neurons participating in this circuit response increased with amplitude (5 μ A: 0.14 ± 0.05, 25 μ A: 0.21 ± 0.08, 50 μ A: 0.52 ± 0.26, and 100 μ A: 0.79 ± 0.21; F amplitude = 175.18, p < 0.001), such that at the highest amplitudes, over half of all recorded cortical units were recruited

    Results, Single pulse electrical stimulation leads to extensive circuit response, paragraph 2; Figure 6E; Figure S7C

    What the paper reports

    • amplitude: multivariate ANOVA, F = 175.18, p < 0.001
    • region: multivariate ANOVA, F = 81.50, p < 0.001
    • polarity: multivariate ANOVA, F = 5.19, p = 0.006 (n = 9 recordings)

    With zetapy 2.7.2 getZeta (the function the code calls) the proportion of cortical circuit responders is 0.125, 0.209, 0.520 and 0.785 and of subcortical ones 0.089, 0.094, 0.164 and 0.360 (7 recordings), above one half for cortex at 50 and 100 uA. The cathodal model gives amplitude F = 89.48 (p = 6.84e-13), region F = 36.59 (p = 1.61e-07) and interaction F = 15.26 (p = 2.72e-04); the Results print F region 40.14 and interaction 15.61 next to F amplitude 175.18, without naming the model. Two values are just outside 10% (cortical 0.125 against 0.14, subcortical 0.089 against 0.08 at 5 uA). The recomputed p-values agree with the paper's cached ones (Pearson r = 0.992, same decision for 0.970 of 7980 tests), and the cached ones give 0.14, 0.21, 0.52, 0.79 and 0.08, 0.09, 0.17, 0.35 on the dataset's units, the paper's values, and amplitude F = 172.99 and region F = 85.99 in the code's full model (paper 175.18 and 81.50). The difference is the resampling of the ZETA test.

    How the test differed from the paper’s method

    • Only the cathodal conditions were tested with ZETA (7980 tests, 9 minutes); the code's full model over three polarities would take about 20 minutes more and is shown with the cached p-values.
    • zetapy 2.7.2 imports scipy.signal.gaussian, which current scipy keeps under scipy.signal.windows; the helper adds it before importing zetapy.
    • As the code (get_zetas), the trials are all those of the run that holds the condition's first trial, with no other filter; each run holds only the condition's single 100 us pulses, except in jlh53concat, where a run number covers 150 trials of monopolar and bipolar pulses. With jlh53concat's own 75 cathodal trials the proportions are 0.121, 0.203, 0.514 and 0.756 (cortical) and 0.095 and 0.339 (subcortical at 5 and 100 uA).
  • Partially reproducedResults, Single pulse electrical stimulation leads to extensive circuit response, paragraph 4; Figure S9B

    Direct response probability is slightly but significantly higher on high than on low cortical-state trials at every cathodal amplitude.

    DRP was slightly higher during high cortical states across stimulation amplitudes

    Results, Single pulse electrical stimulation leads to extensive circuit response, paragraph 4; Figure S9B

    What the paper reports

    • 5 uA: paired Wilcoxon signed-rank, FDR corrected, p < 0.001, FDR-corrected
    • 25 uA: paired Wilcoxon signed-rank, FDR corrected, p < 0.001, FDR-corrected
    • 50 uA: paired Wilcoxon signed-rank, FDR corrected, p = 0.01, FDR-corrected
    • 100 uA: paired Wilcoxon signed-rank, FDR corrected, p = 0.03, FDR-corrected

    Under readings (a, i, x), all good stable units, the whole population and DRP, DRP is higher on high-state trials at every amplitude (0.0124 against 0.0120 at 5 uA, 0.0143 against 0.0126, 0.0159 against 0.0150, 0.0165 against 0.0159 at 100 uA; FDR p 0.0016 or less); the cortical population gives the same (FDR p up to 0.0156). The values are about a quarter of the paper's at 5 uA (0.051 and 0.056) and a sixth at 100 uA (0.104 and 0.107), and no reading reaches them; among the 99 direct responders the values are 0.12 to 0.25 and no difference is significant (FDR p 0.3008 to 0.9463). With nDRP over all units the low state is higher at every amplitude (FDR p 0.0336 or less).

Not reproduced7

  • Not reproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 4; Figure 4Fverifier disagreed at firstrests on an assumptionquestion to the reader

    Direct responders are sparse and mostly lie within the average boundaries of the evoked potential (33.3%, 74.1%, 88.2% and 100% at 5, 25, 50 and 100 uA), with more outside at 5 uA.

    Direct responders were sparse and generally constrained within the average boundaries of the evoked potential (33.3% for 5 μ A, 74.1% for 25 μ A, 88.2% for 50 μ A, and 100% for 100 μ A

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 4; Figure 4F

    With the code's mean evoked-potential boundaries (outside input; -424 to 510, -725 to 605, -812 to 747 and -750 to 888 um at 5 to 100 uA) and the dataset's coordinates, the fractions of RS responders within (reading b, the code's classes) are 0.400, 0.778, 0.882 and 1.000; three are within 10% and 5 uA is the lowest, but 0.400 is outside 0.333 +- 0.033. Over all responders (reading a) they are 0.552, 0.794, 0.889 and 1.000. The 5 uA fraction is outside tolerance under every reading run, and the criterion has no partial clause, so the claim is not reproduced. With the code's coordinates reading (b) gives 0.333, 0.741, 0.882 and 1.000 and reading (a) 0.483, 0.765, 0.889 and 1.000, the paper's and the code's values; the 5 uA difference comes from the contact coordinates of jlh31 and jlh32. Under the Methods' k-means classes the RS fractions are 0.520, 0.774, 0.900 and 1.000. The claim depends on fig2g-extent-increases, which is blocked; it is tested on the boundaries the code computed, which its criterion names as an outside input. Whether such claims should be marked as using outside inputs is asked of a reader.

    How the test differed from the paper’s method

    • The evoked-potential boundaries are the code's EP_distances.pkl (outside input), not computed from the raw AP band, which the dandiset does not hold.
    • Assumed the code's FS and RS classes for reading (b) (see fig5-fsrs-proportion).
    • Signed distances as in fig4h-extent-constant (fitted channel coordinates, the dataset's contact coordinates).
  • Not reproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 4; Figure 4H and 4I

    5 and 100 uA recruit distinct spatial populations of direct responders; the other pairs the paper reports (5-25, 25-50, 25-100 and 50-100 uA) do not differ after correction.

    Further, 5 μ A and 100 μ A recruited distinct spatial populations (pairwise two-sample Kolmogorov-Smirnov tests, corrected for six comparisons using the Holm-Šidák procedure

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 4; Figure 4H and 4I

    What the paper reports

    • 5 against 100 uA: two-sample Kolmogorov-Smirnov, Holm-Sidak corrected, p = 0.016, Holm-Sidak, six comparisons-corrected
    • 50 against 100 uA: two-sample Kolmogorov-Smirnov, Holm-Sidak corrected, p = 0.083, Holm-Sidak, six comparisons-corrected
    • 5 against 25 uA: two-sample Kolmogorov-Smirnov, Holm-Sidak corrected, p = 0.566, Holm-Sidak, six comparisons-corrected
    • 25 against 50 uA: two-sample Kolmogorov-Smirnov, Holm-Sidak corrected, p = 0.566, Holm-Sidak, six comparisons-corrected
    • 25 against 100 uA: two-sample Kolmogorov-Smirnov, Holm-Sidak corrected, p = 0.363, Holm-Sidak, six comparisons-corrected

    With the dataset's coordinates, 5 against 100 uA has KS p = 0.0193 and Holm-Sidak adjusted p = 0.1105, so the two populations do not differ after correction; 5-25, 25-50, 25-100 and 50-100 do not differ (adjusted p 0.5255, 0.8290, 0.7439, 0.8290), and 5-50 has adjusted p 0.1309 (not judged). With the code's coordinates the adjusted p are 0.5662, 0.0164, 0.5662, 0.3625 and 0.8245 for 5-25, 5-100, 25-50, 25-100 and 50-100, which reproduce the claim (for 50-100 the paper prints 0.083, against 0.8245 here). The outcome turns on the stimulation contact coordinates of jlh31 and jlh32, which in the dataset differ from those the paper used.

    How the test differed from the paper’s method

    • Signed distances as in fig4h-extent-constant.
  • Not reproducedResults, Evoked spiking is sparse and increases in density with amplitude, paragraph 6; Figure S3Everifier disagreed at first

    Broader windows (5 and 10 ms) give more responsive units than the 2 ms window (0.092 and 0.175 at 100 uA), and their amplitude dependence is stronger.

    As expected, broader windows increased the proportion of responsive units (0.092 ± 0.01 for the 5 ms window and 0.175 ± 0.01 for the 10 ms window at 100 μ A, FDR-corrected p < 0.001, paired Wilcoxon signed-rank tests

    Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 6; Figure S3E

    What the paper reports

    • paired Wilcoxon signed-rank tests, p < 0.001, FDR-corrected (n = 9 recordings)

    Broader windows give more responders and a stronger amplitude dependence. Under reading (a) the mean proportions at 100 uA are 0.021 (2 ms), 0.058 (5 ms) and 0.106 (10 ms), and both broader windows exceed 2 ms (FDR p = 0.0208; 0.0417 under reading b); from 5 to 100 uA the proportion rises by 0.005, 0.041 and 0.073. The 100 uA proportions are below the paper's 0.092 and 0.175 under both readings (b gives 0.057 and 0.091), and the criterion has no partial clause, so the claim is not reproduced. With seven recordings the smallest exact Wilcoxon p is 0.0156.

    How the test differed from the paper’s method

    • Seven recordings in place of nine.
  • Not reproducedResults, Cell type correlates with direct activation distance distributions, paragraph 1 (the text cites Figure 6I; the panels are Figure 5I and 5J)verifier disagreed at first

    Direct responders stay sparse for both FS and RS units across amplitudes, with a main effect of amplitude but no effect of waveform class or interaction.

    The proportion of direct responders remained sparse across amplitudes for both FS ... and RS neurons ... with a significant main effect of amplitude but no effect of waveform class or interaction

    Results, Cell type correlates with direct activation distance distributions, paragraph 1 (the text cites Figure 6I; the panels are Figure 5I and 5J)

    What the paper reports

    • amplitude: two-way ANOVA, F = 3.04, p = 0.028
    • waveform class: two-way ANOVA, F = 1.33, p = 0.262
    • interaction: two-way ANOVA, F = 1.34, p = 0.21

    No reading gives the four proportions within 10%. Over all FS and RS units (reading a), under the code's classes (616 FS and 1352 RS units), FS 0.0196 and 0.0227 and RS 0.0111 and 0.0221 at 5 and 100 uA. Within the evoked potential (reading b, the code's boundaries as outside input) the closest is reading (c) of fig4i, FS 0.0275 and 0.0324 and RS 0.0069 and 0.0299. Responders are sparse in both classes, and class and interaction have no effect under the code's classes (class p 0.437 over all units; 0.983, 0.608 and 0.821 within, under readings a, b and c of fig4i). The amplitude term, which the paper's F = 3.04 with p = 0.028 implies is categorical, has p = 0.123 over all units and 0.039, 0.005 and 0.051 within under the three readings; with amplitude numeric p is 0.037 or less. Under the Methods' k-means (one run) the class effect is significant (F = 4.85, p = 0.028), and in 4 of 20 k-means runs (that run and 19 from other seeds and initialisations).

    How the test differed from the paper’s method

    • FS and RS classes shown under the code's 0.4 ms duration threshold and the Methods' k-means; the outcome is the same under both.
    • Reading (b) uses the code's evoked-potential boundaries (outside input).
  • Not reproducedResults, Cell type correlates with direct activation distance distributions, paragraph 2; Figure 5I and 5Jverifier disagreed at firstrests on an assumption

    FS responders lie almost entirely within the evoked potential even at low amplitude (75%, 83.3%, 100%, 100%), while only 33% of RS responders do at 5 uA, rising with amplitude.

    FS responders were almost entirely constrained within the evoked-potential boundary, even at low amplitudes (75.0% for 5 μ A; 83.3% for 25 μ A; 100% for 50 μ A; and 100% for 100 μ A). In contrast, only 33.0% of RS responders were within the evoked potential at 5 μ A

    Results, Cell type correlates with direct activation distance distributions, paragraph 2; Figure 5I and 5J

    Under the code's threshold classes and the dataset's coordinates, FS responders lie within the evoked potential at 0.750, 0.833, 1.000 and 1.000 and RS responders at 0.400, 0.778, 0.882 and 1.000 at 5 to 100 uA. FS at 5 and 25 uA and RS at 25 uA are within 10% and FS is at or above RS at every amplitude; RS at 5 uA (0.400) is outside 0.33 +- 0.033. With the code's coordinates RS gives 0.333 and 0.741, and all four match. Under the Methods' k-means FS responders are 2, 2, 4 and 7 units, all within (1.000), outside tolerance. No class choice gives the four values, and the criterion has no partial clause, so the claim is not reproduced. The boundaries are the code's (outside input); the claim depends on fig4f-responders-within-ep.

    How the test differed from the paper’s method

    • The evoked-potential boundaries are the code's EP_distances.pkl (outside input).
    • Assumed the code's FS and RS classes (see fig5j-fs-closer).
    • Signed distances as in fig4h-extent-constant.
  • Not reproducedResults, Single pulse electrical stimulation leads to extensive circuit response, paragraph 1; Figure 6Cverifier disagreed at firstquestion to the reader

    Spiking modulation after a single pulse increases with amplitude for cortical units (12.60, 13.37, 15.11, 20.21) but stays about constant for subcortical units (12.27, 12.60, 12.08, 13.48), differing between them at 50 and 100 uA.

    The magnitude of spiking modulation ( STAR Methods ) increased with amplitude for cortical units ... but remained relatively constant for subcortical units

    Results, Single pulse electrical stimulation leads to extensive circuit response, paragraph 1; Figure 6C

    What the paper reports

    • region: two-way ANOVA, F = 106.28, p < 0.001
    • amplitude: two-way ANOVA, F = 487.36, p < 0.001
    • amplitude x region: two-way ANOVA, F = 90.88, p < 0.001
    • cortical against subcortical, 50 uA and 100 uA: Student's t test, p < 0.001

    Under reading (a), the code's (0 to 500 ms, 10 ms bins), the magnitudes are 12.60, 13.37, 15.11, 20.21 (cortical, 1287 units) and 12.27, 12.60, 12.08, 13.48 (subcortical, 351 units), the paper's values, with region F = 106.27, amplitude F = 487.36 and interaction F = 90.88 (all p < 1e-20), and cortical against subcortical Student p = 4.65e-09 at 50 uA and 5.50e-29 at 100 uA. Under the legend's 5 ms bins (b) the three terms hold (p 4.76e-05 or less) but the 50 uA comparison has p = 1.10e-02. Under the Methods' window (c, 500 to 1000 ms) there is no region effect (p = 0.457), the amplitude term has p = 1.06e-03, and no amplitude comparison has p < 0.001. The criterion asks for every reading and has no partial clause, so the claim is not reproduced.

    How the test differed from the paper’s method

    • Seven recordings with cathodal pulses, where the legend gives nine.
  • Not reproducedResults, Single pulse electrical stimulation leads to extensive circuit response, paragraph 1; Figure 6D; Figure S8A

    At 100 uA higher visual areas, beyond the reach of direct activation, are more modulated than subcortical units (15.82 against 13.48) although both lie at similar distances from the stimulation site.

    Higher visual areas (HVAs) exhibited responses at higher amplitudes (magnitude of modulation = 15.82 ± 9.97) that were significantly different from those of subcortical units (magnitude of modulation = 13.48 ± 7.70) at 100 μ A (Tukey HSD, p = 0.004

    Results, Single pulse electrical stimulation leads to extensive circuit response, paragraph 1; Figure 6D; Figure S8A

    What the paper reports

    • HVA against subcortical modulation, 100 uA: Tukey HSD, p = 0.004
    • distance from stimulation, HVA against subcortical: Mann-Whitney, p = 0.511

    With the legend's groups (HVA LM, AM and PM, 120 units; subcortical CA1, subiculum, dentate gyrus and superior colliculus, 460 units in the seven recordings with cathodal pulses, 273 of them CA1 units that the code's list leaves out) HVA units have a magnitude of 13.91 at 100 uA and subcortical units 13.26, Tukey p = 0.784; the distances do not differ (Mann-Whitney p = 0.54), as in the paper. With the code's subcortical list (351 units) the values are 13.91 and 13.48, Tukey p = 0.902. With the code's broader HVA group (any layered name outside V1, 349 units, of which 91 are superior colliculus or dentate gyrus units) the values are 15.82 and 13.68 and Tukey p = 0.0202, but the distances then differ (p = 5.92e-08). The paper's 15.82 is the code's group; the legend's groups show no difference.

    How the test differed from the paper’s method

    • In the code's variant each unit is put in one group (V1, then HVA, then subcortical); the code lets superior colliculus units fall in both its HVA and subcortical groups.

Not tested in a notebook 21

21 data missing from the dataset

Blocked21

  • BlockedResults, Orthogonal multi-Neuropixels recordings, paragraph 3; Figure 1L· Data missing from the dataset

    Recording amplifier saturation after a 0.2 ms pulse lasts about 0.14, 0.24, 0.26 and 0.30 ms at 5, 25, 50 and 100 uA, increasing with amplitude.

    saturation persisted for approximately 0.14 ± 0.00 ms at 5 μ A, 0.24 ± 0.01 ms at 25 μ A, 0.26 ± 0.01 ms at 50 μ A, and 0.30 ± 0.01 ms at 100 μ A ( Figure 1L , one-way ANOVA, F = 105.62, p < 0.001)

    Results, Orthogonal multi-Neuropixels recordings, paragraph 3; Figure 1L

    What the paper reports

    • one-way ANOVA, F = 105.62, p < 0.001 (n = 25 probes from 9 recordings)

    Saturation is measured on the raw AP band, which the dandiset does not hold (its NWB files have units and trials only, no acquisition). The paper says the raw AP, LFP and NI-DAQ binaries are available on request; no public copy was found on DANDI (searches for Denman and Hickman return only this dandiset) or in the code repository.

  • BlockedResults, Orthogonal multi-Neuropixels recordings, paragraph 3; Figure S1B· Data missing from the dataset

    Saturation duration increases with amplitude, and polarity does not significantly change it.

    despite increases in saturation with amplitude (two-way ANOVA, F = 236.04, p < 0.001), and polarity alone did not significantly change saturation duration ( F = 2.48, p = 0.09)

    Results, Orthogonal multi-Neuropixels recordings, paragraph 3; Figure S1B

    What the paper reports

    • amplitude: two-way ANOVA, F = 236.04, p < 0.001
    • polarity: two-way ANOVA, F = 2.48, p = 0.086 (n = 9 recordings)

    Needs the raw AP band (see fig1l-saturation-amplitude), and depends on that claim.

  • BlockedResults, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 3; Figure 2G· Data missing from the dataset

    The spatial extent of the evoked potential, both superficial and deep to the stimulation source, increases with amplitude (superficial about -424, -683, -742 and -744 um; deep about 510, 591, 662 and 712 um at 5, 25, 50 and 100 uA).

    Both the superficial and deep distance from stimulation increased with amplitude ( Figure 2G , superficial ... p = 0.002, linear mixed-effects model

    Results, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 3; Figure 2G

    What the paper reports

    • amplitude, Results and Figure S2E legend: linear mixed-effects model, random intercept for recording, p = 0.002 (n = 25 probes from 9 recordings)
    • amplitude, Figure 2G legend: linear mixed-effects model, random intercept for recording, p < 0.001 (n = 25 probes from 9 recordings)

    The evoked-potential extents are computed from the raw AP band, which is not in the dandiset (see fig1l-saturation-amplitude). The code repository holds the extents it computed (EP_boundaries_w_array_distances_curated.pkl, EP_distances.pkl), but these are results derived from the missing data, not data, so the claim was not tested on them.

  • BlockedResults, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 3; Figure 2G; Figure S2F· Data missing from the dataset

    The increase in spatial extent with amplitude is sublinear and poorly fit by a linear regression (superficial R2 = 0.02, deep R2 = 0.17), unlike the linear growth of the point-source model.

    the increase in spatial extent of the biologically observed evoked potential was sublinear and poorly fit by linear regression (superficial R 2 = 0.02 and deep R 2 = 0.17)

    Results, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 3; Figure 2G; Figure S2F

    Needs the raw AP band; depends on fig2g-extent-increases.

  • BlockedResults, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 3; Figure 2I; Figure S2E· Data missing from the dataset

    Spatial extents do not differ across polarities.

    Notably, spatial extents did not differ across polarities ( Figures 2J and S2E , linear mixed-effects model, p = 0.151

    Results, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 3; Figure 2I; Figure S2E

    What the paper reports

    • polarity, Results: linear mixed-effects model, p = 0.151
    • polarity, 50 uA, Figure 2I and S2E legends: linear mixed-effects model, random intercept for recording, p = 0.495 (n = 25 probes from 9 recordings)

    Needs the raw AP band; depends on fig2g-extent-increases.

  • BlockedResults, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 4; Figure 2H· Data missing from the dataset

    Probe trajectory affects spatial extent, with larger extents on probe C, which stays in cortex, than on probes A and B, which cross the white matter into subcortical areas.

    we observed probe trajectory had a significant effect on average spatial extent, most notably with an increase in spatial extent for probe C ( Figure 2H , probe A vs. probe C: p < 0.001, probe B vs. probe C: p < 0.001

    Results, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 4; Figure 2H

    What the paper reports

    • probe A against probe C: linear mixed-effects model, p < 0.001
    • probe B against probe C: linear mixed-effects model, p < 0.001 (n = probe A 7, probe B 9, probe C 9 probes from 9 recordings)

    Needs the raw AP band; depends on fig2g-extent-increases.

  • BlockedResults, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 2; Figure S2G and S2H· Data missing from the dataset

    For all polarities, the maximum evoked potential and its area under the curve at 1 ms increase nearly linearly with amplitude.

    For all polarities, both the maximum evoked potential and AUC increased nearly linearly ( Figures S2G and S2H

    Results, Evoked potential spatial extent increases sub-linearly with amplitude, paragraph 2; Figure S2G and S2H

    Needs the raw AP band; depends on fig2g-extent-increases.

  • BlockedResults, Anisotropic volumetric profile of the induced electrical field, paragraph 1; Figure 3C· Data missing from the dataset

    The evoked-potential volume, the convex hull of the six extents on the three probes, tends to increase sublinearly with amplitude (cathodal 0.31, 0.43 and 0.61 mm3 at 25, 50 and 100 uA), with no significant main effect of amplitude or polarity.

    trended to increase sublinearly with stimulation amplitude ( Figure 3C ; cathodal polarity ... 0.31 ± 0.20 mm 3 , −50 μ A: 0.43 ± 0.19 mm 3 , −100 μ A: 0.61 ± 0.22 mm 3 ; linear R 2 = 0.259

    Results, Anisotropic volumetric profile of the induced electrical field, paragraph 1; Figure 3C

    What the paper reports

    • amplitude: two-way ANOVA, F = 1.96, p = 0.16
    • polarity: two-way ANOVA, F = 0.42, p = 0.66 (n = 4 recordings)

    The volumes are built from the six evoked-potential extents, which need the raw AP band; depends on fig2g-extent-increases.

  • BlockedResults, Anisotropic volumetric profile of the induced electrical field, paragraph 2; Figure 3D and 3E· Data missing from the dataset

    Evoked potentials extend broadly within cortex but are compressed ventrally at the corpus callosum; along probe C, which stays in cortex, they extend beyond a fitted isotropic sphere, and along probe B, which crosses into white matter, they are truncated relative to it.

    when probes remained in cortex (probe C, Figure 3E ), responses extended beyond the symmetric fit, whereas when probes traversed into subcortical white matter (probe B, Figure 3E ), responses were truncated relative to the sphere

    Results, Anisotropic volumetric profile of the induced electrical field, paragraph 2; Figure 3D and 3E

    Needs the raw AP band; depends on fig3c-volume.

  • BlockedResults, Anisotropic volumetric profile of the induced electrical field, paragraph 3; Figure 3F and 3G· Data missing from the dataset

    Anisotropy strength averages 0.493 +- 0.288 across amplitudes and polarities and does not vary with amplitude or polarity.

    Across all amplitudes and polarities, anisotropy strength averaged 0.493 ± 0.288 and did not vary systematically with stimulation amplitude or polarity

    Results, Anisotropic volumetric profile of the induced electrical field, paragraph 3; Figure 3F and 3G

    What the paper reports

    • polarity: two-way ANOVA, F = 2.81, p = 0.1
    • amplitude: two-way ANOVA, F = 0.143, p = 0.86 (n = 4 recordings)

    Needs the raw AP band; depends on fig3c-volume.

  • BlockedResults, Anisotropic volumetric profile of the induced electrical field, paragraph 3; Figure 3H· Data missing from the dataset

    The centre of mass of the evoked-potential volume is displaced from the stimulating contact (displacement greater than zero), with a dorsal component of 105 +- 46 um away from the corpus callosum that does not reach significance.

    The displacement magnitude exceeded zero (one-sample t test vs. 0, p = 0.02), indicating anisotropy. The dorsal (superficial) component averaged 105 ± 46 μ m away from the corpus callosum and, while directionally consistent, did not reach significance (one-sample t test vs. 0, p = 0.1095)

    Results, Anisotropic volumetric profile of the induced electrical field, paragraph 3; Figure 3H

    What the paper reports

    • displacement magnitude: one-sample t test against 0, p = 0.02 (n = 4 recordings)
    • dorsal component: one-sample t test against 0, p = 0.1095 (n = 4 recordings)

    Needs the raw AP band; depends on fig3c-volume.

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 2; Figure 7D· Data missing from the dataset

    Mice learned to detect high-contrast static gratings over 1353.8 +- 496.8 rewarded high-contrast trials and 14802 +- 5154 total trials.

    Mice successfully learned to detect high-contrast static gratings (>30% contrast) over an average of 1353.8 ± 496.8 rewarded high-contrast trials and 14802 ± 5154 total trials

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 2; Figure 7D

    The behaviour of the four trained mice is not in the dandiset, which holds only the acute recordings; the paper names only these NWB files as deposited, and its code reads the behaviour from local CSV files. No other dandiset of the lab holds it (DANDI search for Denman and Hickman).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 2; Figure 7E· Data missing from the dataset

    Mice generalized from gratings to the moving-circle stimulus before electrical testing (d' 2.35 +- 0.30), comparable to high-contrast gratings (1.62 +- 0.92).

    These mice generalized detection of static gratings to the moving-circle stimulus before testing electrical stimulation, achieving a sensitivity index of d ′ = 2.35 ± 0.30, comparable to their performance on high-contrast grating stimuli ( d ′ = 1.62 ± 0.92, p = 0.113)

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 2; Figure 7E

    What the paper reports

    • moving circles against gratings: paired t test, t = -2.2, p = 0.113 (n = 4 mice)
    • moving circles against 0: one-sample t test, t = 15.68, p < 0.001
    • gratings against 0: one-sample t test, t = 3.55, p = 0.03

    Behaviour not deposited (see fig7d-learning).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure 7G· Data missing from the dataset

    Mice detect single pulses at similar rates across amplitudes (d' 0.37, 0.32, 0.36 and 0.45 at 5, 25, 50 and 100 uA); amplitude has no significant effect.

    Mice detected single pulses of electrical stimulation at similar rates across amplitudes ( d ′ for 5 μ A: 0.37 ± 0.04, 25 μ A: 0.32 ± 0.14, 50 μ A: 0.36 ± 0.06, and 100 μ A: 0.45 ± 0.14; F amplitude = 1.74, p = 0.21

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure 7G

    What the paper reports

    • one-way ANOVA, F = 1.74, p = 0.21 (n = 4 mice)

    Behaviour not deposited (see fig7d-learning).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure 7H· Data missing from the dataset

    Single pulses are detected above chance (average d' 0.37 +- 0.06).

    and above chance levels (average pulse d ′ = 0.37 ± 0.06; one-sample t test compared to d ′ = 0, p = 0.001

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure 7H

    What the paper reports

    • one-sample t test against 0, t = 12.59, p = 0.001 (n = 4 mice)

    Behaviour not deposited (see fig7d-learning).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure S10C· Data missing from the dataset

    Anodal-leading pulses are marginally more detectable than cathodal-leading pulses (d' 0.44 against 0.33).

    Anodal-leading waveforms were marginally more detectable than cathodal-leading waveforms (anodal d ′ = 0.44 ± 0.06, cathodal d ′ = 0.33 ± 0.09; paired t test, p = 0.028

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure S10C

    What the paper reports

    • paired t test, p = 0.028 (n = 4 mice)

    Behaviour not deposited (see fig7d-learning).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure 7H; Figure S10D· Data missing from the dataset

    High-frequency trains (100 to 300 Hz, 10 to 50 pulses) are detected better than single pulses (d' 0.66 against 0.37), with no significant difference by frequency.

    yielded higher detection sensitivity than single pulses ( d ′ = 0.66 ± 0.11 vs. 0.37 ± 0.06; paired t test, p = 0.019). Within stimulation trains, detection sensitivity did not differ significantly by frequency

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure 7H; Figure S10D

    What the paper reports

    • trains against single pulses: paired t test, t = -4.59, p = 0.019 (n = 4 mice)
    • trains against 0: one-sample t test, t = 12.19, p = 0.001
    • frequency: one-way ANOVA, F = 3.22, p = 0.112

    Behaviour not deposited (see fig7d-learning).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure 7I· Data missing from the dataset

    Detection of visual stimuli increases with contrast for gratings and moving circles, emerges at 4% contrast for gratings and 2% for circles, and single pulses perform like 2 to 4% contrast while trains perform near 8%.

    Perceptual detection in highly trained mice emerged at 4% contrast for static gratings (one-sample t test, p = 0.035) and 2% contrast for moving circles (one-sample t test, p = 0.018, Figure 7I ). Single pulses of electrical stimulation performed comparably to visual contrasts between 2% and 4%, whereas stimulation trains performed near the level of an 8% visual contrast

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 4; Figure 7I

    What the paper reports

    • contrast, gratings: one-way ANOVA, F = 18.28, p < 0.001
    • contrast, moving circles: one-way ANOVA, F = 46.52, p < 0.001
    • gratings at 4% against 0: one-sample t test, p = 0.035
    • circles at 2% against 0: one-sample t test, p = 0.018 (n = 4 mice)

    Behaviour not deposited (see fig7d-learning).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 5; Figure 7J· Data missing from the dataset

    Within sessions, electrical stimulation becomes less detectable across trials (slope -0.004 per trial).

    Within individual sessions, electrical stimulation became less detectable across trials (mixed linear model regression

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 5; Figure 7J

    What the paper reports

    • slope on trial: linear mixed-effects model, z = -15.18, p < 0.001 (n = 4 mice)

    Behaviour not deposited (see fig7d-learning).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 6; Figure S10F· Data missing from the dataset

    Detection of electrical stimulation does not improve across sessions.

    We observed no evidence of such learning: detection performance did not increase across sessions (ordinary least-squares regression, R 2 = 0.01, F = 1.12, slope = −0.004, p slope = 0.29

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 6; Figure S10F

    What the paper reports

    • ordinary least-squares regression, F = 1.12, p = 0.29

    Behaviour not deposited (see fig7d-learning).

  • BlockedResults, Single pulses of electrical stimulation are weakly detected, paragraph 5; Figure S10E (the text cites S10F)· Data missing from the dataset

    Reaction times fall with contrast for both visual classes, while reaction times to electrical stimulation are like those to low-contrast visual stimuli.

    Reaction times decreased with increasing contrast for both classes of visual stimuli, whereas reaction times for electrical stimuli were similar to those for low-contrast visual stimuli

    Results, Single pulses of electrical stimulation are weakly detected, paragraph 5; Figure S10E (the text cites S10F)

    Behaviour not deposited (see fig7d-learning).

Findings about the data

What attempts found about the dataset itself, for its owner and the archive.

  • Data missingThe raw AP band, LFP band and NI-DAQ files are not deposited (the paper says they are available on request), so the evoked-potential analyses of Figures 1 to 3, S1 and S2 cannot be repeated (11 claims blocked), and the comparisons of units with the evoked-potential boundaries (Figures 4F, 4I and 5I) can be made only with the boundaries the paper's code computed, taken here from its repository.
  • Data missingThe behaviour of the four trained mice (Figure 7 and S10) is not deposited on DANDI or elsewhere that the paper names; the paper's code reads it from local CSV files. 10 claims are blocked for this reason.
  • Published code disagreesThe stimulation contact coordinates of jlh31 and jlh32 in the trials table differ from those in the paper's code repository by 206 um (100, 150, -100) and 373 um (370, 0, -50); the other seven recordings agree exactly. With the code's coordinates the distance results of Figure 4H give the paper's statistics and most of its values (not the 50 uA extent, 868 against 965 um, or the 50 against 100 uA KS p, 0.8245 against 0.083); with the dataset's, the 5 against 100 uA KS test is no longer significant after correction.
  • Metadata errorThe file of session jlh49v1, which the paper and its code treat as a separate recording and mouse, is filed under subject jlh48. Its subject age (P381D) is one day more than that of the jlh48 file, which was recorded four months later.
  • Description missingThe electrodes table holds placeholders (location "none" for all 1152 channels, x 5 to 53, y 0, z 3820, impedance 0), so channel CCF coordinates and regions are not given; the paper's signed distances need them and here they are fitted from the units' coordinates. The stimulation contacts have coordinates but no region.
  • Description missingThe units table holds neither the waveform classes (FS, RS, up, axon) nor the stability flags the paper uses to select units; both are recomputed here (the stability filter gives the same number of units per cathodal and anodal condition as the code's intermediate table direct_df.csv, and 1 to 5 more or fewer per bipolar condition).
  • Selection not documentedThe dandiset holds ten recordings; the paper analyses nine. The tenth, jlh41 ("lots of noise issues ... not much analyzed yet here"), has no CCF coordinates or regions, and the dandiset does not say it was left out.
  • OtherCathodal single pulses are in seven of the nine recordings (jlh39 and jlh49v1 hold only bipolar pulses, jlh40 no bipolar), while the legends of the cathodal analyses give n = 9 recordings. jlh33 has 60 trials per condition, where the Methods state 75, and no anodal 50 uA single pulse. Trial and spike times of jlh40 and jlh53concat are negative in part.
  • Published code disagreesThe FS and RS classes of the paper's code are a 0.4 ms threshold on waveform duration (dlab classify_waveform_shape called with kmeans=0), where the Methods describe k-means on three normalised features. The two give 616 and 332 FS units and different outcomes for fig5j-fs-closer and figs6b-fsrs-early-firing.
  • Methods not availableQuestion for the authors. Were the FS and RS classes of Figures 5G to 5J and S6 assigned by k-means (Methods) or by the 0.4 ms duration threshold the code runs? This attempt assumed the threshold, whose classes give the paper's statistics exactly; k-means on the Methods' three features depends on its initialisation here (296 to 470 FS units over 20 runs), and under none of those runs are FS responders significantly closer than RS responders (p 0.087 to 0.401).
  • Published code disagreesFor Figure 6D the code's HVA group takes any CCF name with "layer" outside V1, which adds superior colliculus and dentate gyrus units to LM, AM, PM and other areas (151 good stable units in the nine recordings, 91 in the seven with cathodal pulses), and its subcortical list has no CA1 (315 good stable units); the legend names LM, AM and PM, and CA1, subiculum, dentate gyrus and superior colliculus. The paper's 15.82 comes from the code's groups; with the legend's groups HVA and subcortical modulation do not differ (13.91 against 13.26, Tukey p = 0.784).
  • Published code disagreesThe code's map_region_to_group looks for "posteromedial area", which matches no CCF name in the files ("posteromedial visual area"), so 16 PM units are left out of the cortical circuit responders of Figure 6E.
  • Methods not availableQuestion for the authors. The Methods compute the modulation index over 500 to 1000 ms with 5 ms bins, the Figure 6C legend over a 500 ms window from onset with 5 ms bins, and the code over 0 to 500 ms with 10 ms bins; the printed magnitudes are the code's. Which is meant? Over 500 to 1000 ms there is no cortical against subcortical difference (p = 0.457).
  • Metadata errorIn jlh53concat, a concatenated recording, the trials table's run numbers 0 to 7 are each used by a monopolar and a bipolar run, so a run number does not identify one condition; code that selects trials by run, as the paper's ZETA code does, takes 150 trials of both polarities.
  • OtherIn all nine recordings the population fires less in the 3 ms before the logged pulse onset, down to 0.07 to 0.20 spikes per trial from 2 to 1 ms before it (0.61 in jlh31), against 0.81 to 2.03 from 8 to 7 ms before; the trials table does not say whether the onset is logged after the pulse begins or spikes were removed around it.
  • Methods not availableQuestion for the authors. Was the 4 ms minimum peak width of the Figure 6F peak times meant to apply? The code's configuration sets it, but its conversion to 5 ms bins rounds it to one bin and skips it, and the Methods give no rule. Without it time to peak does not correlate with distance (r = -0.001 and 0.047); with it both correlate (r = 0.131, p = 0.001 and r = -0.115, p = 0.004).

Questions

Where the author or a verifier believed the claims file misread the paper, they asked, citing the paper, and a reader who had not seen the data answered from the paper alone, amending the claims file where it was wrong. The author then ran the claim as answered. Attempts before protocol pilot-11 also asked here for method details the claims file did not give.

  • openasked 2026-10-31about fig4d-latency-trough

    The criterion takes the median trough over the four cathodal amplitudes, as the code does. The Results say "a median value of 1.95 ms across parameters". Does "across parameters" mean the cathodal amplitudes only, or all polarity and amplitude conditions?

    Passage. an early peak immediately after stimulation, followed by a brief pause occurring at a median value of 1.95 ms across parameters and a broader, and sustained peak beginning around 2.5 ms (Figures 4D and S3B) (Results, Evoked spiking is sparse and increases in density with amplitude, paragraph 3)

    Why it matters. Over the four cathodal amplitudes the median trough is 1.95 ms (reproduced); over all 12 polarity and amplitude conditions it is 1.45 ms, outside the 0.195 ms tolerance. The verifier reported that the claims file may misread the paper here.

    No reader has answered it yet.
  • openasked 2026-10-31about fig5bc-amplitude-clusters

    The criterion passes the claim when the fractions match "in some assignment of clusters to the labels", with no condition that a cluster labelled selective for an amplitude responds almost only at it. Is that what the paper's statement, clusters "responding almost only at one amplitude", allows? The same clause is in fig5ef-polarity-clusters.

    Passage. we observed a subset of single units demonstrating amplitude selectivity, responding maximally, and almost exclusively, at specific intermediate amplitudes ... (19.2% for 5 uA, 15.2% for 25 uA, and 7.1% for 100 uA; Figures 5A-5C, iv-vi) (Results, Stimulation amplitude and polarity tune unit-specific recruitment, paragraph 1)

    Why it matters. The only matching assignment for fig5bc labels as 100 uA-selective a cluster whose mean nDRP peaks at 50 uA (0.651 against 0.343 at 100 uA) and as 25 uA-selective one with 0.27 at 5 uA and 0.47 at 25 uA; with labels required to follow the clusters' peaks, no assignment matches and the claim would not be reproduced. For fig5ef the matching assignment labels as bipolar a cluster whose anodal response is 0.79 of its bipolar one; by the half-maximum rule bipolar is 0.205 against 0.277.

    No reader has answered it yet.
  • openasked 2026-10-31about fig6f-no-propagating-wave

    The criterion lists peaks "with prominence 0.2, width at least 4 ms" and the centre of mass over +-6 ms. These are the values in the code's configuration, but the code converts them to 5 ms bins with round(), so the width becomes one bin and is then not applied (the code applies it only above one bin). Is the width filter part of the method to test, or should the criterion follow what the code does?

    Passage. This was confirmed statistically by the absence of a relationship between distance from stimulation and time to peak inhibition (Figure 6F, ii; Pearson r = 0.09) or time to peak rebound excitation (Figure 6F, iii; Pearson r = -0.03) (Results, Single pulse electrical stimulation leads to extensive circuit response, paragraph 3); the Methods give no peak rule.

    Why it matters. With the code's rule neither correlation is significant (r = -0.001, p = 0.987; r = 0.047, p = 0.244), the outcome recorded. With a 4 ms (0.8 bin) minimum width both are (r = 0.131, p = 0.001; r = -0.115, p = 0.004), which would make the claim not reproduced.

    No reader has answered it yet.
  • openasked 2026-10-31about fig4f-responders-within-ep

    The criteria of fig4f-responders-within-ep, fig4i-proportion-within-ep and fig5i-fs-within-ep name the code's EP_distances.pkl (evoked-potential boundaries computed from the raw AP band, which the dandiset does not hold) as an outside input for the author to judge, and the claims depend on fig2g-extent-increases, which is blocked. They are tested here with that file. Should these claims be marked as using outside inputs (`with_inputs`), so that they are counted apart from the claims on the dataset?

    Passage. a responder is within when its signed distance (norm_distance) lies between the mean superficial and mean deep extents of the evoked potential at that amplitude (fig2g-extent-increases; the code reads them from EP_distances.pkl, an outside input for the author to judge) (claims file, criterion of fig4f-responders-within-ep)

    Why it matters. The three outcomes (not reproduced, reproduced, not reproduced) rest on boundaries from the paper's code, not the dandiset. As the claims file stands they are counted with the claims on the dataset.

    No reader has answered it yet.
  • openasked 2026-10-31about fig6c-modulation-cortex

    The criterion asks that the region, amplitude and interaction terms have p < 0.001 under every reading, including reading (c), the Methods' 500 to 1000 ms window. The claims file also says the paper's printed values are reached only with the code's 0 to 500 ms window. Is the test under every reading what the claim intends, or did the Methods mis-state the window, so that the claim should be tested under the code's window alone?

    Passage. the response firing rate was computed from binned spike counts in a 500–1000 ms post-stimulus window using 5 ms bins (STAR Methods, Modulation index calculation); the criterion of fig6c-modulation-cortex reads (a) code, 0 to 500 ms, 10 ms bins, (b) legend, 0 to 500 ms, 5 ms bins, (c) Methods, 500 to 1000 ms, 5 ms bins, and asks for p < 0.001 "under every reading" (claims file).

    Why it matters. Under reading (a) the values and statistics match the paper exactly. Under reading (c) the region term has p = 0.457 and the amplitude term p = 1.06e-03, so the claim is not reproduced as the criterion stands. The verifier (round 3) reported that the outcome turns on a window that does not produce the printed numbers.

    No reader has answered it yet.

Notebooks

Every registered entry that uses this dataset. An entry can hold several notebooks, and can use more than one dataset.

Reproduction Direct electrical stimulation of mouse visual cortex with three Neuropixels probes (Hickman et al. 2026)

Tests the single-unit findings of Hickman, Hughes, Sahai, Miles and Denman (2026) on dandiset 000774, which holds spike times, templates, CCF coordinates and regions of the units of nine acute recordings with three Neuropixels probes around a stimulating electrode in mouse V1, and the table of single electrical pulses. It tests the direct response window, the sparse and amplitude-dependent direct responders, their clustering by amplitude and polarity, the FS and RS differences, the circuit response, cortical state and population dimensionality. The evoked potentials need the raw AP band and the behaviour is not deposited, so those claims are blocked. All files

Verification

  • Round 1
    2026-10-31
    Agreed with 40 of 42 outcomes. Independent agent run by CatalystNeuro (verifier-1; same harness, separate context, did not write the notebook). Reviewed: Round 1. The executed notebook of the author run (commit 432df9787) with the conclusion cells and the summary removed, the helper module and the claims file, with the paper. The verifier read the archive metadata and tested helper functions on synthetic inputs. Its outcomes differ from the record on fig5ef-polarity-clusters and fig5-fsrs-proportion; it asked for fig4f, fig4i, fig5i and reading (b) of fig5-fsrs to be run with the code's EP_distances.pkl.
    • fig5ef-polarity-clusters: the verifier’s outcome was “Reproduced”. Reading a gives n = 220, CCC 0.773, silhouette 0.538. Under the criterion's "some assignment" (clusters 2 and 5 bipolar 0.286, cluster 1 anodal 0.232, cluster 3 cathodal 0.186; each is that cluster's highest polarity) all three fractions are within tolerance. The notebook did no assignment search, although it does one for fig5bc, and its half-maximum rule gives bipolar 0.205 (B1).
    • fig5-fsrs-proportion: the verifier’s outcome was “Not reproduced”. No reading gets the proportions within tolerance (FS at 5 uA 0.006 or 0.020 against 0.032, RS at 100 uA 0.022 against 0.037). The paper's F = 3.04 with p = 0.028 implies 3 df, so amplitude was categorical, and then the amplitude term has p = 0.12. Under the k-means classes the class effect is significant. Reading b, units within the evoked potential, was not run, and the paper's proportions look like that population (B2).
  • Round 2
    2026-10-31
    Agreed with 38 of 42 outcomes. Independent agent run by CatalystNeuro (verifier-2; same harness, separate context, did not write the notebook). Reviewed: Round 2. The full executed notebook of the author-revise-1 run (commit b357179b9) with its conclusions and summary, the README, the attempt record, the helper module and the claims file, with the paper. The verifier disagrees with four outcomes (fig4f, fig5i, figs3e, fig6c: partially reproduced where the criteria have no partial clause) and quoted sentences of the conclusions and the record that the outputs do not support.
    • fig4f-responders-within-ep: the verifier’s outcome was “Not reproduced”. N1. The criterion has no partial clause; RS at 5 uA is 0.400 against 0.333 +- 0.033 (cell 16), so by the rule used for fig5-fsrs this is not reproduced. The FS/RS assumption is stated but the outcome does not depend on it: k-means RS fractions also fail at 5 uA (0.520). The outcome turns on the contact coordinates.
    • figs3e-early-windows: the verifier’s outcome was “Not reproduced”. N1. The criterion has no partial clause; the 100 uA proportions are 0.058 and 0.106 against 0.092 and 0.175 (cell 34). The rule the record used to change fig5-fsrs makes this not reproduced; in round 1 it accepted partial and revises that so one rule is applied throughout.
    • fig5i-fs-within-ep: the verifier’s outcome was “Not reproduced”. N1. The criterion has no partial clause; RS at 5 uA is 0.400 (cell 49). The FS/RS assumption is stated but the outcome does not depend on it: under k-means FS is 1.000 at 5 uA from n = 2. The outcome turns on the contact coordinates.
    • fig6c-modulation-cortex: the verifier’s outcome was “Not reproduced”. N1. The criterion has no partial clause and requires every reading to hold; under reading (c) the region term has p = 0.457 (cell 57). In round 1 it accepted partial and revises that so one rule is applied throughout.
  • Round 3
    2026-10-31
    Agreed with 4 of 4 outcomes. Independent agent run by CatalystNeuro (verifier-3; same harness, separate context, did not write the notebook). Reviewed: Round 3. The four outcome changes proposed by the author-revise-2 run (commit 05897771f) and the analysis behind them in the executed notebook, with the claims file and the paper. Nothing else was reviewed.

Method differences

How the tests differed from the paper’s method or data, and the claims each applies to.

  • Signed distances as in fig4h-extent-constant.fig4h-ks-distinct fig5j-fs-closer fig5i-fs-within-ep
  • Seven recordings with cathodal pulses, where the legend gives nine.figs6b-fsrs-early-firing fig6c-modulation-cortex fig6g-dimensionality-collapse
  • Signed distances as in fig4h-extent-constant (fitted channel coordinates, the dataset's contact coordinates).fig4f-responders-within-ep fig4h-distance-decreases
  • The evoked-potential boundaries are the code's EP_distances.pkl (outside input).fig4i-proportion-within-ep fig5i-fs-within-ep
  • Assumed the code's FS and RS classes (see fig5j-fs-closer).fig5i-fs-within-ep figs6b-fsrs-early-firing
  • Latencies come from all trials of the condition's contact and all good units, as the code (figure4_figure5_figure6.ipynb cell 20), from the seven recordings that hold cathodal pulses.fig4d-latency-trough
  • Seven recordings, where the S3D legend implies nine; two recordings hold no cathodal pulses.figs3cd-gmm-window
  • The crossover is taken as the first latency after the earliest component's mean where the second component's posterior reaches the first's; the Methods say only "crossover points where adjacent component responsibilities were equal" (STAR Methods, Gaussian mixture models).figs3cd-gmm-window
  • The evoked-potential boundaries are the code's EP_distances.pkl (outside input), not computed from the raw AP band, which the dandiset does not hold.fig4f-responders-within-ep
  • Assumed the code's FS and RS classes for reading (b) (see fig5-fsrs-proportion).fig4f-responders-within-ep
  • The signed distance needs CCF coordinates of every channel, which the dandiset does not give; they are fitted per probe by a straight line through the units' coordinates on their channel number. The distances use the dataset's stimulation contact coordinates, which differ from the code's for jlh31 and jlh32.fig4h-extent-constant
  • The stability filter is recomputed from the Methods (the paper's flags are not in the dataset); the time origin of the fit and the set of units the quartiles are taken over are not stated and were chosen as the first spike of the recording and all units (a standard choice). The unit counts per cathodal and anodal condition match the code's; bipolar ones differ by 1 to 5 units.fig4j-drp-increases
  • Seven recordings in place of nine.figs3e-early-windows
  • The contact nearest layer 5 is judged from the region of the nearest unit, as the dataset gives no region for the contacts.fig5ef-polarity-clusters
  • FS and RS classes shown under the code's 0.4 ms duration threshold and the Methods' k-means; the outcome is the same under both.fig5-fsrs-proportion
  • Reading (b) uses the code's evoked-potential boundaries (outside input).fig5-fsrs-proportion
  • Assumed the code's FS and RS classes (a 0.4 ms duration threshold), where the Methods describe k-means. The code is what the paper ran, its classes give the paper's statistics exactly here, they are deterministic, and the k-means classes change with the initialisation (296 to 470 FS units over 20 runs); the k-means reading is shown. Asked of the authors (see findings).fig5j-fs-closer
  • In the code's variant each unit is put in one group (V1, then HVA, then subcortical); the code lets superior colliculus units fall in both its HVA and subcortical groups.fig6d-hva-subcortical
  • Only the cathodal conditions were tested with ZETA (7980 tests, 9 minutes); the code's full model over three polarities would take about 20 minutes more and is shown with the cached p-values.fig6e-circuit-responders
  • zetapy 2.7.2 imports scipy.signal.gaussian, which current scipy keeps under scipy.signal.windows; the helper adds it before importing zetapy.fig6e-circuit-responders
  • As the code (get_zetas), the trials are all those of the run that holds the condition's first trial, with no other filter; each run holds only the condition's single 100 us pulses, except in jlh53concat, where a run number covers 150 trials of monopolar and bipolar pulses. With jlh53concat's own 75 cathodal trials the proportions are 0.121, 0.203, 0.514 and 0.756 (cortical) and 0.095 and 0.339 (subcortical at 5 and 100 uA).fig6e-circuit-responders
  • Assumed the code's peak rule. The criterion lists a minimum width of 4 ms and a centre of mass over +-6 ms, the values in the code's configuration (figure7.ipynb cell 11); the code converts them to 5 ms bins with round(), so the width (0.8 bins, rounded to 1) is not applied, as the code applies it only above 1 bin, and the centre of mass spans +-1 bin. The paper's r values come from this rule and the Methods give none. Applying the width filter changes the outcome; asked of a reader and of the authors.fig6f-no-propagating-wave

Provenance

Attempt record
000774-hickman-2026-01
Paper
Neural population dynamics of direct electrical stimulation of neocortex
Date
2026-10-31
Made by
An agent, run by CatalystNeuro
Dataset
DANDI:000774
Version
0.260520.1753 (doi:10.48324/dandi.000774/0.260520.1753)
Run on
GitHub Actions
Harness
Claude Code
Model
claude-opus-5-5
Protocol
pilot-12

Transcripts of the agent runs

Each run’s prompt, messages and tool calls, with long output cut and the text of the paper left out, and its final report.