An autoregressive model for irregular time series of variable stars