Bayesian analysis of two-regime threshold autoregressive-moving average models with exogenous inputs