In this talk we discuss an extension of standard auto-regressive models to capture heterogeneous behavior in mean level, volatility and multi-modality of the conditional or marginal distributions of time series that are observed in practice in many financial and econometric examples. We will develop new estimation and model selection techniques in these models. The new methods are assessed theoretically, via simulations and a real data analysis.