Your new experience awaits. Try the new design now and help us make it even better

ORIGINAL RESEARCH article

Front. Med., 24 November 2025

Sec. Hematology

Volume 12 - 2025 | https://doi.org/10.3389/fmed.2025.1582277

This article is part of the Research TopicMathematical Approaches in Advancing Medical Science, Physical Fitness, and Clinical SciencesView all 12 articles

Time series forecasting of red blood cell demand in hematology patients using SARIMA and exponential smoothing models: a retrospective analysis in a Chinese tertiary hospital

Jusong LiuJusong Liu1Qiong YuanQiong Yuan1Liang CaiLiang Cai1Zhi CaiZhi Cai2Xuemei Xu
Xuemei Xu1*
  • 1Department of Transfusion, Zigong First People's Hospital, Zigong, Sichuan, China
  • 2Department of Transfusion, The Affiliated Hospital of Southwest Medical University, Luzhou, Sichuan, China

Background: Red blood cells (RBCs) infusion is very important for the treatment of hematology patients, but how to maintain a balanced state between the supply and demand of RBCs is still a major challenge.

Objective: This study aimed to explore the feasibility of seasonal autoregressive integrated moving average (SARIMA) model and exponential smoothing (ES) model in predicting the clinical demand of RBCs for hematology patients each month.

Methods: Our study collected the monthly RBCs usage data of hematology patients from January 2014 to December 2023 to establish the SARIMA model and ES model, respectively. Then, the optimal model was used to forecast the monthly usage of RBCs from January to June 2024, and we subsequently compared the data with actual values to evaluate the prediction effect of the model.

Results: The best fitting SARIMA model was SARIMA (2,1,0)(1,1,1)12, whose R2 = 0.603, MAE = 37.092, MAPE = 13.693, BIC = 7.896. The best fitting ES model was Winters addition model, whose R2 = 0.702, MAE = 32.617, MAPE = 12.138, BIC = 7.485. The mean relative errors of two models were 0.085 and 0.159, respectively. The SARIMA (2,1,0)(1,1,1)12 model performed better in prediction.

Conclusion: Compared with the ES model, the SARIMA model has a smaller mean relative error in predicting RBCs usage in hematology patients. DM test also verify this result. But in the future, more similar research data are needed to make research more convincing.

1 Introduction

The Department of Hematology is a department that diagnoses and treats various blood diseases. Its common disease types are mainly divided into three categories: anemia diseases, tumor diseases and bleeding diseases. Patients with blood diseases often suffer from symptoms such as anemia and bleeding. Therefore, blood transfusion is one of the most commonly used treatments. Timely and effective blood transfusion can not only maintain blood volume, increase hemoglobin concentration, improve the body’s oxygen carrying capacity, and alleviate the condition, but also supplement plasma proteins and other substances, increase the body’s immune ability, and improve the patient’s prognosis. Studies have shown that among all patients who need blood transfusion treatment, more than a quarter are hematology patients, some patients with special blood diseases even need long-term and multiple blood transfusions, this reflects the irreplaceable role of blood transfusion in the treatment of blood diseases (13).

RBCs are the most common blood component in transfusion therapy and play a vital role in the treatment of patients with blood diseases (4). However, due to the short storage time, the collection and storage of RBCs by blood stations and blood transfusion departments are greatly limited, resulting in an imbalance between clinical blood demand and the inventory of RBCs in the blood transfusion department. When the blood transfusion department is short of inventory, patients cannot receive blood transfusion treatment in time, which affects their recovery; when the inventory is in excess, the excess RBCs will lead to waste due to exceeding the storage period (57). Therefore, if the demand for RBCs in the clinic in a certain period of time in future can be accurately predicted, it will provide a certain basis and reference for blood stations to recruit blood donors and the timing of collecting RBCs. This not only ensures that the demand and supply of RBCs are in balance, but also enables blood stations and blood transfusion departments to allocate and manage RBCs more scientifically, ensure the rational use of RBCs, reduce the waste of RBCs, and allow clinical patients to receive blood transfusions in a timely manner, especially fresh blood component treatment, thereby further promoting improvement of the disease (2, 810).

