{"id":"178b355d-4ada-4bf1-8809-fda87b7ca3fe","arxiv_id":"2608.07168","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"A Bayesian bivariate conditional Poisson model of three EPL seasons finds a significant negative dependence between home and away goals and an attendance effect that is positive for home goals only.","lead":"This paper fits a Bayesian model to 1,140 English Premier League matches in which home and away goals are allowed to be correlated. It finds that home goals are linked to fewer away goals, and that crowd size is tied to more home scoring but not more away scoring.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Negative φ in Table 2 is plausibly an artifact of omitted team-strength effects; the central dependence claim needs a robustness check with team effects.","rationale":"The reader's weakest-assumption analysis identifies the exact load-bearing point: the model omits team-specific attack and defense strengths, so the negative dependence estimate is identified only if residual association after three covariates is match-level dynamics. I agree. Section 5 explicitly acknowledges this limitation, but the Section 4 claim of a 'statistically significant negative relationship' is stated without that caveat. The omitted team strength is not a minor nuisance: EPL schedules pair strong and weak teams, and such between-match heterogeneity alone can induce negative residual correlation in a model without team terms. Thus the 95% CI for φ excluding zero is expected even if the true within-match dependence is zero. The proposed simulation test directly checks this: generate data from an independent-Poisson model with team strengths and fit the paper's BCP H→A model; if φ is routinely negative and significant, the original result is not diagnostic of match-level dynamics. The methodology is standard and the posterior predictive checks are appropriate; the concern is about interpretation, not computation. Therefore the verdict should remain CONDITIONAL, requiring either a team-strength robustness check or a toned-down claim. No change to the reader's verdict is needed, so verdict_should_be is UNCHANGED.","tokens_in":12201,"tokens_out":6780,"duration_ms":63592,"concrete_test":"Simulate 1,140 matches from an independent-Poisson model with team-specific attack and defense parameters (Dixon–Coles style, without dependence) and season intercepts calibrated to the observed marginal goal means, then fit the paper's BCP H→A model with the same attendance and foul covariates. Repeat 100 times and record the posterior of φ. If the mean posterior of φ is negative and its 95% interval excludes zero in more than half of replications, the negative dependence reported in Table 2 is quantitatively what one expects from omitted team strength alone, so the central claim should be rephrased as a conditional association rather than a structural match-level dependence.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper's central claim (Section 4) is that the posterior for φ is centered at −0.107 with a 95% CI that excludes zero, indicating a statistically significant negative relationship between home and away goals. However, the BCP regression (Section 3) includes only attendance and fouls as covariates; it has no team-specific attack or defense strengths, as acknowledged in Section 5: 'the current specification does not explicitly account for team-specific attacking and defensive strengths.' In EPL data, strong teams often play at home against weak away teams, producing high home goals and low away goals, while weak home teams against strong away teams produce the opposite pattern. This team-strength mismatch induces a negative association in the residuals that the model attributes to φ. The reported φ is therefore identified only under the strong assumption that the residual negative correlation is match-level dynamics rather than between-match heterogeneity in team quality. Because model selection and the φ estimate come from the same data, the central claim as stated is not fully supported. The omission also affects the attendance coefficient: attendance is almost collinear with season (462 in 2020–21 vs. ~38k in other seasons), so β1,Att may partly absorb season-specific scoring environment shifts, though the authors hedge on causality in Section 4.1. The load-bearing vulnerability is that the negative dependence conclusion may be an artifact of an omitted variable; a refit or simulation with team effects is required to determine whether φ remains negative when team quality is controlled.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a Bayesian bivariate conditional Poisson (BCP) regression model for home and away goal counts in football, using a parameterization that allows both positive and negative dependence between the two counts. The model includes log attendance and fouls suffered as covariates, with separate coefficients for home and away scoring intensities. The authors fit the model to 1,140 English Premier League matches from the 2018–19, 2020–21, and 2023–24 seasons using Hamiltonian Monte Carlo in Stan with empirical Bayes priors. They compare two directional specifications (A→H and H→A) and a baseline independent Poisson model. The preferred H→A specification yields a posterior mean for the dependence parameter φ of −0.107 with a 95% credible interval [−0.147, −0.066], which the paper interprets as statistically significant negative dependence between home and away goals; attendance is positively associated with home goals but not away goals. The paper concludes that BCP models improve predictive performance over the independent Poisson model and better reproduce the observed goal correlation and home-loss proportion.","tokens_in":12513,"tokens_out":4496,"duration_ms":41828,"significance":"If the negative dependence estimate is trustworthy, the paper makes a worthwhile contribution by demonstrating that football goal counts can exhibit negative association after conditioning on few match-level covariates, going beyond the positive-only dependence of standard bivariate Poisson models. The BCP parameterization is analytically tractable (closed-form marginal moments and correlation) and the Bayesian implementation with HMC is clearly described. The paper is also transparent about several limitations, notably the absence of team-specific attack and defense strengths. However, the central empirical claim is conditional on the adequacy of a very simple covariate set, and the acknowledged omission of team strengths is directly load-bearing for the sign and significance of φ. The attendance asymmetry claim is likewise threatened by the near-collinearity of attendance with season and pandemic conditions. These identification concerns, unless addressed with robustness checks, undermine the strength of the paper's headline conclusions.","major_comments":[{"comment":"The model contains no team-specific attack or defense strength terms, as acknowledged in Section 5 ('the current specification does not explicitly account for team-specific attacking and defensive strengths'). In EPL data, strong teams frequently play at home against weak away teams, producing high home goal counts and low away goal counts, while the opposite pairing produces the reverse pattern. This between-match heterogeneity induces a negative association in the residuals even when there is no within-match dependence, so the negative posterior for φ in Table 2 may be an artifact of omitted team quality rather than evidence of match-level negative dependence. The central claim in Section 4 therefore requires a robustness check that adds team strength parameters (e.g., hierarchical attack and defense random effects) and re-examines whether the credible interval for φ still excludes zero.","section":"§3, Eq. (2); §5"},{"comment":"The home-attendance coefficient β1,Att is identified almost entirely by the 2020–21 season: log(Attendance+1) is near zero in that season while taking values around 10.5 in the other two seasons. Consequently β1,Att functions largely as a season contrast rather than a within-season attendance effect, and the seasons also differ in team composition, substitution rules, and the broader scoring environment. The paper itself concedes in Section 4.1 that 'attendance, pandemic conditions, and season are closely intertwined in these data.' The claim that attendance is positively associated with home scoring needs a sensitivity analysis with season-specific intercepts or with the pandemic season excluded before it can be presented as a substantive finding.","section":"§4.1; Table 2"}],"minor_comments":[{"comment":"Please verify the printed formula for Corr(Y1,Y2). As typeset it appears dimensionally inconsistent with the stated covariance λ1λ2(e^φ−1) and variance λ2 + λ1λ2(e^φ−1)^2; the denominator should simplify to 1 + λ1(e^φ−1)^2 rather than the expression involving an exponential of λ1(e^φ−1)^2.","section":"§3, correlation formula"},{"comment":"The statement that 'results are robust to the prior choice as investigated by the authors but not reported here to save space' is not verifiable; please include a prior sensitivity analysis in an appendix or supplementary material, or remove the claim.","section":"§4, model selection"},{"comment":"The Stan implementation is described as available 'upon request'; for a data-analysis paper, providing the code in a public repository and a reproduction script would substantially strengthen reproducibility.","section":"§5, Software and data availability"},{"comment":"There is a minor formatting issue: 'W AIC' should be 'WAIC' in the running text, and the WAIC/LOO-CV comparison should report standard errors for the ELPD differences so that the claimed improvement over the independent Poisson model can be assessed.","section":"§3.3"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of Statistica Applicata and the BCP model is a reasonable extension of existing bivariate Poisson methodology. The main concern is not the model itself but the identification of the negative dependence parameter without team strength controls; the authors acknowledge the omission but do not test its impact. This is fixable within the manuscript's scope by adding a robustness section with team-level random effects. The attendance asymmetry finding is similarly confounded with season effects. I recommend major revision rather than rejection because the methodological contribution is sound and the empirical claims can be made credible with additional sensitivity analyses."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a clean, honest application of an existing BCP parameterization to EPL data, and it makes one genuinely new empirical claim—negative dependence between home and away goals after adjusting for attendance and fouls—plus an attendance asymmetry. The Bayesian machinery is standard and the exposition is careful. But the central φ result is not yet fully supported, because the model has no team attack/defense terms; as the authors admit in Section 5, the specification \"does not explicitly account for team-specific attacking and defensive strengths.\" If strong teams beat weak teams at home, residuals will be negatively correlated even with independent Poisson scoring. With only attendance and fouls in the means, φ can absorb that between-match heterogeneity. So the credible interval excluding zero is real under the model, but the model isn't rich enough to license the match-level \"suppressive effect\" interpretation in Section 4.\n\nThat said, the soft spots are proportionate. The paper is transparent about the limitation, the authors hedge causality on attendance, and the directionality comparison (H→A vs A→H) is a reasonable way to choose a factorization. The empirical Bayes priors use MLEs from the same data, which is a mild circularity but unlikely to drive a φ this far from zero. The bigger issues are the missing team effects and the fact that 2020–21 attendance is near zero, so β_Att is partly a season indicator. Code being \"available upon request\" is a minor irritant, not a fatal flaw—the model is simple enough to reimplement.\n\nCredit where due: the paper clearly distinguishes its BCP parameterization from Petretta et al. (2025), reports convergence diagnostics, and runs posterior predictive checks. It does not oversell the causal reading.\n\nBottom line: the paper is worth a serious referee, but the referee should ask for a refit that includes team strengths (or at least a simulation demonstrating φ stays negative under team heterogeneity). I'd condition acceptance on that. If they can show φ survives, the negative-dependence finding becomes an interesting and citable result.","headline":"The negative dependence estimate is plausible but not yet identified; the paper needs a team-strength robustness check before I'd trust the headline claim.","tokens_in":13066,"tokens_out":1359,"would_cite":false,"duration_ms":12620,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Home and away goals in the English Premier League are negatively associated after accounting for stadium attendance and fouls, and crowd presence is tied to home scoring but not away scoring.","keywords":["Bayesian inference","bivariate conditional Poisson","goal dependence","home advantage","stadium attendance","posterior predictive checks","English Premier League"],"falsifier":"Fit the same BCP regression with team-specific attack and defence random effects (or a ranking-strength covariate) on the same 1,140 matches. If the 95% credible interval for $\\phi$ then straddles zero, the negative dependence conclusion is an artifact of omitted team quality rather than a match-level effect.","tokens_in":12020,"feed_emoji":"⚽","tokens_out":6239,"duration_ms":49499,"temperature":0.7,"pith_summary":"The paper tries to establish that, in the English Premier League, the number of goals scored by the home team and the away team are not independent but negatively associated once stadium attendance and fouls are accounted for. Using a Bayesian bivariate conditional Poisson regression on 1,140 matches across three seasons, the posterior for the dependence parameter $\\phi$ is centered at $-0.107$ with its 95% credible interval entirely below zero. The paper also claims an asymmetry in crowd effects: attendance is positively associated with home scoring but shows no clear association with away scoring. This matters because standard bivariate Poisson models for football scores only permit non-negative correlation, so they cannot represent the tactical suppression of the opponent's scoring that a team protecting a lead may produce.","feed_headline":"Home goals tied to about 10% fewer away goals in the EPL","feed_subtitle":"A Bayesian model of 1,140 EPL matches finds negative goal dependence after crowds and fouls; attendance boosts only home scoring.","key_machinery":"The central object is the reparameterized bivariate conditional Poisson (BCP) distribution, which specifies one goal count marginally as Poisson($\\lambda_1$) and the other conditionally as Poisson($\\mu_2 e^{\\phi y_1}$), with $\\phi \\in \\mathbb{R}$ controlling the sign of dependence. This construction gives a closed-form correlation whose sign follows the sign of $(e^{\\phi} - 1)$, allowing negative association, unlike the classic shared-component bivariate Poisson model, which forces covariance $\\lambda_0 \\geq 0$. The BCP regression places log-linear means on $\\lambda_1$ and $\\lambda_2$ with covariates log(attendance+1) and fouls suffered by each team, and the parameter $\\phi$ carries the paper's main inferential weight.","core_discovery":"On the paper's own terms, the central discovery is that the BCP H→A model fitted to the 2018-19, 2020-21, and 2023-24 Premier League seasons estimates $\\phi = -0.107$ (95% credible interval $[-0.147, -0.066]$), a statistically significant negative dependence between home and away goals after adjusting for log attendance and fouls suffered. Because a one-goal increase in the home score multiplies the conditional expected away goals by $\\exp(-0.107) \\approx 0.899$, the association corresponds to roughly a 10% reduction in the conditional mean away scoring. The same model finds that log attendance has a positive effect on home goals (posterior mean 0.024, 95% CI [0.013, 0.034]) and an effect centered at zero on away goals, confirming an asymmetry consistent with home advantage. The paper interprets the negative dependence as a within-match suppressive effect, such as a team protecting a lead, while noting the estimate is an adjusted association, not a causal claim.","pith_inferences":["If the negative dependence is genuine, real-time prediction models for football should use joint score distributions that permit suppression rather than only positive-copula constructions.","Adding team-specific attack and defence random effects to the same data is a direct test: if the credible interval for $\\phi$ then includes zero, the estimated negative dependence is largely an artifact of team-quality mismatch.","The attendance asymmetry predicts a concrete effect: in the two pandemic-affected seasons with low crowds, the home-goal attendance coefficient should shrink or vanish when estimated season-by-season.","The inconclusive foul effects suggest that granular match-event data, such as expected goals or cards, could separate attacking pressure from refereeing effects in the away-foul-to-home-goals pattern."],"forward_implications":["Both BCP factorizations beat the independent Poisson model on ELPD-LOO (-3552.0 and -3552.1 versus -3565.0), so allowing dependence improves out-of-sample joint prediction.","Conditional ELPD prefers the home-to-away factorization (-1721.8 versus -1817.2), meaning home goals predict away goals better than the reverse, as a predictive factorization rather than a causal claim.","A one-goal increase in home score is associated with about a 10% reduction in the conditional expected number of away goals, holding modeled covariates fixed.","Posterior predictive checks show the BCP models reproduce the observed negative goal correlation and the home-loss proportion, while the independent Poisson model does not.","Attendance is positively associated with home goals (posterior mean 0.024, 95% CI [0.013, 0.034]) but not with away goals (95% CI [-0.011, 0.011])."],"supporting_citations":[{"why":"Supplies the bivariate conditional Poisson construction used as the paper's model family.","marker":"Berkhout and Plug (2004)"},{"why":"Provides the reparameterization in terms of marginal means and the dependence parameter, giving closed-form moments and correlation.","marker":"Piancastelli et al. (2023b)"},{"why":"Establishes the bivariate Poisson benchmark that permits only non-negative correlation, which the BCP model extends.","marker":"Karlis and Ntzoufras (2003)"},{"why":"Introduces a dependence adjustment for low-scoring football matches, the baseline independence-assumption the BCP generalizes.","marker":"Dixon and Coles (1997)"},{"why":"Reports that crowd absence mainly reduces home-team scoring, the asymmetry replicated by the paper's attendance results.","marker":"Piancastelli et al. (2023a)"},{"why":"Presents an alternative bivariate conditional Poisson formulation with intractable marginals, contrasted with the paper's tractable parameterization.","marker":"Petretta et al. (2025)"}],"fun_headline_variants":["Bayesian model: each extra home goal predicts 10% fewer away goals in EPL","Home-away goal link is negative in Premier League, Bayesian study shows","EPL: home goals tied to lower away scoring after crowd, foul adjustment","Goal dependence in EPL is negative: -0.107 adjusted correlation","Attendance boosts home scoring but not away in EPL, model finds"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The model has no team-strength terms, so the negative dependence $\\phi$ is only a structural feature of goal dynamics if the residual association after controlling for attendance and fouls is not actually driven by mismatched team quality.","fun_headline_variants_meta":{"raw":{"variants":["Bayesian model: each extra home goal predicts 10% fewer away goals in EPL","Home-away goal link is negative in Premier League, Bayesian study shows","EPL: home goals tied to lower away scoring after crowd, foul adjustment","Goal dependence in EPL is negative: -0.107 adjusted correlation","Attendance boosts home scoring but not away in EPL, model finds"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000739,"raw_usage":{"total_tokens":3305,"prompt_tokens":953,"completion_tokens":2352,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":569,"completion_tokens_details":{"reasoning_tokens":2252}},"tokens_in":569,"tokens_out":2352,"duration_ms":14354,"temperature":1.0,"reasoning_tokens":2252,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T13:38:50.717055+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the same BCP regression with team-specific attack and defence random effects (or a ranking-strength covariate) on the same 1,140 matches. If the 95% credible interval for $\\phi$ then straddles zero, the negative dependence conclusion is an artifact of omitted team quality rather than a match-level effect.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the bivariate conditional Poisson construction used as the paper's model family."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the bivariate Poisson benchmark that permits only non-negative correlation, which the BCP model extends."},{"cited_title":"J., and Coles, S","cited_arxiv_id":null,"evidence_quote":"Introduces a dependence adjustment for low-scoring football matches, the baseline independence-assumption the BCP generalizes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Presents an alternative bivariate conditional Poisson formulation with intractable marginals, contrasted with the paper's tractable parameterization."}],"review_version":1}