338
Views
17
CrossRef citations to date
0
Altmetric
Original Articles

SEM Modeling with Singular Moment Matrices Part I: ML-Estimation of Time Series

Pages 301-320 | Published online: 20 Sep 2010
 

Abstract

A structural equation model (SEM) with deterministic intercepts is introduced. The Gaussian likelihood function does not contain determinants of sample moment matrices and is thus well-defined for only one statistical unit. The SEM is applied to the dynamic state space model and compared with the Kalman filter (KF) approach. The likelihood of both methods are shown to be equivalent, but for long time series numerical problems occur in the SEM approach, which are traced to the inversion of the latent state covariance matrix. Both approaches are compared on several aspects. The SEM approach is now open for idiographic (N = 1) analysis and estimation of panel data with correlated units.

Notes

1Otherwise there are exact identities in the components of y and the singular normal distribution can be used (Mardia, Kent, and Bibby, Citation1979, p. 41). See also Footnote 4.

2Moreover, the system matrices may depend on lagged measurements Z i = {z i, …, z 0}, and the measurement matrices Hi,d i R i may depend on Z i−1 in order to specify ARCH effects (AutoRegressive Conditional Heteroscedasticity [Bollerslev, Engle, and Nelson, Citation1994]). ARCH models are used in financial econometrics to model the volatility fluctuations of return series of stocks, exchange rates, etc. The so called conditional Gaussian model can be treated by the Kalman filter (Liptser and Shiryayev, Citation2001, Vol. 2), but not by SEM. This is because the SEM approach assumes a joint Gaussian distribution of the manifest and latent states, whereas in the recursive approach only conditional Gaussianity is required.

3SEM with N = 1; SEM states y and x (without index i) must be distinguished from the dynamical states y i , x i in Eqs. (9–10).

4Σ y is positive definite (p.d.), if either (i) Σϵ is p.d. (cf. Magnus and Neudecker, Citation1999, Theorem 22, p. 21) or (ii) if Σζ is p.d. and ΛB 1 has rank K = k(T + 1) (cf. Golub and Van Loan, Citation1996, Theorem 4.2.1, p. 141).

5The prediction errors ν i+1 = z i+1 − E[z i+1|Z i ] are martingale differences, that is, E[νi+1|Z i ] = 0 and thus uncorrelated, since for i > j we have .

For system matrices independent of the measurements they are Gaussian and thus even independent (cf. Liptser and Shiryayev, 2001, Vol. 2, Ch. 13).

6Highest posterior density

Reprints and Corporate Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

To request a reprint or corporate permissions for this article, please click on the relevant link below:

Academic Permissions

Please note: Selecting permissions does not provide access to the full text of the article, please see our help page How do I view content?

Obtain permissions instantly via Rightslink by clicking on the button below:

If you are unable to obtain permissions via Rightslink, please complete and submit this Permissions form. For more information, please visit our Permissions help page.