A Fuzzy ARIMA Correction Model for Transport Volume Forecast

In recent years, few water transportation forecasting studies conduct relative to transportation forecasting. As a neglected area, the inland waterway volume prediction is an important indicator for investment management and government policymaking. Considering the time-series forecasting, some researchers try to narrow the predicted value interval. However, certain limitations detract from their popularity. For instance, if the prediction length is more than ten, the result would not be acceptable. Therefore, we propose a hybrid model that combines both of their unique properties’ advantages to provide more accurate traﬃc volume forecasts. Also, the forecasting process will be more straightforward. The empirical results present the proposed model and improve the long-term predictive accuracy at waterway traﬃc volume.


Introduction
A long-term prediction usually hardly predicts more than ten periods and follows the time pass, and the accumulated error will make the result barely acceptable. On the contrary, how to combine different influence factors and historical time-series data is also an interesting part for researchers. In this case, some scholars propose Fuzzy Time Series (FTS), Autoregressive Integrated Moving Average (ARIMA), Artificial Neural Networks (ANNs), etc., to deal with this problem.
us, this research focused on keeping the long-term prediction result in a narrow interval value range. In other words, the study attempts to combine different methods for exploiting its strong points to increase forecasting accuracy rates. is work will benefit the long-term forecasting area. Usually, the inland waterway volume (VOL) has a long transfer period, and it also contributes to the inland waterway investment management and government policymaking.

Related Works
When the fuzzy prediction methods were first introduced to deal with the forecasting issue, the historical data were linguistic values [1][2][3]. Since then, the researchers have generally accepted the standard form of Fuzzy Time Series (FTS) [4]. After that, numerous research studies in the forecasting area work on the FTS method. Some of them forecasted the pollution concentrations [5]. Duru adopted a fuzzy integrated logical forecasting model for dry bulk shipping index forecasting [6]. Chang and Chen proposed the temperature prediction method based on fuzzy clustering and fuzzy rule interpolation techniques [7]. Chen and Dang mentioned a method to construct with variable spreads [8]. However, both of them used fuzzy rule techniques. Aladag and Basaran proposed an approach that used feedforward neural networks to gain the high-order fuzzy time data [9]. Huang proposed a fuzzy time-series forecasting, which used a multivariable heuristic function to improve forecasting results by integrating various univariate models [10]. Karnik and Mendel designed a type-2 fuzzy logic system to deal with the Mackey-Grass chaotic timeseries forecasting problem [11]. Yu and Huang applied neural networks to fuzzy time series forecasting and proposed bivariate models to improve forecasting quality [12]. Wang and Hsu an improved fuzzy time series model for forecasting short-term time-series data [13]. Lau proposed a method to forecast energy consumption change based on the fuzzy logic method [14].
Another common approach for prediction is Autoregressive Integrated Moving Average (ARIMA) model, which was a typical time-series prediction method. Cullinane applied the Box-Jenkins approach to forecast the shipping market and discuss the different forecasting stages in shipping [15]. Fanoodi et al. predicted the blood platelet demands based on Artificial Neural Networks (ANNs) and ARIMA to reduce the supply chain's uncertainty [16]. Pavlyuk devoted to a multivariate model including ARIMA and Vector Autoregressive Moving Average (VARMA), Error Correction Vector Autoregressive Moving Average (EC-VARMA), Space-time ARMA (STARMA), and Multivariate Autoregressive Space State (MARSS) models to explain short-term traffic flow forecasting [17]. Li and Parson introduced a nonlinear approach in neural networks and evaluated its short and long period forecasting performance relative to ARMA models [18]. Cullinane et al. forecaste the Baltic Freight Index (BFI) and investigate the impact of a change in the BFI composition in 1993, although the results indicated that the modification did not significantly affect the behavior of the BFI [19]. Wong generated the short-and long-term predictions of the Baltic Dry Index (BDI) and used Fuzzy heuristic modeling with Grey System and ARIMA, and the results indicated that the ARIMA has a better performance in the long-term forecast [20].
However, the ARIMA prediction interval was too broad to be precise for the forecasting purpose. In this case, some researchers combine Fuzzy Time Series and ARIMA model to solve this issue. Duru et al. presented a bivariate long-term fuzzy inference system for time-series forecasting tasks in the freight market by combining FTS and ARIMA models [21]. Wong et al. proposed a traditional ARIMA model and the Fuzzy Time Series Method for forecasting the amount of Taiwan export [22]. Torbat et al. proposed a fuzzy autoregressive integrated moving average model as an improved ARIMA version to narrow the forecast interval [23]. Li and Hu raised three examples to test the Neuro-Fuzzy System ARIMA (NFS-ARIMA) model for nonlinear problems forecasting ability [24]. Tseng integrated the time-series ARIMA model and fuzzy regression model as a Fuzzy ARIMA (FARIMA) model to forecast NT dollars' exchange rate of US dollars [25]. Tseng and Tseng propose a Fuzzy Seasonal ARIMA (FSARIMA) forecasting model to predict the total production value of the Taiwan machinery industry [26].
Compared with the previous works, the contributions of this paper are as follows.
rough original multivariate information, transform into integrated information by constructing a fuzzy membership set. erefore, Correction Fuzzy Time Series Autoregressive Integrated Moving Average (CFTSARIMA) method dramatically simplifies calculating and efficiently processing large amounts of data.
Generally, the FTSARIMA and FMLR models provide minimum widths' intervals for a given confidence degree. However, forecasting will significantly deviate from the real data in the long term. In this case, we propose a CFTSARIMA method, and the result could be acceptable. In equation (3), the h value level represents the value range widths' influences. However, when the h value grows, the prediction interval's tendency is not changing, instead the space between Upper Control Limit (UCL) and Lower Control Limit (LCL) is increased. erefore, it was necessary to recalculate a new center of the fuzzy numerical value a i and the width c i which was around the center fuzzy number.
Compared to railways and road deliveries, few research pieces focus on inland waterway transportations; however, precise long-term forecasting could benefit the economic aspect. e example demonstrated provides both the bestand worst-possible situations to decision-maker consideration.

