transformSSM(object, type = c("ldl", "augment"))
'ldl'
performs LDL
decomposition for covariance matrix $H_t$, and
multiplies the observation equation with the
$L_t^{-1}$, so $\epsilon_t^* \sim
N(0,D_t)$. Option
'augment'
adds $\epsilon_t$ to
the state