In a broad sense, a time series is a sequence of values formed at various time points, and is one of the indispensable features of many data (11). It exists in many fields, such as economics, environment, meteorology, humanities, and medicine (1215). The time series model is a tool used to model regression problems. The entire modeling process has only one feature, time, and the prediction process of the time series model is to use the feature of time to predict the continuous value of the corresponding time series object. It mainly studies the changing law of the time series itself. Time series can be divided into different categories according to the number of research objects, continuity, statistical characteristics, and distribution patterns. There are different analysis methods for different time series (1618).

SARIMA model is a time series analysis method that can scientifically extract and analyze the seasonal effect, trend effect and the mutual influence of some confounding factors of the time series. It was first proposed by Box and Jenkins in 1970 and was originally designed for economics. ES model is another method which predicts the future by taking a weighted average of past data in a time series. It is a further development and improvement of the weighted moving average. ES model can integrate the information of time series and perform non-equal weight processing on data of different time periods, it is also easy to operate and adaptable. Based on the above advantages, the ARIMA model and ES model is widely used in the prediction of data such as morbidity, mortality, number of hospitalizations and drug demand in medical fields (1922).

In the past research, Guo (23) used the ARIMA model to predict RBCs demand in pediatric patients. In order to reduce the stability of the platelet supply chain, Bahareh (24) used the ARIMA model to predict the platelet dosage. Volken (25) used the ES model to estimate the trends of whole blood donation and RBCs infusion in a certain area, Han (26) used the ES model to effectively predict the blood components supply in Taiwan’s blood center. While their study applied the models to prediction of blood component dosage and showed good results, none have specifically modeled RBCs demand among hematology patients. Therefore, this study takes the RBCs usage data of the hematology patients as a time series, uses the SARIMA model and ES model to predict their RBCs demand, this is of great significance for them to be able to receive blood transfusion treatment stably.

2 Data and methods

The construction process of the SARIMA model and ES model is shown in Figure 1. The specific process from data statistics and processing to model prediction will be explained in the following content.

Figure 1
Flowchart comparing SARIMA and ES models for RBCs monthly dosage data. Both models begin with data preprocessing. The SARIMA model includes steps: stability test, parameter estimation, model test, model evaluation, and model prediction. The ES model follows with ES model selection, determining weighting coefficients, model test, model evaluation, and model prediction. Both models loop back to model test if evaluation is insufficient.

Figure 1. Model construction flow chart.

2.1 Study sample

The samples included in this study were derived from the information management system of the blood transfusion department of a large general hospital in Sichuan Province. These retrospective transfusion data included the total number of RBCs used by hematology patients each month from January 2014 to June 2024. The types of RBCs included suspended RBCs, washed RBCs, leukocyte-depleted suspended RBCs, irradiated leukocyte-depleted suspended RBCs, and frozen–thawed glycerol-depleted RBCs. The blood types included type A, type B, type O, and type AB. The specific operation process is to use the data from January 2014 to December 2023 as the original sequence for analysis, the data from January 2024 to June 2024 is used to verify the fitting effect of the model. The amount of RBCs is expressed in U. One unit of RBCs (U) is obtained from 200 mL of whole blood. In addition, the specifications of RBCs are mainly divided into three specifications: 1 u, 1.5 u, and 2 u, some newborns may use 0.25 u, therefore, monthly counts are reported as decimals.

2.2 Model construction

2.2.1 SARIMA model

The core idea of SARIMA model is to use the historical time data of the variable to predict itself by describing the relationship between the current value and the historical value. The model must meet the requirements of stationarity. For a non-stationary sequence, the trend component can be eliminated through a certain number of differential operations to make it a stationary sequence, and then model it using autoregression (AR) and moving average (MA) methods.

The parameters of the SARIMA model mainly include the autoregressive order (p), difference order (d), moving average order (q), seasonal autoregressive order (P), seasonal difference order (D), seasonal moving average order (Q), and seasonal cycle length (s), expressed as SARIMA (p, d, q)(P, D, Q)s. For time series analysis, it is very critical to reasonably select the parameters of the SARIMA model, which directly affects the prediction accuracy of the model. Autocorrelation coefficient (ACF) and partial autocorrelation coefficient (PACF) graphs are important tools for selecting SARIMA model parameters. The values of p, q, P, Q are determined by the order corresponding to the tailing or truncation in the ACF and PACF graph, the values of d and D are determined by the number of differences and seasonal differences.

2.2.2 ES model

