However, these associations are not maintained among AHY birds collected in the autumn over a year later (table 3). Precipitation during the wintering season prior to development has a positive association with the allometric intercept (β = 0.03, pd = 90.0%), while the relationship with the allometric slope is still not significant (β = 0.02, pd = 88.7%), and precipitation during the pre-breeding period prior to development is not associated with autumn AHY wing length
Results, Drier wintering seasons and pre-breeding periods prior to development are associated with increases in relative wing length among hatch-year and spring after-hatch-year birds, paragraph 2; Table 3; Table S6
What the paper reports
- autumn AHY preliminary model, previous pre-breeding precipitation: phylogenetic mixed model (MCMCglmm), 95% interval, 95% CI -0 to 0.06
Under reading (a) the claim holds. PrevPrevPrecipWinMo has intercept and slope effects of 0.031 (pd 90.38%, 0.6 Monte Carlo errors above 90%) and 0.027 (pd 88.60%), and in the preliminary model PrevPrecipPreBreedBase (0.027, -0.002 to 0.055) and its tarsus interaction include 0 under both interval readings and all filtering readings. Under (c) the intercept pd is 89.15%, 1.3 errors below 90% (slope 0.020, pd 81.46%). Under (d) chain 4 is stuck (R-hat 1.596), so the four-chain Monte Carlo errors are not valid and the (d) call rests on chains 1 to 3. Pooled, they give an intercept of 0.031 (pd 91.35%) and a slope of 0.030 (pd 91.72%, Monte Carlo error 0.68%, 2.5 errors above 90%; each chain 90.80% to 93.15%, chain 1 under one point from the line), so the slope effect falls on the other side of its line. The round-1 (d) fit gave a slope pd of 93.6%, and the four chains pooled give 93.79%. Signs and means hold under all three readings, and in this run each reading has at most one effect on the wrong side, so it is partially reproduced. All three calls are within a few Monte Carlo errors of the line. Round 1's (d) fit had both effects on the wrong side (intercept pd 89.4%), a case the criterion does not cover; a question asks the reader for it. The paper's intercept pd is 90.0%, on the line. The draws can depend on the CPU, since the model is compiled with -march=native as Script01 does; the run reviewed in round 1 gave other values with the same seed. The re-run of this revision, on an AMD EPYC 7763, gave the same draws as the run reviewed in round 2 (AMD EPYC 9V74).
How the test differed from the paper’s method
- lm_mod.stan run with CmdStan in place of rstan, and MCMCglmm replaced by a Gibbs sampler of the same model, as for the claims above.