Modelling Hepatitis B virus related hospital discharges in Spain: ARIMAX based liver disease forecasting tool for hospital workload and mortality progression

Lesly Acosta, Núria Soldevila, Núria Torner, Ana Avellón, Victoria Hernando, Eva Borràs, Ana Martínez, Carles Pericas, Cristina Rius, Pere Godoy and Angela Domínguez

PLOS ONE, 2026, vol. 21, issue 6, 1-11

Abstract: Introduction: Chronic hepatitis B (CHB) virus infection leads to severe complications, cirrhosis and hepatocellular carcinoma (HCC). The main objective of this study was to develop an ARIMA-based model to forecast the progression of global hospital discharges, due to cirrhosis and HCC related to chronic hepatitis B. Methods: Retrospective observational study of monthly incidence of CHB related hospitalization discharges from 2005 to 2021 in Spain. Data were obtained through the Spanish Minimum dataset of hospital discharge registry (CMBDH) of the Health Ministry. Main diagnosis of CHB, liver cirrhosis and HCC encoded by International Medical codes (ICM–9 and ICM-10) were used. Descriptive and time series analysis was performed with forecasts made for 2022 using seasonal ARIMA and ARIMAX models. Data stationarity was achieved via a square root Box-Cox transformations and differencing. Model selection used was AIC, BIC, MAPE, and forecasting precision. Analysis was performed in R (version 4.5.0). Results: The total number of discharges related to hepatitis B was 6743, 58% were due to HCC and 42% to cirrhosis diagnosis. Median age was 59 years (range: 7 to > 100), being 83.4% men. The global chronic hepatitis B (CHB) related workload values range from 10 to 55 monthly discharges, while hepatitis B related to HCC and cirrhosis range from 4 to 34 and 1–29 discharges, respectively. The best fit and 2022 forecasts found for CHB and HCC time series was obtained with the approximate Gaussian Y-ARIMAX (6,0,0) (0,1,1) _12 model. This model after treating outliers, removes seasonal patterns and captures the series’ autoregressive dynamics with an AR(6) and seasonal MA(1) noise with expression: (1 − ϕ1B − ϕ2B2 − ... − ϕ6B6)(yt −yt−12) = εt + Θ1 εt−12 with ε ~ N(0, σ²) and in the square root scale. Conclusion: Both ARIMA and ARIMAX models play critical roles in forecasting CHB-related HCC and cirrhosis, enabling better disease monitoring, healthcare resources, and intervention assessment. ARIMAX provided more accurate context-aware predictions, making it especially valuable for public health decision-making.

Date: 2026
References: Add references at CitEc
Citations:

Downloads: (external link)
https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0329751 (text/html)
https://journals.plos.org/plosone/article/file?id= ... 29751&type=printable (application/pdf)

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:plo:pone00:0329751

DOI: 10.1371/journal.pone.0329751

Access Statistics for this article

More articles in PLOS ONE from Public Library of Science
Bibliographic data for series maintained by plosone ().