ES model uses the weighted average of the past values of the series to predict the future value. The characteristic of this method is to give greater weight to recent data, while give smaller weights to distant data. Therefore, the predicted value can not only reflect the latest information, but also reflect historical data information, making the prediction results more realistic.

There are several different forms of ES model, including Simple seasonal ES, Winters additive ES and Winters multiplicative ES. It is controlled by three parameters: smoothing coefficient α, trend coefficient γ and seasonal coefficient δ, and their values range from 0 to 1, the closer the value is to 1, the greater the weight it has in prediction. Choosing the appropriate ES type according to whether the sequence is trendy and seasonal, and then perform the model test, compare each evaluation indicator to find the optimal model.

2.3 Model test

Model test mainly includes two parts: model and parameter significance testing. The model and parameters were tested using the t statistic. If the calculated p < 0.05, the model was considered to have passed the significance test of the model and parameters, which meant that the model structure had reached the simplest level and the established model was valid. The Ljung-Box Q test is used to determine whether the residual sequence of the model is a white noise sequence. If the calculated p > 0.05, it can be considered that the residual sequence is a white noise sequence, indicating the model fully extracts the information of the time series. Only models that pass the hypothesis test can be used for time series prediction. If multiple models pass the hypothesis test, the model with a larger R2, smaller mean absolute error (MAE), and smaller mean absolute percentage error (MAPE) is selected as the optimal model. In addition, the Bayesian information criterion (BIC) can be used to assist in model selection. The smaller the BIC value, the better the model fitting effect. The relevant formula is as follows.

R 2 = 1 2 2 ,

where yi represents the actual value, y ̂ i represents the predicted value, y ¯ represents the mean value, the numerator represents the sum of squared errors between the predicted value and the actual value, and the denominator represents the sum of squared errors between the mean value and the actual value. The larger the R2, the better the model fit. When the prediction model does not make any mistakes, R2 is equal to 1.

MAE = 1 m i = 1 m y ̂ i yi , MAPE = 100 % m i = 1 m y ̂ i yi yi ,

where m represents the number of samples. MAE is the average of the absolute errors between the predicted value and the true value, while MAPE is the relative size of the error between the predicted value and the true value. Both MAE and MAPE range from 0 to 1. The smaller their values, the better the model fitting effect.

BIC = k ln ( n ) 2 ln ( L ) ,

Where k represents the number of model parameters, n represents the number of samples, and L represents the likelihood function.

2.4 Model prediction

Both the SARIMA model and the ES model use the mean relative error (MRE) to evaluate the prediction effect of the model.

MRE = y ̂ i yi yi

The smaller MRE, the better the prediction effect of the model.

2.5 Statistical analysis

In this study, α = 0.05 is used as the test level, and IBM SPSS20.0 version software (IBM Corp, Armonk, NY, USA) is used to analyze the data.

3 Results

3.1 Cardinality features

The monthly RBCs usage data of hematology patients from January 2014 to December 2023 were sorted and the original sequence diagram was drawn. The results are shown in Figure 2. It can be seen from the figure that the overall usage of RBCs shows a gradual upward trend. In addition, it is found that the usage is higher in July, August, and September each year, and lower in January, February, and December. The sequence shows obvious periodic and seasonal changes.

Figure 2
Line graph showing dosage levels over time from January 2014 to September 2023. Dosage ranges from 100 to 500, with fluctuations throughout. Peaks occur intermittently, indicating variable dosage trends. Time is on the horizontal axis, dosage on the vertical axis.

Figure 2. Original sequence diagram of RBCs dosage.

3.2 Time series analysis and processing

In order to eliminate the influence of trend and seasonality of the original sequence, it was processed with first-order difference and first-order seasonal difference. The results are shown in Figure 3. By comparison, it is found that after the time series is processed by the first-order difference and the first-order seasonal difference, the upward trend is no longer obvious and to be more stable, each value in the series basically fluctuates around the zero value. Therefore, the seasonal model is used in the subsequent modeling process for research.

Figure 3
Line graph showing dosage variations over time from February 2015 to November 2023. Dosage values fluctuate significantly between two hundred and negative three hundred. Time intervals are marked bi-annually.

Figure 3. Series diagram after first-order difference and first-order seasonal difference processing.

3.3 Model parameter estimation and model testing

3.3.1 SARIMA model

