You signed in with another tab or window. Reload to refresh your session.You signed out in another tab or window. Reload to refresh your session.You switched accounts on another tab or window. Reload to refresh your session.Dismiss alert
Copy file name to clipboardExpand all lines: vignettes/issm_introduction.Rmd
+12-1Lines changed: 12 additions & 1 deletion
Original file line number
Diff line number
Diff line change
@@ -35,7 +35,7 @@ The `tsissm` package implements the linear innovations state space model describ
35
35
36
36
1.**Automatic differentiation and robust inference**:
37
37
Estimation leverages automatic differentiation (autodiff), with multiple sandwich estimators available for standard error calculation. System forecastability and ARMA constraints
38
-
are handled correctly via the [nloptr](https://CRAN.R-project.org/package=nloptr) solver using autodiff-based Jacobians.
38
+
are handled exactly via the [nloptr](https://CRAN.R-project.org/package=nloptr) solver using autodiff-based Jacobians.
39
39
40
40
2.**Flexible error distributions**:
41
41
In addition to Gaussian errors, the model supports heavy-tailed and skewed alternatives, including the Student's $t$ distribution and the Johnson's SU distribution.
@@ -58,6 +58,15 @@ The `tsissm` package implements the linear innovations state space model describ
58
58
59
59
\section{Model Formulation}\label{sec:model}
60
60
61
+
Given an initial state vector of unobserved components (such as level, slope, and seasonality),
62
+
the proposed model evolves the states over time using a linear transition equation,
63
+
incorporating the effect of the most recent observation error. At each time step,
64
+
the observed (Box-Cox transformed) value is modeled as a linear combination of the
65
+
previous state, lagged external regressors, and a normally^[In our formulation we relax this to allow for other choices as well.] distributed
66
+
random error. This structure allows the model to capture complex patterns in the
67
+
data—such as trends, seasonal cycles, and the influence of exogenous variables—while
68
+
dynamically updating its internal state based on new information.
69
+
61
70
Consider the following Single Error Model (SEM) with trigonometric seasonality:
62
71
63
72
\begin{equation}
@@ -72,6 +81,8 @@ where $\lambda$ represents the Box Cox parameter, $\bf{w}$ the observation coeff
72
81
vector, $\bf{x}_t$ the unobserved state vector, and $\bf{c}$ a vector of coefficients
73
82
on the external regressor set $\bf{u}$.^[In the package, it is expected that the regression matrix is already pre-lagged.]
74
83
84
+
85
+
75
86
Define the state vector^[The following equations apply for the case when all components are present (level, slope, seasonal and ARMA).] as:
0 commit comments