{
 "bundle_of": "C4",
 "schema_version": "0.42",
 "paper": {
  "id": "entrylib:2026.000006",
  "version": 1,
  "title": "Staying Positive: Challenges and Solutions in Using Pure Multiplicative ETS Models",
  "licence": "CC-BY-4.0",
  "status": "active",
  "date": "2026-11-12"
 },
 "claims": [
  {
   "id": "C4",
   "statement": "[[D4:Log-Normal error distribution]], [[D5:Gamma error distribution]], and [[D6:Inverse Gaussian error distribution]] distributions defined on $(0, \\infty)$ can be applied directly to the multiplicative error term $(1+\\epsilon_t)$ of conventional [[D1:ETS(M,*,*)]] models without transforming the model into [[D10:Akram's log-ETS model]], provided $E(1+\\epsilon_t)=1$.",
   "about": {
    "subject": "D1",
    "relation": "supports positive error distributions directly on the multiplicative error term",
    "object": null,
    "comparator": "D10",
    "conditions": [
     "unit expectation restriction $E(1 + \\epsilon_t) = 1$"
    ]
   },
   "type": "methodological",
   "role": "principal",
   "novelty": "extension",
   "planning": "exploratory",
   "scope": {
    "conditions": [
     "strictly positive time series modelled within the original $ETS(M,*,*)$ taxonomy"
    ]
   },
   "relations": [
    {
     "type": "extends",
     "target": "10.1111/j.1467-842X.2009.00555.x",
     "id_type": "doi",
     "locator": "Section 2",
     "note": "Implements Log-Normal, Gamma, and Inverse Gaussian distributions directly within original $ETS(M,*,*)$ rather than sidestepping the taxonomy via logarithmic ETS."
    }
   ],
   "support": [
    {
     "kind": "formal_proof",
     "evidence": "E3"
    },
    {
     "kind": "cited",
     "reference": "B4",
     "locator": "Section 2",
     "self": false
    }
   ],
   "assumptions": [
    "C12",
    "C13"
   ],
   "verification": "unchecked",
   "assessments": [
    {
     "by": "author",
     "strength": "high",
     "rationale": "Directly established by deriving valid parameter restrictions and scalability properties on $(0, \\infty)$."
    }
   ],
   "provenance": "author_stated",
   "verifiability": "open",
   "verifiability_flags": []
  },
  {
   "id": "C12",
   "statement": "The multiplicative error term $(1+\\epsilon_t)$ in pure multiplicative ETS models is independent and identically distributed with unit expectation $E(1+\\epsilon_t) = 1$.",
   "type": "formal",
   "role": "supporting",
   "support": [
    {
     "kind": "assumption"
    }
   ],
   "assumption_category": "statistical",
   "check": "unchecked",
   "provenance": "author_stated",
   "verifiability": "not_applicable",
   "verifiability_flags": []
  },
  {
   "id": "C13",
   "statement": "Smoothing parameters $\\alpha, \\beta, \\gamma$ in pure multiplicative ETS models lie within $[0, 1]$ so that state update terms $(1 + g\\epsilon_t)$ remain strictly positive whenever $1 + \\epsilon_t > 0$.",
   "type": "formal",
   "role": "supporting",
   "support": [
    {
     "kind": "assumption"
    }
   ],
   "assumption_category": "methodological",
   "check": "tested",
   "check_evidence": "E1",
   "provenance": "author_stated",
   "verifiability": "not_applicable",
   "verifiability_flags": []
  }
 ],
 "definitions": [
  {
   "id": "D1",
   "term": "ETS(M,*,*)",
   "kind": "method",
   "definition": "Pure multiplicative exponential smoothing state space model family written in logarithmic state space form as $$\\begin{aligned} \\log y_t &= w' \\log(v_{t-l}) + \\log(1 + \\epsilon_t) \\\\ \\log v_t &= F \\log v_{t-l} + \\log(1_k + g \\epsilon_t) \\end{aligned}$$ with conventional multi-step point forecast $\\hat{y}_{t+h} = \\exp\\left((w_1' F_1^{h-1} + w_m' F_m^{\\lceil h/m \\rceil - 1}) \\log v_t\\right)$.",
   "source": "B7"
  },
  {
   "id": "D4",
   "term": "Log-Normal error distribution",
   "kind": "concept",
   "definition": "Probability distribution on $(0, \\infty)$ applied to $1+\\epsilon_t$ with location $\\mu = -\\sigma^2/2$ and scale $\\sigma^2$ so that $E(1+\\epsilon_t) = 1$ and $V(1+\\epsilon_t) = e^{\\sigma^2} - 1$."
  },
  {
   "id": "D5",
   "term": "Gamma error distribution",
   "kind": "concept",
   "definition": "Probability distribution on $(0, \\infty)$ applied to $1+\\epsilon_t$ with shape $\\xi = \\sigma^{-2}$ and scale $\\sigma^2$ so that $E(1+\\epsilon_t) = 1$ and $V(1+\\epsilon_t) = \\sigma^2$."
  },
  {
   "id": "D6",
   "term": "Inverse Gaussian error distribution",
   "kind": "concept",
   "definition": "Probability distribution on $(0, \\infty)$ applied to $1+\\epsilon_t$ with location $\\mu = 1$ and scale $\\sigma^2$ so that $E(1+\\epsilon_t) = 1$ and $V(1+\\epsilon_t) = \\sigma^2$."
  },
  {
   "id": "D10",
   "term": "Akram's log-ETS model",
   "kind": "method",
   "definition": "Additive error ETS model $ETS(A,Ad,A)$ with Normal errors applied to the natural logarithms of strictly positive time series and exponentiated back to the original scale.",
   "source": "B4"
  }
 ],
 "evidence": [
  {
   "id": "E1",
   "profile": "formal",
   "description": "Analytical derivation of recursive conditional expectations and point forecasts for pure multiplicative ETS models in logarithmic state space form (Section 2, Eqs. 2.1-2.11).",
   "foundation": "classical_probability",
   "assumptions": [
    "C12",
    "C13"
   ],
   "steps": [
    {
     "id": "E1.1",
     "statement": "Reformulate pure multiplicative ETS in logarithmic state space form: $\\log y_t = w' \\log(v_{t-l}) + \\log(1 + \\epsilon_t)$ and $\\log v_t = F \\log v_{t-l} + \\log(1_k + g \\epsilon_t)$.",
     "by": "Section 2, Eq. (2.1)",
     "rule": "logarithmic_transformation"
    },
    {
     "id": "E1.2",
     "statement": "Expand $h$-steps-ahead recursion into non-seasonal and seasonal components and take conditional expectation $E(\\log y_{t+h}|t)$.",
     "by": "Section 2, Eqs. (2.5)-(2.6)",
     "rule": "recursive_substitution"
    },
    {
     "id": "E1.3",
     "statement": "Apply Jensen's inequality and non-negativity of $g$ to show point forecast $\\hat{y}_{t+h} \\le \\tilde{y}_{t+h} \\le \\mu_{y,t+h}$, with equality only at $h=1$, $h \\le m$ without trend, or for $ETS(M,N,N)$.",
     "by": "Section 2, Eqs. (2.8)-(2.11)",
     "rule": "jensens_inequality"
    }
   ],
   "availability": {
    "level": "public"
   }
  },
  {
   "id": "E3",
   "profile": "formal",
   "description": "Analytical derivation of unit-mean parameter restrictions, concentrated log-likelihoods, and scale estimators for Normal, Inverse Gaussian, Log-Normal, and Gamma distributions in $ETS(M,*,*)$ models (Section 4.1, Tables 2-3), with author erratum referencing `entrylib:2026.000003@v4#C7` for the Log-Normal scale MLE.",
   "foundation": "mathematical_statistics",
   "assumptions": [
    "C12"
   ],
   "steps": [
    {
     "id": "E3.1",
     "statement": "Impose $E(1 + \\epsilon_t) = 1$ on each distribution: $\\mu = 1$ for Normal and Inverse Gaussian, $\\mu = -\\sigma^2 / 2$ for Log-Normal, and $\\xi = \\sigma^{-2}$ for Gamma.",
     "by": "Section 4.1, Table 2",
     "rule": "moment_constraint"
    },
    {
     "id": "E3.2",
     "statement": "Substitute parameter restrictions into probability density functions and Jacobian terms to obtain concentrated log-likelihood functions.",
     "by": "Section 4.1, Table 3",
     "rule": "likelihood_derivation"
    },
    {
     "id": "E3.3",
     "statement": "Differentiate log-likelihoods with respect to $\\sigma^2$ to obtain closed-form scale MLEs for Normal and Inverse Gaussian, method-of-moments estimator for Gamma, and the corrected Log-Normal scale MLE referenced in `entrylib:2026.000003@v4#C7`.",
     "by": "Section 4.1, pp. 9-10 and entrylib:2026.000003@v4#C7",
     "rule": "maximum_likelihood_estimation"
    }
   ],
   "availability": {
    "level": "public"
   }
  }
 ],
 "references": [
  {
   "id": "B4",
   "target": "10.1111/j.1467-842X.2009.00555.x",
   "id_type": "doi"
  },
  {
   "id": "B7",
   "target": "10.1201/9781003452652",
   "id_type": "doi"
  }
 ]
}