Observe the ACF and PACF graph after the first-order difference and first-order seasonal difference of the time series. The results were shown in Figure 4. The series was subjected to first-order difference and first-order seasonal difference, d = 1, D = 1, ACF was truncated at 1st lag, q = 1, PACF was tailed at 3rd lag, p = 3. Both ACF and PACF are significantly non-zero at lag 12, P and Q = 1, so the alternative fitting model is tentatively determined to be SARIMA (3,1,1)(1,1,1)12.

Figure 4
Two bar graphs represent the number of admissions. The left graph shows the Autocorrelation Function (ACF), and the right graph shows the Partial Autocorrelation Function (PACF), each with different lag orders. Both graphs include coefficients and upper and lower confidence limits marked by lines. The bars fluctuate around zero, indicating varying correlations at different lags.

Figure 4. ACF and PACF after first-order difference and first-order seasonal difference.

The t-test was performed on each variable of the SARIMA (3,1,1)(1,1,1)12 model parameters, and the results are shown in Table 1. It shows that the AR hypothesis test result of the model failed (p > 0.05), so the model needs to be reselected. In order to consider the model more comprehensively, we compared various models by trying from low to high order, then combined various evaluation indicators and hypothesis test result. Finally, we found that only the parameters of the model SARIMA (2,1,0)(1,1,1)12 passed the hypothesis test (Table 1), so the final model selected was SARIMA (2,1,0)(1,1,1)12, its R2 = 0.603, MAE = 37.092, MAPE = 13.693, BIC = 7.896 (Table 2).

Table 1
www.frontiersin.org

Table 1. Parameter estimation and hypothesis testing results of SARIMA alternative models.

Table 2
www.frontiersin.org

Table 2. Evaluation indicators of selected model.

From the ACF and PACF diagrams of the SARIMA (2,1,0)(1,1,1)12 residual sequence (Figure 5), it can be seen that the values of the residual sequence of the model basically fall within the 95% confidence interval (CI). Q-Q plot and histogram show that the residual sequence is basically located on a straight line (Figure 6) and Ljung-Box Q test results show that the Q statistic is 12.684, p > 0.05, indicating that the model residual sequence does not have autocorrelation and partial autocorrelation, meets the randomness assumption, and is a white noise sequence. Therefore, it is believed that the establishment of this model is appropriate and can be used to predict the amount of RBCs.

Figure 5
Two plots comparing ACF and PACF graphs. The left plot shows results for the SARIMA model, and the right plot for the Winter addition model. Both plots display similar patterns with vertical bars representing residual correlations over delay time.

Figure 5. ACF and PACF of model residual sequence.

Figure 6
A side-by-side depiction of a Q-Q plot and a histogram. The Q-Q plot on the left shows data points closely following a straight diagonal line, suggesting normal distribution. The histogram on the right displays a frequency distribution with multiple peaks, centered around zero, indicating the spread of the same data.

Figure 6. Q-Q plot and histogram of residual sequence.

3.3.2 ES model

Different seasonal ES models are established respectively, and the hypothesis test results of each model are shown in Table 3. There were statistical differences in Alpha among the three models (p < 0.05). The simple seasonal model had no statistical significance in Delta (p > 0.05). The Winters additive model had no statistical significance in Gamma and Delta (p > 0.05). The Winters multiplicative model had no statistical significance in Gamma (p > 0.05), but had statistical differences in Delta (p < 0.05). p-value close to 1 indicates that the model has little statistical difference at this parameter level.

Table 3
www.frontiersin.org

Table 3. Statistical analysis of different ES models.

By comparing the R2, MAE, MAPE and BIC values of the three models, it can be found that the Winters additive model is the optimal model (Table 4), the R2 of the model fitting parameters is the largest, the MAE and MAPE are the smallest. The Ljung-Box Q test shows that the statistic is 9.047, p > 0.05, indicating that the residual sequence after data fitting does not have autocorrelation and partial autocorrelation, it is a white noise sequence. From the ACF and PACF diagrams of the residual sequence, it can be seen that each value of the residual sequence falls within the CI (Figure 5), so the model can be used to predict the amount of RBCs.

Table 4
www.frontiersin.org

Table 4. Evaluation indicators of different ES models.

3.4 Model prediction

