Robust Inventory System Optimization Based on Simulation and Multiple Criteria Decision Making

Inventory management in retailers is difficult and complex decision making process which is related to the conflict criteria, also existence of cyclic changes and trend in demand is inevitable in many industries. In this paper, simulation modeling is considered as efficient tool for modeling of retailer multiproduct inventory system. For simulation model optimization, a novel multicriteria and robust surrogate model is designed based on multiple attribute decision making (MADM) method, design of experiments (DOE), and principal component analysis (PCA). This approach as a main contribution of this paper, provides a framework for robust multiple criteria decision making under uncertainty.


Introduction
Supply chain is complicated system that includes many components such as suppliers, factories, distribution centers, and retailers.These components are connected to each other by three streams of financial, information, and material flow.In many references, there is a fundamental hypothesis that financial flow is upstream from customers to suppliers and material flow is downstream from suppliers to customers while information flow is mutual.By the way, sometimes, this hypothesis is neglected about material flow (e.g., reveres logistic).
Material flow plays dominant role in supply chain and is defined as inventory problem.Inventory actually is the bridge that connects material handling and production planning to the supply chain [1].On the other hand, retailers in highly competitive market face dynamic change of demand including seasonal cyclic change and long-term trend.These situations caused multiproduct inventory systems to change into complex multicriteria systems.
Complex systems are challenging in the case of modelling and computation.In most of the problems, modelling of complex systems is very time consuming and also is not error free.Furthermore, modelling of complex systems needs too much computational effort to solve and sometimes they are not solvable in reasonable time.In this situation, importance of data-driven methods emerges [2][3][4].
Although model-based approach (i.e., simulation model) has many advantages, such as what if analysis and ability of scenario generation, data-driven approach is less time consuming and also easy to implement.So, these approaches can be combined to achieve tailored approach which inherits advantages of model-based and data-driven approaches.In this research we used discrete event simulation model and data-driven methods like principal component analysis and multiple attribute decision making to design an efficient framework for robust and multiobjective optimization of retailer inventory system with multiproduct as complex system.
The rest of the paper is structured as follows.Literature review is presented in Section 2. Problem statement is presented in Section 3, proposed framework which contains simulation modeling and surrogate design is described in Section 4, Section 5 provides numerical result, and finally Section 6 is dedicated to conclusion.

Literature Review
As the investigated problem encompasses two issues of simulation optimization and inventory problem, related literature 2 Modelling and Simulation in Engineering is categorized in two separated parts for providing more supportive literature review.

Inventory Problem.
Basic problems of inventory system are studied thoroughly by [5][6][7].Economical order quantity (EOQ) is the simplest model of inventory problem; in this model, demand of each period is constant and time independent.EOQ does not consider lost sale, back order, and other cost for reason of simplification and emphasizes only on holding and ordering cost.Cheng [8] developed a model of inventory system with cost-dependent demand and included production cost in proposed model.Chen et al. [9] studied back order by fuzzy technique and Zhao et al. [10] suggested analytic model with demand according to time series.They concluded that, in time-dependent demand, efficiency of EOQ model increases with shorter lead time and weaker autocorrelation.
Basic models that were offered for inventory problem only have one objective function, including different inventory costs, while advancement in technology and intensification of business competition caused necessity of other criteria.To the best of our knowledge, in recent years, service level has been observed in variety of supply chain and inventory problem as performance criterion.Adding service level criterion, classic definition of inventory problem is changed into multiobjective optimization.In such problem, inventory costs should be minimized while service level should be maximized.Available models for this problem are divided in two groups of deterministic and stochastic.Both of these models can be solved by three approaches of analytic methods, for example, mathematical programming, metaheuristic methods, and simulation.Lau et al. [11] benefited from simulation to compare four inventory management policies with two criteria of cost and service levels.They also surveyed preorder and information sharing impacts on their models.Xu and Zhao [12] used fuzzy rough simulation to optimize multiobjective problem to minimize wasted cost and maximize expected value of revenue.Hnaien et al. [13] Surveyed two-level assembly system with two objectives of service level and maintenance cost.They considered stochastic lead time and applied genetic algorithm to solve this problem.
Although cost and service levels are important criteria for inventory system performance, there are other criteria that should be considered such as amount of systems inventory, which is important factor with significant impact on inventory system behavior.Because of lead time uncertainty, that origins from natural disasters and transportation problems, organizations face delay in delivery, so most of them hold safety stocks.This phenomenon is the main cause that leads to the increase of inventory in hand.As inventory in hand increases, the inventory system is faced with holding cost and other problems like decrease in quality, lack of flexibility, and so on.So, amount of systems inventory can be considered as performance factor and inventory in hand should be minimized as mentioned in the just in time (JIT) philosophy.Purnomo et al. [14] researched about influence of periodic replenishment and continue replenishment inventory policies on supply chain and considered both inventory in hand and work in progress as performance factors.Now in recent years modeling of inventory system as multiobjective and stochastic problem is an interesting area for research.

