In this paper, we propose a novel and efficient two-stage variable selection approach for sparse GLARMA models, which are pervasive modelling discrete-valued time series. Our consists in iteratively combining the estimation of autoregressive moving average (ARMA) coefficients models with regularized methods designed performing regression Generalized Linear Models (GLM). We first establish consi...