The SARIMA (2,1,0)(1,1,1)12 model and the Winters addition model were used to predict the RBCs dosage from January 2024 to June 2024, respectively. The results are shown in Figure 7 and Tables 5, 6. As can be seen from Figure 7, the curve change trend of predicted values and actual values in the fitted graphs of the two models is basically consistent and the empirical coverage is 100%. The MRE of the SARIMA (2,1,0)(1,1,1)12 model and the Winters addition model are 0.085 and 0.159, respectively (Table 5). By comparing the size of MRE, the SARIMA (2,1,0)(1,1,1)12 model has better prediction effect. The Diebold-Mariano (DM) result shows that the p < 0.05, indicating that there is a significant difference in the accuracy of SARIMA prediction and ES prediction. From the perspective of the sum of square errors, the sum of square errors of SARIMA’s prediction is smaller than ES, indicating that ARIMA performs better on this data.

Figure 7
Two line graphs compare forecasting models. The top graph shows the SARIMA model with actual values, fitted values, a ninety-five percent confidence interval, and predicted values. The bottom graph illustrates the Winters addition model with similar data. Both are plotted over time, showing fluctuations between 100 and 500. A vertical line distinguishes historical and forecasted data.

Figure 7. Model prediction fit plot.

Table 5
www.frontiersin.org

Table 5. Prediction of RBCs usage by SARIMA (2,1,0)(1,1,1)12 model and Winters additive model.

Table 6
www.frontiersin.org

Table 6. The Diebold-Mariano test results.

4 Discussion

RBCs are one of the most commonly used blood products in transfusion therapy. The transfusion of RBCs can enhance the body’s oxygen-carrying capacity, it also has better clinical efficacy and a lower risk of viral and bacterial infections (7, 27). However, all blood products, including RBCs, are in short supply across the country and it is even more scarce in first-tier cities due to large population flows and abundant medical resources. In addition, due to the short shelf life of blood products themselves and imperfect storage methods, an imbalance in the supply and demand of blood products often occurs. Patients are often unable to receive blood transfusions in time, which delays the progress of their disease treatment. Therefore, it is of great significance to find a suitable method to scientifically and accurately predict the demand for RBCs and other blood products, it will provide a basis and reference for the blood collection in the future (28, 29).

The SARIMA model and ES model both use time as a single variable and fit the best prediction model through repeated identification and analysis of time series (30, 31). They are two classic statistical time series models (32). In this study, the SARIMA model and ES model were used to statistically analyze and model the monthly RBCs dosage of hematology patients from January 2014 to December 2023. Then used the optimal model SARIMA (2,1,0)(1,1,1)12 and Winters additive model to predict the monthly usage of RBCs from January to June 2024. The results showed that the mean relative error of the SARIMA model (8.5%) was smaller than ES model (15.9%). From this we can see that the prediction effect of the SARIMA model is better. According to Sarvestani’s research results (33), SARIMA can predict the demand for RBCs within 1 year. In Guo’s study (23), they used the SARIMA model to predict the RBC dosage of pediatric patients in a single center in China, with an mean relative error of 6.44% (<10%), which is consistent with our research results. This study also conducted DM test to verify the superiority of the SARIMA model, it is not available in previous studies.

It is worth mentioning that in the process of establishing the SARIMA model, the model selected by the ACF and PACF diagram did not pass the significance test of the model and parameters. According to the operation steps of the relevant research, the optimal fitting model was found by continuously trying from low order to high order. In addition, the model with a smaller BIC value cannot represent the optimal fitting model for the relevant time series prediction. Only when the model pass the significance test, can we ensure that the model extracts sufficient information about the time series, the model structure is as simple as possible and the model establishment is meaningful.

This study also has some limits. First of all, although the sample of this study included data from 10 years from 2014 to 2023, it was only a single-center study. In the future, multi-center research will be needed to obtain data with a larger sample size to make the research results more convincing. Secondly, both the SARIMA model and ES model are built based on existing time series, they do not consider the impact of emergencies such as COVID that occurred in 2020, these mutation points may affect the validity of the model. Thus, medical institutions should constantly update or modify the model according to actual conditions to ensure the accuracy of predictions and provide timely reference for clinical blood use plans.

5 Conclusion

By comparing the mean relative error and DM test results, it can be seen that the prediction effect of the SARIMA model is better and is suitable for today’s blood bank management. Medical decision makers can use models to predict the short-term future RBCs usage to ensure the balance between supply and demand, especially when blood resources are limited, it will be meaningful. In addition, the continuously revised model is also necessary, which is the key to ensuring prediction accuracy.

