Research Article  Open Access
Cem Kocak, "FirstOrder ARMA Type Fuzzy Time Series Method Based on Fuzzy Logic Relation Tables", Mathematical Problems in Engineering, vol. 2013, Article ID 769125, 12 pages, 2013. https://doi.org/10.1155/2013/769125
FirstOrder ARMA Type Fuzzy Time Series Method Based on Fuzzy Logic Relation Tables
Abstract
Fuzzy time series approaches have an important deficiency according to classical time series approaches. This deficiency comes from the fact that all of the fuzzy time series models developed in the literature use autoregressive (AR) variables, without any studies that also make use of moving averages (MAs) variables with the exception of only one study (Egrioglu et al. (2013)). In order to eliminate this deficiency, it is necessary to have many of daily life time series be expressed with Autoregressive Moving Averages (ARMAs) models that are based not only on the lagged values of the time series (AR variables) but also on the lagged values of the error series (MA variables). To that end, a new firstorder fuzzy ARMA(1,1) time series forecasting method solution algorithm based on fuzzy logic group relation tables has been developed. The new method proposed has been compared against some methods in the literature by applying them on Istanbul Stock Exchange national 100 index (IMKB) and Gold Prices time series in regards to forecasting performance.
1. Introduction
Fuzzy time series approaches not necessitating many of the limitations seen in classical time series approaches such as linearity, stationarity, and number of observations have increased the interest towards these approaches. Fuzzy time series concept first mentioned in the literature by Song and Chissom [1] was based on the fuzzy set theory of Zadeh [2]. Song and Chissom [3, 4] divided fuzzy time series into two groups, namely, time variant and time invariant. A vast majority of the studies in the literature are methods proposed for solving time invariant fuzzy time series. Because models are significantly effective on the forecasting performances during the determination of fuzzy relations stage, different approaches have been proposed in the literature. In the studies of Song and Chissom [1, 3, 4], relations are determined with complex matrix operations. In order to eliminate this complexity, a new firstorder fuzzy time series model has been proposed where fuzzy logic group relation tables are used in Chen’s [5] study with simplified operations not necessitating complex matrix processes. This approach of Chen [5] is used in many studies due to its positive effect on forecasting performance. Therefore, Chen [6] developed a new approach by using the fuzzy logic relation tables also in highorder fuzzy time series models. Because the methods proposed in the studies of Chen [5, 6] necessitate obtainment of many fuzzy logic group relation tables, they require numerous operations. Thus, studies where fuzzy relations are determined with artificial neural networks are commonly seen. Some Studies where artificial neural networks are used for determining fuzzy relations may be listed as the studies of Huarng and Yu [7], Aladag et al. [8], Yu and Huarng [9], and Yolcu et al. [10].
In majority of the studies in the literature, interval lengths are specified intuitively. In his study, Huarng [11] has proposed two separate approaches based on average and distribution to specify optimal interval length. The optimal interval lengths determined by the approach of Huarng [11] may be obtained to be very large values. Thus, Egrioglu et al. [12] proposed approaches based on the optimisation of interval length. Different from these studies, Huarng and Yu [13] proposed an approach based on ratio with interval length varying exponentially instead of determining a fixed interval length in the solution of firstorder fuzzy time series. Determining the ratio in the study of Huarng and Yu [13] requires many complex calculations. Therefore, Yolcu et al. [14] proposed a new approach that improves the approach of Huarng and Yu [13], based on the optimisation of ratio.
There are 3 most commonly used models in the analysis of single variable time series in classical time series approach. These are autoregressive (AR), moving averages (MAs), and mixed autoregressive moving averages (ARMAs) models. However, fuzzy time series methods developed in the literature focus on the AR model of the classical time series theory, without any study conducted on the utilisation of MA and ARMA models except Uslu et al. [15], Alpaslan et al. [16], and Aladag et al. [17]. These studies [15–17] on the other hand are methods that have been proposed for the solution of seasonal time series. No study has been made on the inclusion of error variable to the model for nonseasonal fuzzy time series with the exception of Egrioglu et al. [18] study. This study [18] is the first study developed as first order fuzzy ARMA type model based on particle swarm optimization for the solution of nonseasonal time series in the literature. All of the models in the literature have mentioned issues such as the use of universe of discourse partitioning, membership order, model order, and artificial intelligence approaches. However, the fuzzy time series models proposed in the literature only including AR variables which may lead to a model specification error. In many modelling of the real life time series, MA variables are also required. In this sense, use of only AR variables for the solution of fuzzy time series requiring also MA variables for modelling becomes insufficient regarding the forecasting performance.
For the purpose of eliminating the adverse effects mentioned, a solution algorithm for a new firstorder fuzzy ARMA(1,1) time series forecast model where fuzzy relations are determined based on fuzzy logic group relations has been proposed in this study. The logic for determining fuzzy relations in the proposed method is an approach similar to that of the study of Chen [5], aiming to show that the forecasting performance can be significantly improved when the model specification error in the method of Chen [5], accepted as a fundamental approach in literature, is eliminated. For many real life time series, performance can be increased through the use of highorder fuzzy time series models due to the ability to realise solution with more information. However, because the proposed fuzzy ARMA(1,1) model uses a second variable (error variable), it utilises more information to the extent it eliminates the model error, and although just a first order model, it occurs to have a better forecasting performance than that of high order fuzzy time series methods.
In this paper, the time series which were used in application are examined according to longrange dependence. Time series can be classified short range or longrange dependence time series. Many time series have longrange dependence. The longrange dependence time series are forecasted different methods like autoregressive fractionally integrated moving average (ARFIMA). ARFIMA models have fractionally differencing parameter. First studies were concerned with estimation of fractional differencing parameter in fractional white noise processes. statistic was proposed in Hurst [19]. The other important studies about fractional differenced processes are Li and Zhao [20, 21], Li [22], Stanley et al. [23], Werner [24], Beran [25], Ivanov et al. [26], Podobnik et al. [27], Zevallos and Palma [28], and Bhansali and Kokoszka [29].
In the second section of the study, basic definitions regarding fuzzy time series have been provided. In the third section, the proposed fuzzy time series forecasting model has been defined and the solution algorithm has been provided. In the fourth section, the proposed method has been applied to Istanbul Stock Exchange (IMKB) national 100 index time series and gold prices in 2009 taken from Turkish Republic Central Bank (TCMB) website, comparing it to some other methods in the literature regarding forecasting performance. An indepth comparison has been made in this section. In the fifth section, the study has been summed up by discussing the results obtained.
2. Definition of Fuzzy Time Series
Fuzzy time series concepts and definitions have been developed in accordance to the lagged variables of times series (AR, autoregressive) in all studies conducted in the literature. Main time series definitions developed using AR variables are listed below.
Definition 1. Let , a subset of real numbers, be the universe of discourse on which fuzzy sets are defined. If is a collection of then is called a fuzzy time series defined on [1, 3, 4].
Definition 2. Let us consider the fuzzy relation between , , and . For any value, if is independent from , then . In this case, is called the time invariant fuzzy time series, while otherwise called as time variant fuzzy time series [1].
Definition 3. If the fuzzy time series is only affected by one lagged fuzzy time series, then the fuzzy relation between and is expressed as This is called as a firstorder fuzzy time series forecasting model. Then this relation can be expressed as The “” operator in (2) had been determined as the maxmin operator by Song and Chissom [1, 3, 4].
Definition 4. If fuzzy time series is affected by the lagged fuzzy time series of , then the fuzzy relation between fuzzy time series and fuzzy time series may be expressed as and is called the th order fuzzy time series forecasting model [6].
3. The Proposed Method
Based on the definition of fuzzy AR(1) given in Definition 3, the definition of the main fuzzy time series to represent the fuzzy ARMA(1,1) model is expressed as follows.
Definition 5. Let be a fuzzy time series and let be the fuzzy error series obtained from fuzzy time series. If is affected by one lagged and one lagged fuzzy time series, then the relationship can be expressed as This is called as firstorder fuzzy autoregressive moving averages (ARMA(1,1)) time series forecasting model [18].
In this study, an algorithm has been proposed for solving the ARMA(1,1) fuzzy time series forecasting model defined in (4). In the algorithm proposed, initially the AR(1) fuzzy time series model defined in (1) is estimated. Later on, errors are calculated by taking the differences between the observed values of the times series and the forecasts obtained through the solution of fuzzy AR(1). By using these errors, the fuzzy ARMA(1,1) model defined in (4) is estimated. The algorithm of the proposed approach is given below.
Algorithm 6 6. The proposed method’s algorithm.
Step 1 (the universe of discourse and subintervals are defined). The beginning and the ending points of the universe of discourse for time series are determined. Then is divided into subintervals according to appropriate interval length. Definition of interval length is up to the researcher. It should not be forgotten that the interval length to be determined affects the number of subintervals. If the smallest value of the time series is taken as , the largest value as , and two arbitrary values as and , the universal set may be defined as the closed interval of subintervals determined for are the subintervals of the universal set , which is defined as For example, for and when is selected as [40, 100] and interval length is selected as 10, subintervals are specified as , , , , , and .
Step 2. For the time series, fuzzy sets are defined according to the universal set and the divisions of (). These fuzzy sets are expressed as For , For example, according to (7) and (8), the fuzzy set of can be expressed as .
Step 3 (observations are fuzzified). Subintervals () where each observation occurs are defined. Then the fuzzy set where the defined subinterval has the highest membership value is determined. The fuzzy value of the observation is this fuzzy set defined.
Step 4. For the purpose of determining fuzzy relations, fuzzy logic relations are identified and a fuzzy logic group relation table is formed.
For example, where the fuzzy logic relations are as , , , according to Chen method [5], the fuzzy logic group relation for fuzzy value is as . In the proposed method, this relation occurs to be as . Thus, a little improvement has been realised in the fuzzy AR(1) model of Chen [5] with the proposed method.
Step 5 (fuzzy forecasts are obtained). As and , 3 possible situations regarding forecast obtainment are as follows.
Situation 1. If a relation of is valid on the fuzzy group relation table where affects only pieces of , then the fuzzy forecast is . For example, if the group relation for is as that this relation is repeated a few times in the time series, then the fuzzy forecast is specified as .
Situation 2. If a relation of is valid on the fuzzy group relation table where affects pieces of , pieces and pieces , then the fuzzy forecast is , comprising of pieces of fuzzy values. For example; if the group relation for is as , then the fuzzy forecast is specified as .
Situation 3. If on the fuzzy group relation table, the fuzzy forecast occurs to be . For example, if the group relation for is , then the fuzzy forecast is .
Step 6 (defuzzification process is executed). In this step, centralisation method is used. When the fuzzy forecast for Situation 1 and Situation 3 defined in Step 5 is , then the defuzzy forecast should be the middle point of the subinterval that has the highest membership value within the fuzzy set . For Situation 2, the defuzzy forecast is calculated with the weighted average formula below, by using the middle points () of the intervals that have the highest membership value of each of the fuzzy sets. Consider On an example time series, when solutions are made according to the proposed method from Step 1 to Step 6 as per the determined subintervals of , , and , the solutions displayed in Table 1 are obtained.
Step 7. Errors are calculated by taking the differences between the observed time series values and the defuzzified forecast values obtained in Step 6. Real values of the time series are and the defuzzified forecast values obtained in Step 6 are ; the error series is calculated as follows:
Step 8. For the errors, the universe of discourse set is defined as and subintervals are defined as .
The same partition of the universe of discourse processes done in Step 1 is made for the error series.
Step 9. Fuzzy sets based on the universal set () and partitions () are defined for the errors. The fuzzy sets are expressed as For , For example, according to (11) and (12), the fuzzy set of can be specified as .
Step 10 (error series is fuzzified). Subintervals () for each observation are determined. Then the fuzzy set where the determined subinterval has the highest membership value is defined. The fuzzy value of the observation is this fuzzy set .
Step 11 (fuzzy relations are determined and fuzzy logic group relation table is formed). The fuzzy value of the fuzzy relations time series and the fuzzy value of the error series are determined by taking the fuzzy values into consideration together. For the one lagged fuzzy value of the th observation being , one lagged error value being , and the fuzzy value being , the fuzzy relation given in (4) occurs to be . Thus, the fuzzy values are formed of ordered pairs, and a relation in the manner of one lagged time series and error affecting the time series is mentionable. For example, when the fuzzy logic relations are as , , , the fuzzy logic group relation for fuzzy value occurs to be .
Step 12 (fuzzy forecasts are obtained). As ve and , 3 possible situations regarding forecast obtainment are as follows.
Situation 1. If a relation of is valid on the fuzzy group relation table where affects only pieces of , then the fuzzy forecast is . For example, if the group relation for is as that this relation is repeated within a few times in the time series, then the fuzzy forecast is determined as .
Situation 2. If a relation of is valid on the fuzzy group relation table where affects pieces of , pieces and pieces , then the fuzzy forecast is , comprising of pieces of fuzzy values. For example, if the group relation for is as , then the fuzzy forecast is determined as .
Situation 3. If on the fuzzy group relation table, the fuzzy forecast occurs to be . For example, if the group relation for is , then the fuzzy forecast is .
Step 13 (defuzzification process is made). In this step, centralisation method is used. When the fuzzy forecast for Situation 1 and Situation 3 defined in Step 5 is , then the defuzzy forecast should be the middle point of the subinterval that has the highest membership value within the fuzzy set . For Situation 2, the defuzzy forecast is calculated with the weighted average formula given in (9).
On an example time series, when solutions are made according to the proposed method as per the determined time series subintervals of , , and and error series subintervals of , , and , solutions displayed in Table 2 are obtained.
4. Application
The performance indicators of the root mean square error (RMSE), mean average percentage error (MAPE), and direction accuracy (DA) values used for comparison of the results obtained are as follows: Data is divided into two within the applications, assigning the first part as training set and the second part as test set obtained through taking the last observations into consideration of which number was predetermined. By looking up the fuzzy relation table obtained for the training set as per the Steps 5 and 12 of the proposed method, the fuzzy forecast of the test set and from there the defuzzy forecasts of the test set are calculated by utilizing Steps 6 and 13. After, the RMSE, MAPE, and DA values for the test set are calculated. Therefore, the future performances of the methods are determined. As the forecasts with the lowest RMSE value calculated for the test set provide the best result of the used method, the future performances of the forecasts are obtained along with the aid of MAPE and DA values.
Solution of vast majority of some fuzzy forecasting methods in the literature is realised according to the specified number of fuzzy sets, and some are realised according to interval lengths. For the purpose of maintaining consistency during the comparison of forecast performances, the interval length to be used at fuzzification stage is determined as to have the number of fuzzy sets 5 as the lowest and 35 as the highest for each application data and for all methods. Therefore, in case the universal set division is realised in accordance to interval length, the interval lengths to be tried have been specified by calculation with the formula below: The operations below have been realised when making solution via the proposed method.(i)The last number of data has been specified as test set, aiming to increase the future performances.(ii)During the fuzzification of time series stage of the proposed method, different time series interval lengths for the division of universal set have been tried. Data have been solved from Step 1 to Step 6 according to these intervals lengths. And so, a lot of forecasts have been obtained. The test set forecast with the smallest RMSE value among these forecasts has been determined as the best result of the fuzzy AR(1) model.(iii)The error value of the first observation of data has been assumed as 0, while the error values of other observations have been calculated with the formula (9) by using the training set and test set forecasts obtained through the best result of the fuzzy AR(1) model. Thus the error series has been obtained.(iv)Different time series interval lengths and different error series interval lengths have been tried by solving data from Step 1 to Step 13. Among these trials, the test set forecast with the smallest RMSE value has been determined as the best result of the fuzzy ARMA(1,1) model.
For the purpose of comparing the proposed method with the other fuzzy time series methods in the literature, 2 different data sets comprising of less observations (smaller sample size) and more observations (larger sample size) have been used. One of these data sets is the IMKB time series seen in Figure 1 comprising of 53 observations between the dates of 01.10.2010 and 23.12.2010. The second data set is the gold prices time series seen in Figure 2 comprising of 248 observations received from Turkish Republic Central Bank (TCMB) website between the dates of 02.01.2009 and 31.12.2009.
In the solution of the IMKB data given in Figure 1 via the methods in the literature by taking the last 7 observations and the last 15 observations as test set, one has the following.(i)During the fuzzification stage of the Song and Chissom [1] firstorder fuzzy time series method, among the 31 different results obtained from increasing the fuzzy set number between 5 and 35, the test set forecast that has the minimum RMSE value has been determined as the best result of the method from this results.(ii)For each of the 2nd, 3rd, 4th, and 5thorder models of Chen [5] firstorder fuzzy time series forecast method, Chen [6] high order fuzzy time series forecast method, and Aladag et al. [8] high order fuzzy time series forecast method, the RMSE values have been found for different intervals lengths being increased 100 units between 300 and 2300. The test set forecasts that have the minimum RMSE values among these 21 trials have been determined as the best results of the methods.(iii)By using the optimal interval lengths calculated with the distributionbased approach of Huarng [11] and the averagebased approaches, solution has been executed as per the firstorder fuzzy time series forecast method of Chen [5]. Thus, the best results of the test set from the distribution based approach and the averagebased approaches have been obtained via a single trial.(iv)In the application of the ratiobased approach of Huarng and Yu [13], alpha parameter has been taken as 0.50, obtaining the best result for the test set of the method in one trial.
In the solution of the IMKB time series seen in Figure 1 via the proposed method by taking the number of test sets 7 and 15, one has the following.(i)The division of universal set have been taken different values as increasing the interval length between 300 and 2300 by 100 units and different forecasts have been obtained by solving from Step 1 to Step 6. The test set forecast with the smallest RMSE value among these forecasts has been determined as the best result of the fuzzy AR(1) model. The best fuzzy AR(1) results have been obtained when the interval length is 300 for 7 as the number of test sets and when the interval length is 900 for 15 as the number of test sets.(ii)The error value of the first observation of data has been assumed as 0, while the error values of other observations have been calculated via the formula (9) by utilising data and forecasts obtained through the best result of the fuzzy AR(1) model. Thus the error series has been obtained for 53 observations.(iii)Different trials have been made by increasing the interval length between 300 and 2300 by 100 units for the time series and by increasing the interval length between 300 and 2100 by 100 units for the error series by solving from Step 1 to Step 13. Among these trials, the test set forecast with the smallest RMSE value has been determined as the best result of the fuzzy ARMA(1,1) model.
During the application of the proposed method and the methods in the literature on IMKB time series, the parameters with which the forecasts with the best test set performance for 7 and 15 numbers of test sets occurred to be(i)for the application of Song and Chissom [1] method, when the number of fuzzy sets is 9 for 7 and 20 for 15,(ii)for the application of Chen [5] method, when the interval length is 300 for 7 and 900 for 15,(iii)for the application of distributionbased Huarng [11] approach, when the interval length is 1000, and for the application of average based approach, when the interval length is 200,(iv)for the application of ratiobased Huarng and Yu [13] approach, when the sample percentile alpha = 0.5,(v)for the application of Chen [6] method, when the interval length is 2200 in 3rdorder model for 7 and 1400 in 2ndorder model for 15,(vi)for the application of Aladag et al. [8] method, when the interval length is 600 on 2nddegree model and unit number of artificial neural network hidden layers is 5 for 7, and when the interval length is 1500 on 2nddegree model and unit number of artificial neural network hidden layers is 6 for 15,(vii)for the application of the proposed fuzzy ARMA(1,1) method, when the interval length of time series is 2200 and the interval length of error series is 1400 for 7, and when the interval length of time series is 1400 and the interval length of error series is 400 for 15.
Best forecasts and forecast performances of all methods in result of IMKB time series solution for 7 observation test set are summarised in Table 3.
 
