Skip to main content

A comparative study on the prediction of the BP artificial neural network model and the ARIMA model in the incidence of AIDS



As a kind of widely distributed disease in China, acquired immune deficiency syndrome (AIDS) has been quickly growing each year, become a serious problem and caused serious damage to the life and health of people and the social events of China and the world because of its high fatality rate. It has been much concerned by all aspects of society. Therefore, developing early warning technology and finding the trend of early development are of quite significance to prevent and control human immunodeficiency virus (HIV)/AIDS. This study aimed to explore a suitable model for the morbidity of AIDS in China and establish a professional and feasible disease prediction model for the prevention and control works of AIDS.


At present, the traditional linear model is still utilized by most scholars to predict the incidence of HIV/AIDS. In addition, some scholars may attempt to use the nonlinear prediction model. Both prediction models showed good fitting and prediction effects. In China, the incidence of AIDS presents linear and nonlinear characteristics. In this research, the nonlinear back propagation artificial neural network (BP-ANN) model and the typical auto-regressive integrated moving average (ARIMA) linear model were applied to predict the incidence of HIV/AIDS and compare their fitting effects.


Both models were capable of predicting the expected cases of AIDS. It was seen that ARIMA and BP-ANN models could be used to forecast the monthly incidence of HIV/AIDS, but the fitting and forecasting effects of the nonlinear BP neural network model were better than those of the traditional linear ARIMA model.


In summary, it was further concluded that the BP-ANN model was a suitable way to monitor and predict the change trend and morbidity of AIDS in China.

Peer Review reports


Human Immunodeficiency Virus (HIV) is a deadly virus weakening and attacking the immunity system, which can induce Acquired Immune Deficiency Syndrome (AIDS) that is recognized as one of notifiable communicable diseases around the world [1]. During the last decades, AIDS has been seen as an epidemic that becomes a serious public health problem and social event all over the world, causes serious damage to the life and health of people and affects all aspects of society. In the global context, 36.9 million people were carried with HIV, and 0.94 million people died of HIV-associated diseases by the end of 2017 [2]. Since 1998, the number of provinces affected by HIV/AIDS has reached 31, which still sees a rapid increase in China [3]. The epidemic of AIDS/HIV has been worsened to pose serious threats to public health. Each year, it seems that new infection cases are increasing in China [4, 5]. In 2015, about 571,000 people (15 years old and above) were infected with HIV [6].

Therefore, it is a must to prevent and control the prevalence of AIDS in China. A number of policies on the prevention and control works of HIV disease have been issued by the government. In order to supervise the spread of HIV/AIDS, the National Notifiable Disease Surveillance System was organized in 1995, and the surveillance data for primarily affected populations was collected [7, 8]. Since 2004, this system has been applied to monitor the prevalence of HIV and HIV-related behaviors [9].

Over the past few years, mathematical models have been used to successfully predict the incidence of HIV/AIDS. In the 1980s, the model suggested by the Joint United Nations Programme on HIV/AIDS (UNAIDS) was adopted to forecast HIV-infected patients in many countries so as to identify the growing trend of the disease. The methods are the Workbook Method [10], Estimation and Projection Package (EPP) method [11], Spectrum AIDS Impact Model [12] as well as Asian Epidemic Model (AEM) [13]. Due to the changing incidence of AIDS, it is necessary to think through its influence factors. In these models, adequate indicators are required to fit in different estimation and prediction curves about the epidemic situation of HIV/AIDS. Otherwise, the results will greatly deviate from the actual situation. The features of four models are as follows, Workbook, the parameters required are some relatively fixed demographic indicators, including local adult population, gender composition, base of various high-risk groups and high and low values of infection rates, base of various low-risk groups and high and low values of infection rates, etc. [10]. Spectrum AIDS Impact Model, HIV-infected people receive Antiretroviral Therapy (ART) to extend their survival time. The change in survival time will affect the prediction results of SPECTRUM [14]. EPP, the number of people receiving treatment has increased with the promotion and use of condoms. The improvement of treatment methods and other prevention and control works have reduced the quality and representativeness of monitored data, which exerts a direct influence on EPP’s estimation and prediction of epidemic situations [15]. AEM, its monitoring indicators have a large number of difficult items. Monitoring data has high-quality requirements. Only on the premise of sufficient data and quality assurance can appropriate model parameters be obtained. Then, predictions can be made. Otherwise, major mistakes are easy to make [16].