Data availability statement

The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author.

Ethics statement

The study was approved by the ethical committee of the Zigong First People’s Hospital, which abides by the revised World Declaration of Helsinki on ethical principles for medical research involving human subjects.

Author contributions

JL: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Validation, Writing – original draft, Writing – review & editing. QY: Data curation, Formal analysis, Investigation, Methodology, Software, Supervision, Validation, Writing – review & editing. LC: Conceptualization, Data curation, Investigation, Methodology, Writing – original draft, Writing – review & editing. ZC: Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Software, Writing – original draft, Writing – review & editing. XX: Conceptualization, Data curation, Investigation, Methodology, Resources, Supervision, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The authors declare that no Gen AI was used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmed.2025.1582277/full#supplementary-material

References

1. Rehni, AK, Cho, S, Quero, HN, Shukla, V, Zhang, Z, Dong, C, et al. Red blood cell microparticles limit hematoma growth in intracerebral hemorrhage. Stroke. (2022) 53:3182–91. doi: 10.1161/strokeaha.122.039641

PubMed Abstract | Crossref Full Text | Google Scholar

2. Bowman, Z, Fei, N, Ahn, J, Wen, S, Cumpston, A, Shah, N, et al. Single versus double-unit transfusion: safety and efficacy for patients with hematologic malignancies. Eur J Haematol. (2019) 102:383–8. doi: 10.1111/ejh.13211

PubMed Abstract | Crossref Full Text | Google Scholar

3. Candelario, N, and Klein, C. Megaloblastic anemia due to severe vitamin B(12) deficiency. Cleve Clin J Med. (2022) 89:8–9. doi: 10.3949/ccjm.89a.21041

PubMed Abstract | Crossref Full Text | Google Scholar

4. Hadjesfandiari, N, Serrano, K, Levin, E, Johal, P, Feenstra, S, Shih, AW, et al. Effect of modern infusion pumps on RBC quality. Transfusion. (2022) 62:797–808. doi: 10.1111/trf.16833

PubMed Abstract | Crossref Full Text | Google Scholar

5. Feyisa, T, Tesfaye Kiya, G, and Maleko, WA. Assessment of recipients' characteristics, transfusion appropriateness, and utilization pattern of blood and blood products in Jimma medical center, Jimma, Ethiopia. PLoS One. (2021) 16:e0250623. doi: 10.1371/journal.pone.0250623

PubMed Abstract | Crossref Full Text | Google Scholar

6. Dzik, WS, Healy, B, Brunker, P, Ruby, K, Collins, J, Paik, HI, et al. Platelet transfusion in critical care: a new method to analyze transfusion practice based on decision time intervals. Transfusion. (2023) 63:1661–76. doi: 10.1111/trf.17508

PubMed Abstract | Crossref Full Text | Google Scholar

7. Harvey, AR, Basavaraju, SV, Chung, KW, and Kuehnert, MJ. Transfusion-related adverse reactions reported to the National Healthcare Safety Network hemovigilance module, United States, 2010 to 2012. Transfusion. (2015) 55:709–18. doi: 10.1111/trf.12918

PubMed Abstract | Crossref Full Text | Google Scholar

8. Berra, L, Pinciroli, R, Stowell, CP, Wang, L, Yu, B, Fernandez, BO, et al. Autologous transfusion of stored red blood cells increases pulmonary artery pressure. Am J Respir Crit Care Med. (2014) 190:800–7. doi: 10.1164/rccm.201405-0850OC

PubMed Abstract | Crossref Full Text | Google Scholar

9. Carmel, R. Anemia and aging: an overview of clinical, diagnostic and biological issues. Blood Rev. (2001) 15:9–18. doi: 10.1054/blre.2001.0146

PubMed Abstract | Crossref Full Text | Google Scholar

10. Pomper, GJ, Pham, HP, Maracaja, DLV, and Fadeyi, EA. Emergency release transfusion practices provide an enduring reminder of festina lente-to make haste slowly. Am J Clin Pathol. (2022) 158:445–6. doi: 10.1093/ajcp/aqac085

PubMed Abstract | Crossref Full Text | Google Scholar

