Pith. sign in

REVIEW 4 major objections 6 minor 57 references

Bayesian Multilevel Bivariate Spatial Modelling of Italian School Data

T0 review · 4 major / 6 minor · reviewed 2026-08-11 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read The paper claims that Italian municipal test scores depend on school infrastructure and broadband, but that unmeasured province-level spatial effects—spanning roughly 40 points—remain necessary to explain the geography of outcomes.

desk verdict A solid applied spatial statistics paper: the central finding that school infrastructure effects don't explain away spatial divides is plausible, but the specific effect sizes rest on unaddressed aggregation and model-selection choices. read the letter →

arxiv 2412.17710 v1 pith:ROLGXUQL submitted 2024-12-23 stat.ME stat.AP

Leonardo Cefalo , Alessio Pollice , Virgilio Gómez-Rubio This is my paper · ORCID
classification stat.MEstat.AP MSC 62F1562M3062P25
keywords ItalianschooldataInvalsiscoresbivariateICARspatialconfoundingSpatial+2.0INLAinfrastructureeffectsterritorialdisparities
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper tries to establish that after taking school infrastructure into account—centrality, ultra-broadband, and urban transport access—a spatially structured latent effect at the province level is still required to explain municipal Invalsi scores in Italian and mathematics. It builds a multilevel bivariate ICAR model estimated with INLA and applies Spatial+ 2.0 deconfounding to stop covariates and spatial effects from competing. In the preferred model, central municipalities gain about 2.5 points in mathematics, full broadband coverage about 3.3 points, and full urban transport access 2.5–2.8 points, while province-level spatial effects span roughly 40 points. If right, this means infrastructure investment alone cannot close territorial divides; unmeasured place-level factors remain decisive.

What carries the argument

The engine is the bivariate Intrinsic Conditional Autoregressive (ICAR) latent effect, defined on a graph whose nodes are macro-areas (provinces or infrastructural catchment areas), with precision matrix $\Lambda \otimes R$ scaled to remove graph-induced effects. It is estimated by INLA. To separate covariate and spatial effects, the paper uses Spatial+ 2.0: each covariate's macro-area mean is projected onto the graph Laplacian eigenvectors, the last K eigenvectors (chosen per covariate by WAIC or by shrinking Moran's I) are declared the spatial component, and only the nonspatial remainder enters the regression.

What would settle it

Compare leave-one-province-out predictive scores for the S+(2) model against a comparable model without any spatial effect: if the spatial model does not clearly predict better for all held-out provinces, the conclusion that spatially structured latent effects are necessary collapses.

Watch

Extended reading notes

Core claim

The authors find that student outcomes in the second year of Italian high school are significantly associated with the infrastructural endowment of municipalities, yet the association is only part of the story. In their preferred S+(2) model—a bivariate ICAR spatial effect at province level with Spatial+ 2.0 correction—expected advantages are about 2.5 points (mathematics) for central municipalities, 3.3 points (mathematics) for full ultra-broadband coverage, and 2.5 to 2.8 points for full urban transport access. Between-macro-area spatial effects are large and strongly correlated across subjects (posterior median correlation about 0.98), spanning roughly 40 points from the most to least advantaged provinces. The authors conclude that spatially structured latent effects are still necessary to explain different outcomes across municipalities.

Load-bearing premise

The results rest on two premises: that the Spatial+ 2.0 eigenvector cutoff removes exactly the spatial signal from each covariate, and that the 873 municipalities with at least two high schools represent all Italian municipalities.

Editorial extensions

If this is right

  • A municipality classified as central is expected to score about 2.5 points higher in mathematics than an intermediate one, even after removing spatial trends.
  • Full ultra-broadband coverage of a municipality's schools is tied to an expected gain of more than 3 points in mathematics, with a weaker and less certain effect on Italian.
  • Full urban-transport school access is worth about 2.5–2.8 points in both subjects.
  • Province-level spatial effects in the S+(2) model span roughly 40 points, with Calabria among the most disadvantaged and Lombardia the most advantaged, so infrastructure alone does not close territorial gaps.
  • The ICAR model is preferred over models without spatial effects, with independent ICAR effects, or with unstructured random effects by WAIC, DIC, and LPML.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The roughly 40-point spatial spread compared with few-point infrastructure coefficients suggests unobserved structural factors—labour markets, governance, teaching quality—may dominate test-score gaps; infrastructure policy alone would likely shift scores only modestly.
  • A natural test is to extend the same model to different school years and grades; if the province-level ICAR effects are stable across cohorts, the unmeasured place effect is a durable structural feature rather than a cohort-specific artifact.
  • Applying the deconfounding at the individual student level, if Invalsi microdata were available, could reveal whether municipal averaging masks within-municipality composition effects and whether the skew-normal error for Italian scores is truly needed.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The paper analyzes municipality-level average Invalsi scores in Mathematics and Italian for 873 Italian municipalities, using a bivariate ICAR spatial effect defined at the province or infrastructural-catchment level, together with covariates measuring centrality, ultra-broadband availability, and urban transport access. The authors apply a variant of the Spatial+ deconfounding method (Spatial+ 2.0), estimate the models with R-INLA using variational Bayes and simplified Laplace approximations, and compare a large set of model formulations. They report that infrastructure covariates are associated with higher scores and that a spatially structured latent effect remains necessary after deconfounding, with province-level spatial effects spanning roughly 40 points.