Also known as the historical extension prediction method, the time series prediction method is a kind of historical data extension prediction that is a method of extrapolating and predicting the development trend of things, which can be reflected by time series. More common traditional time series prediction methods include the Auto-Regressive Integrated Moving Average (ARIMA) model, exponential smoothing method, etc., among which ARIMA is the most representative. Considered as one of the major ways to make time series analysis, the ARIMA model involves the changes of trends, random interference and periodic variations and the invariance of other related random variables during time series analysis. Earnest et al. believed that the ARIMA model was quite easy and fast to set related parameters on the prediction of communicable diseases [17]. The establishment of the ARIMA model requires collecting relevant historical data, processing data in advance according to its stability requirements, drawing the diagram of autocorrelation coefficients and partial correlation coefficients to determine the optimal model and finally use it to predict the development trend. Nowadays, ARIMA is used to estimate the mortality of influenza, malaria and other infectious diseases.

In most cases, nonlinear structures are adopted during time series analysis as adequate results cannot be obtained from linear models. In many domains, the Artificial Neural Network (ANN) is applied due to its possibility of getting over the limitations of linear models [18] and analyzing the strongly-coupled and highly-nonlinear correlations between multiple input and output variables. In nonlinear artificial neural network models, particularly the Back Propagation Artificial Neural Network (BP-ANN), the BP-ANN model can improve prediction accuracy close to various functions of arbitrary nonlinear structures [19], and accommodate more multidimensional inputs to improve the accuracy of predictions because of its inherent self-learning property, simple structure and strong simulation ability.

The data of AIDS incidence in China has shown a coexistence of linearity and nonlinearity. In this paper, it was suggested that the nonlinear relationships should exist for the monthly morbidity of AIDS while accuracy relations should not be extracted from the linear model. Two models, namely ARIMA and BP-ANN, were established to forecast the morbidity of HIV/AIDS during the period of 2007–2016. By comparison, the future growing trend of HIV/AIDS was described for early detection and warning.


ARIMA model

As a common linear model in time series analysis, the ARIMA model is usually constructed as ARIMA (p, d, q) (P, D, Q) S, p, d, q, P, D, Q and S refer to autoregressive order, number of difference, moving average order, seasonal autoregressive order, number of seasonal difference, seasonal moving average order and time-series of cyclical pattern respectively. Graphs of Auto-Correlation Function (ACF) and Partial Auto-Correlation Function (PACF) were utilized to determine the ARIMA model [20]. The construction of an optimal model needed to think about minimum Bayesian Information Criterions (BIC) and stable multi-correlation coefficient, statistically significant parameter estimates and residuals as white noise. The ARIMA model was constructed through former forecasting errors and past series values, and developed according to the following procedures: Diagnostic checking, estimation and identification. During the identification process, the ACF and PACF of transformed information would determine seasonal and non-seasonal orders. Conditional least-squares modes were used to estimate parameters. During the diagnosis process, white noise tests were conducted to verify the adequacy of the model in the series and check whether residuals were independently and positively distributed. In this way, a few ARIMA models would be possibly identified [21]. Finally, a suitable model would be selected to forecast morbidity.

BP-ANN model

As one of artificial intelligence (AI) technologies, ANN has been generally applied to fit in nonlinear models with the capability of recognizing the principles of accurate forecasting and offering help to make decisions [22]. A large number of connected nonlinear units are contained in the ANN model for data storage self-learning process [23]. Among ANN models, the BP-ANN model is a type of multi-layered feed forward neural network.

As a system with learning ability, ANN can develop knowledge so as to exceed the original knowledge level of designers. Its learning and training methods can be divided into two types: One is supervised or tutored learning in which given sample criteria are used for classification or imitation; the other is unsupervised or untutored learning in which only learning styles or certain rules are set and the specific learning content varies with the environment (namely the situation of input signal) of the system that can automatically find the characteristics and regularity of the environment.

ANN is an implicit mathematical processing method and a typical black-box modeling tool. In general, it is only necessary to give the input and output data of the modeling object instead of knowing its structure, parameters and dynamic characteristics. Through the training of information samples, the neural network can have the brain’s ability of memory and recognition. Without any prior formulas or modeling, the ANN can self-learn, obtain the mapping relationship between input and output from existing data, store the mapping relationship in each neuron in the form of multigroup weights and thresholds to constitute network knowledge, and use it to predict similar factors. Neural network models are widely used in signal processing, pattern recognition, control, analysis and prediction and other aspects because of their nonlinear characteristics, numerous parallel distribution structures as well as learning and inductive ability.