The best situation. 
When Table 3 is analyzed, it is seen in result of the solution of IMKB time series for 7 observation test set that the proposed method produced the best forecasting performance with a minimum RMSE value of 606.07, minimum MAPE value of 0.762%, and maximum direction accuracy of 83.33%. The graphs of the last 7 observations of IMKB time series along with the 7observation test set forecasts obtained with the proposed method are shown together in Figure 3.
Best forecasts and forecast performances of all methods in result of IMKB time series solution for 15observation test set are summarised in Table 4.
 
The best result. 
When Table 4 is analyzed, it is seen that the proposed method produced the best forecasting performance with a minimum RMSE value of 865.28, minimum MAPE value of 1.029%, and maximum direction accuracy of 71.43% in result of the solution of IMKB time series for 15observation test set. The graphs of the last 15 observations of IMKB time series along with the 15observation test set forecasts obtained with the proposed method are shown together in Figure 4.
In result of the solutions of IMKB time series, it has been observed that the proposed method significantly increased the future forecasting performance compared to other methods. Also in the graphs within Figures 3 and 4, the results of the proposed method are seen to be considerably similar to the test set values.
In the solution of the gold prices data given in Figure 2 via the methods in the literature by taking the last 30 observations and the last 45 observations as test set, one has the following.(i)Gold prices solution of Song and Chissom [1], Huarng [11], and Huarng and Yu [13] methods has been conducted just as previously done on the abovementioned IMKB time series.(ii)For each of the 2nd, 3rd, 4th, and 5thorder models of Chen [5] firstorder fuzzy time series forecast method, Chen [6] highorder fuzzy time series forecast method, and Aladag et al. [8] high order fuzzy time series forecast method, the RMSE values have been found for different lengths being increased 100 units between 500 and 3500. The test set forecasts that have the minimum RMSE values among these 31 trials have been determined as the best results of the methods.
In the solution of the gold prices data given in Figure 2 via the proposed method by taking the last 30 observations and the last 45 observations as test set, one has the following.(i)During the fuzzification of time series of the proposed method, interval lengths of have been tried as increasing the interval length between 500 and 3500 by 100 units by solving from Steps 1 and 6. The test set forecast with the smallest RMSE value among these forecasts has been determined as the best result of the fuzzy AR(1) model.(ii)The error series have been obtained via the formula (9) for 248 observations with the same calculation made previously in the application of IMKB data. (iii)Different trials have been made by increasing the interval length between 500 and 3500 by 100 units for the time series and by increasing the interval length between 100 and 1100 by 50 units for the error series by solving from Step 1 to Step 13. Among these trials, the test set forecast with the smallest RMSE value has been determined as the best result of the fuzzy ARMA(1,1) model.
During the application of the proposed method and the methods in the literature on gold prices time series, the parameters with which the forecasts with the best test set performance for 30 and 45 test sets occurred to be the following:(i)for the application of Song and Chissom [1] method, when the number of fuzzy sets is 10 for both 30 and 45 test sets,(ii)for the application of Chen [5] method, when the interval length is 600 for 30 and 1700 for 45,(iii)for the application of distributionbased Huarng approach [11], when the interval length is 400, and for the application of average based approach, when the interval length is 200,(iv)for the application of ratiobased Huarng and Yu [13] approach, when the sample percentile alpha = 0.5,(v)for the application of Chen method [6], when the interval length is 1900 in 5th order model for 30 and 800 in 3rdorder model for 45,(vi)for the application of Aladag et al. [8], when the interval length is 600 on 5thdegree model and unit number of artificial neural network hidden layers is 5 for 30, and when the interval length is 800 on 3rd degree model and unit number of artificial neural network hidden layers is 4 for 45,(vii)for the application of the proposed fuzzy ARMA(1,1) method, when the interval length of time series is 800 and the interval length of error series is 2500 for 30, and when the interval length of time series is 900 and the interval length of error series is 1000 for 45.
Best forecasts and forecast performances of all methods in result of gold prices time series solution for 30observation test set are summarised in Table 5.
 