11. Anwar, MY, Lewnard, JA, Parikh, S, and Pitzer, VE. Time series analysis of malaria in Afghanistan: using ARIMA models to predict future trends in incidence. Malar J. (2016) 15:566. doi: 10.1186/s12936-016-1602-1

PubMed Abstract | Crossref Full Text | Google Scholar

12. Ali, H, Patel, P, Dahiya, DS, Gangwani, MK, Basuli, D, and Mohan, BP. Prediction of early-onset colorectal cancer mortality rates in the United States using machine learning. Cancer Med. (2024) 13:e6880. doi: 10.1002/cam4.6880

PubMed Abstract | Crossref Full Text | Google Scholar

13. Asante, DO, Walker, AN, Seidu, TA, Kpogo, SA, and Zou, J. Hypertension and diabetes in Akatsi South District, Ghana: modeling and forecasting. Biomed Res Int. (2022) 2022:9690964. doi: 10.1155/2022/9690964

PubMed Abstract | Crossref Full Text | Google Scholar

14. Asghar, ZB, Wankhade, P, Bell, F, Sanderson, K, Hird, K, Phung, VH, et al. Trends, variations and prediction of staff sickness absence rates among NHS ambulance services in England: a time series study. BMJ Open. (2021) 11:e053885. doi: 10.1136/bmjopen-2021-053885

PubMed Abstract | Crossref Full Text | Google Scholar

15. Barnes, SJ. Stuck in the past or living in the present? Temporal focus and the spread of COVID-19. Soc Sci Med. (2021) 280:114057. doi: 10.1016/j.socscimed.2021.114057

PubMed Abstract | Crossref Full Text | Google Scholar

16. Beard, E, Marsden, J, Brown, J, Tombor, I, Stapleton, J, Michie, S, et al. Understanding and using time series analyses in addiction research. Addiction. (2019) 114:1866–84. doi: 10.1111/add.14643

PubMed Abstract | Crossref Full Text | Google Scholar

17. Benoni, R, Sartorello, A, Uliana, M, Solomon, H, Bertolino, A, Pedot, A, et al. Epidemiological factors affecting outpatient department service utilization and hospitalization in patients with diabetes: a time-series analysis from an Ethiopian hospital between 2018 and 2021. J Glob Health. (2022) 12:04087. doi: 10.7189/jogh.12.04087

PubMed Abstract | Crossref Full Text | Google Scholar

18. Zrieq, R, Kamel, S, Boubaker, S, Algahtani, FD, Alzain, MA, Alshammari, F, et al. Time-series analysis and healthcare implications of COVID-19 pandemic in Saudi Arabia. Healthcare (Basel). (2022) 10:1874. doi: 10.3390/healthcare10101874

PubMed Abstract | Crossref Full Text | Google Scholar

19. Zukaib, U, Maray, M, Mustafa, S, Haq, NU, Khan, AUR, and Rehman, F. Impact of COVID-19 lockdown on air quality analyzed through machine learning techniques. PeerJ Comput Sci. (2023) 9:e1270. doi: 10.7717/peerj-cs.1270

PubMed Abstract | Crossref Full Text | Google Scholar

20. Xie, Y, Shi, D, Wang, X, Guan, Y, Wu, W, and Wang, Y. Prevalence trend and burden of neglected parasitic diseases in China from 1990 to 2019: findings from global burden of disease study. Front Public Health. (2023) 11:1077723. doi: 10.3389/fpubh.2023.1077723

PubMed Abstract | Crossref Full Text | Google Scholar

21. Truong, D, and Truong, MD. Projecting daily travel behavior by distance during the pandemic and the spread of COVID-19 infections - are we in a closed loop scenario? Transp Res Interdis Perspect. (2021) 9:100283. doi: 10.1016/j.trip.2020.100283

PubMed Abstract | Crossref Full Text | Google Scholar

22. de Andrade, D, Coomber, K, Livingston, M, Taylor, N, Moayeri, F, Miller, PG, et al. The impact of late-night alcohol restrictions on ambulance call-outs in entertainment precincts. Drug Alcohol Rev. (2021) 40:708–16. doi: 10.1111/dar.13308

Crossref Full Text | Google Scholar

23. Guo, K, Song, S, Qiu, L, Wang, X, and Ma, S. Prediction of red blood cell demand for pediatric patients using a time-series model: a single-center study in China. Front Med. (2022) 9:706284. doi: 10.3389/fmed.2022.706284