Three layers of the BP-ANN model architecture were contained in the paper, which including the input layer, the hidden layer, and the output layer. Each layer has at least one neuron, which connects to neurons in different layers. The classic structure is shown in Fig. 1. This structure is simple, clarity and can enable each neuron to establish a suitable linear or non-linear relationship between input and output, while without limiting in output between − 1 and 1. The core of BP-ANN is each neuron in the input layer as an independent variable; the hidden layer is responsible for internal operations (imitating the human brain), especially non-linear operations; each neuron in the output layer represents a dependent variable. The calculation of BP-ANN is to find the minimum value of the error function.

Fig. 1

BP-ANN classic structure diagram

Model validation and statistical comparisons

Based on the same training set, ARIMA and BP-ANN models were subsequently established to forecast exclusively experimental information. The validity of these models was evaluated by cross validation. Mean Absolute Error (MAE), Mean Square Error (MSE) and Mean Absolute Percentage Error (MAPE) were used to make a statistical comparison of forecast and real morbidity.

Information analysis based on computer software

The ARIMA model was analyzed by using software SPSS26 and Eviews6.0. Neural Network Toolbox in Matlab 2019 was used to evaluate the BP-ANN model. All the analysis results showed significant differences, namely P < 0.05.

Data sources

According to the report on statutory infectious diseases in China, the monthly data about China’s AIDS cases reported from January 2004 to December 2016 was collected as the original data to establish the models, to predict the incidence of AIDS in 2017. Compare forecast incidence of AIDS and actual incidence of 2017, to verify the model fitting effect.

In ARIMA model, The monthly incidence of AIDS in China from 2004 to 2016 was modeled, and predicting the monthly incidence in 2017. The actual value of monthly incidence in 2017 was used as a reference to verify the model. In BP-ANN model, the set of information was classified into three subsets, namely training, validation and test sets. In the training set, the incidence data of the past three years was used to predict the incidence of the fourth year in validation set. The incidence rate in January of t1-t3 years was used to estimate that in January of the t4 year; the incidence rate in February of t1-t3 years was used to estimate that in February of the t4 year, and so on. Then, the incidence rate in the same month of t2-t4 years was used to predict that in the same month of the t5 year, the same month of t3-t5 years was used to predict the incidence rate of the same month of the t6 year, in turns. Finally, the data of 2017 was selected as the test set to verify network performance. All incidence data were numbered in chronological order, with P1, P2 and P13 representing respectively the incidence data in January 2004, February 2004 and January 2005, and so on. The specific data diversity is presented in the following Table 1. Such data diversity method could be adopted to better learn and train network models, and avoid overlearning and overfitting [24].

Table 1 Three date set in BP-ANN


Features of time series analysis in the report rate of AIDS

According to the surveillance data from January 2004 to December 2016, the figure of monthly incidence rates showed a trend of sharp increase from 2010 to 2016 the peak incidence existed in 2012 (Fig. 2). Table 2 showed the average of monthly morbidity of AIDS at the period between 2004 and 2016. The annual incidence rate was between 0.2648 and 4.0211 per 100,000 people from 2004 to 2016. Figure 3 shown that the monthly incidence of AIDS in China was cyclical. The lowest point was generally between January and February of each year, and the highest point was generally from July to December of each year. In summary, the monthly incidence of AIDS in China during the 13-year period from 2004 to 2016 had been cyclical and increasing year by year.

Fig 2

The yearly incidence of AIDS/HIV in China from 2004 to 2016

Table 2 The average of yearly Incidence and growth rate of HIV/AIDS in China, 2004–2016
Fig. 3

Chinese AIDS monthly incidence from 2004 to 2016

ARIMA model

Model identification

The time series from January 2004 to December 2016 were used to establish the model for the morbidity of AIDS, which were not stationary owing to seasonality. After the natural logarithmic transformation was performed, one general difference, one seasonal difference, time plots after transformation are shown in Fig. 4. Transformed time series appeared to be quite stationary.

Fig. 4

AIDS monthly incidence transformation: natural logarithm, difference, seasonal difference

ACF and PACF were used to describe the characteristics of series, select models and determine the order of key points. ACF was utilized to explain the correlation of several adjacent data as the coefficient of the relationship between series and their own historical or stagnant series. When the lag = 1, it is the first-order autocorrelation coefficient (p = 1), which indicates that there is a correlation between adjoining points; lag = 2 means the second-order autocorrelation coefficient (p = 2), which indicates that two adjoining points are also closely related, but generally the autocorrelation coefficient in ACF does not exceed 2. The ACF in Fig. 5-a shows that the autoregressive value broke through the wireframe of confidence interval only when lag = 1, indicating that the series had a high correlation within the first order. PACF was to test whether the partial correlation coefficient of each order was statistically significant one by one from higher to lower order until the first one was significant. The order of coefficients of PACF determines the highest order in the model. As shown in the PACF diagram (Fig. 5-b), the partial regression coefficient exceeded the confidence interval when lag = 1 and 2, indicating the feasibility of modeling within two orders. Therefore, this study considered that the partial regression coefficient decrease sharply after lag = 1, so neglected lag = 2.

