<?xml version="1.0" encoding="UTF-8"?><!DOCTYPE article  PUBLIC "-//NLM//DTD Journal Publishing DTD v3.0 20080202//EN" "http://dtd.nlm.nih.gov/publishing/3.0/journalpublishing3.dtd"><article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" dtd-version="3.0" xml:lang="en" article-type="research article"><front><journal-meta><journal-id journal-id-type="publisher-id">AJIBM</journal-id><journal-title-group><journal-title>American Journal of Industrial and Business Management</journal-title></journal-title-group><issn pub-type="epub">2164-5167</issn><publisher><publisher-name>Scientific Research Publishing</publisher-name></publisher></journal-meta><article-meta><article-id pub-id-type="doi">10.4236/ajibm.2021.1110063</article-id><article-id pub-id-type="publisher-id">AJIBM-112873</article-id><article-categories><subj-group subj-group-type="heading"><subject>Articles</subject></subj-group><subj-group subj-group-type="Discipline-v2"><subject>Business&amp;Economics</subject></subj-group></article-categories><title-group><article-title>
 
 
  An Analysis and Validation of Effectiveness on Demand Forecasting Considering the Planned Order Effect&lt;sup&gt;*&lt;/sup&gt;
 
</article-title></title-group><contrib-group><contrib contrib-type="author" xlink:type="simple"><name name-style="western"><surname>Young</surname><given-names>Rhee</given-names></name><xref ref-type="aff" rid="aff1"><sub>1</sub></xref><xref ref-type="corresp" rid="cor1"><sup>*</sup></xref></contrib></contrib-group><aff id="aff1"><label>1</label><addr-line>Department of Industrial Engineering, Keimyung University, Daegu, South Korea</addr-line></aff><pub-date pub-type="epub"><day>15</day><month>10</month><year>2021</year></pub-date><volume>11</volume><issue>10</issue><fpage>1036</fpage><lpage>1051</lpage><history><date date-type="received"><day>14,</day>	<month>September</month>	<year>2021</year></date><date date-type="rev-recd"><day>26,</day>	<month>October</month>	<year>2021</year>	</date><date date-type="accepted"><day>29,</day>	<month>October</month>	<year>2021</year></date></history><permissions><copyright-statement>&#169; Copyright  2014 by authors and Scientific Research Publishing Inc. </copyright-statement><copyright-year>2014</copyright-year><license><license-p>This work is licensed under the Creative Commons Attribution International License (CC BY). http://creativecommons.org/licenses/by/4.0/</license-p></license></permissions><abstract><p>
 
 
  The accuracy of forecasting is guaranteed for applying the time related forecasting method under normal circumstances. However, man-made situations, such as planned order can lead to completely different forecasting result
  s
  .
   The research to enhance the accuracy of forecasting is meaningful in that it can be improved by controlling the details of planned order. In this study, an effect of planned order is introduced and the method to enhance the accuracy of demand forecasting in the planned order system is proposed. The validity of the system is examined by comparing the effectiveness of demand forecasting under the planned order system.
 
