from sklearn.model_selection import cross_val_score

# Test model on multiple train/test splits
scores = cross_val_score(model, X, y, cv=5)  # 5-fold cross-validation
print(f"Average score: {scores.mean():.2f} (+/- {scores.std()*2:.2f})")