Fig. 5

ACF and PACF graphs of AIDS monthly incidence

As displayed from Fig. 5a and b, the model was initially determined as ARIMA(p, d, q) × (P, D, Q) s (General Multiplicative Seasonal Model). Since one general difference (d = 1) and one seasonal difference (D = 1) were performed in data pre-processing, ARIMA(p,1,q) × (p,1,q) 12 models with all order combinations for all autocorrelation delay coefficients p ≤ 1(P ≤ 1) and moving average delay coefficients q ≤ 1(Q ≤ 1) were selected as primary models.

All primary models were used to simulate and model the monthly incidence of AIDS. The statistics, BIC and parameter estimates of the models obtained are shown in Table 3. The table selected stationary R-squared and BIC with the relatively smallest value, and the model whose residual was white noise was the optimal one. According to the minimum BIC = -6.091 and white noise test for residual errors, Ljung-Box Q [18] =13.909, P > 0.05, which indicated that goodness-of-fit considered ARIMA (0,1,1) × (0,1,1)12 as the most suitable model.

Table 3 Parameter estimation and model verification of ARIMA model

Forecast analysis with ARIMA

ARIMA(0,1,1) × (0,1,1)12 was used to predict the monthly incidence of ADIS from January to December 2017. The results are shown in Table 4. It can be seen from Fig. 6 that the change trend of monthly incidence fitted by this model was basically consistent with original data, and the fitting effect was satisfactory. With the extension of prediction time, 95% confidence interval of predicted value would widen and the accuracy of predictions saw a gradual decline, which was consistent with the conclusion of Xiao-Mei M [25] and Li-Ping R [26].

Table 4 The predictive monthly incidence of AIDS in 2017 based by ARIMA(0,1,1) × (0,1,1)12
Fig. 6

Comparison of ARIMA model prediction and the actual incidence

BP-ANN model

The set of information was divided into training, test and validation data sets in the ARIMA model. The BP-ANN model was established by Matlab 2019 to predict the incidence of HIV/AIDS in China in 2017.

Network architecture

The BP-ANN modeling process has the following three steps:

1) original data was divided into three data sets, namely training, validation and test sets. The training set was used to train models and select the optimal network; the verification set was utilized to monitor the entire network training process; the test set was applied to verify the performance of the selected optimal network model. In network training, training and validation sets are usually selected to enter the network alternately in order to avoid overfitting, which means that established network models explain not only the variation of the observed population but also the fluctuations and errors of individual samples in the population [24].

2) After centralized training, repeated learning, forward and backward propagation of information, and continuous adjustment of network weights, the mean square error (MSE) of validation set would be minimized or reach a predetermined number of iterations [27, 28].

3) As a set of data coming from the same whole with training and verification sets and failing to enter network training, the test set can be used to evaluate established network models to obtain objective and extrapolative effective results.

In the training set of this model, the incidence data of the past three years was used to learn the incidence of the fourth year, such as the incidence rate in January of t1-t3 years was used to estimate that in January of the t4 year, then the incidence rate in the same month of t2-t4 years was used to predict that in the same month of the t5 year, in turn. With such data diversity method could be adopted to better learn and train network models, and avoid overlearning and overfitting.

After dividing data into three sets, network parameters are set up, such as number of network layers, nodes and iterations, the allowable error, and the learning algorithm used.

After the data set has been partitioned, the number of network layers, number of neural nodes, number of iterations, allowable error, learning algorithms and other network parameters of the model should be set before starting training.

To determine the number of network layers. A study by Robert Hecht-Nielsen in 1989 has shown that the feedforward network of a hidden layer can map continuous functions within all closed intervals [29]. A three-layer BP network model can complete any mapping from n to m dimensions. More than two hidden layers should only be considered when learning discontinuous functions. As long as the number of nodes in the hidden layer can be reasonably selected, the BP network model of a hidden layer has also strong nonlinear mapping capability, fast training speed, and good convergence ability. Hence, a three-layer BP network model was selected and a hidden layer was adopted in this study.

