Introduction to
(Demand) Forecasting
Topics
• Introduction to (demand) forecasting
• Overview of forecasting methods
• A generic approach to quantitative forecasting
• Time series-based forecasting
• Building causal models through multiple linear
regression
• Confidence Intervals and their application in forecasting
Forecasting
• The process of predicting the values of a certain quantity,
Q, over a certain time horizon, T, based on past trends
and/or a number of relevant factors.
• Some forecasted quantities in manufacturing
– demand
– equipment and employee availability
– technological forecasts
– economic forecasts (e.g., inflation rates, exchange rates, housing
starts, etc.)
• The time horizon depends on
– the nature of the forecasted quantity
– the intended use of the forecast
Forecasting future demand
• Demand forecasting is based on:
– extrapolating to the future past trends observed in the company
sales;
– understanding the impact of various factors on the company future
sales:
• market data
• strategic plans of the company
• technology trends
• social/economic/political factors
• environmental factors
• etc
• Remark: The longer the forecasting horizon, the more
crucial the impact of the factors listed above.
Demand Patterns
• The observed demand is the cumulative result of:
– systematic variation, due to a number of identified factors, and
– a random component, incorporating all the remaining unaccounted
effects.
• Patterns of systematic variation
– seasonal: cyclical patterns related to the calendar (e.g., holidays,
weather)
– cyclical: patterns related to changes of the market size, due to, e.g.,
economics and politics
– business: patterns related to changes in the company market share,
due to e.g., marketing activity and competition
– product life cycle: patterns reflecting changes to the product life
The problem of demand forecasting
– Identify and characterize the systematic
variation, as a set of trends.
– Characterize the variability in the demand.
Forecasting Methods
• Qualitative (Subjective): Incorporate factors like
the forecaster’s intuition, emotions, personal
experience, and value system.
• These methods include:
– Jury of executive opinion
– Sales force composites
– Delphi method
– Consumer market surveys
Forecasting Methods (cont.)
• Quantitative (Objective): Employ one or more
mathematical models that rely on historical data
and/or causal/indicator variables to forecast
demand.
• Major methods include:
– time series methods: F(t+1) = f (D(t), D(t-1), …)
– causal models: F(t+1) = f(X1(t), X2(t), …)
Selecting a Forecasting Method
• It should be based on the following considerations:
– Forecasting horizon (validity of extrapolating past data)
– Availability and quality of data
– Lead Times (time pressures)
– Cost of forecasting (understanding the value of
forecasting accuracy)
– Forecasting flexibility (amenability of the model to
revision; quite often, a trade-off between filtering out
noise and the ability of the model to respond to abrupt
and/or drastic changes)
Implementing Quantitative Forecasting
Determine Method
•Time Series
•Causal Model
Collect data:
Fit an analytical model - Determine
functional form
to the data: - Estimate parameters
F(t+1) = f(X1, X2,…) - Validate
Use the model for
forecasting future
demand
Update Model
Parameters
Monitor error:
e(t+1) = D(t+1)-F(t+1)
Yes Model No
Valid?
Time Series-based Forecasting
Basic Model:
D(i), i 1,...,t
Time Series ˆ
D(t ), 1,2,...
Model
Historical Forecasts
Data
Remark: The exact model to be used depends on the expected /
observed trends in the data.
Cases typically considered:
• Constant mean series
• Series with linear trend
• Series with seasonalities (and possibly a linear trend)
A constant mean series
14.00
12.00
10.00
8.00
Series1
6.00
4.00
2.00
0.00
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20
The above data points have been sampled from a normal distribution with a
mean value equal to 10.0 and a variance equal to 4.0.
Forecasting constant mean series:
The Moving Average model
The presumed model for the observed data:
D(t ) D e(t )
where
D is the constant mean of the series and
e(t ) is normally 2
distributed with zero mean and some unknown
variance
Then, under a Moving Average of Order N model, denoted as MA(N),
the estimate of D returned at period t, is equal to:
N 1
ˆ (t ) 1
D
N
D(t i)
i 0
The forecasting error
• The forecasting error
N 1
ˆ (t ) D(t 1) 1
(t 1) D
N
D(t i) D(t 1)
i 0
• Also
1 N 1 1
E[ (t 1)] E[ D(t i )] E[ D(t 1)] ( ND ) D 0
N i 0 N
1 N 1 1 1 2
Var[ (t 1)] 2 Var[ D(t i )] Var[ D(t 1)] 2 N (1 )
2 2
N i 0 N N
Forecasting error (cont.)
• (t+1) is normally distributed with the mean and
variance computed in the previous slide.
ˆ
• D (t ) D
follows a normal distribution with zero
mean and variance 2/N.
Selecting an appropriate order N
• Smaller values of N provide more flexibility.
• Larger values of N provide more accuracy (c.f., the formula for the
variance of the forecasting error).
• Hence, the more stable (stationary) the process, the larger the N.
• In practice, N is selected through trial and error, such that it
minimizes one of the following criteria:
t
1
i) MAD (t )
t N iN
1| (t ) |
t
ii) MSD (t ) 1
t N i N 1
( (t )) 2
iii) 1 t
(t )
MAPE(t )
t N i N 1 D(t )
Demonstrating the impact of N on
the model performance
25.00
20.00
15.00
Series1
Series2
Series3
10.00
5.00
0.00
1
4
7
10
13
16
19
22
25
28
31
34
37
40
• blue series: the original data series, distributed according to N(10,4) for the first 20
points, and N(20,4) for the last 20 points.
• magenta series: the predictions of the MA(5) forecasting model.
• yellow series: the predictions of the MA(10) forecasting model.
• Remark: the MA(5) model adjusts faster to the experienced jump of the data mean
value, but the mean estimates that it provides under stationary operation are less accurate
than those provided by the MA(10) model.
Forecasting constant mean series:
The Simple Exponential Smoothing model
• The presumed demand model:
D(t ) D e(t )
where D is an unknown constant and e(t ) is normally distributed
with zero mean and an unknown variance 2 .
ˆ
• The forecast D (t ) , at the end of period t:
ˆ ˆ ˆ ˆ
D (t ) aD(t ) (1 a) D (t 1) D (t 1) a[D(t ) D (t 1)]
where (0,1) is known as the “smoothing constant”.
• Remark: The updating equation constitutes a correction of the
previous estimate in the direction suggested by the forecasting error,
ˆ
D(t ) D (t 1)
Expanding the Model Recursion
ˆ ˆ
D (t ) aD(t ) (1 a) D (t 1)
ˆ
aD(t) a(1 a)D(t 1) (1 a) 2 D (t 2)
.................................................................................................
t 1
ˆ
a (1 a) i D(t i ) (1 a ) t D (0)
i 0
Implications
1. The model considers all the past observations and the
ˆ
initializing value D (0) in the determination of the
ˆ
estimate D (t ) .
2. The weight of the various data points decreases
exponentially with their age.
3. As 1, the model places more emphasis on the most
recent observations.
4. As t,
ˆ (t )] D and ˆ (t )] 2 2
E[ D Var[ D(t 1) D
2a
ˆ (0)
The impact of and of D on
the model performance
25.00
20.00
15.00 Series1
Series2
Series3
10.00 Series4
5.00
0.00
1
4
7
10
13
16
19
22
25
28
31
34
37
40
• dark blue series: the original data series, distributed according to N(10,4) for the first 20
points, and N(20,4) for the last 20 points.
• magenta series: the predictions of the ES(0.2) model initialized at the value of 10.0.
• yellow series: the predictions of the ES(0.2) model initialized as 0.0.
• light blue series: the predictions of the ES(0.8) model initialized at 10.0.
• Remark: the ES(0.8) model adjusts faster to the jump of the series mean value, but the
estimates that it provides under stationary operation are less accurate than those provided by
the ES(0.2) model. Also, notice that the effect of the initial value is only transient.
The inadequacy of SES and MA
models for data with linear trends
12
10
8
Dt
6 SES(0.5)
SES(1.0)
4
2
0
1 2 3 4 5 6 7 8 9 10
• blue series: a deterministic data series increasing linearly with a slope of 1.0.
• magenta series: the predictions obtained from the SES(0.5) model initialized at
the exact value of 1.0.
• yellow series: the predictions obtained from the SES(1.0) model initialized at
the exact value of 1.0.
• Remark: Both models under-estimate the actual values, with the most inert
model SES(0.5) under-estimating the most. This should be expected since both of
these models (as well as any MA model) essentially average the past
observations. Therefore, neither the MA nor the SES model are appropriate for
forecasting a data series with a linear trend in it.
Forecasting series with linear trend:
The Double Exponential Smoothing Model
The presumed data model:
D(t ) I T t e(t )
where
I is the model intercept, i.e., the unknown mean value for t=0,
T is the model trend, i.e., the mean increase per unit of time, and
e(t ) is normally distributed with zero mean and some unknown
variance 2
The Double Exponential Smoothing
Model (cont.)
The model forecasts at period t for periods t+, =1,2,…, are
given by:
ˆ ˆ ˆ
D (t ) I (t ) T (t )
ˆ ˆ
with the quantities I (t ) and T (t ) obtained through the following
recursions:
ˆ ˆ ˆ
I (t ) a D(t ) (1 a )[ I (t 1) T (t 1)]
ˆ ˆ ˆ ˆ
T (t ) b [ I (t ) I (t 1)] (1 b ) T (t 1)
The parameters a and btake values in the interval (0,1) and are the
ˆ ˆ
model smoothing constants, while the values I (0) and T (0) are the
initializing values.
The Double Exponential Smoothing
Model (cont.)
• The smoothing constants are chosen by trial and error, using the
MAD, MSD and/or MAPE indices.
ˆ ˆ
• For t, I (t ) I and T (t ) T
• The variance of the forecasting error, , can be estimated as a
2
function of the noise variance 2 through techniques similar to
those used in the case of the Simple Exp. Smoothing model, but in
practice, it is frequently approximated by
2 1.25 MAD (t )
ˆ
where
MAD(t ) g (t ) (1 g )MAD(t 1)
for some appropriately selected smoothing constant g(0,1) or by
2 MSD (t )
ˆ
DES Example
12
10
8
Dt
6 DES(T0=1)
DES(T0=0)
4
2
0
1 2 3 4 5 6 7 8 9 10
• blue series: a deterministic data series increasing linearly with a slope of 1.0.
• magenta series: the predictions obtained from the DES(0.5;0.2) model initialized
at the exact value of 1.0.
• yellow series: the predictions obtained from the DES(0.5;0.2) model initialized
at the value of 0.0.
• Remark: In the absence of variability in the original data, the first model is
completely accurate (the blue and the magenta series overlap completely), while
the second model overcomes the deficiency of the wrong initial estimate and
eventually converges to the correct values.
Time Series-based Forecasting:
Accommodating seasonal behavior
The data demonstrate a periodic behavior (and maybe some
additional linear trend).
Example: Consider the following data, describing a quarterly
demand over the last 3 years, in 1000’s:
Year 1 Year 2 Year 3
Spring 90 115 120
Summer 180 230 290
Fall 70 85 105
Winter 60 70 100
Total 400 500 615
Seasonal Indices
Plotting the demand data:
350
300
250
200
Series1
150
100
50
0
0 2 4 6 8 10 12 14
Remarks:
• At each cycle, the demand of a particular season is a fairly stable percentage of
the total demand over the cycle.
• Hence, the ratio of a seasonal demand to the average seasonal demand of the
corresponding cycle will be fairly constant.
• This ratio is characterized as the corresponding seasonal index.
A forecasting methodology
Forecasts for the seasonal demand for subsequent years can be obtained by:
i. estimating the seasonal indices corresponding to the various seasons in the
cycle;
ii. estimating the average seasonal demand for the considered cycle (using, for
instance, a forecasting model for a series with constant mean or linear trend,
depending on the situation);
iii. adjusting the average seasonal demand by multiplying it with the
corresponding seasonal index.
Example (cont.):
Year 1 Year 2 Year 3 SI(1) SI(2) SI(3) SI
Spring 90 115 120 0.9 0.92 0.78 0.87
Summer 180 230 290 1.8 1.84 1.88 1.84
Fall 70 85 105 0.7 0.68 0.68 0.69
Winter 60 70 100 0.6 0.56 0.65 0.6
Total 400 500 615 4 4 4 4
Average 100 125 153.75
Winter’s Method for Seasonal
Forecasting
The presumed model for the observed data:
D(t ) ( I T t ) c( t 1) m odN 1 e(t )
where
• N denotes the number of seasons in a cycle;
• ci, i=1,2,…N, is the seasonal index for the i-th season in the cycle;
• I is the intercept for the de-seasonalized series obtained by dividing the
original demand series with the corresponding seasonal indices;
• T is the trend of the de-seasonalized series;
• e(t) is normally distributed with zero mean and some unknown variance 2
Winter’s Method for Seasonal Forecasting
(cont.)
The model forecasts at period t for periods t+, 1,2,…, are given by:
ˆ ˆ ˆ
D (t ) [ I (t ) T (t ) ] c(t 1) modN 1 (t )
ˆ
where the quantities I (t ) , T (t ) and ci (t ), i 1,..., N ,are obtained from the
ˆ ˆ ˆ
following recursions, performed in the indicated sequence:
ˆ D(t ) ˆ ˆ
I (t ) : a (1 a )[ I (t 1) T (t 1)]
c( t 1) m odN 1 (t 1)
ˆ
ˆ ˆ ˆ ˆ
T (t ) : b [ I (t ) I (t 1)] (1 b ) T (t 1)
D(t )
c(t 1) m odN 1 (t ) : g
ˆ (1 g ) c(t 1) m odN 1 (t 1)
ˆ
ˆ(t )
I
ci (t ) : ci (t 1), i (t 1) mod N 1
ˆ ˆ
The parameters ,b,gtake values in the interval (0,1) and are the model smoothing
constants, while I (0), T (0) and ci (0), i 1,..., N , are the initializing values.
ˆ ˆ ˆ
Causal Models:
Multiple Linear Regression
• The basic model:
D b0 b1 X 1 ... bk X k e
where
• Xi, i=1,…,k, are the model independent variables (otherwise known as the
explanatory variables);
• bi, i=0,…,k, are unknown model parameters;
• e is the a random variable following a normal distribution with zero mean and
some unknown variance 2.
• D follows a normal distribution N ( D , ) where
2
D b0 b1 X 1 ... bk X k
• We need to estimate and 2 from a set of n observations
{ D j ; X 1 j , X 2 j ,..., X kj , j 1,..., n}
Estimating the parameters bi
• The observed data satisfy the following equation:
D1 1 X 11 ... X k1 b0 e1
D 1 X ... X k 2 b1 e2
2 12
... ... ... ... ... ... ...
Dn 1 X 1n ... X kn bk en
or in a more concise form
d X b e
• The vector
e d X b
denotes the difference between the actual observations and the corresponding
ˆ
mean values, and therefore, b is selected such that it minimizes the Euclidean
ˆ
norm of the resulting vector e d X b .
ˆ
ˆ ˆ
• The minimizing value for b is equal to b ( X T X ) 1 X T d
• The necessary and sufficient condition for the existence of ( X T X ) 1 is that the
columns of matrix X are linearly independent.
Characterizing the model variance
• An unbiased estimate of 2is given by
SSE
MSE (Mean Squared Error)
n k 1
where
ˆ ˆ
SSE eT e (d X b)T (d X b)
ˆ ˆ (Sum of Squared Errors)
• The quantity SSE/2 follows a Chi-square distribution with n-k-1 degrees of
freedom.
• Given a point x0T=(1,x10,…,xk0), an unbiased estimator of D ( x0 )is given by
ˆ ˆ ˆ ˆ
D ( x0 ) b0 b1 x10 ... bk xk 0
• This estimator is normally distributed with mean D ( x0 ) and variance 2 x0 ( X T X )1 x0
T
ˆ
• The random variable D ( x0 ) can function also as an estimator for any single
ˆ
observation D(x0). The resulting error D ( x0 ) D( x0 ) will have zero mean and
variance 2[1 x0 ( X T X )1 x0 ]
T
Assessing the goodness of fit
• A rigorous characterization of the quality of the resulting approximation can be
obtained through Analysis of Variance, that can be traced in any introductory book
on statistics.
•A more empirical test considers the coefficient of multiple determination
SSR
R
2
SYY
where n
ˆ
SSR bT ( X T d ) nd 2 ( D j d ) 2
ˆ
n
j 1
1
d Dj
n j 1
and
SYY SSE SSR
• Remark: A natural way to interpret R2 is as the fraction of the variability in the
observed data interpreted by the model over the total variability in this data.
Multiple Linear Regression and
Time Series-based forecasting
• The model needs to be linear with respect to the parameters bi but not the
explanatory variables Xi. Hence, the factor multiplying the parameter bi can be any
function fi of the underlying explanatory variables.
• When the only explanatory variable is just the time variable t, the resulting multiple
linear regression model essentially supports time-series analysis.
• The above approach for time-series analysis enables the study of more complex
dependencies on time than those addressed by the moving average and exponential
smoothing models.
• The integration of a new observation in multiple linear regression models is much
more cumbersome than the updating performed by the moving average and
exponential smoothing models (although there is an incremental linear regression
model that alleviates this problem).
Confidence Intervals
• Given a random variable X and p(0,1), a p100% confidence
interval (CI) for it is an interval [a,b] such that
P(a X b) p
• Confidence intervals are used in:
i. monitoring the performance of the applied forecasting model;
ii. adjusting an obtained forecast in order to achieve a certain
performance level
• The necessary confidence intervals are obtained by exploiting
the statistics for the forecasting error, derived in the previous
slides.
Variance estimation and the t distribution
• The variance of the forecasting error is a function of the unknown variance,
2, of the model disturbance, e.
• E.g., in the case of multiple linear regression, the variance of the forecasting
ˆ
error D ( x0 ) D( x0 ) is equal to 2[1 x0T ( X T X )1 x0 ] .
• Hence, one cannot take advantage directly of the normality of the forecasting
error in order to build the sought confidence intervals.
• This problem can be circumvented by exploiting the fact that the quantity
SSE/2 follows a Chi-square distribution with n-k-1 degrees of freedom.
Then, the quantity
ˆ
[ D ( x0 ) D( x0 )] 1 x0 ( X T X ) 1 x0 ˆ
D ( x0 ) D( x0 )
T
T
SSE 2
MSE [1 x0 ( X T X ) 1 x0 ]
T
n k 1
follows a t distribution with n-k-1 degrees of freedom.
• For large samples, T can also be approximated by a standardized normal
distribution.
Adjusting the forecasted demand in
order to achieve a target service level p
Letting y denote the required adjustment, we essentially need to solve the following
equation:
ˆ
P( D( x0 ) D ( x0 ) y) p
ˆ
D( x0 ) D ( x0 ) y
P( ) p
1 1
MSE[1 x0 ( X X ) x0 ] MSE[1 x0 ( X X ) x0 ]
T T T T
y
t p ,nk 1
MSE[1 x0 ( X T X ) 1 x0 ]
T
y t p ,nk 1 MSE[1 x0 ( X T X ) 1 x0 ]
T
Remark: The two-sided confidence interval that is necessary for monitoring the
model performance can be obtained through a straightforward modification of the
above reasoning.