bivariate time series of intensities and angles. It allows us to segment cylindrical time series
according to a finite number of regimes that represent the conditional distributions of the
data under specific environmental conditions. The model parsimoniously accommodates for
circular–linear correlation, multimodality, skewness, and temporal autocorrelation. A
computationally efficient expectation–maximization algorithm is described to estimate the …