Simulation Optimization.
One of the well-known simulation methods is discrete event simulation that is based on stochastic processes and could be efficient tool to capture stochastic behaviour of different systems.While discrete event simulation has several advantages, it is not optimization tool individually [15].By the way, because of its flexibility, simulation can be coupled with other techniques such as metaheuristic algorithms or stochastic methods.This synthesizes makes the powerful and advantageous approach of simulation optimization with vast area of research.Simulation optimization is powerful arsenal for optimization of complex systems such as military, aerospace, and supply chain [16].
To the best of our knowledge, there are three main optimization techniques that were reported as suitable techniques for simulation optimization.These considered methods are metaheuristic optimization, stochastic approximation (SA) methods, and surrogate models.
Fu [17] extensively described role of applied methods in simulation optimization and also surveyed techniques employed in optimization package of simulation software.Wang [18] used hybrid approach including genetic algorithm and artificial neural network for simulation optimization and Keskin et al. [19] applied discrete event simulation and scatter search algorithm for optimization of inventory system and vendor selection.Mazhari et al. [20] developed a simulation optimization framework based on hybrid simulation model (system dynamic and agent based model) and metaheuristic algorithm.Also, Duan and Liao [21] applied metaheuristic approach for developing simulation optimization framework in order to optimize replenishment policy of inventory system in capacitated supply chain.
Although using metaheuristic optimization algorithms is straight forward approach for simulation optimization, it is time consuming and needs high level of computational effort.So, it is inefficient in case of simulation optimization with more than one objective function.
While metaheuristic algorithms use stochastic searching methods, SA is based on gradient search.Because of noisy situation of observations, SA algorithms consider expected value of objective function.SA family includes attractive methods because their convergence is guaranteed theoretically.Simultaneous perturbation stochastic approximation (SPSA) is noteworthy algorithm of SA family.The theoretical aspects of SPSA are deeply described by Spall's [22][23][24][25] proposed simulation optimization framework for inventory control in supply chain based on SPSA.
In contrary to two former methods, surrogate modelling is postprocessing method, so it is less time consuming.In surrogate modelling, the main idea is to fit single surface to the decision space and use this surface instead of simulation model for optimization.In this area, response surface methodology (RSM) [26] and supervised learning methods (e.g., artificial neural network or support vector machine) are considerable.For instance, Can and Heavey [27] applied artificial neural network to develop surrogate model for discrete event simulation.Azadeh et al. [28] used artificial neural network for the designing of simulation optimization framework and they applied proposed framework for optimization of waiting time in tandem queue systems.Bornatico et al. [29] proposed a surrogate model based on redial basis function for simulation optimization of energy systems.Wan et al. [30] designed simulation optimization framework using least square support vector machine (LSSVM) for optimization of inventory level in three-stage supply chain.They also showed that proposed framework leads to better solution with less number of simulation runs in comparison with SPSA algorithm.Surrogate modelling is less time consuming in comparison to metaheuristic or SPSA approaches, but this approach loses accuracy in multiple objectives problem solving.
Although a bunch of papers published in simulation optimization area, to the best of our knowledge a tiny number of them are dedicated to the multiobjective optimization [31] and robustness [32,33].In this case, using metaheuristic and SPSA approaches is very time consuming and it is not economic for optimization of simulation model with accurate details.On the other hand, all of the reviewed approaches lose their accuracy when there are multiple objectives.With these considerations, this paper purposes a framework for optimization of detailed simulation model of inventory system with multiple objectives.Proposed framework is less time consuming in comparison with metaheuristic or SPSA approaches while it provides robust and accurate solutions.So, the proposed framework is relatively new and contribution of this research entails threefold as follows.
(i) We modelled cyclic and long-term demand based on nonparametric time series modelling for more realistic consideration.
(ii) We proposed surrogate model for robust and multiobjective optimization of multiproduct inventory system based on discrete event simulation, full factorial design of experiments (DOE), and multiple attribute decision making (MADM) technique.
(iii) Due to the stochastic nature of objective function, we employed principal component analysis (PCA) as statistical method to improve MADM performance.

