We develop threshold models that allow volatilities and copula functions or their association parameters to change across time. The number and location of the thresholds is assumed unknown. We use a Markov chain Monte Carlo strategy combined with Laplace estimates that evaluate the required...