holdout function

Computes indexes for holdout data split into training and test sets.