Problem Statement
The problem is concerning retailer who sales office furniture and facility.The retailer sales four products, respectively, , , , and ; the aim is the optimization of inventory system according to information which is adapted from local business.Key features of retailer products from inventory view point are as shown in Table 1.
In this table, second column gives average demand of each product type in a year, third column provides holding cost of each product in a planning period, fourth column is dedicated to ordering cost of each type of products, and finally numbers of fifth column are cost of lost sales which are incurred to retailers when they cannot satisfy the demand of customers for each type of products.
The fundamental assumption that should be considered in this problem is as follows.
(1) Order cost for each type of products includes transportation and order registration cost.
(2) There is no backlog inventory so inventory level is nonnegative all the time.
(3) Profit of each product is considered as lost sale cost because unavailable products incur lost profits that are interpreted as cost of lost sale.
(4) According to the adapted information, these products have five years life cycle and then will be substituted with new products.
(5) Planning periods for system under study are as long as 20 days.
The notations that will be used to describe the problems are as follows.
Indices.Consider the following: In ( 1), the objective is minimizing the total cost including costs that depend on reordering, handling, and lost sales.Lost sales not only incur excess cost but also decrease retailer credit.So, ( 2) is considered to maximize service level independently.In (2), the objective function increases, where   is 1 for th demand in th period and such situation is possible if    is greater than   .In fact, (2) causes increase in inventory level while (1)-(3) causes decrease in inventory level.State of inventory level depends on number of orders in each planning period (  ), quantity of orders (  ), and quantity of demands (  ) while state of inventory position depends on inventory level and lead time (), so (3) is responsible for minimizing average of inventory position including inventory in hand (   ) and on-order inventory.Average of inventory position should be minimized in order to improve flexibility of retailer and approach to the just in time (JIT).
Demand of each product follows different pattern with both long-term and cyclic trend.So, mentioned objective functions are considered individually for each product type.

Proposed Framework
In the defined problem, demand of products exposes highly dynamic pattern and as time passes, demand and its variation increase; hence, multiresolution method is employed for demand modelling.Also, three different policies for inventory control are considered, which are reordered based on fixed quantity (FQ), fixed interval (FI), and demand forecasting (DF).Simulation of developed model is implemented in Arena 13.5.Optimization of simulation model is performed by surrogate model that is based on full factorial design of experiment (DOE).For construction of decision space, DOE factors include inventory policy, reorder point, and lead time with three levels for each of them.So, there are 3 3 = 27 different combinations of decision variables to form feasible scenarios.Ranking of produced scenarios is accomplished by MADM technique.For ranking of scenarios, three objective function values are considered (i.e., cost, service level, and average of inventory position) and, in addition, robustness of service level against demand fluctuation is considered.Also PCA is applied for more realistic weighting of objective function values based on their statistical influence on improvement of other objectives.Finally, interacting plot is employed for sensitivity analysis and investigation of solutions in detail.

Simulation Modelling.
Simulation modelling of problem consists of two parts which are modelling of demand and modelling of inventory policies.In this paper, demand is nonhomogeneous Poisson process and three different inventory policies based on continues reviewing, periodic reviewing, and periodic reviewing with forecasting of future demand are considered.

