Saved in:
Bibliographic Details
Main Authors: Binks, Rachel L., Heaps, Sarah E., Panagiotopoulou, Mariella, Wang, Yujiang, Wilkinson, Darren J.
Format: Preprint
Published: 2023
Subjects:
Online Access:https://arxiv.org/abs/2307.05708
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866910725916590080
author Binks, Rachel L.
Heaps, Sarah E.
Panagiotopoulou, Mariella
Wang, Yujiang
Wilkinson, Darren J.
author_facet Binks, Rachel L.
Heaps, Sarah E.
Panagiotopoulou, Mariella
Wang, Yujiang
Wilkinson, Darren J.
contents Vector autoregressions (VARs) are a widely used tool for modelling multivariate time-series. It is common to assume a VAR is stationary; this can be enforced by imposing the stationarity condition which restricts the parameter space of the autoregressive coefficients to the stationary region. However, implementing this constraint is difficult due to the complex geometry of the stationary region. Fortunately, recent work has provided a solution for autoregressions of fixed order $p$ based on a reparameterization in terms of a set of interpretable and unconstrained transformed partial autocorrelation matrices. In this work, focus is placed on the difficult problem of allowing $p$ to be unknown, developing a prior and computational inference that takes full account of order uncertainty. Specifically, the multiplicative gamma process is used to build a prior which encourages increasing shrinkage of the partial autocorrelations with increasing lag. Identifying the lag beyond which the partial autocorrelations become equal to zero then determines $p$. Based on classic time-series theory, a principled choice of truncation criterion identifies whether a partial autocorrelation matrix is effectively zero. Posterior inference utilizes Hamiltonian Monte Carlo via Stan. The work is illustrated in a substantive application to neural activity data to investigate ultradian brain rhythms.
format Preprint
id arxiv_https___arxiv_org_abs_2307_05708
institution arXiv
publishDate 2023
record_format arxiv
spellingShingle Bayesian inference on the order of stationary vector autoregressions
Binks, Rachel L.
Heaps, Sarah E.
Panagiotopoulou, Mariella
Wang, Yujiang
Wilkinson, Darren J.
Methodology
62F15, 62M10
Vector autoregressions (VARs) are a widely used tool for modelling multivariate time-series. It is common to assume a VAR is stationary; this can be enforced by imposing the stationarity condition which restricts the parameter space of the autoregressive coefficients to the stationary region. However, implementing this constraint is difficult due to the complex geometry of the stationary region. Fortunately, recent work has provided a solution for autoregressions of fixed order $p$ based on a reparameterization in terms of a set of interpretable and unconstrained transformed partial autocorrelation matrices. In this work, focus is placed on the difficult problem of allowing $p$ to be unknown, developing a prior and computational inference that takes full account of order uncertainty. Specifically, the multiplicative gamma process is used to build a prior which encourages increasing shrinkage of the partial autocorrelations with increasing lag. Identifying the lag beyond which the partial autocorrelations become equal to zero then determines $p$. Based on classic time-series theory, a principled choice of truncation criterion identifies whether a partial autocorrelation matrix is effectively zero. Posterior inference utilizes Hamiltonian Monte Carlo via Stan. The work is illustrated in a substantive application to neural activity data to investigate ultradian brain rhythms.
title Bayesian inference on the order of stationary vector autoregressions
topic Methodology
62F15, 62M10
url https://arxiv.org/abs/2307.05708