To determine the number of neurons (also called nodes) in each layer. The number of nodes in the input layer and the output layer is generally determined according to the data characteristics of the study. In this study, according to the data diversity and the predicted monthly incidence rate, the number of nodes in the input layer is 3, and the number of nodes in the output layer is 1. The number of nodes in the hidden layer has a certain influence on the performance of the neural network model. Too few neural nodes will cause small learning capacity, and failure to completely learn samples and laws of sample storage; Too many neural nodes will cause the network to be bloated, so that the learning speed may slow down and the irregular parts (such as white noise) of sample data may be stored into the network, resulting in poor network performance and generalization ability. At present, the number of nodes in the BP-ANN hidden layer is almost calculated and estimated by the empirical formula. Based on the literature review, this study uses four formulas and two empirical formulas that are the most commonly used to infer the approximate number of neural nodes in the hidden layer, and the formulas are as follows:

$$ \left(\mathrm{Kuarycki}\right)\mathrm{m}=3\mathrm{N} $$
$$ (Maren)\mathrm{m}=\mathrm{M}\left(\mathrm{N}+1\right) $$
$$ Lippmann:m=\sqrt{M}\times N $$
$$ \left( Hecht- Nielsen\right)\mathrm{m}=2\mathrm{N}+1 $$
$$ {\displaystyle \begin{array}{l} Empirical\ Formulas:m=\sqrt{M\times N}+a\\ {}m={\log}_2M\end{array}} $$

where M represents the number of input layer nodes; N represents the number of output layer nodes; “m” represents the number of hidden layer nodes; “a” is the regulation constant with values between 1 and 10. In this study, the number of nodes in hidden layer ranges from 3 to 12.

Select the learning algorithms and structures, initialize the model. Matlab provides 10 (a total of 11) BP neural network model learning algorithms, including Levenberg-Marquardt algorithm (Train-lm), One Step Secant (OSS) algorithm (Trainoss), conjugate direction algorithm (T-trainscg), Polak-Ribiere algorithm (Traincgp), Fletcher-Reeves algorithm (Traincgf), resilient BP algorithm (Trainrp), self-adaptive learning rate algorithm (Traingda and Traingdx), gradient descent with momentum (Traingdm) and batch gradient descent training function (Traingd). Among them, the Levenberg-Marquardt algorithm, the L-M algorithm for short, is the most widely used nonlinear least square algorithm at present because of its fast convergence speed.

In this study, three years of data were randomly selected from the data set of the monthly incidence of AIDS. After the normalization of data by the PRESTD function, estimated from 3 to 12 nodes in the hidden layer and above 11 algorithms were used to combine into the neural network models of 110 structures. Small sample data was input, and each structure was iterated 20 times to calculate their MSEs respectively. The smaller MSE was, the better the fitting effect of the network model would be and the closer the neural network prediction would be to the real value. The combination of the structure and algorithm of the minimum MSE was shown in Table 5. The combined BP neural network model with the minimum MSE = 0.001863 was the optimal model, and the optimal learning algorithm was the L-M algorithm.

Table 5 MSE of 11 BP-ANN algorithms base on 3–12 neurons in the hidden layer

Forecast analysis with BP-ANN

The BP neural network fitting curve for the incidence of HIV/AIDS in 2017 was obtained by inputting the test set into the trained BP-ANN and using the stored black box to operate network models (Fig. 7 and Table 6 of fitted value). By comparing the predicted value with the actual incidence, the fitted value of the BP-ANN model was very close to the actual monthly incidence of AIDS.

Fig. 7

Comparison of BP-ANN model prediction and the actual incidence

Table 6 The predictive monthly incidence of AIDS in 2017 based by BP-ANN

Comparative analysis

This study mainly compared and evaluated the prediction effects of the ARIMA time series model and BP-ANN model of the following three error evaluation indicators. In Table 7, the observed values were compared with the predicted ones in a point-to-point manner. The modeled MSE, MAE and MAPE in the ARIMA model were 0.0020, 0.0301 and 22.4638 respectively. However, three residuals in the BP-ANN model were 0.0019, 0.0129 and 1.2139 respectively.

Table 7 Comparison of the fitting and prediction performance of the two models

When the morbidity of HIV/AIDS from 2004 to 2016 was set as the original data, models were established to forecast the morbidity of AIDS in 2017. The predicted incidence of AIDS was compared with the actual incidence of AIDS in 2017 so as to verify the fitting effects of models. Ultimately, the ARIMA (0,1,1) (0,1,1)12 structure was considered to be the most suitable time series model with white noise testing LB [18] = 13.909, P > 0.05, which meant that the model was effective. In the model, error parameters were MSE = 0.0020, MAE = 0.0301 and MAPE = 22.4638. The selected BP neural network model was seen as the optimal one with the L-M algorithm. In the model, MSE iterated 16 times was 0.0019, MAE was 0.0129 and MAPE was 1.2139. The fitting error of the BP-ANN model was significantly smaller than that of the ARIMA model while its forecasting accuracy was higher than that of the ARIMA model [30,31,32]. It was seen that the BP-ANN model was more effective in predicting the morbidity of AIDS in China.