Modelling of Demand.
For customers demand modeling, multiresolution method is applied.Kuhl and Wilson [34] developed this method for simulation of nonhomogeneous Poisson process with trend and cyclic changes.This method estimates mean intensity function and the nonparametric nature of this method is one of the most important advantages in comparison with other methods.So, it is independent of statistical parameters and applicable in variety of problems.Furthermore, multiresolution method can support combination of multiple cyclic changes simultaneously.Another advantage of multiresolution is its capability in the modelling of nonsymmetric cyclic pattern.As in our case, demand has nonsymmetric pattern with cyclic changes and longterm trend, and multiresolution approach is a reasonable choice.More theoretical and application of used method are provided in [34].

Modelling of Inventory Policies.
As the main effort of this paper is inventory system optimization, modelling of inventory policies plays an important role in this problem.In this problem, optimization is manipulated by selection of appropriate inventory policy and configuration of its parameters to the way that leads to the optimal state of inventory system.In the inventory management, three approaches are common strategies which are fixed order quantity, fixed time interval, and forecasting methods [35].In the first strategy, inventory level should be reviewed continuously until it reaches below predetermined quantity (reorder point) then order would be organized with fixed quantity of inventory.In the second strategy, reviewing period is a fixed time interval, but quantity order is variable for each order that is based on consumption rate.While first strategy needs more effort for continuously reviewing of inventory level, the second strategy is easier to handle, but the risk of shortage in fixed interval strategy is more in comparison to fixed order quantity.So, due to the mitigation of shortage risk in fixed interval strategy, order quantity is slightly more than fixed order quantity [35].
In the third strategy, reviewing period is fixed as second policy, but reorder quantity is based on forecasting of future demands.As the mentioned strategies are fundamental in inventory management literature and are also common among retailers of office furniture, in this study three policies based on fixed order quantity, fixed time interval, and demand forecasting are developed as follows.
(1) Continuous reviewing with economic quantity order.
(2) Periodic reviewing with order quantity based on demand confidence interval during the consumption period.
(3) Periodic reviewing with order quantity based on forecasting of future demands.
Policy 1 (Fixed Quantity).Based on this policy, each demand will be satisfied if there is sufficient inventory in hand.
After satisfaction of each demand, inventory in hand will be checked to see if the inventory level reaches to reorder point (  ).If inventory level has reached to reorder point, economic order quantity (  ) would be organized, otherwise system waits for next demand.Economic order quantity is derived by Wilson formula [36].If there is not sufficient inventory to satisfy arrived demand, quantity of demand is considered as lost sale and lost profit treated as cost.
If no order has been organized, system reorder inventory otherwise waits for arriving of organized order according to adjusted lead time ().Logic of this policy is visible in Figure 1(a) and is labelled as FQ policy.
Policy 2 (Fixed Interval).In this policy, criterion for reordering is fixed time interval that is known as planning period.After this period, inventory level would be examined and reorder will be organized on condition that inventory level has reached to reorder point (  ).For more realistic consideration, order quantity (  ) is calculated based on demand cumulative distribution function in planning period.For example, in the inventory system, that is planned for 10% lost sale with demand which is distributed based on exponential distribution function, it should be ordered as much as cumulative probability of exponential distribution equals to 0.9.In this policy, inventory level would be updated after lead time ().Logic of this policy is shown in Figure 1(b) and is labelled as FI policy.
Policy 3 (Demand Forecasting).This policy is similar to policy 2 with some differences.In policy 2, probability distribution function of demand in planning period is estimated by historical data and then reorder is organized based on service level (cumulative probability of demand satisfaction).But, in this policy after fixed interval, forecasting of future demand will be performed.If inventory level has reached to reorder point (  ), order would be organized and inventory level would be updated after lead time (), otherwise only information would be updated.If forecasted demand is less than inventory capacity, order quantity (  ) would be as much as forecasted quantity otherwise inventory order will be as much as inventory capacity for each product.As demands follow time series with autocorrelation, forecasting is implemented by exponential smoothing method [37].
In this procedure,   is defined as demand estimation of th planning period.This term is calculated according to real demand of period  indicated by   and estimation of demand in previous period which is indicated by  −1 plus trend adjustment value of previous period which is indicated by  −1 .Impact of current real demand is considered by parameter  that is defined between zero and one.Trend adjustment value is derived based on   and  −1 plus trend adjustment value of previous period with consideration of  as weight parameter which is defined between zero and one.Equations related to   and   are expressed in ( 4) and ( 5), respectively.Finally, demand estimation plus trend adjustment form of future demand forecasting is indicated by  +1 which is described in (6).Logic of this policy is shown in Figure 1(c) and is labelled as DF policy.Consider this research, there are 3 factors), using full factorial design is preferable.On the other hand, because nonlinear behaviour of response variables (i.e., objective functions) centre point is considered.Each factor has high level "+, " low level "−, " and centre point "0." Factors and levels configurations are shown in Table 2 while Table 3 shows experimental design.
In both of these tables, minus sign is representative of low level of related factors while plus sign is indicator of high level of factors and zero determines centre point.For example, in the case of lead time (), 3 is low level, 7 is high level, and 5 is centre point.Lead time is considered with days as measuring unit and begins from time order is organized until the arrival of inventory.Reorder point is labelled as   and is intended to satisfy product demand during lead time ().Eventually, inventory policy is labelled as IP and refers to inventory policies that were described in Section 4.1.2.The mentioned design includes 27 runs and there are four response variables: cost, service level, average of inventory position, and robustness.In fact, first three response variables are expected values of defined objective functions (i.e., resp., (1), (2), and (3)), but robustness is a defined criterion which guarantees that optimal solution remains valid in situation where demand variation and intensification increase.For evaluation of robustness, demands are multiplied by random variable with normal distribution ( = 1,  = 1.5).This random variable is also restricted to positive values to prevent negative demands.Expected value of 4000 generated numbers with mentioned constraint is equal to 1.48; this result implies that multiplied random number increases both demand variation and intensification, simultaneously.Then, the percent of decrease in service level is considered as robustness for all 27 scenarios produced by DOE.It is obvious that a lower decrease in service level is desirable.