The best result. 
When Table 5 is observed, it is seen in result of the solution of gold prices time series for 30 observation test set that the proposed method produced the best forecasting performance with a minimum RMSE value of 707.71, minimum MAPE value of 1.028%, and maximum direction accuracy of 62.07%. The graphs of the last 30 observations of gold prices time series along with the 30observation test set forecasts obtained with the proposed method are shown together in Figure 5.
Best forecasts and forecast performances of all methods in result of gold prices time series solution for 45observation test set are summarised in Table 6. Furthermore, the graphs of the last 45 observations of gold prices time series along with the 45observation test set forecasts obtained with the proposed method are shown together in Figure 6.
 
The best result. 
When Table 6 is evaluated, it is seen in result of the solution of gold prices time series for 45 observation test set that the proposed method produced the best forecasting performance with a minimum RMSE value of 719.69, minimum MAPE value of 1.072%.
In result of the solutions of gold prices time series, it has been observed that the proposed method significantly increased the future forecasting performance compared to other methods. Also in the graphs within Figures 5 and 6, the results of the proposed method are seen to be considerably similar to the test set values.
5. Discussion and Conclusion
MA variables are not included in the fuzzy time series forecast models proposed in the literature whereas real life time series are also influenced from MA variables in addition to AR variables. Therefore, redefining fuzzy times series methods as models including also MA variables are a more realistic act. In this study, a solution algorithm for a new firstorder fuzzy ARMA(1,1) time series forecast model containing not only AR but also MA variables is proposed based on group relation tables. The method proposed is a basic algorithm similar to Chen [5] approach, aiming at eliminating the model specification error formed due to the exclusion of MA variables. In conclusion of the applications, it has been observed that the proposed method has higher forecasting performance than many of the fuzzy time series forecasting methods commonly used in the literature. It is an important finding that although the proposed method is a basic method based on group relation tables, it may have a higher forecasting performance even than the high order fuzzy time series methods based on artificial neural networks. Therefore, it is obvious that forecasting performance is going to increase significantly when fuzzy ARMA models are developed where fuzzy relations are specified with artificial neural networks and artificial intelligence methods or where membership values are used for specification of fuzzy relations. Thus, the proposed method may be provided with a more systematic structure and higher forecasting performance with improvements that may be done on various stages of the method during future studies. Moreover, there is no linear model assumption in the proposed fuzzy time series method. Thus, the proposed method and other fuzzy time series methods can be applied to nonlinear time series. In the application, proposed method is applied to short range dependent time series. In the future studies, it will be researched about the performance of the proposed method for long range dependent time series.
Appendix
The test was applied to two time series which are used in the application. test was applied by using “FinMetrics module of SPLUS package program”. The obtained results are shown in Table 7. It is obtained that both of time series have short range dependence.

