Every study, searchable

Foster and Rahmstorf (2026) remove El Niño, volcanic and solar influences from five global temperature records with a regression fitted to those records. They then report high confidence that warming accelerated, using critical values that treat the adjusted series as raw data. We ran their whole procedure, adjustment and tests together, on 2000 records per dataset, simulated under their own null: a straight-line trend since 1970 plus stationary noise. With noise fitted under that null, their changepoint test at nominal level 0.05 rejected 0.059 to 0.07 of null records, and the observed statistics had full-pipeline p of 0.0005 to 0.0025 in every record. Their quadratic test fared worse, with full-pipeline p up to 0.054. Under long-memory noise the changepoint evidence weakened, to p of at most 0.04, but held. The acceleration finding rests on the changepoint test.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: A numpy port of Foster and Rahmstorf's 2026 adjustment for El Niño, volcanic and solar influences reproduces their published adjusted NASA, NOAA, HadCRU, Berkeley and ERA5 series to within 0.0001 °C, and the changepoint test of Beaulieu et al. (2024) on the adjusted annual means for 1970-2024 gives maximum slope-change statistics from 4.91 to 6.68 with the change in 2013 or 2014 and nominal AR(1) Monte Carlo p-values from 0.00014 to 0.0033, against 0.041 to 0.19 for the unadjusted series.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: When 2,000 records per dataset are simulated under Foster and Rahmstorf's null hypothesis (a linear trend since 1970, their fitted exogenous terms, and ARMA(1,1) noise fitted under the null) and passed through their whole adjustment and changepoint test, the nominal 5% test rejects 5.9% to 7.0% of null records, and the observed statistics have full-pipeline p-values from 0.0005 to 0.0025 in all five records, so the changepoint finding is robust by the pre-registered rule.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Under the same full-pipeline null, Foster and Rahmstorf's quadratic test with the ARMA(1,1) correction of Foster and Rahmstorf (2011) rejects 5.8% to 8.0% of null records at nominal 5% on monthly values and 7.2% to 9.2% on annual means, and the observed quadratic terms have full-pipeline p-values from 0.0015 to 0.054 (monthly) and 0.0005 to 0.032 (annual), below 0.02 in only three records for monthly values and two for annual means, so the reported confidence above 98% in all five records does not hold and the finding is mixed by the pre-registered rule.

1 more claim

Why did United States heart-disease deaths increase between 2018 and 2024? An updated analysis of final national death counts and published population estimates averages all six orders of a demographic decomposition. Known-age deaths increase by 28135. Population size contributes 25981.4 deaths, age composition 26512.2, and age-specific rates -24358.5. Demographic contributions outweigh the negative rate contribution. These are arithmetic contributions conditional on the source denominators, which change estimation methods during the period, rather than causal effects or individual mortality risks.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In the bundled CDC WONDER D158 final snapshot for United States residents with known age, 2018-2024 heart-disease deaths increase by 28,135. Averaging all six orders of the population-size, age-share-vector, and age-specific-rate-vector decomposition yields contributions of +25,981.4, +26,512.2, and -24,358.5 deaths, respectively (display-rounded), so positive demographic contributions outweigh the negative rate contribution. The statement is conditional on the published population denominators and specified age bins, and is arithmetic rather than causal or an individual clinical-risk estimate.

How sensitive is a focused meta-analysis of walking and jogging for depressive symptoms? A registered audit of a fixed public extraction pools 21 direct-comparison trials and 994 participants. Post-treatment Hedges g is -0.762, with a confidence interval from -1.203 to -0.321. Every leave-one-trial-out interval excludes zero, but the usual-care-only interval and the between-study prediction interval include zero. No included trial has an overall low-risk-of-bias rating. These conditional computations support neither a uniform treatment effect nor a treatment recommendation.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In the fixed Noetel et al. public extraction, the preregistered post-treatment direct comparison of walking/jogging with its specified active controls pools 21 trials and 994 participants at Hedges g=-0.762 (95% modified Hartung-Knapp confidence interval -1.203 to -0.321), all 21 single-trial-omission intervals remain below zero, but the overall prediction interval (-2.613 to 1.088) and usual-care-only confidence interval (-1.650 to 0.237) include zero and the overall-low-risk-of-bias subset has no trials, conditional on unchanged source outcomes and ratings.