Weight Assignment
Based on PCA Method.PCA is the abbreviation of "principal component analysis" and is known as powerful data-driven method successfully applied in several areas [2].PCA analysis can determine key variables that govern most of the variability.This technique helps to remove less important variables to reduce dimension of data sets.PCA considers two criteria for each variable in data sets: the first is variability and the second is correlation with other variables.So, input for this method consists of variance covariance matrix [39,40].Although other methods like entropy is applicable in multiple attribute decision making problem, PCA advantage is considerable because variables correlations are not ignored.As our problem is constructed of highly correlated variables, PCA is the preferable choice.
For correlation justification, remember that improvement of inventory position criterion causes negative influence on shortage cost or improvement of service level criterion causes negative influence on inventory in hand and holding cost.In these circumstances, there is no reason for applying variation based methods like entropy, so in this paper PCA is applied for weight assignment.These weights are used by MADM technique (described in Section 4.2.3) to improve quality of ranking.Because of relative utility of scenarios, linear scaleless method employed to scaleless data obtained from simulation.
(2) Calculate sample mean vector;  = ( 1 ,  2 , . . .,   )  and covariance matrix based on ( 7) and ( 8), respectively, (3)  1 /√  is diagonal  ×  matrix, th diagonal element is 1/√  , and sample correlation matrix () is calculated as mentioned in the following equation: (5) Select sufficient number of principal components.In this problem, first,  components whose cumulative percentage of their eigenvalues (  ) is greater than 95% are selected.Cumulative percentage is calculated based on the following equation: (6) Finally, weights vector that is indicated by  is derived from the mathematical relation in ( 13) while   is equal to   /, if the entire elements of PC  are positive, otherwise −  / if entire component are negative In other cases, positive or negative   should be defined in the order that final weight vector is positive.For more detailed discussion about   , refer to [41].Finally, for each product, there are four  parameters that determine relative importance for the following criteria: cost, service level, average of inventory position, and robustness; each element of  vector is divided by the sum of vector elements for normalization.

