stable_bvar()
detects and discards all posterior draws of an
bayesianVARs_bvar
object that do not fulfill the stability condition:
A VAR(p) model is considered as stable only if the eigenvalues of the
companion form matrix lie inside the unit circle.
Arguments
- object
A
bayesianVARs_bvar
object obtained viabvar()
.- quiet
logical indicating whether informative output should be omitted.
Examples
# Access a subset of the usmacro_growth dataset
data <- usmacro_growth[,c("GDPC1", "CPIAUCSL", "FEDFUNDS")]
# Estimate a model
mod <- bvar(data, sv_keep = "all", quiet = TRUE)
# Discard "unstable" draws
stable_mod <- stable_bvar(mod)
#>
#> Original 'bayesianVARs_bvar' object consists of 1000 posterior draws.
#>
#> Detected 475 unstable draws.
#>
#> Remaining draws: 525 !