Significance. If the results are valid, the paper provides a useful applied contribution to the spatial statistics and education-policy literature, and it demonstrates a workflow for multivariate areal modeling with INLA, including scaling of ICAR priors on disconnected graphs and application of a recent Spatial+ variant. The manuscript is transparent in reporting many model comparisons, honest about effects that become nonsignificant after deconfounding (e.g., broadband on Italian scores), and detailed about computational settings, which aids reproducibility. However, the headline covariate effects are small (about 2.5-3.3 points) relative to the residual variability (residual standard deviations around 10 points), so the central empirical claims are sensitive to modeling assumptions that receive limited scrutiny. The main value is applied rather than methodological, since the method itself is largely drawn from prior work by Urdangarin et al. (2024).

major comments (4)
  1. [Section 2.1, Eqs. (1)-(2), Table 3] The response variable is a municipality-level average of student-level Invalsi scores, but the likelihood assigns a single error variance omega1 or omega2 to every municipality. For an average of n_i students, the sampling variance is approximately sigma^2_within / n_i, and n_i varies widely across municipalities and is likely correlated with population size, centrality, broadband access, and transport access. Ignoring this heteroscedasticity can bias posterior means of the covariate effects and mis-calibrate the credible intervals in Table 3. Since the reported effects are on the order of 2.5-3.3 points while the residual standard deviations are about 10 points (Table 4), even a moderate re-weighting could change which infrastructure coefficients have intervals excluding zero. The authors should report n_i, use a weighted likelihood or a measurement-error component for the aggregation step, or at least provide a sensitivity analysis that re-weights observations by n_i.
  2. [Section 2, Abstract] Only 873 of 7904 Italian municipalities are analyzed, and the paper acknowledges this in Section 2 but does not model the selection mechanism or assess representativeness. The inclusion criterion (at least two high schools) excludes small and rural municipalities, which are exactly the peripheral areas central to the paper's policy conclusions, and Trentino-Alto Adige is excluded due to missing auxiliary data. Without a missing-data model or a comparison of included and excluded municipalities on key covariates, the abstract's claim to study 'all Italian municipalities' is too strong, and the estimated associations may not generalize to the municipalities most relevant to the inner-areas narrative.
  3. [Section 3.2.1, Appendix A, Table 2] The Spatial+2.0 deconfounding relies on the assumption that the spatial component of each covariate is exactly the part spanned by the last K eigenvectors of the graph Laplacian, with K selected by WAIC or by Moran's I on the same dataset. The model selection evidence for the preferred deconfounding pattern is very weak: in Table 2, the WAIC differences among the base model, S+(1), and S+(2) are less than one unit (13379.149, 13378.651, 13378.303), which is far below the scale at which WAIC differences are usually considered meaningful. The authors should provide a sensitivity analysis across a range of K values, or a simulation-based check of the eigenvector-cutoff assumption, before presenting the deconfounded coefficients in Table 3 as the central estimates.
  4. [Section 2.2, Table 3] The ultra-broadband activation status is imputed as zero for all schools not listed in the activation plan. If the plan does not list every school (as the paper's own wording suggests), this imputation introduces measurement error in a key covariate whose effect is a headline result (about 3.3 points for Mathematics in S+(2)). The authors should justify the completeness of the plan or examine the sensitivity of the broadband coefficient to alternative missing-status assumptions.
minor comments (6)
  1. [Table 1 caption] The caption reads 'T able 1' with an erroneous space; please correct the typo.
  2. [Section 4, paragraph 1] The text says 'in the next session 5' but should refer to 'Section 5'; the same issue appears in the following paragraph ('discussed in the next session 5').
  3. [References] Some references are incomplete or have formatting issues: Dupont et al. (2023) is listed as an arXiv preprint without an arXiv identifier or journal information, and the Lamouroux et al. entry repeats the arXiv number in an inconsistent way. Please harmonize the reference list.
  4. [Section 4.1, Table 2] The text states that the province-level model is 'overall preferable' based on WAIC, DIC, and LPML, but the differences among Base, S+(1), and S+(2) are smaller than one unit on all criteria. Please soften this claim or report the uncertainty in the model comparison.
  5. [Reproducibility] The paper describes running R-INLA on a single core with internal optimization disabled, which is commendable, but no statement about data and code availability is included. Please add a data-availability statement or a link to a repository with the analysis code.
  6. [Section 5, Table 3] The paper describes the BB Activation effect on Italian scores as 'possibly not significant' after deconfounding; this is accurate since the 90% interval in S+(2) is (-0.024, 4.291), but the same phrasing should be used consistently in the abstract and conclusion, where the language is more definitive about significant infrastructure associations.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the reported associations are estimated from a standard hierarchical spatial model with externally sourced data and external methodology, not reduced to fitted inputs or self-cited uniqueness claims.

full rationale

The paper's central derivation is a Bayesian bivariate ICAR regression of municipality-level Invalsi averages on exogenous infrastructural covariates. The outcome data come from the Invalsi census survey, and the covariates come from ISTAT, Infratel Italia, and the Ministry of Education; these are external sources. The ICAR prior, bivariate extension, scaling for disconnected graphs, and the Spatial+ 2.0 eigenvector decomposition are all cited from independent methodological literature (Besag, York and Mollie; Mardia; Freni-Sterrantino, Ventrucci and Rue; Dupont, Wood and Augustin; Urdangarin et al.; Lamouroux et al.), not from a self-cited uniqueness theorem or a result that presupposes the paper's conclusions. The choices of K eigenvectors by WAIC or Moran's I, and the skew-normal likelihood chosen after residual inspection, are empirical model-building decisions; they affect estimates but are not parameters fitted to the target quantities and then relabelled as predictions. There is no equation in the paper whose output is definitionally equal to an input. The self-citations that exist are to the SchoolDataIT data-retrieval package and to INLAMSA/INLAMSM software tooling; these are code and data-access tools, not load-bearing mathematical assumptions, and the underlying data and estimation routines are public and externally reproducible. Accordingly, no circular step satisfying the required quote-and-reduction standard is present.

Assumptions & free parameters 3 free parameters · 5 assumptions · 0 invented entities

The paper introduces no new physical or statistical entities. Its load-bearing assumptions are the Spatial+2.0 eigenvector decomposition, the likelihood choices, and the treatment of missing municipalities and broadband data. The main free parameters are tuning choices for the deconfounding step and the skewness prior, which are not part of a derivation but are selected from the data.

free parameters (3)
  • PC prior rate lambda for skewness = 4
    Chosen by the authors based on empirical considerations balancing model complexity and fit (Section 3).
  • Number of eigenvectors removed K per covariate and component (S+(2) province level) = e.g., Central 5, Peripheral 4, BB 5, Urban 0 on the continent; 1 for each covariate on Sicily and Sardinia
    Selected by minimizing WAIC or shrinking standardized Moran's I below 1.96 on the same dataset; affects covariate estimates (Appendix A, Table A1).
  • Moran's I threshold for eigenvector removal = standardized I < 1.96 (95th percentile of normal)
    Used in the alternative strategy to choose K; data-dependent (Section 2.4).
assumptions (5)
  • domain assumption Municipality-level Invalsi averages follow a Gaussian (math) or skew-normal (Italian) likelihood conditional on covariates and random effects
    Used in Eq. 2; the skew-normal choice for Italian is motivated by observed residual skewness from the same data.
  • domain assumption Spatial confounding is fully captured by the low-frequency eigenvectors of the graph Laplacian; removing the last K eigenvectors leaves the nonspatial covariate component
    Spatial+2.0 methodology from Urdangarin et al. 2024, applied in Section 3.2.
  • standard math The bivariate ICAR prior with sum-to-zero constraints within connected components is a valid prior for the spatial random effect
    Uses Mardia 1988 and Besag-Kooperberg 1995; the improper prior is handled by constraints and scaling.
  • domain assumption Municipalities with at least two high schools (873 of 7904) provide an ignorable selection mechanism for the associations estimated
    Invalsi privacy restriction described in Section 2.1; no missing-data model is used.
  • ad hoc to paper Schools not listed in the ultra-broadband activation plan had no broadband (imputed zero)
    Section 2.2; this imputation may misclassify schools and affect the broadband effect estimate.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Bayesian Multilevel Bivariate Spatial Modelling of Italian School Data." pith.science (2026). https://pith.science/paper/ROLGXUQL

@misc{pith2026241217710,
  author       = {Pith},
  title        = {Pith review of: Bayesian Multilevel Bivariate Spatial Modelling of Italian School Data},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/ROLGXUQL}},
  note         = {Machine review of arXiv:2412.17710}
}
read the original abstract