VIKOR Method.
This method was developed by Opricovic and is the abbreviation of expression that in Serbian means "multiple criteria optimization and compromising solutions" [42].VIKOR method is based on compromising solutions which is the result of Yu [43] and Zeleny [44] studies.This approach considers closeness to ideal solution.
In comparison with other MADM methods, VIKOR is able to present compromising solutions and substitutes, these solutions with best one in the case of necessity [45].VIKOR method is applicable for problems with multiple discrete criteria.Generally, this method can be considered as ranking tool for scenarios which are generated by DOE.Ranking is based on L-P metric function according to (14) with definition of   as the value of th criterion for th scenario.In construction of L-P metric function,  *  indicates best value of all scenarios in the view point of th criterion and  −  represents the worst value of all scenarios with consideration of th criterion.  is the weight parameter for th criterion and is derived by PCA method.Consider VIKOR method only involves  1 and  ∞ which are known as   and   , respectively. and  index refer to scenarios so both of these functions are regarded for each scenario.VIKOR steps are as follows.
(1) Calculate  −  and  *  for entire of criteria.(2) Calculate   and   for entire of scenarios.
(3) With the assumption that  * and  * are the smallest values among all of   and   , respectively, and also  − and  − are the greatest values in comparison with all of   and   , respectively,   is calculated for th scenario based on (15).In this equation, V is strategic weight that determines the importance of individual criterion (  ) against group importance of criteria (  ) for each scenario.Consider (4) Scenarios should be sorted in descending order based on , , and  specifically that result in three different ranking lists.
(5) According to , the scenario which is indicated by   is the best scenario, if two conditions, namely, (C1) and (C2) are satisfied.
(C1).Acceptable Advantage.This condition implies that the second best solution (  ) should be far enough from the best solution (  ) to be accepted as unique solution.The mathematical equation is defined in (16), considering  as the total number of scenarios.Consider (C2).Acceptable Stability.This condition validates   as the best stable solution if   is also the best solution based on  or .
In case that one or both of the mentioned conditions are not satisfied, compromised solutions are involved.VIKOR method suggests that in such case compromising solutions have equal value to the best solution and could be considered interchangeably by decision makers.
If (C1) is not satisfied, scenarios   ,   , . . .,   are considered.  is the best solution and   , . . .,   are compromising solutions with priority of . on the other hand,   ,   are considered if (C2) is not satisfied. superscript is determined by ( 17)

Numerical Result
Simulation process was executed on laptop with 1.8 GH CPU and 4 GB of RAM.Running of each replication with animation and maximum speed takes 7 minutes and 35 seconds.
As simulation of each scenario contains 10 replications and there are 27 scenarios, the total simulation time will be 2046 minutes and 36 seconds.

Multiresolution Result for Demand Modelling.
Result of demand modelling is presented for product  because this analysis for other products is similar.For this problem, polynomial functions are applied for yearly and monthly demand rates.Equations ( 18) and ( 19) are estimated by nonlinear regression model.In these functions   and   are time variables and 95% confidence interval is presented in Table 4 for each coefficient of (18).In this table, coefficients of (18) are visible in the first column while lower bund and upper bund of each coefficient are presented in the second and third columns, respectively.There is no interval including zero, so it is concluded that coefficient estimation is valid.
Estimation error for (  ) is reported based on the sum of square errors (SSE) which is 1.504 × 10 −6 .Consider Monthly demand rate function that is indicated by (  ) and its coefficients confidence interval are presented in Table 5. SEE is equal to 0.0005825.Based on local business information, total demand of product  is estimated to be about 6822 units in five years.According to the information about product , Figure 2 shows final multiresolution result for this product.