In Fig. 8, the BP-ANN model had a fit value closer to the true value compared with the ARIMA model. Both prediction methods could be adopted to predict the incidence of AIDS in China. In terms of prediction accuracy, the BP-ANN model would be more suitable. The BP-ANNmodel could better improve forecasting duration than the ARIMA model. In this study, both methods just took into account the temporal variations of time series. However, the BP-ANN model was a nonlinear model, whose prediction accuracy could be enhanced by adjusting more dimensional inputs and development space was larger than that of the ARIMA model.

Fig. 8

Two kinds of models to predict the monthly incidence of AIDS from January 2017 to April 2018 compared with the actual monthly incidence


Monitoring the prevalence of infectious illnesses is of great importance for conventional health education. The prediction of anticipated AIDS cases will not only detect outburst conditions or report the possibility of outburst cases, but also help decision-makers to know about possible future change trends and past and present data [33].

Both ARIMA and BP-ANN models were based on the time series data prediction method with which the time series was extrapolated to the future through special development principles. In the model, morbidity could be predicted as special risk factors were not involved. Without complex transformations or additional alternative variables, auto-correlation, seasonal variations and secular change trends in the ARIMA model could be simply managed through seasonal functions, moving average, auto-regression and difference. As long as the suitable model was established, it would be possible to predict anticipated cases at a given time interval in the future [34].

Both models were capable of predicting the expected cases of AIDS. It was seen that both ARIMA and BP-ANN models could be used to predict the monthly incidence of HIV/AIDS, but the fitting and forecasting effects of the nonlinear BP-ANN model were superior to those of the traditional linear ARIMA model. First, the modeling method of the BP-ANNmodel was simpler than that of the ARIMA model, while it was unnecessary to set up a complicated mathematical model or understand its mathematical structure and the correlation between variables. Second, the ANN was able to compute and deal with data spontaneously through a number of simple units. It was much better to fulfill the works that were involved with pattern recognition. The professional idea was compared with traditional statistics to significantly improve the precision accuracy in neural networks. The ARIMA model might be more suitable for making short-term forecast analysis because of a gradual decline in its long-term prediction effect. As a whole, the nonlinear BP-ANN model forecasting the morbidity of AIDS in China was the most appropriate way for complicated dynamic and nonlinear systems [35]. Therefore, multi-dimensional inputs in the BP neural network would be gradually improved to find out the best model and accurately make predictions. It will be very promising in future [36].


In summary, an agreement was further reached that the BP-ANN model was a suitable way to monitor and predict the change trend and morbidity of AIDS in China. According to the prediction results, more health investments would be made during outburst periods while fewer investments would be made during low-risk periods, which thus improved intervention effect and source scheduling.


Several limitations still exist in this study. First of all, time series analysis was carried out without considering the factors affecting the incidence of AIDS, such as production methods, social environment, epidemic variations and humanities.

Secondly, the research objects were required to remain relatively constant in prediction models during the whole process. Meanwhile, diversified infection channels and disease prevalence would be generated for a variety of people under distinct living conditions. In local places, it was necessary to re-learn and train prediction according to local conditions.

Thirdly, the BP-ANN model under black-box testing would affect the possibility of extrapolation beyond its training information and the fulfillment of subjective initiatives by operators during the process of BP-ANN analysis.

Availability of data and materials

The data that support the findings of this study are available from China’s Statutory Infectious Disease Report of National Health Commission of the People’s Republic of China, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available.



Human Immunodeficiency Virus


Acquired Immune Deficiency Syndrome


Joint United Nations Programme on HIV/AIDS


Estimation and Projection Package


Asian Epidemic Model


Auto-Regressive Integrated Moving Average


Artificial Neural Network


Back Propagation Artificial Neural Network


Auto-Correlation Function


Partial Auto-Correlation Function


Bayesian Information Criterions


Artificial Intelligence


Mean Absolute Percentage Error


Mean Absolute Error