This paper studies the relationship between the student's abilities in the second year of high school and the infrastructural endowment in all Italian municipalities, using spatial Bayesian modelling. Municipal student scores are obtained by averaging standardized and spatially homogeneous indicators of student outcomes provided by the Invalsi Institute for two subjects, Italian and Mathematics. Given the nature of the data, we employ a multilevel regression model assuming a bivariate Intrinsic Conditionally Autoregressive (ICAR) latent effect to explain the spatial variability and account for the correlation between the two subjects. Bayesian model estimation is obtained by the Integrated Nested Laplace Approximation (INLA), implemented in the \texttt{R-INLA} package. We find that alongside a significant association with the current state of school infrastructure and facilities, spatially structured latent effects are still necessary to explain the different student outcomes across municipalities.

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

57 extracted references · 37 canonical work pages

  1. [1]

    \ Vittadini, G

    Agasisti APACrefauthors Agasisti, T. \ Vittadini, G. APACrefauthors \ 2012 . Regional economic disparities as determinants of student's achievement in Italy Regional economic disparities as determinants of student's achievement in Italy . Research in Applied Economics 4 2 33, APACrefDOI doi:http://dx.doi.org/10.5296/rae.v4i2.1316 APACrefDOI

    work page doi:10.5296/rae.v4i2.1316 2012

  2. [2]

    \ Capitanio, A

    SN APACrefauthors Azzalini, A. \ Capitanio, A. APACrefauthors \ 2014 . The Skew-Normal and Related Families The skew-normal and related families \ ( 3). Cambridge University Press

    work page 2014

  3. [3]

    , Rue, H

    INLArev APACrefauthors Bakka, H. , Rue, H. , Fuglstad, G A. , Riebler, A. , Bolin, D. , Illian, J. Lindgren, F. APACrefauthors \ 2018 . Spatial modeling with R-INLA : A review Spatial modeling with R-INLA : A review . WIREs Computational Statistics 10 6 e1443, APACrefDOI doi:https://doi.org/10.1002/wics.1443 APACrefDOI

    doi:10.1002/wics.1443 2018

  4. [4]

    , Treves, A

    WB APACrefauthors Barrett, P. , Treves, A. , Shmis, T. Ambasz, D. APACrefauthors \ 2019 . The impact of school infrastructure on learning: A synthesis of the evidence. The impact of school infrastructure on learning: A synthesis of the evidence. World Bank Publications . APACrefURL http://hdl.handle.net/10986/30920 APACrefURL

    work page 2019

  5. [5]

    APACrefauthors \ 1974

    CAR APACrefauthors Besag, J. APACrefauthors \ 1974 . Spatial interaction and the statistical analysis of lattice systems Spatial interaction and the statistical analysis of lattice systems . Journal of the Royal Statistical Society: Series B (Methodological) 36 2 192--225, APACrefDOI doi:10.1111/j.2517-6161.1974.tb00999.x APACrefDOI

    arXiv 1974

  6. [6]

    \ Kooperberg, C

    ICAR APACrefauthors Besag, J. \ Kooperberg, C. APACrefauthors \ 1995 12 . On conditional and intrinsic autoregressions On conditional and intrinsic autoregressions . Biometrika 82 4 733-746, APACrefDOI doi:https://doi.org/10.1093/biomet/82.4.733 APACrefDOI

    doi:10.1093/biomet/82.4.733 1995

  7. [7]

    , York, J

    BYM APACrefauthors Besag, J. , York, J. Molli \'e , A. APACrefauthors \ 1991 . Bayesian image restoration, with two applications in spatial statistics Bayesian image restoration, with two applications in spatial statistics . Annals of the institute of statistical mathematics 43 1--20, APACrefDOI doi:https://doi.org/10.1007/BF00116466 APACrefDOI

    doi:10.1007/bf00116466 1991

  8. [8]

    , Altman, M

    spdep APACrefauthors Bivand, R. , Altman, M. , Anselin, L. , Assun c \ a o, R. , Berke, O. , Bernat, A. Blanchet, G. APACrefauthors \ 2017 . Package ‘spdep’: Spatial dependence: Weighting schemes, statistics, R package version. Package ‘spdep’: Spatial dependence: Weighting schemes, statistics, R package version. APACrefURL https://cran.r-project.org/pack...

    work page 2017