Simulation Parameters and Configuration.
For simulation of problem, it is necessary to configure inventory policy parameters and variables.We have three different policies, namely, FQ, FI, and DF.Variables of FQ policy are reorder point (  ), lead time (), and economic quantity orders (  ).First two variables are configured based on Table 2 and third variable is derived from Wilson formula, according to average Orders quantity is unknown for FI policy and is calculated based on probability distribution function which is obtained by demand simulation of each product for 20 days' time interval.Distribution fitting is performed by input analyser of Arena software and graphical result for product  is shown in Figure 3. Results for other products are reported in Table 6.
In this table, fitted probability distribution function for each product is presented in second column while  value of fitting which is greater than 0.15 for fitted distribution grantees goodness of fitting in 95% confidence level.Finally, according to fitted distribution functions, for FI policy, proper ordering amount of each product is presented in the fourth column.
DF policy has two unknown parameters, namely,  and .The first parameter is weight of current real demand in (4) and second one is the weight of current trend in (5).These parameters are adjusted by simulation of five years demand.Demand for the entire four products simultaneously is considered and an average of lost sale percentage is considered as response variable.Different values of response variables are reported in Table 7.According to the different values of  and , the best values for  and  are 0.9 and 0.1, respectively, which results in minimum lost percentage in average.

Weight Assignment.
As mentioned before, PCA is applied for weight assignment.For this purpose, DOE result is needed in the form of decision matrix that is shown in Table 8.In this table, four criteria: cost, average of inventory position, service level, and robustness are evaluated for each scenario of product .Decision matrix should be scaleless that is performed by linear method.Then, the result of eigenvalues is calculated and shown in Table 9.Also result of principal components is presented in Table 10.
As it is visible in Table 9, cumulative proportion of the first and second components is 0.966 which is greater than 0.95, so the first two components are considered as principal components.Their relative weights are 0.78 and 0.186, respectively.Final weights are normalized and results are reported in Table 11.

Ranking and Selection of Scenarios.
Weights of criteria obtained by PCA method in addition to decision matrix are input of VIKOR method., , and  values are calculated based on ( 14) and ( 15) while strategic weight indicated by V is equal to 0.5, so individual and group utility have same   Scenario number 17: inventory policy is FI,   is equal to 15 products, and lead time is 7 days.

Sensitivity Analysis.
As mentioned in the previous section, all scenarios related to product  prefer DF inventory policy, so it is concluded that DF policy is unique optimal policy while lead time includes values of 3, 5, and 7 in selected scenarios.This situation implies that product  is insensitive to the delivery time.It can be wise to decide lead time about 5 days, but there is no need for strict control on lead time while   should be controlled strictly to prevent inventory in hand from reaching lower than 20 products.
In the case of product , (C1) is satisfied so there is a unique scenario while there are four compromising solution for product .In this situation, based on VIKOR method, compromising solutions have same value while interaction plot can be applied for investigation of existing difference between compromising solutions.As mentioned before, the best scenario suggests FI inventory policy for this product but in compromising solutions there are both FI and DF policies.Lead time varies from 3 to 7 days.So, this confusing situation is resolved by deeper analysis and applying interaction plot.Figure 4 shows interaction plot for DOE factors and average of inventory position as response variables.
Figure 4 shows that with either DF or FI policy, there is no significant difference among different lead times in the view point of inventory position, but DF leads to less average of inventory position.With   equaling 15, there is less sensitivity to different lead times.So, configuring   which is equal to 15 products, the average of inventory position criterion will be robust against lead time variations.
Figure 5 regards cost criterion and shows that FI and DF policies are equal and make least sensitivity related to lead time.On the other hand, three plots of   against lead time are parallel.The parallel plots are interpreted as lack of relation between   and lead time.With   equaling 15, increase in inventory cost could be mitigated when lead time sets to 3 days.
In Figure 6, with 3 days for lead time, there is no significant difference between inventory policies in the view point of service level.  and lead time are independent and the best service level is related to FI policy, which is also very sensitive to   .As   increases, service level will grow.According to the interaction plot analyses, optimal decision for product  is FI inventory policy, 3 days lead time, and 15 or more products for   .
Decision about product  is simple.Based on the information obtained from solutions, it is recommended to employ DF inventory policy with   being equal to 15 products and 3 days for lead time.