Mean Square Error


  1. 1.

    Gottlieb MS, Schroff R, Schanker HM, et al. Pneumocystis carinii pneumonia and mucosal candidiasis in previously healthy homosexual men: evidence of a new acquired cellular immunodeficiency [J]. N Engl J Med. 1981;305(24):1425–31.

    CAS  Article  Google Scholar 

  2. 2.

    World Health Organization. Global Health Observatory data repository [DB/OL]. Accessed 24 June 2020.

  3. 3.

    Huang M B, Ye L, Liang B Y, et al. Characterizing the HIV/AIDS Epidemic in the United States and China [J]. Int J Environ Res Public Health, 2015, 13(1):30. doi:

  4. 4.

    WHO. Evaluation Report on China HIV/AIDS Epidemic 2011[J]. Chin J AIDS STD. 2012;18(01):1–5.

    Google Scholar 

  5. 5.

    Murray CJ, Ortblad KF, Guinovart C, et al. Global, regional, and national incidence and mortality for HIV, tuberculosis, and malaria during 1990-2013: A systematic analysis for the global burden of disease study 2013.[J]. Lancet. 2014;384(9947):1005–70.

    Article  Google Scholar 

  6. 6.

    Wang LY, Qin QQ, Ding ZW, et al. Current Situation of AIDS epidemic in China [J]. Chin J AIDS STD. 2017;23(04):330–3.

    Google Scholar 

  7. 7.

    Wu Z, Wang Y, Mao Y, et al. The integration of multiple HIV/AIDS projects into a coordinated national programme in China [J]. Bull World Health Organ. 2011;89(3):227.

    Article  Google Scholar 

  8. 8.

    Liu E, Rou K, Mcgoogan JM, et al. Factors associated with mortality of HIV-positive clients receiving methadone maintenance treatment in China.[J]. J Infect Dis. 2013;208(3):442–53.

    CAS  Article  Google Scholar 

  9. 9.

    Wang L, Guo W, Li D, et al. HIV epidemic among drug users in China: 1995–2011[J]. Addiction. 2015;110(Suppl 1(S1)):20.

    Article  Google Scholar 

  10. 10.

    Walker N, Stover J, Stanecki K, et al. The workbook approach to making estimates and projecting future scenarios of HIV/AIDS in countries with low level and concentrated epidemics [J]. Bri J Venereal Dis. 2004;80(suppl 1):i10.

    Google Scholar 

  11. 11.

    Brown T, Le B, Eaton JW, et al. Improvements in prevalence trend fitting and incidence estimation in EPP 2013[J]. Aids. 2014;28(4):S415–25.

    Article  Google Scholar 

  12. 12.

    Stover J, Mckinnon R, Winfrey B. Spectrum: a model platform for linking maternal and child survival interventions with AIDS, family planning and demographic, projections [J]. Int J Epidemiol. 2010;39(Suppl 1):i7.

    Article  Google Scholar 

  13. 13.

    Lim SH, Cheung DH, Guadamuz TE, et al. Latent class analysis of substance use among men who have sex with men in Malaysia: findings from the Asian internet MSM sex survey [J]. Drug Alcohol Depend. 2015;151:31–7.

    Article  Google Scholar 

  14. 14.

    Stover J. Projecting the demographic impact of AIDS and the number of people in need of treatment: updates to the Spectrum projection package [J]. Sexually Transmitted Infections. 2006;82(suppl_3):iii45–50.

    Article  Google Scholar 

  15. 15.

    Tuhuma T, Gideon K, Japhet K, et al. Estimating and projecting HIV prevalence and AIDS deaths in Tanzania using antenatal surveillance data [J]. BMC Public Health. 2006;6(1):120.

    Article  Google Scholar 

  16. 16.

    Sharma SK, Kadhiravan T. Management of the Patient with HIV Disease[J]. Disease A Month. 2008;54(3):162–95.

  17. 17.

    Earnest A, Tan SB, Wildersmith A, et al. Comparing statistical models to predict dengue fever notifications.[J]. Computational and Mathematical Methods in Medicine, ,2012,(2012-3-8). 2012;2012(1):758674.

    Google Scholar 

  18. 18.

    Yolcu U, Egrioglu E, Aladag CH. A new linear & nonlinear artificial neural network model for time series forecasting [J]. Decis Support Syst. 2013;54(3):1340–7.

    Article  Google Scholar 

  19. 19.

    Hornik K, Stinchcombe M, White H. Multilayer feedforward networks are universal approximators.[J]. Neural Netw. 1989;2(5):359–66.

    Article  Google Scholar 

  20. 20.

    Yu HK, Na-Young K, Soon KS, et al. Forecasting the number of human immunodeficiency virus infections in the Korean population using the autoregressive integrated moving average model [J]. Osong Public Health & Research Perspectives. 2013;4(6):358–62.

    Article  Google Scholar 

  21. 21.

    Wang T, Zhou Y, Wang L, et al. Using an autoregressive integrated moving average model to predict the incidence of hemorrhagic fever with renal syndrome in Zibo, China, 2004-2014[J]. Jpn J Infect Dis. 2015;69(4):279–84.

    Article  Google Scholar 

  22. 22.

    Guan P, Huang DS, Zhou BS. Forecasting model for the incidence of hepatitis a based on artificial neural network.[J]. World J Gastroenterol. 2004;10(24):3579–82.

    Article  Google Scholar 

  23. 23.

    Connor JT, Martin RD, Atlas LE. Recurrent neural networks and robust time series prediction[J]. IEEE Trans Neural Netw. 1994;5(2):240–54.

  24. 24.

    Mcguire V, Nelson LM, Koepsell TD, et al. Assessment of occupational exposures in community-based CASE-control studies [j]. Annu Rev Public Health. 1998;19(1):35–53.

    CAS  Article  Google Scholar 

  25. 25.

    Xiao M, Xu Q X, Shi ZX, et al. Application of ARIMA model in predicting monthly incidence of syphilis[J]. Journal of Xi'an Jiaotong University (Medical Sciences). 2018;39(1):131–134,152.

  26. 26.

    Li PR , Ming L , Xin GY , et al. The prediction of the Japanese encephalitis invasion based on the ARIMA model in Guizhou in 2017[J]. Modern Preventive Medicine. 2018;45(08):1349–53.

  27. 27.

    Roman J, Jameel A. Backpropagation and Recurrent Neural Networks in Financial Analysis of Multiple Stock Market Returns [C]// System Sciences, 1996, Proceedings of the Twenty-Ninth Hawaii International Conference on. IEEE, vol. 2; 1996. p. 454–60.

    Google Scholar 

  28. 28.

    Chua CG, Goh ATC. A hybrid Bayesian back-propagation neural network approach to multivariate modelling [J]. Int J Numerical Analytical Methods Geomechanics. 2003;27(27):651–67.

    Article  Google Scholar 

  29. 29.

    Robert HN. Theory of the backpropagation neural network[J]. Proc. 1989 IEEE IJCNN. 1989;1:593–605.

  30. 30.

    Rathnayaka RMKT, Seneviratna D, Jian GW, et al. A hybrid statistical approach for stock market forecasting based on artificial neural network and ARIMA time series models[C].//2015 International Conference on Behavioral, Economic and Socio-cultural Computing (BESC). IEEE. 2015:54–60.

  31. 31.

    Chuang FK, Hung CY, Chang CY, et al. Deploying Arima and artificial neural networks models to predict energy consumption in Taiwan [J]. Sens Lett. 2013.

  32. 32.

    Lewis CD. Industrial and business forecasting methods : a practical guide to exponential smoothing and curve fitting [M]; 1982.

    Google Scholar 

  33. 33.

    Fan YG, Wang J, Su H, et al. Prediction on the number of HIV with models of ARIMA and GM(1,1)[J]. Chin J Control Prev. 2012;12:1100–3.

    Google Scholar 

  34. 34.

    Luo J, Yang S, Zhang Q, Wang L. ARIMA model of time series for forecasting epidemic situation of AIDS [J]. Chongqing Med. 2012;13:1255–9.

    Google Scholar 

  35. 35.

    Jain A, Srinivasulu S. Development of effective and Efficientra- infall -runoff models using integration of deterministic,real coded genetic algorithms, and artificial neural networktechniques. Water Resour Res. 2004;40:W04302.

    Article  Google Scholar 

  36. 36.

    Ran L, Ma N. Comparison of four AIDS epidemic estimation and models [J]. Chin J AIDS STD. 2012;5:347–50.

    Google Scholar 

Download references


We wish to thank the timely help given by The Proofreading Team in improving the language of this article.


This work was supported by the National Natural Science Foundation of China [No.81260450].

Author information




ZL contributed significantly to analysis and manuscript preparation; YL contributed to the conception of the study and helped perform the analysis with constructive discussions. All authors have read and approved the manuscript.

Corresponding author

Correspondence to Yanning Li.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Li, Z., Li, Y. A comparative study on the prediction of the BP artificial neural network model and the ARIMA model in the incidence of AIDS. BMC Med Inform Decis Mak 20, 143 (2020).

Download citation


  • AIDS
  • Prediction
  • BP artificial neural network model
  • ARIMA model