Could unobserved cognitive-test scores overturn the low sensitivity of a question about worsening memory or confusion? We preregistered an analysis of public National Health and Nutrition Examination Survey data and assigned missing low-performance labels in the most and least favorable possible ways. Among 3469 examined older adults with a recorded answer, 457 lacked a Digit Symbol Substitution Test score. Observed-score sensitivity was 24.15%. Allowing arbitrary missing-score labels gives a sensitivity range of 19.27% to 33.78%; the upper endpoint has a confidence interval from 30.75% to 36.82%. Missing scores alone therefore cannot make the question identify most low performers under the prespecified definition. This extends an established low-sensitivity result with a missingness robustness analysis; it is not a dementia diagnostic study.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Among examined NHANES 2011–2014 adults aged 60 or older with recorded MCQ084 responses, unrestricted low-performance labels for missing DSST scores at the fixed cutoff of 40 yield survey-weighted SCD sensitivity bounds of 0.192660 to 0.337814, with a 95% confidence interval of 0.307465–0.368164 for the upper endpoint, below 0.5.

A regulated system nearing instability fluctuates more slowly, raising its lag-1 autocorrelation ("critical slowing down"). In a registered analysis of public sit-to-stand recordings from 27 adults, standing raised the autocorrelation of beat-to-beat systolic pressure by 0.0431 (CI 0.0066 to 0.0761). The rise was not larger with diabetes (probability of a larger rise 0.505, CI 0.2802 to 0.7308), and postural sway was not coupled to pressure beyond chance (median surrogate z -0.028). Standing shifts pressure dynamics toward slower fluctuations, but in this small sample the shift does not track diabetic status.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In 27 adults from the PhysioNet Cerebral Vasoregulation in Diabetes sit-to-stand recordings, standing with eyes open raised the lag-1 autocorrelation of linearly detrended beat-to-beat systolic blood pressure relative to the preceding 240 s of sitting, by a Hodges-Lehmann shift of 0.043 (bootstrap 95% CI 0.007 to 0.076), with 19 of 27 participants rising (one-sided Wilcoxon signed-rank p = 0.012; Holm-adjusted across two primary tests p = 0.025).
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In the same recordings, the standing-induced rise in systolic-pressure lag-1 autocorrelation was not larger in 13 adults with type 2 diabetes than in 14 controls: the probability that a diabetic participant's rise exceeds a control's was 0.505 (bootstrap 95% CI 0.28 to 0.73; one-sided Mann-Whitney p = 0.49), which rules out only differences larger than that interval's upper end, since the test had 80% power only at a standardised difference of about 1.0.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: During eyes-open standing, magnitude-squared coherence in 0.05-0.15 Hz between anteroposterior whole-body sway acceleration (horizontal ground-reaction force over body weight) and beat-to-beat systolic pressure did not exceed phase-randomised surrogates: median surrogate z-score -0.03 across 26 participants, 2 of whom exceeded 1.645 (one-sided Wilcoxon signed-rank p = 0.54).

Primes other than $2$ fall into two classes modulo $4$, congruent to $1$ or to $3$, equally common in the long run, yet the class of $3$ is almost always ahead: Chebyshev's bias. Following the race exactly at every integer $x \le 10^8$, the class of $1$ first leads at x = 26861, and leads at only 30624 of the hundred million integers, in 128 stretches, none after x = 12382326. Weighting each $x$ by $1/x$, the class of $1$ leads on a share 0.00040587 of the race, about a tenth of the limiting share Rubinstein and Sarnak (1994) computed under two standard hypotheses, so the race is still short of its limiting behavior this far out. Modulo $3$, primes congruent to $2$ lead those congruent to $1$ strictly at every $x$ from $2$ to $10^8$, consistent with the first reversal Bays and Hudson (1978) found near $6.1 \times 10^{11}$.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: The least x at which the primes up to x congruent to 1 modulo 4 outnumber those congruent to 3 modulo 4 is 26861.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Of the integers x from 1 to 100000000, the primes congruent to 1 modulo 4 lead at exactly 30624, in 128 separate stretches, the last ending at x = 12382326, and the two classes are tied at exactly 3866 integers, the last being x = 12424002.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Weighting each integer x from 1 to 100000000 by 1/x, the primes congruent to 1 modulo 4 lead on a share 0.00040587 of the race, about a tenth of the limiting share of about 0.0041 implied by Rubinstein and Sarnak's logarithmic density of 0.9959 for the other side.