Conclusion
In this paper, a simulation optimization framework is proposed for robust optimization of retailer inventory system.This framework is less time consuming in comparison with metaheuristic or SPSA approach, it also provides robust solution for multiple objective functions.In this research also trend and cyclic change of customers demand are considered for more realistic view of problem.Proposed framework consists of discrete event simulation and surrogate model.In this framework, full factorial design of experiment is employed for producing decision space.Based on the three decision factors, 27 different scenarios are produced and simulated.In simulation modelling, multiresolution method is applied for simulation of demand with highly dynamic pattern.Because there are multiple criteria for inventory system performance, VIKOR method is employed as multicriteria decision making technique.In addition, robustness is considered as performance criterion and PCA is used for improving performance of VIKOR method.Finally, developed framework gave the best ranked and compromising solutions, so interaction plot applied for more investigation and sensitivity analysis of obtained solutions.
Due to the novelty of the proposed framework, it is recommended that future studies encompass optimization of inventory system of integrated supply chain.Also improvement of developed framework can be considered as future study.

Figure 2 :
Figure 2: Comparison between real and estimated demand.

Figure 3 :
Figure 3: Distribution function of customer demand for reorder interval of product .

level in 𝑡th planning period; 𝐼 𝑝 𝑡 : inventory position in 𝑡th planning period; 𝑑 𝑖𝑡 : quantity of 𝑖th demand in 𝑡th planning period;
: number of reorder in th planning period;   : 1 if   is less than    , 0 otherwise;   : quantity of th order in th planning period;   ,  ℎ , and   are different costs of inventory system according to Table1.   is inventory level and refers to physical quantity of inventory which is available in retailer while    is the reorder point for organizing of new order.In other words, if inventory level reaches   or less, a new order would be organized with quantity of   .Each organized order reaches the retailer and increases inventory level after passing of lead time ().Considering assumptions and described notations, the following equations are the main objective functions of defined problem: is the position of inventory and includes quantity of on-order inventory in addition to inventory level in planning period.  is binary variable which is one if inventory level is greater than the arrived demand.So, if   is one,   can be satisfied and otherwise
nested design.Specific design varies for each problem and should be selected based on problem condition.Our problem consists of tree effective factors, namely, inventory policy, reorder point, and lead time, these factors have nonlinear effect on objective function values.As the aim of DOE in this research is producing decision space with few factors (in

( 4 )
Solve the following equation, where   is a  ×  identity matrix:

Table 6 :
Estimation of distribution functions of demand and reorders quantity.

Table 7 :
Average of lost sale percentage for different values of  and .

Table 8 :
Decision matrix for product .

Table 9 :
Result of eigenvalues.

Table 10 :
Result of principle components.

Table 11 :
Weights of criteria for different products.Inventory policy is DF with   being equal to 30 products and 7 days lead time.According to decision matrix which is shown in Table8, the performance of this scenario for product  results in 38.3944 products as average of inventory position and expected cost is equal to 35269609 for product life cycle.73% of customer demand will be satisfied and, in the case of increasing demand intensity and variation, 2% decrease is expected in service level, so it will change into 71%.Analyses for other products are similar, so only final results are presented in this paper.
Although the best scenario is number 15, VIKOR method implies that compromising solutions have the same values.
Scenario number 6: inventory policy is DF,   is 15 products, and lead time is 5 days.Scenario number 3: inventory policy is DF,   is 15 products, and lead time is 3 days.Scenario number 7: inventory policy is FI,   is 15 products, and lead time is 7 days.Scenario number 4: inventory policy is FI,   is 15 products, and lead time is 5 days.Product .Scenario number 3 has the best rank for this product.This scenario consists of DF inventory policy,   is 15 products, and lead time is 3 days.This scenario leads to 36.3896 products as average of inventory position with 62648337 expected cost and 61% service level with 1% robustness.Compromising solutions are as follows.Scenario number 6: this scenario is mentioned before as a compromising solution for product C.