References
 Q. Song and B. S. Chissom, “Fuzzy time series and its models,” Fuzzy Sets and Systems, vol. 54, no. 3, pp. 269–277, 1993. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 L. A. Zadeh, “Fuzzy sets,” Information and Computation, vol. 8, pp. 338–353, 1965. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 Q. Song and B. S. Chissom, “Forecasting enrollments with fuzzy time series. Part I,” Fuzzy Sets and Systems, vol. 54, pp. 1–10, 1993. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 Q. Song and B. S. Chissom, “Forecasting enrollments with fuzzy time series. Part II,” Fuzzy Sets and Systems, vol. 62, no. 1, pp. 1–8, 1994. View at: Google Scholar
 S.M. Chen, “Forecasting enrollments based on fuzzy time series,” Fuzzy Sets and Systems, vol. 81, no. 3, pp. 311–319, 1996. View at: Google Scholar
 S.M. Chen, “Forecasting enrollments based on highorder fuzzy time series,” Cybernetics and Systems, vol. 33, no. 1, pp. 1–16, 2002. View at: Publisher Site  Google Scholar
 K. Huarng and H.K. Yu, “The application of neural networks to forecast fuzzy time series,” Physica A, vol. 363, no. 2, pp. 481–491, 2006. View at: Publisher Site  Google Scholar
 C. H. Aladag, M. A. Basaran, E. Egrioglu, U. Yolcu, and V. R. Uslu, “Forecasting in high order fuzzy times series by using neural networks to define fuzzy relations,” Expert Systems with Applications, vol. 36, no. 3, pp. 4228–4231, 2009. View at: Publisher Site  Google Scholar
 T. H.K. Yu and K.H. Huarng, “A neural networkbased fuzzy time series model to improve forecasting,” Expert Systems with Applications, vol. 37, no. 4, pp. 3366–3372, 2010. View at: Publisher Site  Google Scholar
 U. Yolcu, C. H. Aladag, E. Egrioglu, and V. R. Uslu, “Time series forecasting with a novel fuzzy time series approach: an example for Istanbul stock market,” Journal of Statistical Computation and Simulation, vol. 83, no. 4, pp. 599–612, 2013. View at: Google Scholar
 K. Huarng, “Effective lengths of intervals to improve forecasting in fuzzy time series,” Fuzzy Sets and Systems, vol. 123, no. 3, pp. 387–394, 2001. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 E. Egrioglu, C. H. Aladag, U. Yolcu, V. R. Uslu, and M. A. Basaran, “Finding an optimal interval length in high order fuzzy time series,” Expert Systems with Applications, vol. 37, no. 7, pp. 5052–5055, 2010. View at: Publisher Site  Google Scholar
 K. Huarng and T. H.K. Yu, “Ratiobased lengths of intervals to improve fuzzy time series forecasting,” IEEE Transactions on Systems, Man, and Cybernetics B, vol. 36, no. 2, pp. 328–340, 2006. View at: Publisher Site  Google Scholar
 U. Yolcu, E. Egrioglu, V. R. Uslu, M. A. Basaran, and C. H. Aladag, “A new approach for determining the length of intervals for fuzzy time series,” Applied Soft Computing Journal, vol. 9, no. 2, pp. 647–651, 2009. View at: Publisher Site  Google Scholar
 V. R. Uslu, C. H. Aladag, U. Yolcu, and E. Egrioğlu, “A new hybrid approach for forecasting a seasonal fuzzy time series,” in Proceedings of the International Symposium Computing Science and Engineering Proceeding Book, pp. 1152–1158, 2010. View at: Google Scholar
 F. Alpaslan, O. Cagcag, C. H. Aladag, U. Yolcu, and E. Egrioglu, “A novel seasonal fuzzy time series method,” Hacettepe Journal of Mathematics and Statistics, vol. 41, no. 3, pp. 375–385, 2012. View at: Google Scholar  MathSciNet
 S. Aladag, C. H. Aladag, T. Mentes, and E. Egrioglu, “A new seasonal fuzzy time series method based on the multiplicative neuron model and SARIMA,” Hacettepe Journal of Mathematics and Statistics, vol. 41, no. 3, pp. 337–345, 2012. View at: Google Scholar  MathSciNet
 E. Egrioglu, U. Yolcu, C. H. Aladag, and C. Kocak, “An ARMA type fuzzy time series forecasting method based on particle swarm optimization,” Mathematical Problems in Engineering, vol. 2013, Article ID 935815, 12 pages, 2013. View at: Publisher Site  Google Scholar
 H. E. Hurst, “Longterm storage capacity of reservoirs,” Transactions of the American Society of Civil Engineers, vol. 116, pp. 770–779, 1951. View at: Google Scholar
 M. Li and W. Zhao, “On $1/f$ noise,” Mathematical Problems in Engineering, vol. 2012, Article ID 673648, 22 pages, 2012. View at: Google Scholar  Zentralblatt MATH  MathSciNet
 M. Li and W. Zhao, “Visiting power laws in cyberphysical networking systems,” Mathematical Problems in Engineering, vol. 2012, Article ID 302786, 13 pages, 2012. View at: Publisher Site  Google Scholar
 M. Li, “Fractal time series: a tutorial review,” Mathematical Problems in Engineering, vol. 2010, Article ID 157264, 26 pages, 2010. View at: Publisher Site  Google Scholar  Zentralblatt MATH  MathSciNet
 H. E. Stanley, S. V. Buldyrev, A. L. Goldberger, S. Havlin, C.K. Peng, and M. Simons, “Longrange powerlaw correlations in condensed matter physics and biophysics,” Physica A, vol. 200, no. 1 4, pp. 4–24, 1993. View at: Google Scholar
 G. Werner, “Fractals in the nervous system: conceptual implications for theoretical neuroscience,” Frontiers in Fractal Physiology, vol. 1, article 15, 28 pages, 2010. View at: Publisher Site  Google Scholar
 J. Beran, “Statistical methods for data with longrange dependence,” Statisticla Science, vol. 7, no. 4, pp. 404–416, 1992. View at: Google Scholar
 P. C. Ivanov, L. A. Nunes Amaral, A. L. Goldberger et al., “From 1/f noise to multifractal cascades in heartbeat dynamics,” Chaos, vol. 11, no. 3, pp. 641–652, 2001. View at: Publisher Site  Google Scholar
 B. Podobnik, D. Horvatic, A. L. Ng, H. E. Stanley, and P. Ch. Ivanov, “Modeling longrange crosscorrelations in twocomponent ARFIMA and FIARCH processes,” Physica A, vol. 387, no. 15, pp. 3954–3959, 2008. View at: Publisher Site  Google Scholar  MathSciNet
 M. Zevallos and W. Palma, “Minimum distance estimation of ARFIMA processes,” Computational Statistics & Data Analysis, vol. 58, pp. 242–256, 2013. View at: Publisher Site  Google Scholar  MathSciNet
 R. J. Bhansali and P. S. Kokoszka, “Prediction of longmemory time series: a tutorial review,” in Processes With LongRange Correlations, vol. 621 of Lecture Notes in Physics, pp. 3–21, 2003. View at: Google Scholar
Copyright
Copyright © 2013 Cem Kocak. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.