nest entry (middle; 6/21 OT-PVN neurons, p=0.01; 33/126 PVN neurons, p=0.0001
Fig. 2i; Extended Data Fig. 7a (middle); Source Data Fig. 2i
What the paper reports
- OT-PVN units: one-sided Fisher's exact test, p = 0.01 (n = 21 OT-PVN units)
- PVN units: one-sided Fisher's exact test, p = 10⁻⁴ (n = 126 PVN units)
42 of 194 units (21.6%) are activated during nest entry, outside the tolerance of both readings (pooled 23.9 to 29.2%, PVN only 23.6 to 28.8%); the Fisher test against zero gives p = 1.9e-14, and 7 units are suppressed. By mouse, V116 has 37 of 70 (52.9%), V110 1 of 31 and V111 4 of 93. Without the 40 s cap, as in the code, 2 of 194 (1.0%, p = 0.25). Without the 20 units above 1,048,500 spikes, 42 of 174 (24.1%), inside the tolerance of both readings, so this variant would reproduce the fraction. With q = 0.01 and 0.1, 14.4% and 28.4%. The paper's statement separates OT-PVN (6 of 21) and PVN (33 of 126) units; the OT-PVN part cannot be tested. The clock shared by spikes and annotations could not be established in any session. Poisson units with each unit's spike count and no response give 2 of 194 (1.0%) on the same episodes and test. Real episodes without spikes are left out of the mean in 116 units (median 17% of their episodes); counting them as 0, as the shuffles do, gives 31 of 194 (16.0%, p = 1.3e-10). Counting V111 once gives 38 to 40 of 132 (28.8 to 30.3%; 28.8% is inside the pooled reading's tolerance, so this variant would reproduce the fraction); without V116 03102018, whose clock is not known, 25 of 168 (14.9%, p = 1.1e-08). The baselines of 3 to 20 entries per session overlap an earlier nest dwell or a shepherding episode; the largest count is in V110 03032018 (20 of 72, all an earlier dwell) and the largest share in V116 03092018 (16 of 29, all shepherding).
How the test differed from the paper’s method
- Units from 4 virgins and 7 sessions (194 units), against the paper's 147 units from more animals; which sessions the paper used is not stated.
- OT-PVN and PVN units are pooled and the false discovery rate is applied over each session's units, as OT identity is not in the dandiset; the Methods apply it over OT-PVN or PVN units of a day.
- The baseline ends 0.5 s before the episode, from the code; the Methods do not give it.
- V111's three sessions probably hold the same clusters (spike counts correlated 0.93 to 0.97) and are counted as separate units, as the paper counts units per day.
- As in the code, real episodes with no spike in either window are left out of a unit's mean, while shuffled ones count as 0; counting both as 0 is shown as a variant, as is counting V111 once.