# Facebook Prophet: Multiplicative Seasonality

from fbprophet import Prophet import pandas as pd import numpy as np from matplotlib import pyplot as plt import logging logging.getLogger('fbprophet').setLevel(logging.ERROR) import warnings warnings.filterwarnings("ignore")

library(prophet) library(ggplot2) # Required for Nextjournal plotting

By default Prophet fits additive seasonalities, meaning the effect of the seasonality is added to the trend to get the forecast. This time series of the number of air passengers is an example of when additive seasonality does not work:

df <- read.csv(example_wp_log_peyton_manning.csv) m <- prophet(df) future <- make_future_dataframe(m, 50, freq = 'm') forecast <- predict(m, future) plot(m, forecast)

df = pd.read_csv(example_wp_log_peyton_manning.csv) m = Prophet() m.fit(df) future = m.make_future_dataframe(50, freq='MS') forecast = m.predict(future) fig = m.plot(forecast) fig

This time series has a clear yearly cycle, but the seasonality in the forecast is too large at the start of the time series and too small at the end. In this time series, the seasonality is not a constant additive factor as assumed by Prophet, rather it grows with the trend. This is multiplicative seasonality.

Prophet can model multiplicative seasonality by setting `seasonality_mode='multiplicative'`

in the input arguments:

m <- prophet(df, seasonality.mode = 'multiplicative') forecast <- predict(m, future) plot(m, forecast)

m = Prophet(seasonality_mode='multiplicative') m.fit(df) forecast = m.predict(future) fig = m.plot(forecast) fig

The components figure will now show the seasonality as a percent of the trend:

prophet_plot_components(m, forecast)

fig = m.plot_components(forecast) fig

With `seasonality_mode='multiplicative'`

, holiday effects will also be modeled as multiplicative. Any added seasonalities or extra regressors will by default use whatever `seasonality_mode`

is set to, but can be overriden by specifying `mode='additive'`

or `mode='multiplicative'`

as an argument when adding the seasonality or regressor.

For example, this block sets the built-in seasonalities to multiplicative, but includes an additive quarterly seasonality and an additive regressor:

m <- prophet(seasonality.mode = 'multiplicative') m <- add_seasonality(m, 'quarterly', period = 91.25, fourier.order = 8, mode = 'additive') m <- add_regressor(m, 'regressor', mode = 'additive')

m = Prophet(seasonality_mode='multiplicative') m.add_seasonality('quarterly', period=91.25, fourier_order=8, mode='additive') m.add_regressor('regressor', mode='additive')

Additive and multiplicative extra regressors will show up in separate panels on the components plot.