cvts function

Cross validation for time series