Modeling Regime Shifts in Multiple Time Series