Forecasting Population of Madiun Regency Using ARIMA Method CAUCHY –Jurnal Matematika Murni dan Aplikasi Volume 7(3) (2022), Pages 420-431 p-ISSN: 2086-0382; e-ISSN: 2477-3344 Submitted: May 26, 2022 Reviewed: July 19, 2022 Accepted: July 25, 2022 DOI: http://dx.doi.org/10.18860/ca.v7i3.16156 Forecasting Population of Madiun Regency Using ARIMA Method Yuniar Farida*, Mayandah Farmita, Nurissaidah Ulinnuha, Dian Yuliati Department of Mathematics, Faculty of Science and Technology, Sunan Ampel State Islamic University Surabaya, Indonesia Email: yuniar_farida@uinsby.ac.id ABSTRACT The high population growth of the Madiun Regency can cause population density that can have implications for other problems, both in terms of social, economic, welfare, security, land availability, availability of clean water, and food needs. This study aims to predict the population growth of Madiun Regency using the ARIMA method. The ARIMA method is popular for forecasting time series data, which is reliable because the calculation process is done gradually. The ARIMA method has three models, namely AR (Autoregressive), MA (Moving Average), and ARMA (Autoregressive Moving Average). This study uses annual population data of Madiun Regency from 1983 to 2021 and produces an ARIMA forecasting model (0,2,1) with a MAPE value of 8.42%. This study also showed that from 2022 to 2024 is predicted to increase by 17947 people or 2.39%. The results of this study are expected to be used as information from the Madiun Regency government in anticipating the emergence of problems caused by the population level of Madiun Regency in the future. Keywords: ARIMA; forecasting; population; time series analysis INTRODUCTION One of the most important problems globally is the high population growth in developing countries [1], [2]. Indonesia is listed as one of the five most populous countries in the world. Indonesia ranks fourth after China, India, and the United States and is the most populous Asian continent [3]. According to the 2020 Population Census, conducted in 2020, the population of Indonesia reached 270.20 million. Indonesia has a land area of 1.9 million km2 and a population density of 141 people per km2, with an average annual population growth rate of 1.25% between 2010 and 2020 [4]. In Indonesia, Java ranks first as the most populous island, and East Java is the second-most populous province after West Java, with 40.67 million people. [4]. East Java Province consists of several cities and regencies, one of which is Madiun Regency, whose population growth rate ranks fifth in the 2020 period and experiences population growth every year. In 2015 the population was 676,087 people, while its development in 2019 was 749,070 [5]. Population growth in Madiun Regency is affected by the high birth rate. In 2020, the birth rate in Madiun Regency will be 3928, with a population growth rate of 0.92% [5]. High population growth can cause various problems, such as regional spatial problems, housing, employment, education, economy, and security. In addition, it can also http://dx.doi.org/10.18860/ca.v7i3.16156 mailto:yuniar_farida@uinsby.ac.id Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 421 cause problems in social aspects, welfare, availability of clean water, food needs, and can cause environmental damage [1], [6], [7]. Population density, which can cause problems, needs to be anticipated by making predictions so that various handling strategies can be carried out. Several statistical mathematical-based forecasting models can predict it, including exponential smoothing, moving average, and ARIMA models (Box-Jenkins). Other forecasting models are based on artificial intelligence, such as neural networks, genetic algorithmic, simulated annealing, and classification [8]. Several previous studies by Xu et al., predicted Beijing's main area population using the Long Short Term Memory (LSTM) Model with MAPE 4.35% [9]. Next, a study that indicated the population of residents in East Kalimantan using Exponential Smoothing by Pakpahan, Basani, and Hariani that yielded a MAPE of 14.81% [10]. On the other, there are studies related to the prediction of population prediction using the ARIMA method conducted by Mardiyah et al. in Pasuruan City [11], Nyoni, Mutongi, and Munyaradzi in the Gambia [12], and Nyoni in Zimbabwe obtained MAPE < 3.94% [13]. Based on the previous studies, the authors are interested in the ARIMA method, which results in an excellent level of accuracy in some cases of population forecasting. The ARIMA method is flexible and straightforward in an application, and accurate prediction results for the short term, but the forecasting accuracy for long-term forecasting is not good and will usually tend to be flat for a long time [14], [15]. In applying practice and forecasting the population, ARIMA is also widely used in various case studies, including research by Alabdulrazzaq related to predicting the spread of COVID-19 with MAPE 4.2% [16]. Other research by Swaraj about COVID-19 predictions in India with MAPE 4.7% [17]. Additional research by Guha and Bandyopadhyay predicts the price of gold with a MAPE of 3.25% [18]. Another study by Banerjee related forecasting on the Indian stock market with MAPE 3.33% [19]. Then there was research by Grigonytė and Butkevičiūtė about predicting wind speed in Latvia with MAPE 1% [20]. Based on the explanation above, this study used the ARIMA method to predict the population of the Madiun Regency. This research is expected to provide information for the Madiun Regency Government to take policy steps to minimize and reduce risk due to the high rate of population growth of Madiun Regency. METHODS The Data The data used in this study is data on residents of Madiun Regency from 1983 to 2021, which is taken from the Central Bureau of Statistics of Madiun Regency, from the website https://bit.ly/PendudukKabMadiun [21]. Table 1. The population of Madiun Regency No. Year Population 1. 1983 638586 2. 1984 627467 ⋮ ⋮ ⋮ 38 2020 744350 39 2021 750143 https://bit.ly/PendudukKabMadiun Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 422 ARIMA Box and Jenkins first developed the ARIMA model in the 1970s [22]. The ARIMA is one of the econometric methods used to predict univariate time-series data. Box and Jenkins state that this model does not use independent variables but instead utilizes the information in the circuit to generate pre-predicted values. Therefore, the ARIMA model requires an autocorrelation process in the series. Autocorrelation is the correlation between two observations at different points in a time series. In other words, time series data is self-correlated. Time series models in the ARIMA method include autoregressive (AR), moving average (MA), and autoregressive moving average (ARMA) [23], [24]. The analysis with ARIMA Box Jenkins begins by creating a series of plot periods and plotting the ACF to determine whether the data is mean-stationary or variance- stationary. Differentiation must be done if the data are not stationary to the mean. Otherwise, if the data is not stationary to the variance, a Box-Cox transformation is performed. Repeat the process for the data stationer. After getting stationary data, the next step is to predict the data from ARIMA based on the ACF and PACF plots. Then use Ljung-Box to test the parameters of the test model as well as test the residual hypothesis, which is residual white noise. It can be concluded, there are several stages of forecasting in ARIMA, namely model identification, parameter estimation, diagnostic testing, and prediction. Model Identification When identifying the model on the ARIMA method, the data used must meet the stationary or stability requirements. If the data does not meet the stationery requirements, the data must be stationary for the variance and average (mean) [25]. The transformation equation is as follows [26]. 𝑇(𝑍𝑡) ′ = 𝑍𝑡 𝜆 𝜆 (1) Where: 𝑇(𝑍𝑡) : transformed data value 𝑍𝑡 : i th time data value 𝜆 : the estimated value of transformation parameters The transformed data is determined by the lambda value. For example, the following table shows some commonly used 𝜆 values and associated transformations. Table 2. 𝜆 value and transformation 𝝀 value Transformation -1.0 1 𝑍𝑡⁄ -0.5 1 √𝑍𝑡⁄ 0.0 ln𝑍𝑡 0.5 √𝑍𝑡 1.0 𝑍𝑡 For time-series data that have not satisfied the stationarity of the average, the data must be processed differentially to find the difference between one data and the previous data in sequence. The differencing equation is as follows. 𝑍𝑡 ′ = 𝑍𝑡 − 𝑍𝑡−1 (2) Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 423 𝑍𝑡 ′ is the differentiated data value, where 𝑍𝑡 is the i th time data value. If the data is already stationary, a tentative model of ARIMA (p, d, q) is obtained. Annotation p is a lag that exceeds the significance limit on the Partial Autocorrelation Function (PACF) plot graph, d is the level of differencing performed, q is the lag that crosses the significance limit of the Autocorrelation Function (ACF) plot. Autoregressive is a model in which a dependent variable is influenced by the value of the dependent variable itself because the data used is single. In general, AR is p ordo, with the form 𝐴𝑅(𝑝) as follows [27]. 𝑋𝑡 = ∅0 + ∅1𝑋𝑡−1 + ∅2𝑋𝑡−2 + ⋯+ ∅𝑖𝑋𝑡−𝑖 + 𝛼𝑡 (3) Where: 𝑋𝑡 : time series t 𝑋𝑡−𝑖 : time series t-i 𝛼𝑡 : time error value t ∅0 : constant ∅𝑖 : coefficients of autoregressive Moving Average is a model that measures autocorrelation between the error or residual values. The MA is generally ornate q, with the following form of 𝑀𝐴(𝑞) [28]. 𝑋𝑡 = 𝑒𝑡 − 𝜃1𝛼𝑡−1 − 𝜃2𝛼𝑡−2 − ⋯− 𝜃𝑖𝛼𝑡−𝑖 (4) Where: 𝛼𝑡−𝑖 : time error value t-i 𝜃𝑖 : coefficient of moving average Autoregressive Moving Average or ARMA (p, q), with the following general equations [29]. 𝑋𝑡 = ∅0 + ∅1𝑋𝑡−1 + ⋯+ ∅𝑖𝑋𝑡−𝑖 + 𝛼𝑡 − 𝜃1𝛼𝑡−1 − ⋯− 𝜃𝑖𝛼𝑡−𝑖 (5) Where: 𝑋𝑡 : stationary time series Autoregressive Integrated Moving Average data used must be stationary. ARIMA's general statement is as follows [30]. ∅0(𝐵)(1 − 𝐵) 𝑑𝑍𝑡 = 𝜃0 + 𝜃𝑞(𝐵)𝑎𝑡 (6) Where: ∅0 : Autoregressive process 𝜃𝑞 : moving average process (1 − 𝐵)𝑑 : differentiating operator 𝑑 : differencing parameter 𝐵 : step-back operator 𝑍𝑡 : deviations from the average process Parameter Estimation Tentative model determination requires several estimation stages through model feasibility tests to find the best model. The significance test hypothesis is as follows [29]. 𝐻0:∅ = 0 (indicates parameters are not yet significant) 𝐻1:∅ ≠ 0 (Shows the parameters are significant) Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 424 𝑡𝑐𝑜𝑢𝑛𝑡 = 𝜃 𝑆𝐸(𝜃𝑗) (7) Where: 𝜃 : estimation of autoregressive model parameters and moving averages 𝑆𝐸(𝜃𝑗) : standard errors Diagnostic Test Diagnostic tests are used to determine whether or not the model is the best. A good model, where the residual results of the white noise assumption test using the Ljung-Box test are as follows [6], [29]. 𝑄 = 𝑛(𝑛 + 2)∑ �̂�𝑘 2 (𝑛 − 𝑘) 𝑖 𝑘 (8) Where: �̂�𝑘 : lag autocorrelation value k Q : Ljung-Box test 𝑘 : lag time Prediction Accuracy Value The results produced by the ARIMA model are measured in terms of forecast accuracy. Each method has a MAPE (Mean Absolute Percentage Error) error value that can be used to calculate the error value with the following formula [16], [31]. 𝑀𝐴𝑃𝐸 = ∑ |𝑃𝐸𝑡| 𝑛 𝑡=1 𝑛 (9) With 𝑃𝐸𝑡 = 𝑒𝑡 𝑍𝑡 × 100 (8) Where: 𝑃𝐸𝑡 : percentage of errors at t 𝑒𝑡 : t-time error value 𝑍𝑡 : actual data of t-time The quality of the prediction can be shown by the MAPE value, which can be interpreted into four categories, namely excellent (MAPE < 10%), good (MAPE 11% - 20%), good enough (MAPE 21% - 50%), and not good (MAPE > 50%). RESULTS AND DISCUSSION Based on Table 1, a time series plot is performed to determine the ARIMA model and identify the stationarity of the data. Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 425 Figure 1. Plot data on the population of Madiun Regency Based on Figure 1, the plot data shows an uptrend (positive). The data is not stationary because in 2019 there was an increase seen from the previous year's difference of 67,676 people. If there is no increase or decrease invariance and average, the data is stationary. Figure 2. Plot Box-Cox transformation Figure 2 shows that the lambda value is equal to 1, the data can be said to be stationary in variance. Stationary data on the average can be seen from the ACF plot and time series plot. The data does not yet have a fixed pattern. Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 426 Figure 3. ACF plot of Madiun Regency population From the plot Figure 3, it appears that the lag-lag is falling slowly. The plot time series data also does not have a fixed pattern, so the data is not stationary against the average. As a result, it is necessary to do a further transformation process through differencing so that the data is stationary. Figure 4. ACF plot after differencing Figure 4 shows that the data is stationary against the mean. If the data is stationary, the next step is to plot the autocorrelation function (PACF). Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 427 Figure 5. PACF plot after differencing Based on Figures 4 and 5 shows that the plot does not have an autocorrelation on the model, so the values 𝑀𝐴(𝑞) = 0 and 𝐴𝑅(𝑝) = 0, then obtained the tentative model ARIMA (0,2,0). The model is a random walk where the autocorrelation coefficient is equal to 1, so the tentative models of ARIMA are ARIMA models (1,2,0), (0,2,1), and (1,2,1). A significance test and a residual white noise test were carried out to choose the model used in the prediction. Test the significance of the parameters by knowing the p- value. If the p-value is less than 0.05, then the model is significant. The results of the ARIMA model's tentative significance test (1,2,0), (0,2,1), and (1,2,1) are as follows. Table 3. Significance Test Results Model Parameters Coef SE Coef T-Value P-Value ARIMA (1,2,0) AR (1) -0.580 0.139 -4.19 0.000 ARIMA (0,2,1) MA (1) 0.950 0.131 7.26 0.000 ARIMA (1,2,1) AR (1) -0.221 0.180 -1.23 0.228 MA (1) 0.946 0.122 7.73 0.000 Table 3 shows that the ARIMA models (1,2,1) are not significant because p-values are more than 0.05, ARIMA models (1,2,0) and (0,2,1) are significant because p-values are less than 0.05. After conducting a parameter significance test, it is necessary to perform a residual white noise test to determine which model to use for prediction in performing residual tests using Ljung-Box. If the p-value is more than 0.05, the model meets the white noise requirement. Ljung-Box test results for ARIMA models (1,2,0), (0,2,1), and ARIMA models (1,2,1). Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 428 Table 4. Ljung-Box Test Results Model Lag Chi-square DF P-Value ARIMA (1,1,1) 12 5.34 10 0.867 24 11.17 22 0.972 36 23.15 34 0.920 ARIMA (2,1,1) 12 1.83 10 0.998 24 7.23 24 0.999 36 11.94 36 1.000 ARIMA (0,1,2) 12 0.53 9 1.000 24 6.15 21 0.999 36 11.41 33 1.000 Table 4, after the Ljung-Box test, shows that the ARIMA model (1,2,0), (0,2,1), and ARIMA model (1,2,1) are white noise because the p-value is more than 0.05. After performing the white noise test, determine the MAPE value of the ARIMA provisional model. Table 5. MAPE value Model MAPE ARIMA (1,2,0) 12.45 ARIMA (0,2,1) 8.42% ARIMA (1,2,1) 8.92% Based on Table 3, Table 4, and Table 5 of the ARIMA model, whose parameters are significant and meet the assumption of white noise, the ARIMA model (0,2,1) has the smallest MAPE value. After choosing the best model, then predicted the number of residents of Madiun Regency. Figure 6. Actual Data Plot and Forecasting Figure 6 showed that there is not much difference between the actual data and the forecast and obtained MAPE < 10%, which means that the prediction model is already very good. The prediction of the number of residents of Madiun regency for the next three years (2022 to 2024) is presented in table 6 below: Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 429 Table 6. Population Prediction Year Population 2022 758561 2023 767346 2024 776508 Table 6 shows that the Madiun Regency population from 2022 to 2024 is predicted to increase by 17947 people or 2.39%. The results of these predictions show an uptrend every year. Population growth, if followed by an increase in the quality of human resources, will become a regional potential for development. On the other hand, if the increase in population is not accompanied by good quality human resources, it will become a burden for regional development [32]. Population data prediction is needed in the planning and evaluating of human- oriented development as the primary target because the population is both an object and a subject of the action. The object’s function means the population as a target, and the people carry out the development mark. The function of the issue means that the people are the sole actor in action. The two functions are expected to go hand in hand and line integrally [32]. CONCLUSIONS Based on the research results on the prediction of the population of Madiun Regency using the ARIMA method, it can be concluded that the best model for predicting the number of residents of Madiun Regency is ARIMA (0,2,1) with a MAPE of 8.42%. The predicted number of residents of Madiun Regency in 2022 amounted to 758561 people. The selection of the ARIMA method in this research for forecasting the number of residents in Madiun is very appropriate because it produces an error value of less than 10%. This method can be applied to similar case studies, especially the case of forecasting the number of residents in other areas. REFERENCES [1] C. Christiani, P. Tedjo, and B. Martono, “Analisis Dampak Kepadatan Penduduk Terhadap Kualitas Hidup Masyarakat Provinsi Jawa Tengah,” Ilmiah, vol. 3, no. 1, pp. 102–114, 2014. [2] J. Dai and S. Chen, “The application of ARIMA model in forecasting population data,” J. Phys. Conf. Ser., vol. 1324, no. 1, 2019, doi: 10.1088/1742-6596/1324/1/012100. [3] F. Fejriani, M. Hendrawansyah, L. Muharni, S. F. Handayani, and Syaharuddin, “Forecasting Peningkatan Jumlah Penduduk Berdasarkan Jenis Kelamin menggunakan Metode Arima,” J. Kajian, Penelit. dan Pengemb. Pendidik., vol. 8, no. 1 April, pp. 27–36, 2020. [4] Direktorat Statistik Kependudukan dan Ketenagakerjaan, Potret Sensus Penduduk 2020 Menuju Satu Data Kependudukan Indonesia. Jakarta: BPS RI, 2021. [5] B. K. Madiun, Kabupaten Madiun dalam Angka 2021. Madiun: BPS Kabupaten Madiun, 2021. [6] Haslina, Hasmah, K. W. Fitriani, M. Asbar, and Asrirawan, “Penerapan Metode ARIMA (Autoregressive Integrated Moving Average) Box Jenkins Untuk Memprediksi Pertambahan Jumlah Penduduk Transmigran (Jawa dan Bali) di Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 430 Kecamatan Sukamaju, Kabupaten Luwu Utara Propinsi Sulawesi Selatan,” Dinamika, vol. 9, no. 1, pp. 55–67, 2018. [7] H. Yoshikura, “Negative impacts of large population size and high population density on the progress of measles elimination,” Jpn. J. Infect. Dis., vol. 65, no. 5, pp. 450–454, 2012, doi: 10.7883/yoken.65.450. [8] N. L. A. K. Yuniastari and I. W. W. Wirawan, “Peramalan Permintaan Produk Perak Menggunakan Metode Simple Moving Average Dan Single Exponential Smoothing,” Sist. dan Inform., vol. 9, no. 1, pp. 97–106, 2016. [9] Z. Xu, J. Li, Z. Lv, Y. Wang, L. Fu, and X. Wang, “A graph spatial-temporal model for predicting population density of key areas,” Comput. Electr. Eng., vol. 93, no. May, p. 107235, 2021, doi: 10.1016/j.compeleceng.2021.107235. [10] H. S. Pakpahan, Y. Basani, and R. R. Hariani, “Prediksi Jumlah Penduduk Miskin Kalimantan Timur Menggunakan Single dan Double Exponential Smoothing,” Inform. Mulawarman J. Ilm. Ilmu Komput., vol. 15, no. 1, pp. 47–51, 2020. [11] I. Mardiyah, W. D. Utami, D. C. R. Novitasari, M. Hafiyusholeh, and D. Sulistiyawati, “ANALISIS PREDIKSI JUMLAH PENDUDUK DI KOTA PASURUAN MENGGUNAKAN METODE ARIMA,” Ilmu Mat. dan Terap., vol. 15, no. 3, pp. 525–534, 2021. [12] T. Nyoni, C. Mutongi, and N. Munyaradzi, “Population dynamics in Gambia: An ARIMA approach,” Munich Pers. RePEc Arch., 2019. [13] T. Nyoni, “The population question in Zimbabwe : reliable projections from the Box - jenkins ARIMA approach,” Munich Pers. RePEc Arch., pp. 0–15, 2019. [14] M. S. K. Abhilash, A. Thakur, D. Gupta, and B. Sreevidya, “Time series analysis of air pollution in bengaluru using ARIMA model,” Adv. Intell. Syst. Comput., vol. 696, pp. 413–426, 2018, doi: 10.1007/978-981-10-7386-1_36. [15] Nurviana, R. P. Sari, U. Nabilla, and T. Talib, “Forecasting Rice Paddy Production in Aceh Using ARIMA and Exponential Smoothing Models,” CAUCHY, vol. 7, no. 2, pp. 281–292, 2022. [16] H. Alabdulrazzaq, M. N. Alenezi, Y. Rawajfih, B. A. Alghannam, A. A. Al-Hassan, and F. S. Al-Anzi, “On the accuracy of ARIMA based prediction of COVID-19 spread,” Results Phys., vol. 27, p. 104509, 2021, doi: 10.1016/j.rinp.2021.104509. [17] A. Swaraj, K. Verma, A. Kaur, G. Singh, A. Kumar, and L. Melo de Sales, “Implementation of stacking based ARIMA model for prediction of Covid-19 cases in India,” J. Biomed. Inform., vol. 121, no. August 2020, p. 103887, 2021, doi: 10.1016/j.jbi.2021.103887. [18] B. Guha and G. Bandyopadhyay, “Gold Price Forecasting Using ARIMA Model,” J. Adv. Manag. Sci., no. March, pp. 117–121, 2016, doi: 10.12720/joams.4.2.117-121. [19] D. Banerjee, “Forecasting of Indian stock market using time-series ARIMA model,” Int. Conf. Bus. Inf. Manag. ICBIM 2014, pp. 131–135, 2014, doi: 10.1109/ICBIM.2014.6970973. [20] E. Grigonytė and E. Butkevičiūtė, “Short-term wind speed forecasting using ARIMA model,” Energetika, vol. 62, no. 1–2, pp. 45–55, 2016, doi: 10.6001/energetika.v62i1-2.3313. [21] BPS Kabupaten Madiun, Kabupaten Madiun Dalam Angka Madiun Regency In Figures 2021. 2021. [22] S. Ozturk and F. Ozturk, “Forecasting Energy Consumption of Turkey by Arima Model,” J. Asian Sci. Res., vol. 8, no. 2, pp. 52–60, 2018, doi: 10.18488/journal.2.2018.82.52.60. [23] J. Sun, “Forecasting COVID-19 pandemic in Alberta, Canada using modified ARIMA models,” Comput. Methods Programs Biomed. Updat., vol. 1, no. September, p. Population Forecasting of Madiun Regency Using ARIMA Method Yuniar Farida 431 100029, 2021, doi: 10.1016/j.cmpbup.2021.100029. [24] C. B. A. Satrio, W. Darmawan, B. U. Nadia, and N. Hanafiah, “Time series analysis and forecasting of coronavirus disease in Indonesia using ARIMA model and PROPHET,” Procedia Comput. Sci., vol. 179, no. 2020, pp. 524–532, 2021, doi: 10.1016/j.procs.2021.01.036. [25] D. Didiharyono and M. Syukri, “Forecasting with arima model in anticipating open unemployment rates in south sulawesi,” Int. J. Sci. Technol. Res., vol. 9, no. 3, pp. 3838–3841, 2020. [26] D. S. Domingos, J. F. L. de Oliveira, and P. S. G. de Mattos Neto, “An intelligent hybridization of ARIMA with machine learning models for time series forecasting,” Knowledge-Based Syst., vol. 175, pp. 72–86, 2019, doi: 10.1016/j.knosys.2019.03.011. [27] F. A. Chyon, M. N. H. Suman, M. R. I. Fahim, and M. S. Ahmmed, “Time Series Analysis and Predicting COVID-19 Affected Patients by ARIMA Model Using Machine Learning,” J. Virol. Methods, vol. 301, no. December 2021, p. 114433, 2021, doi: 10.1016/j.jviromet.2021.114433. [28] N. Ulinnuha and Y. Farida, “Prediksi Cuaca Kota Surabaya Menggunakan Autoregressive Integrated Moving Average (Arima) Box Jenkins dan Kalman Filter,” J. Mat. “MANTIK,” vol. 4, no. 1, pp. 59–67, 2018, doi: 10.15642/mantik.2018.4.1.59-67. [29] T. Yunita, “Peramalan Jumlah Penggunaan Kuota Internet Menggunakan Metode Autoregressive Integrated Moving Average ( ARIMA ),” J. Math. Theory Appl., vol. 1, no. 2, pp. 16–22, 2019. [30] L. Wulandari, Y. Farida, A. Fanani, and M. Syai’in, “Optimization of Autoregressive Integrated Moving Average (ARIMA) for Forecasting Indonesia Sharia Stock of Index (ISSI) using Kalman Filter,” pp. 295–303, 2020, doi: 10.5220/0008906902950303. [31] M. B. S. Junianto, “Fuzzy Inference System Mamdani dan the Mean Absolute Percentage Error (MAPE) untuk Prediksi Permintaan Dompet Pulsa pada XL Axiata Depok,” J. Inform. Univ. Pamulang, vol. 2, no. 2, p. 97, 2017, doi: 10.32493/informatika.v2i2.1511. [32] P. K. Madiun, Data Demografi, Ekonomi dan Sosial Budaya Kota Madiun 2017. Madiun: Pemerintah Kota Madiun, 2017.