Problem Statement
In this section, we first briefly review the coupled Fuzzy Clustering, Fuzzy Multiple Linear Regression method, and Fuzzy Time Series Autoregressive Integrated Moving Average method for long-term forecasting. However, the results diverge from observing data for solving the issue. We apply the Correction Fuzzy Time Series Autoregressive Integrated Moving Average (CFTSARIMA) method to solve the projection bias problems specifically. Table 1 summarizes the code of nomenclature for fuzzy ARIMA prediction and the correspondence between the symbol and the definition. [7,27]). Different time-series groups have been chosen to participate in the prediction operation. For reducing the noise in the forecasting information, a fuzzy clustering method is useful.

Fuzzy Clustering (See
To measure the cluster similarity, the Membership formula for different fuzzy sets is as follows. When the predictors are positively correlated with the sample of eigenvalues, r ij � (x ij − x j min )/(x j max − x j min ).
When the predictors are negatively correlated with the sample of eigenvalues, in which x ij are fuzzy sets, where x i , i � 1, 2, . . . , j, is the ith influence factor and j is the number of variable collected data sets. Let elements x i and x j belong to set r ij . en, the correlation coefficient ρ i is the similarity between two fuzzy clusters y ij and x ij which are computed through where y ij , i � 1, 2, . . . , j, is the ith predicted output, j is the number of data sets, And y i and x i are the average value of y i and x i , respectively. In Table 2, the time series are divided into different groups, and Table 3 selects the high similarity data from fuzzy clusters and compares them to y. [8,28]). e general form of a fuzzy multiple linear regression model can be expressed as follows:

Fuzzy Multiple Linear Regression (See
2 Mathematical Problems in Engineering where y i , i � 1, 2, . . . , j, is the ith predicted output, j is the number of datasets, and x j , j � 0, 1, 2, . . . , j, is the jth variable of the collected data. e number of independent variables A j is the fuzzy coefficient of the jth independent variable. e criterion of minimizing the total vagueness S is defined as the sum of individual spreads of the fuzzy parameters of the model: e membership value of each observation y i is considering the condition that should be greater than an imposed threshold h, h ∈ [0, 1]. e choice of the h level value influences the widths c j of the fuzzy parameters. is criterion simply expresses the fact that the fuzzy output of the model should "cover" all the data points y 1 , y 2 , . . . , y k to a certain h level. a j denotes a center of a fuzzy parameter, and c j shows the fuzziness of its parameter. e fuzzy coefficients can be expressed as In the following, take A 1 , A 2 , . . . , A j in equation (2), and the result is as follows: [23,25,26,29]). e basis of this technique is the Autoregressive Integrated Moving Average ARIMA (p, d, q). p is the order of autocorrelation, d is the number of differencing to achieve stationarity, and q is the order of the moving average. e mathematic formula (5) is given as

FTSARIMA (See
where y t is a time series, and let y denote the differenced (stationarized) version of x, e.g., y t � X t − X t− 1 . K is a constant, the random variables ε t and ε t− i are error/disturbance terms, and φ i and θ i present correlative parameters. e formula gives the mathematic representation of an ARMA (p, d, q) model. e index t refers to the number of nonfuzzy data used in constructing the model. Transfer the fuzzy regression parameters problem as a linear programming problem [30].
e fuzzy model interval is as  Series Table 3: e values of y corresponding to fuzzy base clustering ρ i .
where a i is the center of the fuzzy numbers and c i presents the extension value. Despite the advantages of the fuzzy time series autoregressive integrated moving average models (FTSARIMA), forecasted intervals are too wide when data are fluctuations or outliers in the datasets. Besides, the results also include the maximum and minimum value of y s . e result is an interval between y s(UCL) and y s(LCL) , obtained by the FTSARIMA and a confidence degree of a.
3.4. Combined Model. As mentioned above, the weight factors could participate in the combined model and the variance analysis, which presents the impact of different factors through data analysis. Variance decomposition here is used as a weight to help combine the FMLR and FTSARIMA results. e new upper bound is y ucl and the new lower bound is y lcl : where α is the FMLR's weight and β is the FTSARIMA's weight. rough equation (3), Subject to where α ′ � (α 1 ′ , α 2 ′ , . . . , α j ′ ) and c ′ � (c 1 ′ , c 2 ′ , . . . , c j ′ ) are vectors of unknown variables. An interval with new width is obtained by the FTSARIMA and FMLR. is combined model considers the history data and other influence factors.
e results of this model for long-term forecasting are more accurate than the others.

The Empirical Results and Validation
In general, the Mississippi shipping volume y i appears to be influenced mainly by the local urbanization rate x i1 , gross domestic product x i2 , rate of industrialization x i3 , proportion of railway transportation x i4 , proportion of road transportation x i5 , proportion of waterway transportation x i6 , and rate of agriculture x i7 . e information used consists of 47 annual observations of Mississippi River shipping volume from 1970 to 2016, 30 observations are initially used to formulate the model, and the last 17 observations are used to verify the results (see Supplementary Material for data analysis). e research step is shown in Figure 1. From Figure 1, there were four steps to process the FTSARIMA correction model, and the details are as follows.

U.S. Mississippi Inland Waterway Influencing Factors' Fuzzy Clustering.
e purpose of using the fuzzy clustering method is to reduce the computational complexity and disturbance term, and a simple algorithm is introduced.
en, an equation is proposed to generate fuzzy membership matrices, which are the basis for the final clustering results. Finally, the closer the relationship between the sample and the factor, the greater the influence of the predictive factor on the phenomenon. erefore, we can consider removing some effect factors. e similarity of fuzzy clusters is computed through equation (1). Let x ij be fuzzy sets defined on x i and x j with membership functions r ij . Let an element x i ∈ X and x j ∈ X belong to set x ij .
Four different groups are selected and divided from the similarity coefficient which is greater than or equal to 0.9 (see Table 4).
Group 1: x 1 Group 2: x 2 , x 3 , and x 7 Group 3: x 4 and x 5 Group 4: x 6 rough the analysis of pertinence between y and x 1 , specially x 4 , x 5 , and x 6 are less correlated with the y set (see Table 5). erefore, group1 and group2 are retained to rearrange as time series g 1 and g 2 .

Fuzzy Multiple Linear Regression.
is fuzzy multiple linear regression model (FMLR) is applied to the Mississippi River's shipping volume. e input data are obtained from the fuzzy clusters and the output data (see Figure 2). In the column x i , two clustering categories are indicated: group 1 and group 2. From these data, the fuzzy linear system is y i � A 0 + A 1 g 1 + A 2 g 2 and the fitting model for the data is given. e results of fuzzy parameters A 1 and A 2 are given as (a 1 � 0.24, c 1 � 0; a 2 � 3.44, c 2 � 0.91), where h � 0, a j denotes a center of a fuzzy parameter, and c j shows the fuzziness of its parameter: y r � 2.23 + <0.24, 0.00 > g 1 + <3.44, 0.91 > g 2 .
In this case, the fuzziness of A 0 is 2.23. A 1 and A 2 are positive values and rely on the correlations variables g 1 and 4 Mathematical Problems in Engineering e conventional regression model's confidential interval seems to estimate the upper and lower observation errors' limits. Setting (a 1 , a 2 ) � (0.24, 3.44) and h � 0. e following linear interval model has shown its result in Figure 2.
Following the FMLR 17 predictions, see Figure 3, the actual value is lower than estimated, which means the FMLR method for shipping volume forecasting still needs improvement.

Fuzzy Time Series ARIMA Building the ARIMA Model.
e time-series data were preprocessed using the first-order regular differencing to stabilize the variance and remove the growth trend. e derived model is ARIMA (1, 1, 1), and the equation is FML result upper control limit y set lower control limit y set ARIMA prediction Fuzzy time series ARIMA prediction FSARIMA result upper control limit y set lower control limit y set

Mathematical Problems in Engineering
After determining the minimal fuzziness equation (3), the FTSARIMA model is Setting (a 1 , a 2 ) � (0.536, 1) and h � 0, the linear interval model is obtained, and the results are demonstrated in Figure 4.
e actual values locate under the prediction interval value set (see Figure 5) below even the most pessimistic forecasts. To improve the accuracy of model prediction in the long term, the influence of the other factors x i is necessary to be discussed.

Correction Fuzzy Time Series ARIMA Model.
From FMLR and FTSARIMA, two different estimated equations (10) and (13) are obtained. To restructure y ucl and y lcl , the variance decomposition is introduced to assign a weight to them.
Summing up the results of the variance decomposition, the 10th period (see Table 6) could be considered as the weight value that 0.88 belong to y s and (1-0.88), is the allocation to y r . According to the FMLR model and FTSARIMA method, there are two different upper control limit (UCL) and lower control limit (LCL) lines. Based on the table weights 0.88 and 0.12, respectively, the new series are as follows: in which y s is the result of the FTSARIMA model and y r are the FMLR's model result. From equation (13), we can obtain equation (14), and Figure 6 reveals the tendency: By substituting (a 1 , a 2 ) � (0.536, 1), y ucl , and y lcl into equation (9), the new c value has been determined, and the following linear interval model is given in equation (14) and plotted in Figure 7: Usually, three metrics are used to evaluate prediction intervals: coverage rate, calibration, and sharpness [31,32]. e coverage refers to the statistical consistency between the forecasts and the observations, and it measures how many observations are inside the prediction interval. e properties of sharpness and resolution refer to the concentration of the predictive distribution, or how wide and variable are the intervals, and refer uniquely to the forecasts [33].    Compared to FMLR, FTSARIMA, and CFTSARIMA with 17 length predictions, only the CFTSARIMA conclusion has covered half of the observed values, and other prediction intervals are much higher than observed values. On the contrary, CFTSARIMA's prediction intervals are wide-ranging than FMLR and FTSARIMA, and the results are narrow than the traditional ARIMA method.

Results and Discussion
eses examples' empirical results expose that the possible interval of the fuzzy ARIMA is narrower than the interval of classical ARIMA. e fuzzy ARIMA usually is indicated as formulating the model. e output is fuzziness, leading to the assumption of white noise; in this case, FTSARIMA requires fewer observations than ARIMA.
In Table 7, UCL and LCL denote the upper control limit and the lower control limit and PV indicates the predicted values from different forecasting methods. e fuzzy linear regression (FMLR) model is based on the probability distribution of statistics. e relationship between them can only manage the input and output information; therefore, massive information is required. is method is conducive to observing information with trend growth, and the short-term prediction accuracy is only acceptable. On the contrary, long-term forecasts still need to be advanced. Table 7 is compared with four different results. e FTSARIMA is learned from the Fuzzy Linear Regression method and the ARIMA model; compared to both of them, the prediction interval has significantly narrowed. A four-or five-period forecasting result can be acceptable. If the prediction period is more than five, the actual value will be under the prediction interval; neglected other influence factors may cause this.
According to Figure 7 and Table 7, CFTSARIMA showed that the prediction interval contains the most observed value than FMLR and FTSARIMA. Although this method has also enlarged the range of forecasting value, the LCL line keeps a steady trend, and the prediction interval is still narrow than others.
To compare the proposed method's performance with the existing techniques, we apply the proposed approach to long-time forecasting by using 30 observations and 7 influence indexes. Table 7 shows a comparison of the upper control limit (UCL), lower control limit (LCL), and predicted values (PV) from different methods. e proposed method comes with prediction intervals that have significantly narrowed, and the predicted values are more precise than the existing methods. In other words, the proposed Mathematical Problems in Engineering method produces better forecasting results than the strategies presented in [8,23,25,26,28,29] for long-time forecasting.

Conclusions
In this study, the primary objective is to find an accurate prediction method for long-term waterway traffic volume prediction, and the proposed algorithm consists of four parts.
In the first part, a fuzzy cluster reduces the algorithm design difficulty, and nonessential double-counting raises the counting yield. After cluster, seven influence factors are classified into 4 groups, and only two of them satisfy the criteria.
According to the fuzzy clustering, two group influence factors are converted to two different time series sets. A fuzzy multiple linear regression is used to estimate the upper bound y r(UCL) and lower bound y r(LCL) .
In the third part, a typically fuzzy time series ARIMA method is adopted to define the center of fuzzy number a i and the width around the center of the fuzzy number c i . After that, FTSARIMA prediction intervals of y s(UCL) and y s(LCL) could be obtained.
In the final part, by utilizing the variance decomposition to identify the weight coefficient, then combine them into the FMLR interval value and FTSARIMA interval value to obtain a new result. Repetition of the FTSARIMA is carried out, using new y LCL and y UCL , to obtain the CFTSARIMA model. When the influence factors participate in forecasting processing, a long-term prediction becomes possible.
For future work, the prediction could improve in the following aspects: e fuzzy clustering still needs to increase the method's robustness to nonstationarity and concept-drift in this case, and incremental learning may be helpful.
Although most of the actual observations in the longterm prediction have entered the prediction range, there is still a gap compared with the accurate prediction. Consider combining the Markov chain in fuzzy clustering to improve accuracy.
Supposing the influence factors have nonlinearity by employing the adaptive filtering algorithm combined with the fuzzy method to forecast the chaotic time series may solve this issue.

Data Availability
e data used to support this findings of the study are included within the supplementary information file.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of the paper.

Acknowledgments
is work was substantially supported by the Research and Demonstration of Regulation Technology of the "Golden Waterway" of the Yangtze River Projects (no. 2016YFC0402103) and a project sponsored by the Yangtze River Survey, Planning, Design and Research Co., Ltd..

Supplementary Materials
In the Supplementary Materials, most data came from "WATERBORNE COMMERCE OF THE United States" and the first line was given the timeline, from the year 1970 to 2016. e first part from lines 2 to 10 refers to the Mississippi River (river system), coal, petro and petro products, chemicals, crude materials, manufactured goods, food and farm, and other. e second part from lines 11 to 19 refers to the Mississippi River (mainstream), coal, petro and petro products, chemicals, crude materials, manufactured goods, food and farm, and other. Both the first and second part unit of measurement is one hundred million tons. e third part from lines 20 to 26 refers to population (100 million), urban population (100 million), urbanization rate (percentage), GDP (USD 1 billion), per capita GDP (USD 10 million), US manufacturing ($1 billion), and US rate of industrialization (percentage). (Supplementary  Materials)