h inline">\(\epsilon_t\)</span> being the
noise at time <span class="math inline">\(t\)</span>, the model for outcome <span class="math inline">\(y_t\)</span> is defined by
<span class="math display">\[
y_t = \mu + \theta_1 \epsilon_{t-1} + \dotsb + \theta_Q \epsilon_{t-Q}
+ \epsilon_t,
\]</span>
with the noise term <span class="math inline">\(\epsilon_t\)</span> for outcome <span class="math inline">\(y_t\)</span> modeled as
normal,
<span class="math display">\[
\epsilon_t \sim \textsf{normal}(0,\sigma).
\