Show all 57 references
  1. [9]

    \ Cameletti, M

    Blangiardo APACrefauthors Blangiardo, M. \ Cameletti, M. APACrefauthors \ 2015 . Spatial and spatio-temporal Bayesian models with R-INLA Spatial and spatio-temporal Bayesian models with R-INLA . John Wiley & Sons

    2015

  2. [10]

    , Checchi, D

    Bratti APACrefauthors Bratti, M. , Checchi, D. Filippin, A. APACrefauthors \ 2007 . Territorial differences in Italian students' mathematical competencies: evidence from Pisa 2003 Territorial differences in Italian students' mathematical competencies: evidence from pisa 2003 ....

    2007

  3. [11]

    , Gazzano, L

    BDI APACrefauthors Bucci, M. , Gazzano, L. , Gennari, E. , Grompone, A. , Ivaldi, G. , Messina, G. Ziglio, G. APACrefauthors \ 2023 . Per chi suona la campan (ell) a? La dotazione di infrastrutture scolastiche in Italia (For whom the bell tolls? The availability of school infr...

    2023

  4. [12]

    , Pollice, A

    SchoolDataIT APACrefauthors Cefalo, L. , Pollice, A. Maranzano, P. APACrefauthors \ 2024 . SchoolDataIT : Retrieve, Harmonise and Map Open Data regarding the Italian School System. SchoolDataIT : Retrieve, harmonise and map open data regarding the Italian school system. APACre...

    2024

  5. [13]

    \ Ord, J.K

    Cliff_Ord APACrefauthors Cliff, A.D. \ Ord, J.K. APACrefauthors \ 1981 . Spatial Processes: Models and Applications Spatial processes: Models and applications . Pion, London

    1981

  6. [14]

    , Marques, I

    DupontArXiv APACrefauthors Dupont, E. , Marques, I. Kneib, T. APACrefauthors \ 2023 . Demystifying spatial confounding. Demystifying spatial confounding

    2023

  7. [15]

    , Wood, S.N

    Dupont APACrefauthors Dupont, E. , Wood, S.N. Augustin, N.H. APACrefauthors \ 2022 . Spatial+: a novel approach to spatial confounding Spatial+: a novel approach to spatial confounding . Biometrics 78 4 1279--1290, APACrefDOI doi:https://doi.org/10.1111/biom.13656 APACrefDOI

    2022 doi

  8. [16]

    , Ventrucci, M

    Freni APACrefauthors Freni-Sterrantino, A. , Ventrucci, M. Rue, H. APACrefauthors \ 2018 . A note on intrinsic conditional autoregressive models for disconnected graphs A note on intrinsic conditional autoregressive models for disconnected graphs . Spatial and spatio-temporal ...

    2018 doi

  9. [17]

    \ Eddy, W.F

    Geisser APACrefauthors Geisser, S. \ Eddy, W.F. APACrefauthors \ 1979 . A predictive approach to model selection A predictive approach to model selection . Journal of the American Statistical Association 74 365 153--160, APACrefDOI doi:https://doi.org/10.1080/01621459.1979.104...

    1979

  10. [18]

    , Dey, D.K

    Gelfand APACrefauthors Gelfand, A.E. , Dey, D.K. Chang, H. APACrefauthors \ 1992 . Model determination using predictive distributions with implementation via sampling-based methods Model determination using predictive distributions with implementation via sampling-based method...

    1992

  11. [19]

    \ Vounatsou, P

    PCAR_Gelfand APACrefauthors Gelfand, A.E. \ Vounatsou, P. APACrefauthors \ 2003 . Proper multivariate conditional autoregressive models for spatial data analysis Proper multivariate conditional autoregressive models for spatial data analysis . Biostatistics 4 1 11--15, APACref...

    2003 doi

  12. [20]

    , Carlin, J.B

    Gelman APACrefauthors Gelman, A. , Carlin, J.B. , Stern, H.S. Rubin, D.B. APACrefauthors \ 2004 . Bayesian Data Analysis Bayesian data analysis \ ( 2nd \ ). Chapman and Hall/CRC

    2004

  13. [21]

    , Hwang, J

    GelmanWAIC APACrefauthors Gelman, A. , Hwang, J. Vehtari, A. APACrefauthors \ 2014 . Understanding predictive information criteria for Bayesian models Understanding predictive information criteria for Bayesian models . Statistics and Computing 24 6 997--1016, APACrefDOI doi:10...

    2014 doi

  14. [22]

    , Benadusi, L

    Giancola APACrefauthors Giancola, O. , Benadusi, L. Fornari, R. APACrefauthors \ 2010 . Cos \` vicine, cos \` lontane. La questione dell’equit \`a scolastica nelle regioni italiane Cos \` vicine, cos \` lontane. la questione dell’equit \`a scolastica nelle regioni italiane . S...

    2010

  15. [23]

    \ Salmieri, L

    UniromaWP1 APACrefauthors Giancola, O. \ Salmieri, L. APACrefauthors \ 2020 . Family Background, School-Track and Macro-Area: the Complex Chains of Education Inequalities in Italy. Family background, school-track and macro-area: the complex chains of education inequalities in ...

    2020

  16. [24]

    APACrefauthors \ 2020

    INLAbook APACrefauthors Gómez-Rubio, V. APACrefauthors \ 2020 . Bayesian Inference with INLA Bayesian inference with INLA . Chapman and Hall/CRC

    2020

  17. [25]

    , Schrödle, B

    CPOINLA APACrefauthors Held, L. , Schrödle, B. Rue, H. APACrefauthors \ 2010 . Posterior and Cross-validatory Predictive Checks: A Comparison of MCMC and INLA Posterior and cross-validatory predictive checks: A comparison of MCMC and INLA . T. Kneib\ G. Tutz\ ( ), Statistical ...

    2010

  18. [26]

    , Carlin, B.P

    Hodges2003 APACrefauthors Hodges, J.S. , Carlin, B.P. Fan, Q. APACrefauthors \ 2003 . On the precision of the conditionally autoregressive prior in spatial models On the precision of the conditionally autoregressive prior in spatial models . Biometrics 59 2 317--322, APACrefDO...

    2003

  19. [27]

    \ Reich, B.J

    Hodges APACrefauthors Hodges, J.S. \ Reich, B.J. APACrefauthors \ 2010 . Adding Spatially-Correlated Errors Can Mess Up the Fixed Effect You Love Adding spatially-correlated errors can mess up the fixed effect you love . The American Statistician 64 4 325--334, APACrefDOI doi:...

    2010 arXiv

  20. [28]

    Schools Dashboard, part of the Ultra-Broadband Activation plan

    BB APACrefauthors Infratel Italia APACrefauthors \ 2024 . Schools Dashboard, part of the Ultra-Broadband Activation plan. Schools dashboard, part of the ultra-broadband activation plan. APACrefURL https://bandaultralarga.italia.it/scuole-voucher/progetto-scuole/ APACrefURL las...

    2024

  21. [29]

    Municipality Standardized Assessment Data (Invalsi Censuary Survey)

    Invalsi_IS APACrefauthors INVALSI - National Institute for the Evaluation of the Education System APACrefauthors \ 2024 . Municipality Standardized Assessment Data (Invalsi Censuary Survey). Invalsi Statistical Services. Municipality standardized assessment data (invalsi censu...

    2024

  22. [30]

    La geografia delle aree interne nel 2020: vasti territori tra potenzialità e debolezze

    InnerAreas APACrefauthors ISTAT - Italian National Institute of Statistics APACrefauthors \ 2022 . La geografia delle aree interne nel 2020: vasti territori tra potenzialità e debolezze. La geografia delle aree interne nel 2020: vasti territori tra potenzialità e debolezze. AP...

    2022

  23. [31]

    Portale Unico dei Dati sulla scuola (Unique School Data Portal)

    MIUR APACrefauthors Italian Ministry of Education, University and Research APACrefauthors \ 2024 . Portale Unico dei Dati sulla scuola (Unique School Data Portal). Portale unico dei dati sulla scuola (unique school data portal). APACrefURL https://dati.istruzione.it/opendata/ ...

    2024

  24. [32]

    Law 176 of October 25th 2007

    InvalsiLaw APACrefauthors Italian Official Journal APACrefauthors \ 2007 . Law 176 of October 25th 2007. Law 176 of october 25th 2007. APACrefURL https://www.gazzettaufficiale.it/atto/serie_generale/caricaDettaglioAtto/originario?atto.dataPubblicazioneGazzetta=2007-10-26&atto....

    2007

  25. [33]

    \ Calder, C.A

    Khan APACrefauthors Khan, K. \ Calder, C.A. APACrefauthors \ 2022 . Restricted spatial regression methods: Implications for inference Restricted spatial regression methods: Implications for inference . Journal of the American Statistical Association 117 537 482--494, APACrefDO...

    2022

  26. [34]

    , Geffroy, A

    Lamouroux APACrefauthors Lamouroux, J. , Geffroy, A. , Leblond, S. , Meyer, C. Albert, I. APACrefauthors \ 2024 . Addressing Spatial Confounding in geostatistical regression models: An R-INLA approach Addressing spatial confounding in geostatistical regression models: An R-INL...

    2024 doi arXiv

  27. [35]

    APACrefauthors \ 1988

    Mardia APACrefauthors Mardia, K. APACrefauthors \ 1988 . Multi-dimensional multivariate Gaussian Markov random fields with application to image processing Multi-dimensional multivariate Gaussian Markov random fields with application to image processing . Journal of Multivariat...

    1988 doi

  28. [36]

    APACrefauthors \ 2020

    Invalsi2020 APACrefauthors Martini, A. APACrefauthors \ 2020 . Il Divario Nord-Sud nei Risultati delle Prove INVALSI. Il divario nord-sud nei risultati delle prove invalsi. Invalsi Working Paper n. 52

    2020

  29. [37]

    \ Mignani, S

    Matteucci APACrefauthors Matteucci, M. \ Mignani, S. APACrefauthors \ 2014 . Exploring regional differences in the reading competencies of Italian students Exploring regional differences in the reading competencies of Italian students . Evaluation review 38 3 251--290, APACref...

    2014 doi

  30. [38]

    , Schmidt, A.M

    Nobre APACrefauthors Nobre, W.S. , Schmidt, A.M. Pereira, J.B. APACrefauthors \ 2021 . On the effects of spatial confounding in hierarchical models On the effects of spatial confounding in hierarchical models . International Statistical Review 89 2 302--322, APACrefDOI doi:10....

    2021 doi

  31. [39]

    PISA Data Analysis Manual: SPSS, Second Edition Pisa data analysis manual: Spss, second edition

    PISA APACrefauthors OECD APACrefauthors \ 2009 . PISA Data Analysis Manual: SPSS, Second Edition Pisa data analysis manual: Spss, second edition . Paris OECD Publishing

    2009

  32. [40]

    , Gómez-Rubio, V

    INLAMSM2 APACrefauthors Palmí-Perales, F. , Gómez-Rubio, V. , Bivand, R.S. , Cameletti, M. Rue, H. APACrefauthors \ 2023 . Bayesian Inference for Multivariate Spatial Models with INLA . Bayesian inference for multivariate spatial models with INLA . The R Journal

    2023

  33. [41]

    , Gómez-Rubio, V

    INLAMSM APACrefauthors Palmí-Perales, F. , Gómez-Rubio, V. Martinez-Beneito, M.A. APACrefauthors \ 2021 . Bayesian Multivariate Spatial Models for Lattice Data with INLA Bayesian multivariate spatial models for lattice data with INLA . Journal of Statistical Software 98 2 1--2...

    2021 doi

  34. [42]

    APACrefauthors \ 1990

    CPO APACrefauthors Pettit, L.I. APACrefauthors \ 1990 . The Conditional Predictive Ordinate for the Normal Distribution The conditional predictive ordinate for the normal distribution . Journal of the Royal Statistical Society. Series B (Methodological) 52 1 175--184, APACrefD...

    1990

  35. [43]

    , Hodges, J.S

    RHZ APACrefauthors Reich, B.J. , Hodges, J.S. Zadnik, V. APACrefauthors \ 2006 . Effects of residual smoothing on the posterior of the fixed effects in disease-mapping models Effects of residual smoothing on the posterior of the fixed effects in disease-mapping models . Biomet...

    2006

  36. [44]

    , Martino, S

    INLA APACrefauthors Rue, H. , Martino, S. Chopin, N. APACrefauthors \ 2009 . Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace...

    2009

  37. [45]

    APACrefauthors \ 1986

    Silverman APACrefauthors Silverman, B.W. APACrefauthors \ 1986 . Density Estimation for Statistics and Data Analysis Density estimation for statistics and data analysis . London Chapman and Hall

    1986

  38. [46]

    , Rue, H

    PCprior APACrefauthors Simpson, D. , Rue, H. , Riebler, A. , Martins, T.G. Sørbye, S.H. APACrefauthors \ 2017 . Penalising Model Component Complexity: A Principled, Practical Approach to Constructing Priors Penalising model component complexity: A principled, practical approac...

    2017 doi

  39. [47]

    , Best, N.G

    DIC APACrefauthors Spiegelhalter, D.J. , Best, N.G. , Carlin, B.P. Van Der Linde, A. APACrefauthors \ 2002 10 . Bayesian measures of model complexity and fit Bayesian measures of model complexity and fit . Journal of the Royal Statistical Society: Series B (Statistical Methodo...

    2002

  40. [48]

    \ Rue, H

    Sorbye APACrefauthors Sørbye, S. \ Rue, H. APACrefauthors \ 2014 . Scaling intrinsic Gaussian Markov random field priors in spatial modelling Scaling intrinsic Gaussian Markov random field priors in spatial modelling . Spatial Statistics 8 39--51, APACrefDOI doi:10.1016/j.spas...

    2014 doi

  41. [49]

    , Goicoa, T

    Urdangarin24 APACrefauthors Urdangarin, A. , Goicoa, T. , Kneib, T. Ugarte, M.D. APACrefauthors \ 2024 . A simplified spatial+ approach to mitigate spatial confounding in multivariate spatial areal models A simplified spatial+ approach to mitigate spatial confounding in multiv...

    2024

  42. [50]

    , Goicoa, T

    Urdangarin23 APACrefauthors Urdangarin, A. , Goicoa, T. Ugarte, M.D. APACrefauthors \ 2023 . Evaluating recent methods to overcome spatial confounding Evaluating recent methods to overcome spatial confounding . Revista Matemática Complutense 36 2 333--360, APACrefDOI doi:10.10...

    2023 doi

  43. [51]

    , Krainski, E

    VBINLA APACrefauthors Van Niekerk, J. , Krainski, E. , Rustand, D. Rue, H. APACrefauthors \ 2023 . A new avenue for Bayesian inference with INLA A new avenue for Bayesian inference with INLA . Computational Statistics and Data Analysis 181 , APACrefDOI doi:10.1016/j.csda.2023....

    2023

  44. [52]

    \ Rue, H

    SNprior APACrefauthors Van Niekerk, J. \ Rue, H. APACrefauthors \ 2021 . Skewed Probit Regression - Identifiability, Contraction and Reformulation Skewed probit regression - identifiability, contraction and reformulation . Revstat - Statistical Journal 19 1--22, APACrefDOI doi...

    2021 doi

  45. [53]

    \ Rue, H

    VB APACrefauthors Van Niekerk, J. \ Rue, H. APACrefauthors \ 2024 . Low-rank variational Bayes correction to the Laplace method Low-rank variational Bayes correction to the Laplace method . The Journal of Machine Learning Research 25 62 1--25, APACrefURL http://jmlr.org/papers...

    2024

  46. [54]

    , Yue, Y.R

    Wang APACrefauthors Wang, X. , Yue, Y.R. Faraway, J.J. APACrefauthors \ 2018 . Bayesian Regression Modeling with INLA Bayesian regression modeling with INLA . Chapman and Hall/CRC

    2018

  47. [55]

    APACrefauthors \ 2013

    WAIC APACrefauthors Watanabe, S. APACrefauthors \ 2013 . A widely applicable Bayesian information criterion A widely applicable Bayesian information criterion . The Journal of Machine Learning Research 14 1 867--897, APACrefURL http://jmlr.org/papers/v14/watanabe13a.html APACrefURL

    2013

  48. [56]

    APACrefauthors \ 2011

    ggplot APACrefauthors Wickham, H. APACrefauthors \ 2011 . ggplot2 ggplot2 . Wiley interdisciplinary reviews: computational statistics 3 2 180--185, APACrefURL https://cran.r-project.org/web/packages/ggplot2/index.html APACrefURL R package version 3.5.1

    2011

  49. [57]

    write newline

    " write newline " cite write " FUNCTION editor.postfix editor num.names #1 > "( )" "( )" if FUNCTION editor.trans.postfix editor num.names #1 > "( )" "( )" if FUNCTION trans.postfix translator num.names #1 > "( )" "( )" if FUNCTION authors.editors.reflist.apa5 'field := 'dot :...

Pith tools

Reviewed August 11, 2026 · model on record in the stance chip above.