2 more claims

How does transmission heterogeneity govern whether an emerging pathogen dies out or causes an epidemic? In negative binomial branching processes at reproduction number 2.5, increasing overdispersion from the Poisson limit to geometric variance to superspreading raises the extinction probability from 0.107355 to 0.4 to 0.860685. Epidemic certainty at 95 percent requires 2 introductions under Poisson transmission, 4 under geometric variance, and 20 under superspreading. Otter-Dwass evaluation establishes that solitary index cases account for 0.721943 extinction probability, or 83.88 percent of extinctions. Cluster size summation through 500 cases recovers the extinction probability with an absolute residual of 0, confirmed across 50000 simulations. Superspreading erects an extinction barrier against solitary introductions while concentrating epidemic risk in explosive clusters.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In negative binomial branching processes with reproduction number R0 = 2.50, increasing transmission overdispersion from the Poisson limit (k -> inf) to geometric variance (k = 1.00) to SARS-like superspreading (k = 0.10) increases the single-case extinction probability from 0.107355 to 0.400000 to 0.860685, expanding the critical introduction threshold for 95% epidemic certainty from 2 to 4 to 20 independent introductions.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: For SARS-like superspreading dynamics (R0 = 2.50, k = 0.10), closed-form Otter-Dwass evaluation establishes that solitary index cases without secondary spread account for 0.721943 extinction probability (83.88% of all extinctions), and cumulative summation over finite cluster sizes s from 1 through 500 recovers the complete extinction probability with an absolute residual under 1.0e-6, verified across 50000 independent branching process simulations.

What are the exact finite-sample distributions of longest increasing subsequence lengths in random permutations, and how rapidly do their moments approach the limiting Tracy-Widom law? We compute the exact distribution of longest increasing subsequence lengths for every sample size up to 50 using the Robinson-Schensted-Knuth correspondence and the hook-length formula. At the largest sample size, the exact expected length is 11.309389 with variance 1.94134, corresponding to a normalized centered mean of -1.475863 and normalized variance of 0.526961. Across all 50 sample sizes, the expected length remains strictly below two times the square root of the sample size, with maximum ratio 0.799695. Exhaustive enumeration via patience sorting independently confirms every distribution count across 409113 permutations.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: For uniformly random permutations of size n from 1 through 50, the exact distribution of longest increasing subsequence lengths computed via the Robinson-Schensted-Knuth correspondence yields an expected length of 11.309389 at n=50, with all expected lengths strictly bounded below 2*sqrt(n).
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: The benchmark independently verifies the exact distribution of longest increasing subsequence lengths against exhaustive enumeration via patience sorting for all 409113 permutations across n from 1 through 9, and confirms the Robinson-Schensted-Knuth sum-of-squares identity across all 50 sample sizes.

Adding floating-point numbers one after another loses a little at each step, and how much depends on the order of the additions. For sums of up to a million doubles, drawn uniformly from $[0, 1)$ and from $[-1, 1)$, this study measures how far four common methods land from the correctly rounded sum, in units in the last place (ulps). Adding left to right, the error averages 196.86 ulps for a million positive values and grows with their number $n$ as about $n^{\alpha}$ with fitted $\alpha$ = 0.545, close to the square root that independent rounding errors predict. Pairwise summation stays within 2 ulps at every length. Kahan's and Neumaier's compensated sums return the correctly rounded sum in all 400 draws, including mixed-sign sums that nearly cancel. Every value is reproducible exactly, on any machine and Python version.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Adding 1000000 doubles drawn uniformly from [0, 1) from left to right lands a mean of 196.86 units in the last place (at most 589) from the correctly rounded sum over 50 draws, and the mean error grows as about n to the power 0.545 for n from 1000 to 1000000, near the square-root growth that independent rounding errors predict.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Pairwise summation of the same positive draws stays within 2 units in the last place of the correctly rounded sum in all 200 draws for n from 1000 to 1000000, with a mean error of 0.4 units at n = 1000000.
• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Kahan's and Neumaier's compensated sums equal the correctly rounded sum in all 400 draws, positive and mixed-sign, including mixed-sign sums of 1000000 values whose median condition number is 991.

1 more claim