PubMed Abstract | Crossref Full Text | Google Scholar

24. Fanoodi, B, Malmir, B, and Jahantigh, FF. Reducing demand uncertainty in the platelet supply chain through artificial neural networks and ARIMA models. Comput Biol Med. (2019) 113:103415. doi: 10.1016/j.compbiomed.2019.103415

PubMed Abstract | Crossref Full Text | Google Scholar

25. Volken, T, Buser, A, Castelli, D, Fontana, S, Frey, BM, Rüsges-Wolter, I, et al. Red blood cell use in Switzerland: trends and demographic challenges. Blood Transfus. (2018) 16:73–82. doi: 10.2450/2016.0079-16

PubMed Abstract | Crossref Full Text | Google Scholar

26. Shih, H, and Rajendran, S. Comparison of time series methods and machine learning algorithms for forecasting Taiwan blood services foundation's blood supply. J Healthc Eng. (2019) 2019:6123745. doi: 10.1155/2019/6123745

PubMed Abstract | Crossref Full Text | Google Scholar

27. Cap, AP, Beckett, A, Benov, A, Borgman, M, Chen, J, Corley, JB, et al. Whole blood transfusion. Mil Med. (2018) 183:44–51. doi: 10.1093/milmed/usy120

PubMed Abstract | Crossref Full Text | Google Scholar

28. Muir, K, Firoz, A, Kashem, MA, Shigemura, N, and Toyoda, Y. Red blood cell transfusion prior to lung transplantation: impact on patient outcomes. ASAIO J. (2023) 69:625–31. doi: 10.1097/mat.0000000000001898

PubMed Abstract | Crossref Full Text | Google Scholar

29. Nellis, ME, Goel, R, and Karam, O. Transfusion management in pediatric oncology patients. Hematol Oncol Clin North Am. (2019) 33:903–13. doi: 10.1016/j.hoc.2019.05.011

PubMed Abstract | Crossref Full Text | Google Scholar

30. Zheng, W, Li, H, Yang, X, Wang, L, Shi, Y, Shan, H, et al. Trends and prediction in the incidence rate of hepatitis C in Shandong Province in China from 2004 to 2030. Prev Med. (2023) 177:107749. doi: 10.1016/j.ypmed.2023.107749

PubMed Abstract | Crossref Full Text | Google Scholar

31. Zhou, L, Zhao, P, Wu, D, Cheng, C, and Huang, H. Time series model for forecasting the number of new admission inpatients. BMC Med Inform Decis Mak. (2018) 18:39. doi: 10.1186/s12911-018-0616-8

PubMed Abstract | Crossref Full Text | Google Scholar

32. Yang, W, Su, A, and Ding, L. Application of exponential smoothing method and SARIMA model in predicting the number of admissions in a third-class hospital in Zhejiang Province. BMC Public Health. (2023) 23:2309. doi: 10.1186/s12889-023-17218-x

PubMed Abstract | Crossref Full Text | Google Scholar

33. Sarvestani, SE, Hatam, N, Seif, M, Kasraian, L, Lari, FS, and Bayati, M. Forecasting blood demand for different blood groups in shiraz using auto regressive integrated moving average (ARIMA) and artificial neural network (ANN) and a hybrid approaches. Sci Rep. (2022) 12:22031. doi: 10.1038/s41598-022-26461-y

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: SARIMA model, ES model, red blood cells usage, model prediction, hematology patients

Citation: Liu J, Yuan Q, Cai L, Cai Z and Xu X (2025) Time series forecasting of red blood cell demand in hematology patients using SARIMA and exponential smoothing models: a retrospective analysis in a Chinese tertiary hospital. Front. Med. 12:1582277. doi: 10.3389/fmed.2025.1582277

Received: 24 February 2025; Revised: 08 October 2025; Accepted: 07 November 2025;
Published: 24 November 2025.

Edited by:

Eleni Gavriilaki, Aristotle University of Thessaloniki, Greece

Reviewed by:

Alfry Aristo Jansen Sinlae, Universitas Katolik Widya Mandira, Indonesia
Michael Makris, AHEPA University Hospital, Greece

Copyright © 2025 Liu, Yuan, Cai, Cai and Xu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Xuemei Xu, eG14dTEwMjNAMTYzLmNvbQ==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.