</p></abstract><kwd-group><kwd>Demand Forecasting</kwd><kwd> Validity</kwd><kwd> Planned Order Effect</kwd></kwd-group></article-meta></front><body><sec id="s1"><title>1. Introduction</title><p>Product differentiation is one of the business strategies of companies to increase competitiveness by emphasizing that their products are superior to their competitors. However, securing competitiveness through product differentiation or price differentiation has reached its limit. And now service level is the main issue for securing a competitive advantage in the market as technology and price competitiveness enter the leveling stage. In particular, post-sales service management is the face of a company because it occurs at the point of contact with customers. Therefore, companies are placing more weight on consumers’ loyalty to companies and products, and companies that cannot respond to this cannot survive. In fact, there are a lot of companies around the world that fail due to their excellent product competitiveness but lack of after-sales ability, and in the case of the United States and some European countries, if the after sale service are delayed by more than two days, the penalty is aggravated, which is deepening the concerns of companies.</p><p>Service parts in the manufacturing industry have various characteristics, and as many as tens of thousands of parts are required for the final product. Demand forecasting is almost impossible to meet the demand since demand for most parts is unspecified demand caused by accidents or breakdowns. However, prompt service to customers and customer expectations are continuously increasing, and the product storage period is continuously increasing due to the service warranty period. For this reason, excess inventory and shortage of inventory occur frequently, and enormous logistics expenses are wasted in order to cope with urgent demand. Therefore, it is becoming an issue of efficient service parts management in most companies because it is necessary to consider the two goals of reducing inventory cost as well as improving customer service level. Although the level of awareness of the importance of demand forecasting is increasing, few companies invest in developing their own demand forecasting models and training forecasting experts to enhance demand forecasting accuracy.</p><p>Demand forecasting has an important meaning in that it provides the basis for various decision making in the process of creating value added. First, if the company’s management policies including planned production can be properly implemented according to changes in the environment, the material supply and demand plan is regularly implemented and it is possible to respond to changes in customer demand (Riezebos &amp; Zhu, 2014) (Farhangmehr &amp; Brito, 1997). Second, it is possible to efficiently manage the excess inventory cost. Third, planned production is possible, so budgeting can be established well, production costs can be reduced. In addition to product related matters, product price, which is the value that a company provides to customers, has been used most frequently and widely as a means to induce consumer reaction due to its ability to change in a relatively short period of time. Manufacturing companies are implementing a PO (Planned Order) system for giving price discount to customer for the purpose of customer service and increasing sales (Kolter, 2003). In other words, it is a system that supplies service parts at low prices for a certain period of time, creates consumer purchasing power, and increases sales. The reason for implementing the PO system positively is that first, it can reduce the increase in purchase avoidance due to the economic recession, and second, it can prevent the decrease in demand due to seasonal factors and insufficient demand by new products. Recently, regardless of the external environment, the PO system is being adopted throughout the year in terms of sales strategy. As consumers perceive various reasons or motives for a price discount, consumers’ reactions also change (Lichtenstein, Burton, &amp; O’Hara, 1989). Subsequent research has expanded the research area to see what kind of discriminatory responses consumers show when discount reasons such as bulk purchases and inventory disposal are presented in addition to price discounts (Sheng, Parker, &amp; Nakamoto, 2007) (Bobinski, Cox, &amp; Cox, 1996). However, although studies related to consumer response to price discounts are active, studies related to the future order quantity due to the PO system are not considered. The time series data on the order demand for parts is based on four variables, such as trend, cycle, seasonal and irregular, and has data made up of two or more combinations. Under the PO system, the quantity of orders increases in demand, causing statistically outliers intentionally. If demand forecasting is performed with data containing these outliers, the forecast is distorted and the accuracy is lowered. This is because the existing method has a drawback in that the outlier cannot reflected be properly. However, a new method that reflects the effect of PO on outliers is needed, since the inclusion of outliers due to the PO system is a strategy to create planned demand. The inventory management and the opportunity cost management would be easy since the accuracy of forecasting is guaranteed by the demand forecasting method considering the effect of the PO. As a result, it increases the availability of parts and customer service satisfaction. In this study, alternatives are suggested and analyzed by comparing the accuracy of demand forecasting by the PO effect with the accuracy of the ARIMA model. To the best of my knowledge, these kinds of model have not as yet been studied in the open literature.</p><p>The paper is organized as follows. We review related work briefly in Section 2. The demand forecasting model related to service parts is presented in Section 3. Section 4 introduces the PO effect model considering the PO effect. In Section 5, the validity and comparison of the model is examined against the model without PO. Finally, Section 6 gives concluding remarks.</p></sec><sec id="s2"><title>2. Related Studies</title><p>Time series data analysis is a dynamic analysis method that attempts to predict the future by analyzing changes in data from the past to the present using time as an independent variable. The characteristics of time series data are divided into four components, such as trend, cycle, season, and irregularity. Letting Z<sub>t</sub> be the observed time series data at time t, Z<sub>t</sub> is generally expressed as data in multiplicative model or additive model involving the four components. If seasonality exists before implementing time series analysis, time series decomposition, winters model, and ARIMA can be used for demand forecasting. On the other hand, seasonality or trend does not exist, demand forecasting can be made with a simple exponential smoothing method or a moving average method. If there is no seasonality but a trend exists, trend analysis, double exponential smoothing, or ARIMA is used. As one of quantity indices for determining the existence seasonality and of a trend, it is possible to calculate the auto correlation coefficient, which is a correlation coefficient between the average deviations of all possible pairs of data. The estimated auto correlation coefficient means the degree of correlation between time series data, and it can be said to be a statistical relationship expressing a direction and intensity between a pair of observation. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x3.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x4.png" xlink:type="simple"/></inline-formula> mean the observation at time t and the auto correlation coefficient between observations separated by k time lags in the time series respectively, and is expressed as follows.</p><disp-formula id="scirp.112873-formula9"><label>(1)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x5.png"  xlink:type="simple"/></disp-formula><p>In (1), k means a time differences separated by k, and is mainly used to calculate a pair of autocorrelation coefficients.</p><sec id="s2_1"><title>2.1. Box-Jenkins Model (Box &amp; Jenkins, 1976)</title><p>The study on time series analysis started under the assumption that the decomposition method can decompose each of the 4 components of the time series data. However, the reality is that these 4 components can be easily decomposed by the decomposition method. Accordingly, a new approach to time series analysis was devised, and a statistical theoretic system based on ARIMA model is established by generalizing this method. The Box-Jenkins method implements three steps: Model Identification, Parameter Estimation, and Model Diagnostic Checking to find out which probabilistic properties the time series data have and which model is appropriate. Model identification is to find the most appropriate model by statistical comparison after selecting several models that are considered appropriate for the observed time series data. The principle of parsimony should be adhered to in selecting a model. This is to express the observed data as simply as possible by selecting the model with the smallest number of parameters while appropriately expressing the observed data. The parameter estimation and the fitness of good test of the model are performed through statistical procedures, and if assumptions and conditions are not satisfied, return to the model identification step and repeat the above three steps until an appropriate model is found. The assumption of stationarity in stochastic processes and an understanding of various types of autocorrelation functions help to extend the Box-Jenkins model.</p><p>Letting the time series data at an arbitrary time point t is expressed as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x6.png" xlink:type="simple"/></inline-formula>, it can be regarded as a random variable obtained from the population. In general, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x7.png" xlink:type="simple"/></inline-formula>which is a set of random variables, is called a stochastic process, and is also called a time series process. In statistical estimation and testing for time series models, one of the crucial assumptions for simplifying analysis is stationarity. In a stationary time series, the mean, (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x8.png" xlink:type="simple"/></inline-formula>) and variance, (<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x9.png" xlink:type="simple"/></inline-formula>) of time series data generated in a stochastic process are constant over time. And auto covariance and auto correlation functions are time-dependent on time lag k, not dependent on time t.</p><p>Correlation functions used in the model diagnosis step include ACF (autocorrelation function), PACF (part autocorrelation function), and SACF (sample autocorrelation function). When two time points with a time difference k are t and t + k, the expected value of the product of the sample values of the stochastic process <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x10.png" xlink:type="simple"/></inline-formula> is called the autocorrelation function of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x11.png" xlink:type="simple"/></inline-formula>.</p><p>If the covariance of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x12.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x13.png" xlink:type="simple"/></inline-formula> is expressed as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x14.png" xlink:type="simple"/></inline-formula>. <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x15.png" xlink:type="simple"/></inline-formula>is called an auto-covariance function since <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x16.png" xlink:type="simple"/></inline-formula> is a function of k. Therefore, ACF, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x17.png" xlink:type="simple"/></inline-formula>is expressed as (2), where<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x18.png" xlink:type="simple"/></inline-formula>.</p><disp-formula id="scirp.112873-formula10"><label>(2)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x19.png"  xlink:type="simple"/></disp-formula><p>The PACF (partial autocorrelation function) is an autocorrelation function just only for <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x20.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x21.png" xlink:type="simple"/></inline-formula>, and is expressed as a conditional auto-correlation function as shown in (3), after removing the mutual linear dependence of<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x22.png" xlink:type="simple"/></inline-formula>, the random variables of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x23.png" xlink:type="simple"/></inline-formula> from time <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x20.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x21.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x22.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x23.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x24.png" xlink:type="simple"/></inline-formula> to time t +k − 1.</p><disp-formula id="scirp.112873-formula11"><label>(3)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x25.png"  xlink:type="simple"/></disp-formula><p>As another approach, the PACF, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x26.png" xlink:type="simple"/></inline-formula>is a function of time differencek and can be obtained by considering the regression model. That is, if the dependent variable <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x27.png" xlink:type="simple"/></inline-formula> is regression-analyzed against k random variables, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x28.png" xlink:type="simple"/></inline-formula>, n a normal process with a mean of 0, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x29.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x30.png" xlink:type="simple"/></inline-formula> is the PACF. If <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x31.png" xlink:type="simple"/></inline-formula> is given for the stationary time series<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x32.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x26.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x27.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x28.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x29.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x30.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x31.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x32.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x33.png" xlink:type="simple"/></inline-formula>becomes the value of the PACF in the solution of the equation with the time difference k.</p><p>The SACF (sample autocorrelation function) can be used by substituting <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x34.png" xlink:type="simple"/></inline-formula> instead of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x35.png" xlink:type="simple"/></inline-formula> to estimate<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x34.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x35.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x36.png" xlink:type="simple"/></inline-formula>. The estimated SACF is applied by a reduced model using an algorithm instead of a complex determinant as in (4).</p><disp-formula id="scirp.112873-formula12"><label>(4)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x37.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x38.png" xlink:type="simple"/></inline-formula></p></sec><sec id="s2_2"><title>2.2. ARIMA Model</title><p>The ARIMA model is a generalization of the ARMA (autoregressive moving average) and predicts a given time series based on its own past values. It can be used for any non-seasonal series of numbers that exhibits patterns and is not a series of random events. The ARIMA model is suitable for short-term demand forecasting since it gives more weight to observations close to the most recent time. The ARMA model is a combined model that has both the parameters of the AR (autoregressive) model and the MA (moving average) model.</p><p>The basic model applied to time series analysis is the AR model, and the AR(1) model with 1 autoregressive parameter is defined as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x39.png" xlink:type="simple"/></inline-formula> where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x40.png" xlink:type="simple"/></inline-formula> is a random variable in the stationary time series, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x41.png" xlink:type="simple"/></inline-formula>is autoregressive parameter of an order 1, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x42.png" xlink:type="simple"/></inline-formula> is the probability error in time series data at timet. The AR(2) model has 2 parameters of order 1 and order 2 such as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x43.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x44.png" xlink:type="simple"/></inline-formula>, expressed as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x39.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x40.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x41.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x42.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x43.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x44.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x45.png" xlink:type="simple"/></inline-formula> by generalizing this concept, an AR(p) model with p parameters can be expressed as (5).</p><disp-formula id="scirp.112873-formula13"><label>(5)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x46.png"  xlink:type="simple"/></disp-formula><p>The second type of Box-Jenkins model is the MA (moving average) model. The basic MA(1) model with 1 parameter is expressed as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x47.png" xlink:type="simple"/></inline-formula>. Here, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x48.png" xlink:type="simple"/></inline-formula>is the fitted value fitted to the series value<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x49.png" xlink:type="simple"/></inline-formula>, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x50.png" xlink:type="simple"/></inline-formula> is called the MA parameter of an order 1. The MA model is almost identical to the AR model in appearance, but the meaning is completely different. In other words, the MA parameter is related only to the probability errors, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x51.png" xlink:type="simple"/></inline-formula>, that the series values occurred before time t. The MA model is also generalized like the AR model, and the general MA(p) model with p parameters is expressed as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x52.png" xlink:type="simple"/></inline-formula>. The time series value <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x47.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x48.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x49.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x50.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x51.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x52.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x53.png" xlink:type="simple"/></inline-formula> in the general MA model is expressed as the sum of the p probability errors at the previous p time period and the probability errors at the current time. Therefore, the p-order MA(p) model can be represented as (6).</p><disp-formula id="scirp.112873-formula14"><label>(6)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x54.png"  xlink:type="simple"/></disp-formula><p>Since it is represented as a model of order q with only one MA parameter of order q, it has fewer than q MA parameters and order terms lower than q are excluded. The MA model is derived simply as a weighted average of the probability errors of the previous shifted time as time t increases.</p><p>The ARMA (p, q) (autoregressive moving average) model has both p AR parameters and q MA parameters, and its equation is expressed as (7).</p><disp-formula id="scirp.112873-formula15"><label>(7)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x55.png"  xlink:type="simple"/></disp-formula><p>When a stationary time series is given, the criteria for determining whether the time series data to be analyzed are suitable for an AR model or a MA model depends on the value of the ACF and the value of the PACF mentioned in section 2.1. The relationships between ACF and PACF of each AR process, MA process, and ARMA process as follows; the ACF value is exponentially decreases as the time difference k increases in the AR (p) process. And the ACF value remains 0 at time difference q + 1 in the MA (q) process. On the other hand, The PACF value is exponentially decreases as the time difference k increases in the MA (q) process. And the PACF value remains 0 at time difference p + 1 in the AR (p) process. Therefore, in the ARMA (p, q) process, both the ACF value and the PACF value decrease as the time difference k increases. Since the Box-Jenkins identification method does not identify the orders of the ARMA model, it is possible by applying a statistical estimation method.</p><p>The ARIMA (p, d, q) model is a model in which the concept of differencing is inserted into the ARMA (p, q) model, where differencing refers to calculating the difference between consecutively measured observations. The nth differencing order calculates the difference between every nth observation. Determination of the differencing order, d helps to make non-stationary time series into a stationary time series. When the order of auto-regression is p, the order of the moving average is q, and the differencing order is d, the ARIMA (p, d, q) model can be expressed as (8).</p><disp-formula id="scirp.112873-formula16"><label>(8)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x56.png"  xlink:type="simple"/></disp-formula><p>where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x57.png" xlink:type="simple"/></inline-formula> and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x57.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x58.png" xlink:type="simple"/></inline-formula></p><p>In (8), B is a backshift operator and is defined as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x59.png" xlink:type="simple"/></inline-formula>. In the ARIMA (p, d, q) model, if p = 0, d = 1, q = 0, it is called a random walk, and it becomes <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x60.png" xlink:type="simple"/></inline-formula> or<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x59.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x60.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x61.png" xlink:type="simple"/></inline-formula>.</p><p>In the analysis of time series data having seasonality, if only the term corresponding to seasonality is required for the model, a model is constructed using a partial model. When AR operator or MA operator is expressed in the form of product of non-seasonal operator and seasonal operator, multiplicative model (ARIMA (p, d, q) &#215; (P, D, Q) 12) is used. If the order of the time series model is (p, d, q), the seasonal period is s, and the order of the seasonal time series model is (P, D, Q), the multiplicative seasonal ARIMA model is expressed as (9).</p><disp-formula id="scirp.112873-formula17"><label>(9)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x62.png"  xlink:type="simple"/></disp-formula><p>In (9), <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x63.png" xlink:type="simple"/></inline-formula>and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x64.png" xlink:type="simple"/></inline-formula> represent seasonal auto-regression and seasonal moving average polynomials respectively. Accordingly, autoregressive and moving average polynomials are expressed as <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x65.png" xlink:type="simple"/></inline-formula> and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x63.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x64.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x65.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x66.png" xlink:type="simple"/></inline-formula>.</p></sec></sec><sec id="s3"><title>3. Demand Forecasting Using ARIMA Model</title><p>The PO system has been mainly applied as a marketing strategy for price discounts as customer service and sales increase by manufacturers. It cannot be viewed as a time series data because the demand increases abnormally since the PO system is a strategy to create customer orders at a specific point in time.</p><p>In this study, the demand is forecasted by applying the actual service parts orders for 6 years of Daedong Engineering in Daegu as a time series data. Thousands of service parts are forecasted per month, and demand forecasting is simply based on experience in consideration of past demand data, seasonal characteristics, and PO implementation. By grouping tens of thousands of parts, five representative parts are selected and applied to demand forecasting analysis. The monthly data of service parts for 70 months from January 2013 to October 2018 are used for the demand forecasting. The order quantity of filter parts including PO is shown in <xref ref-type="fig" rid="fig1">Figure 1</xref>.</p><p>In addition to the Filter part shown as an example, the rest of the parts also show similar demand trends. In <xref ref-type="fig" rid="fig1">Figure 1</xref>, it can be inferred that outliers and slope changes may exist in time series data. The first step in model identification to find an appropriate model is to determine the differencing order of ARIMA (p, d, q). Characteristics of time series data can be identified through ACF estimation and PACF estimation and PACF is used to determine the differencing order of the AR model, and ACF is used to determine the differencing order of the MA model.</p><p>The autocorrelation function and the partial autocorrelation function of time series data for the Filter part are shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>. In order to determine the</p><fig id="fig1"  position="float"><label><xref ref-type="fig" rid="fig1">Figure 1</xref></label><caption><title> Filter order quantity</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2122170x67.png"/></fig><fig id="fig2"  position="float"><label><xref ref-type="fig" rid="fig2">Figure 2</xref></label><caption><title> ACF and PACF</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2122170x68.png"/></fig><p>presence of seasonality, a sharp slope change such as the cut shape of the ACF cannot be found at the periodic time difference is 12, and there is no appearance where the signs are all positive and the peaks are repeated. From these observations, it can be said that there is no seasonal effect in Filter part. As shown in <xref ref-type="fig" rid="fig2">Figure 2</xref>, the estimation of the correlation function is set at the 5% significant level, it can be seen that the time difference is 3 at the point outside the confidence limit line indicated by the dotted line. Therefore, the first attempts to identify the model are AR (3) and MA (3). The general formula of the AR (3) model is<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x69.png" xlink:type="simple"/></inline-formula>.</p><p><xref ref-type="table" rid="table1">Table 1</xref> shows the estimation results of parameters on AR (3). It satisfies the stationary condition since the absolute values of the estimated autocorrelation coefficients are<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x70.png" xlink:type="simple"/></inline-formula>, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x71.png" xlink:type="simple"/></inline-formula>, and<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x72.png" xlink:type="simple"/></inline-formula>. The auto-regression coefficient, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x70.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x71.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x72.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x73.png" xlink:type="simple"/></inline-formula>should be included in this model.</p><p>After identifying the model, the residual ACF and the residual PACF are used to diagnose the model. Model selection criteria include<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x74.png" xlink:type="simple"/></inline-formula>, residual autocorrelation function test statistic, Ljung-Box test statistic, and residual chart. The residual autocorrelation function is calculated based on the residuals obtained by estimating the AR (1) model. In this case, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x75.png" xlink:type="simple"/></inline-formula>is an estimate of the<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x74.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x75.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x76.png" xlink:type="simple"/></inline-formula>, unobservable white noise. These white noises are assumed to be statistically independent.</p><p>As seen in <xref ref-type="fig" rid="fig3">Figure 3</xref>, none of the values of the residual ACF protrudes outside the confidence limit with a 5% significant level. The absence of a residual PACF value protruding outside the confidence interval indicates that the addition of an autoregressive factor is no longer necessary.</p></sec><sec id="s4"><title>4. Demand Forecasting with PO Effect</title><p>The research for demand forecasting is to find a way to reduce the difference between the actual value and the forecast value, and forecasting accuracy depends on how the forecasting effect is set and reflected in the model. Since the PO system is a policy in which a manufacturer strategically creates a customer’s order quantity at a specific point in time, the application of the existing normal</p><table-wrap id="table1" ><label><xref ref-type="table" rid="table1">Table 1</xref></label><caption><title> AR (3) parameter estimation</title></caption><table><tbody><thead><tr><th align="center" valign="middle"  colspan="5"  >Parameter estimation</th></tr></thead><tr><td align="center" valign="middle" >type</td><td align="center" valign="middle" >coefficient</td><td align="center" valign="middle" >SE coefficient</td><td align="center" valign="middle" >T</td><td align="center" valign="middle" >P</td></tr><tr><td align="center" valign="middle" >AR (1)</td><td align="center" valign="middle" >−0.0010</td><td align="center" valign="middle" >0.1179</td><td align="center" valign="middle" >−0.01</td><td align="center" valign="middle" >0.993</td></tr><tr><td align="center" valign="middle" >AR (2)</td><td align="center" valign="middle" >−0.1696</td><td align="center" valign="middle" >0.1168</td><td align="center" valign="middle" >−1.45</td><td align="center" valign="middle" >0.151</td></tr><tr><td align="center" valign="middle" >AR (3)</td><td align="center" valign="middle" >−0.2903</td><td align="center" valign="middle" >0.1191</td><td align="center" valign="middle" >−2.44</td><td align="center" valign="middle" >0.017</td></tr><tr><td align="center" valign="middle" >constant</td><td align="center" valign="middle" >4920.1</td><td align="center" valign="middle" >293.7</td><td align="center" valign="middle" >16.75</td><td align="center" valign="middle" >0.000</td></tr><tr><td align="center" valign="middle" >mean</td><td align="center" valign="middle" >3367.6</td><td align="center" valign="middle" >210</td><td align="center" valign="middle" ></td><td align="center" valign="middle" ></td></tr><tr><td align="center" valign="middle" >No. of obs.</td><td align="center" valign="middle"  colspan="4"  >70</td></tr><tr><td align="center" valign="middle"  colspan="5"  >residual SE (sample error) = 3,982,968.3 MS (mean square) = 6,034,800.0 DF(degree of freedom) = 66.0</td></tr></tbody></table></table-wrap><fig id="fig3"  position="float"><label><xref ref-type="fig" rid="fig3">Figure 3</xref></label><caption><title> Residual ACF and Residual PACF</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2122170x77.png"/></fig><p>demand forecasting method raises questions about the accuracy of demand forecasting. In this study, the demand forecasting method that reflects the PO effect is proposed by modifying the time series analysis method under the PO system. In other words, if PO occurs at<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x78.png" xlink:type="simple"/></inline-formula>, the PO effect is calculated at each time and reflected in the forecasting amount calculation. By doing so, it is expected that the accuracy of demand forecasting will be improved by the PO effect.</p><p>In this study, although the outlier limit is usually defined as the inter-quartile range, the average monthly order quantity <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x79.png" xlink:type="simple"/></inline-formula> and the standard deviation <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x79.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x80.png" xlink:type="simple"/></inline-formula> for a certain period(window) of n months are obtained to set the outlier limit using the time series data<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x79.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x81.png" xlink:type="simple"/></inline-formula>. For simplicity, the values outside of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x79.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x82.png" xlink:type="simple"/></inline-formula> are considered as outliers. The average <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x79.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x83.png" xlink:type="simple"/></inline-formula> of the remaining data in the period is calculated after removing the data identified as outliers. And as shown in (10), the value obtained by dividing the difference between the actual value <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x79.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x84.png" xlink:type="simple"/></inline-formula> and the average monthly order quantity by <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x79.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x80.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x81.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x82.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x83.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x84.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x85.png" xlink:type="simple"/></inline-formula> is defined as the PO effect.</p><disp-formula id="scirp.112873-formula18"><label>(10)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x86.png"  xlink:type="simple"/></disp-formula><p>As shown in (10), the PO effect can take various modifications, such as a conservative approach and a sensitive approach, considering the characteristics of time series data. In other words, various alternatives can be suggested by adjusting the outlier limit or by changing the time period n. The adjusted demand forecast <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x87.png" xlink:type="simple"/></inline-formula> reflecting the PO effect is shown as (11), where <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x87.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x88.png" xlink:type="simple"/></inline-formula> represents the demand forecast at time t.</p><disp-formula id="scirp.112873-formula19"><label>(11)</label><graphic position="anchor" xlink:href="http://html.scirp.org/file/2-2122170x89.png"  xlink:type="simple"/></disp-formula><p>For the example given in Section 3, the process of finding outliers is carried out by setting the PO period n to be 12-month (1 yr.) and 60-month (5 yrs.). In <xref ref-type="fig" rid="fig4">Figure 4</xref>, the time series data for two cases of 1 year (upper) and 5 years (lower) is shown, the outliers are the values outside of <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x90.png" xlink:type="simple"/></inline-formula> and are indicated by circles.</p><fig id="fig4"  position="float"><label><xref ref-type="fig" rid="fig4">Figure 4</xref></label><caption><title> Filter outliers</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2122170x91.png"/></fig><p>In order to apply the PO effect proposed in (10) and (11) to time series data, the 12-month average and the correction value replacing the PO effect are applied after removing the order quantity at the time of PO. The time series data to which the PO effect is applied is shown in <xref ref-type="fig" rid="fig5">Figure 5</xref>.</p><p>As seen in <xref ref-type="fig" rid="fig5">Figure 5</xref>, the monthly time series order quantity has a different scale and trend compared to <xref ref-type="fig" rid="fig1">Figure 1</xref>, which shows the time series data before the revision. In the next step, ACF and PACF for model identification are examined. (2) and (3) were used to calculate ACF and PACF for each time difference, and <xref ref-type="fig" rid="fig6">Figure 6</xref> shows ACF and PACF in relation to the time difference.</p><p>As seen in <xref ref-type="fig" rid="fig6">Figure 6</xref>, ACF and PACF for each time difference are indicated by a bar chart, and the confidence limit applied with a 5% significant level is indicated by a dotted line. It can be said that there is no seasonal effect since ACF is not beyond the confidence limit at time difference 12. Therefore, the forecasting order quantity is obtained using AR (1) and MA (1) according to the principle of parsimony. And the modified forecasting amount is obtained by applying the PO effect.</p><p>In this study, from the perspective of a conservative approach, the PO effect reflecting the average of 60-month is performed with the same procedures as the one presented before. After removing the order quantity at the time of PO, the 60-month average and the correction value substituted for the PO effect are applied. The time series data to which the PO effect is applied is shown in <xref ref-type="fig" rid="fig7">Figure 7</xref>.</p><p>Comparing <xref ref-type="fig" rid="fig5">Figure 5</xref> and <xref ref-type="fig" rid="fig7">Figure 7</xref>, even if the period for reflecting the PO effect in the demand is different, the overall ordering quantity trend is similar, but a slight difference can be found.</p><p>ACF and PACF were examined for the purpose of model identification. As shown in <xref ref-type="fig" rid="fig8">Figure 8</xref>, since the time difference 12 does not deviate from the confidence limit, it can be said that there is no trend and seasonal effect as in the</p><fig id="fig5"  position="float"><label><xref ref-type="fig" rid="fig5">Figure 5</xref></label><caption><title> Ordering quantity with 12-month PO period</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2122170x92.png"/></fig><fig id="fig6"  position="float"><label><xref ref-type="fig" rid="fig6">Figure 6</xref></label><caption><title> ACF and PACF with PO effect</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2122170x93.png"/></fig><fig id="fig7"  position="float"><label><xref ref-type="fig" rid="fig7">Figure 7</xref></label><caption><title> Ordering quantity with 60-month PO period</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2122170x94.png"/></fig><fig id="fig8"  position="float"><label><xref ref-type="fig" rid="fig8">Figure 8</xref></label><caption><title> ACF and PACF with 60-month PO period</title></caption><graphic mimetype="image"   position="float"  xlink:type="simple"  xlink:href="http://html.scirp.org/file/2-2122170x95.png"/></fig><p>original time series data. The result is the same as the result in <xref ref-type="fig" rid="fig6">Figure 6</xref>. However, a closer look reveals a difference between the two figures. There is a difference in the PACF value at each time difference, which can be said to have different influences in determining the order of the AR model.</p></sec><sec id="s5"><title>5. Comparisons and Validations</title><p>After comparative analysis of various suggested models, model selection criteria include<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x96.png" xlink:type="simple"/></inline-formula>, residual autocorrelation function test statistic, and residual chart. After selecting the model, the accuracy of forecasting is performed based on the error between the estimated value and the actual value. Among the various scales for accuracy comparison, in this study, well known RMSE (root mean square error) and MAPE (mean absolute percentage error) are applied. RMSE is a frequently used measure of the differences between an estimator and the values observed.</p><p>Letting k be the time difference, <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x97.png" xlink:type="simple"/></inline-formula>be the order quantity, and <inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x98.png" xlink:type="simple"/></inline-formula> be the demand forecast quantity, RMSE is expressed as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x99.png" xlink:type="simple"/></inline-formula>, and MAPE as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x97.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x98.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x99.png" xlink:type="simple"/></inline-formula><inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x100.png" xlink:type="simple"/></inline-formula>.</p><p>In this study, the accuracy of demand forecasting is reviewed for disk, gasket, cable and other general items in addition to the filter. As shown in <xref ref-type="table" rid="table2">Table 2</xref>, the forecasts obtained with PO show better results than those without PO by examining RMSE and MAPE, which are indicators of the accuracy of the forecasts. In addition, although the policy is normally implemented as a PO period of 12-month, a policy with 60-month PO period was proposed and compared with the results of 12-month PO period. In comparisons of the forecasting accuracy results by 12-month PO period and the results of 60-month PO period, the superiority of the policy could not be decided. This is because the characteristic of ARIMA is that the more the data goes into the past, the less influence to the forecast. The fact that the results of the 12-month PO period are somewhat superior can be a result of reflecting this. As forecasting accuracy indicators, there are MAE, MSE, and MPE and etc. in addition to RMSE and MAPE, but they were not investigated in this study. However, it is expected that the analysis results for other forecasting accuracy indicators would be similar.</p><table-wrap id="table2" ><label><xref ref-type="table" rid="table2">Table 2</xref></label><caption><title> RMSE and MAPE (%) comparison</title></caption><table><tbody><thead><tr><th align="center" valign="middle" ></th><th align="center" valign="middle" >Parts</th><th align="center" valign="middle" >Without PO</th><th align="center" valign="middle" >12-month PO period</th><th align="center" valign="middle" >60-month PO period</th></tr></thead><tr><td align="center" valign="middle"  rowspan="5"  >RMSE</td><td align="center" valign="middle" >Filter</td><td align="center" valign="middle" >129.27</td><td align="center" valign="middle" >89.34</td><td align="center" valign="middle" >96.39</td></tr><tr><td align="center" valign="middle" >Disk</td><td align="center" valign="middle" >138.43</td><td align="center" valign="middle" >99.27</td><td align="center" valign="middle" >93.26</td></tr><tr><td align="center" valign="middle" >General</td><td align="center" valign="middle" >134.27</td><td align="center" valign="middle" >87.48</td><td align="center" valign="middle" >109.43</td></tr><tr><td align="center" valign="middle" >Gasket</td><td align="center" valign="middle" >138.35</td><td align="center" valign="middle" >112.34</td><td align="center" valign="middle" >108.37</td></tr><tr><td align="center" valign="middle" >Cable</td><td align="center" valign="middle" >143.86</td><td align="center" valign="middle" >109.26</td><td align="center" valign="middle" >120.74</td></tr><tr><td align="center" valign="middle"  rowspan="5"  >MAPE (%)</td><td align="center" valign="middle" >Filter</td><td align="center" valign="middle" >62.21</td><td align="center" valign="middle" >40.97</td><td align="center" valign="middle" >46.31</td></tr><tr><td align="center" valign="middle" >Disk</td><td align="center" valign="middle" >98.65</td><td align="center" valign="middle" >70.80</td><td align="center" valign="middle" >72.13</td></tr><tr><td align="center" valign="middle" >General</td><td align="center" valign="middle" >109.16</td><td align="center" valign="middle" >92.15</td><td align="center" valign="middle" >110.40</td></tr><tr><td align="center" valign="middle" >Gasket</td><td align="center" valign="middle" >121.94</td><td align="center" valign="middle" >99.47</td><td align="center" valign="middle" >92.04</td></tr><tr><td align="center" valign="middle" >Cable</td><td align="center" valign="middle" >128.27</td><td align="center" valign="middle" >76.23</td><td align="center" valign="middle" >80.81</td></tr></tbody></table></table-wrap></sec><sec id="s6"><title>6. Conclusion</title><p>Customer service, inventory costs, and company credit are interrelated and demand forecasting is at the center of them. The purpose of this study is to propose a method to improve the accuracy of demand forecasting in the presence of a PO system. The demand can be forecasted by applying ARIMA, but a relatively large error occurs in the forecast amount when PO occurs. Therefore, these shortcomings make it inconvenient to apply ARIMA directly to time series data with planned uncertainty.</p><p>In this study, the demand forecasting method reflecting the PO effect was proposed as one of the methods to improve the forecasting accuracy under the PO system. The time at which the outlier occurred was regarded as the time at which the PO occurred, and the criterion for the outlier was set as<inline-formula><inline-graphic xlink:href="http://html.scirp.org/file/2-2122170x101.png" xlink:type="simple"/></inline-formula>. It is believed that the PO effect can be further enhanced by lowering the threshold for outliers, and the accuracy of the forecasting can also be improved. On the other hand, the PO effect can vary with the PO period. Although 12 months is mainly applied as a PO period, a 60-month PO period is also considered with a conservative approach in this paper. The comparison of the results by these two PO periods cannot be definitive because it depends on the characteristics of the data, but the short-term PO effect seems better. This is because the characteristics of ARIMA have less effect on forecasting as the data goes into the past. The extension of this study can be possible in two directions: lowering the criteria for outliers and adjusting the PO period.</p></sec><sec id="s7"><title>Conflicts of Interest</title><p>The author declares no conflicts of interest regarding the publication of this paper.</p></sec><sec id="s8"><title>Cite this paper</title><p>Rhee, Y. (2021). An Analysis and Validation of Effectiveness on Demand Forecasting Considering the Planned Order Effect. American Journal of Industrial and Business Management, 11, 1036-1051. https://doi.org/10.4236/ajibm.2021.1110063</p></sec><sec id="s9"><title>NOTES</title></sec></body><back><ref-list><title>References</title><ref id="scirp.112873-ref1"><label>1</label><mixed-citation publication-type="other" xlink:type="simple">Sheng, S., Parker, A.M., &amp; Nakamoto, K. (2007). The Effects of Price Discount and Product Complementarity on Consumer Evaluations of Bundle Components. Journal of Marketing Theory and Practice, 15, 53-64. https://doi.org/10.2753/MTP1069-6679150104</mixed-citation></ref><ref id="scirp.112873-ref2"><label>2</label><mixed-citation publication-type="other" xlink:type="simple">Riezebos, J., &amp; Zhu, S.X. (2014). MRP Planned Orders in a Multiple-Supplier Environment with Differing Lead Times. Production and Operation Management, 24, 833-895. https://doi.org/10.1111/poms.12318</mixed-citation></ref><ref id="scirp.112873-ref3"><label>3</label><mixed-citation publication-type="other" xlink:type="simple">Lichtenstein, D.R., Burton, S., &amp; O’Hara, B.S. (1989). Marketplace Attributions and Consumer Evaluations of Discount Claims. Psychology &amp; Marketing, 6, 163-180. https://doi.org/10.1002/mar.4220060302</mixed-citation></ref><ref id="scirp.112873-ref4"><label>4</label><mixed-citation publication-type="other" xlink:type="simple">Kolter, P. (2003). Marketing Management (11th ed.). Prentice Hall.</mixed-citation></ref><ref id="scirp.112873-ref5"><label>5</label><mixed-citation publication-type="journal" xlink:type="simple"><name name-style="western"><surname>Farhangmehr</surname><given-names> M.</given-names></name>,<name name-style="western"><surname> &amp; Brito</surname><given-names> P. </given-names></name>,<etal>et al</etal>. (<year>1997</year>)<article-title>. Sales Promotions: Prescriptive Definitions and the Managers View</article-title><source> Proceedings of the 1997 World Marketing Congress</source><volume> 1</volume>,<fpage> 45</fpage>-<lpage>49</lpage>.<pub-id pub-id-type="doi"></pub-id></mixed-citation></ref><ref id="scirp.112873-ref6"><label>6</label><mixed-citation publication-type="other" xlink:type="simple">Box, G., &amp; Jenkins, G. (1976). Time Series Analysis: Forecasting Control (2nd ed.). San Francisco: Holden Day.</mixed-citation></ref><ref id="scirp.112873-ref7"><label>7</label><mixed-citation publication-type="other" xlink:type="simple">Bobinski, G., Cox, D., &amp; Cox, A. (1996). Retail ‘Sale’ Advertising, Perceived Retailer Credibility and Price Rationale. Journal of Retailing, 72, 291-306. https://doi.org/10.1016/S0022-4359(96)90031-7</mixed-citation></ref></ref-list></back></article>