How often do stars like the Sun release superflares of $10^{34}$ erg or more? A recent Kepler analysis found once per century, far more often than earlier ones. We tested this with TESS: a pre-registered search of the two-minute light curves of 32908 stars selected as that analysis selected its Sun-like stars, with flares of known energy injected into every light curve to measure what it finds. It counted 78 superflares where the once-per-century rate predicts 23.4: 3.33 times that rate (interval 2.63 to 4.15), one per 33 years per star, far above the earlier rates. Every planned subset and sensitivity analysis gives a ratio above one. An exploratory emulation of the Kepler search suggests its correction for missed flares, measured with the profile it injected, understates the rate of flares shaped as assumed here by a factor of 1.79 between $10^{34}$ and $10^{35}$ erg; its published rates don't follow from its own figures, so how much of the gap is Kepler's remains open. All of it rests on the stars selected and the flare energies and shapes assumed.

• Published.Published: In 94,747 TESS 2-minute light curves of 32,908 stars of 5000 to 6500 K with absolute G magnitudes of 4 to 6 and RUWE below 1.4 (5,605 star-years), the pre-registered search, with its SAP test corrected after registration to require the uncorrected flux to rise, counted 78 flares of bolometric energy 10^34 erg or more where the rate Vasilyev et al. (2024) measured with Kepler predicts 23.4 under the search's measured completeness and the assumed flare durations: 3.33 times that rate (95% CI 2.63 to 4.15), about once per 33 years per star; the registered code counted 79 (3.37 times, 2.67 to 4.20).
• Published.Published: The same TESS data are incompatible with the earlier Kepler rate of one superflare of 10^34 erg or more per 2,500 years on Sun-like stars, which under the same completeness predicts 1.04 counted flares where 78 were found (Poisson probability below 10^-100).
• Published.Published: Among the 18,203 stars of 5500 to 6000 K whose median TESS variability range is below 0.3%, the 50 counted superflares are 3.63 times (95% CI 2.70 to 4.79) the 13.8 that the rate Vasilyev et al. (2024) give for their Sun-like stars predicts.

5 more claims

How much does scheduling a daily batch by historical grid intensity help when computing capacity and idle power are explicit? We evaluate seasonal clock-time profiles trained on 2023 public Great Britain data, with equal work and equal energy in every comparison. At 50% utilization, assigned carbon falls by 11.63% in 2024 and 12.81% in 2025 relative to uniform processing. Setting idle power to 50% of maximum power reduces these advantages to 3.88% and 4.27%. Some dates have higher assigned carbon after shifting. These results describe a fixed historical workload benchmark and an accounting signal; they do not measure avoided emissions.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: In the frozen Great Britain daily-batch benchmark, seasonal profiles trained on 2023 data reduce assigned carbon versus uniform processing by 11.63% and 12.81% at 50% utilization with zero idle power in 2024 and 2025, but by 3.88% and 4.27% when idle power is 50% of maximum; at 87.5% utilization with zero idle power the reductions are 3.37% and 3.64%, while the half-utilization policy increases assigned carbon on 16 of 364 and 10 of 365 evaluation dates, respectively, conditional on the specified processor model and estimated average intensity data.

Do historical electricity schedules help when jobs arrive throughout the day and cannot wait indefinitely? We evaluated fixed seasonal profiles and recent-history profiles on public Great Britain grid records. For the primary workload, seasonal scheduling lowered mean assigned carbon intensity by 11.24% in 2024 and 12.01% in 2025 relative to immediate execution. Yet 17.61% and 15.55% of jobs received higher intensity. Recent-minus-seasonal mean differences were -0.21 and -0.07 g CO2/kWh. These are retrospective accounting results for a synthetic workload, not measured emissions avoided.

• Published, Reproduced, Reviewed.Published, Reproduced, Reviewed: Conditional on the frozen Great Britain national intensity snapshot and a uniform half-hour arrival workload, three-hour jobs allowed up to six hours of start delay have 11.24% and 12.01% lower mean assigned intensity under 2023-trained seasonal profiles than immediate execution across 17,424 and 17,472 jobs in 2024 and 2025, respectively, while 17.61% and 15.55% have higher intensity and a lagged recent-history profile changes the seasonal mean by -0.21 and -0.07 g CO2/kWh, without establishing causal emissions avoided or future workload performance.