engine.py 35.7 KB
Newer Older
wxchan's avatar
wxchan committed
1
# coding: utf-8
2
"""Library with training routines of LightGBM."""
3

4
import copy
5
import json
6
import warnings
7
from collections import OrderedDict, defaultdict
wxchan's avatar
wxchan committed
8
from operator import attrgetter
9
from pathlib import Path
10
from typing import Any, Callable, Dict, Iterable, List, Optional, Tuple, Union
11

wxchan's avatar
wxchan committed
12
import numpy as np
13

wxchan's avatar
wxchan committed
14
from . import callback
15
16
17
from .basic import (
    Booster,
    Dataset,
18
    LGBMDeprecationWarning,
19
20
21
22
23
24
25
26
27
28
29
30
    LightGBMError,
    _choose_param_value,
    _ConfigAliases,
    _InnerPredictor,
    _LGBM_BoosterEvalMethodResultType,
    _LGBM_BoosterEvalMethodResultWithStandardDeviationType,
    _LGBM_CategoricalFeatureConfiguration,
    _LGBM_CustomObjectiveFunction,
    _LGBM_EvalFunctionResultType,
    _LGBM_FeatureNameConfiguration,
    _log_warning,
)
31
from .compat import SKLEARN_INSTALLED, _LGBMBaseCrossValidator, _LGBMGroupKFold, _LGBMStratifiedKFold
wxchan's avatar
wxchan committed
32

33
__all__ = [
34
35
36
    "cv",
    "CVBooster",
    "train",
37
38
39
]


40
41
42
43
44
45
46
_LGBM_CustomMetricFunction = Union[
    Callable[
        [np.ndarray, Dataset],
        _LGBM_EvalFunctionResultType,
    ],
    Callable[
        [np.ndarray, Dataset],
47
        List[_LGBM_EvalFunctionResultType],
48
    ],
49
]
wxchan's avatar
wxchan committed
50

51
52
_LGBM_PreprocFunction = Callable[
    [Dataset, Dataset, Dict[str, Any]],
53
    Tuple[Dataset, Dataset, Dict[str, Any]],
54
55
]

56

57
58
59
60
61
62
63
64
65
def _emit_dataset_kwarg_warning(calling_function: str, argname: str) -> None:
    msg = (
        f"Argument '{argname}' to {calling_function}() is deprecated and will be removed in "
        f"a future release. Set '{argname}' when calling lightgbm.Dataset() instead. "
        "See https://github.com/microsoft/LightGBM/issues/6435."
    )
    warnings.warn(msg, category=LGBMDeprecationWarning, stacklevel=2)


66
67
68
69
70
71
72
73
def train(
    params: Dict[str, Any],
    train_set: Dataset,
    num_boost_round: int = 100,
    valid_sets: Optional[List[Dataset]] = None,
    valid_names: Optional[List[str]] = None,
    feval: Optional[Union[_LGBM_CustomMetricFunction, List[_LGBM_CustomMetricFunction]]] = None,
    init_model: Optional[Union[str, Path, Booster]] = None,
74
75
    feature_name: _LGBM_FeatureNameConfiguration = "auto",
    categorical_feature: _LGBM_CategoricalFeatureConfiguration = "auto",
76
    keep_training_booster: bool = False,
77
    callbacks: Optional[List[Callable]] = None,
78
) -> Booster:
79
    """Perform the training with given parameters.
wxchan's avatar
wxchan committed
80
81
82
83

    Parameters
    ----------
    params : dict
84
85
        Parameters for training. Values passed through ``params`` take precedence over those
        supplied via arguments.
Guolin Ke's avatar
Guolin Ke committed
86
    train_set : Dataset
87
88
        Data to be trained on.
    num_boost_round : int, optional (default=100)
wxchan's avatar
wxchan committed
89
        Number of boosting iterations.
90
    valid_sets : list of Dataset, or None, optional (default=None)
91
        List of data to be evaluated on during training.
92
    valid_names : list of str, or None, optional (default=None)
93
        Names of ``valid_sets``.
94
    feval : callable, list of callable, or None, optional (default=None)
wxchan's avatar
wxchan committed
95
        Customized evaluation function.
Akshita Dixit's avatar
Akshita Dixit committed
96
        Each evaluation function should accept two parameters: preds, eval_data,
97
        and return (eval_name, eval_result, is_higher_better) or list of such tuples.
98

99
            preds : numpy 1-D array or numpy 2-D array (for multi-class task)
100
                The predicted values.
101
                For multi-class task, preds are numpy 2-D array of shape = [n_samples, n_classes].
102
                If custom objective function is used, predicted values are returned before any transformation,
103
                e.g. they are raw margin instead of probability of positive class for binary task in this case.
Akshita Dixit's avatar
Akshita Dixit committed
104
            eval_data : Dataset
105
                A ``Dataset`` to evaluate.
106
            eval_name : str
107
                The name of evaluation function (without whitespaces).
108
109
110
111
112
            eval_result : float
                The eval result.
            is_higher_better : bool
                Is eval result higher better, e.g. AUC is ``is_higher_better``.

113
114
        To ignore the default metric corresponding to the used objective,
        set the ``metric`` parameter to the string ``"None"`` in ``params``.
115
    init_model : str, pathlib.Path, Booster or None, optional (default=None)
116
        Filename of LightGBM model or Booster instance used for continue training.
117
    feature_name : list of str, or 'auto', optional (default="auto")
118
        **Deprecated.** Set ``feature_name`` on ``train_set`` instead.
119
120
        Feature names.
        If 'auto' and data is pandas DataFrame, data columns names are used.
121
    categorical_feature : list of str or int, or 'auto', optional (default="auto")
122
        **Deprecated.** Set ``categorical_feature`` on ``train_set`` instead.
123
124
        Categorical features.
        If list of int, interpreted as indices.
125
        If list of str, interpreted as feature names (need to specify ``feature_name`` as well).
126
        If 'auto' and data is pandas DataFrame, pandas unordered categorical columns are used.
127
        All values in categorical features will be cast to int32 and thus should be less than int32 max value (2147483647).
128
        Large values could be memory consuming. Consider using consecutive integers starting from zero.
129
        All negative values in categorical features will be treated as missing values.
130
        The output cannot be monotonically constrained with respect to a categorical feature.
131
        Floating point numbers in categorical features will be rounded towards 0.
132
133
134
    keep_training_booster : bool, optional (default=False)
        Whether the returned Booster will be used to keep training.
        If False, the returned value will be converted into _InnerPredictor before returning.
135
        This means you won't be able to use ``eval``, ``eval_train`` or ``eval_valid`` methods of the returned Booster.
136
137
        When your model is very large and cause the memory error,
        you can try to set this param to ``True`` to avoid the model conversion performed during the internal call of ``model_to_string``.
138
        You can still use _InnerPredictor as ``init_model`` for future continue training.
139
    callbacks : list of callable, or None, optional (default=None)
140
        List of callback functions that are applied at each iteration.
141
        See Callbacks in Python API for more information.
wxchan's avatar
wxchan committed
142

143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
    Note
    ----
    A custom objective function can be provided for the ``objective`` parameter.
    It should accept two parameters: preds, train_data and return (grad, hess).

        preds : numpy 1-D array or numpy 2-D array (for multi-class task)
            The predicted values.
            Predicted values are returned before any transformation,
            e.g. they are raw margin instead of probability of positive class for binary task.
        train_data : Dataset
            The training dataset.
        grad : numpy 1-D array or numpy 2-D array (for multi-class task)
            The value of the first order derivative (gradient) of the loss
            with respect to the elements of preds for each sample point.
        hess : numpy 1-D array or numpy 2-D array (for multi-class task)
            The value of the second order derivative (Hessian) of the loss
            with respect to the elements of preds for each sample point.

    For multi-class task, preds are numpy 2-D array of shape = [n_samples, n_classes],
    and grad and hess should be returned in the same format.

wxchan's avatar
wxchan committed
164
165
    Returns
    -------
166
167
    booster : Booster
        The trained Booster model.
wxchan's avatar
wxchan committed
168
    """
169
170
171
172
173
174
175
176
177
178
179
180
181
182
    if not isinstance(train_set, Dataset):
        raise TypeError(f"train() only accepts Dataset object, train_set has type '{type(train_set).__name__}'.")

    if num_boost_round <= 0:
        raise ValueError(f"num_boost_round must be greater than 0. Got {num_boost_round}.")

    if isinstance(valid_sets, list):
        for i, valid_item in enumerate(valid_sets):
            if not isinstance(valid_item, Dataset):
                raise TypeError(
                    "Every item in valid_sets must be a Dataset object. "
                    f"Item {i} has type '{type(valid_item).__name__}'."
                )

183
184
185
186
187
188
189
    # raise deprecation warnings if necessary
    # ref: https://github.com/microsoft/LightGBM/issues/6435
    if categorical_feature != "auto":
        _emit_dataset_kwarg_warning("train", "categorical_feature")
    if feature_name != "auto":
        _emit_dataset_kwarg_warning("train", "feature_name")

190
    # create predictor first
191
    params = copy.deepcopy(params)
192
    params = _choose_param_value(
193
        main_param_name="objective",
194
        params=params,
195
        default_value=None,
196
    )
197
    fobj: Optional[_LGBM_CustomObjectiveFunction] = None
198
199
    if callable(params["objective"]):
        fobj = params["objective"]
200
        params["objective"] = "none"
201
    for alias in _ConfigAliases.get("num_iterations"):
202
        if alias in params:
203
            num_boost_round = params.pop(alias)
204
            _log_warning(f"Found `{alias}` in params. Will use it instead of argument")
205
    params["num_iterations"] = num_boost_round
206
207
208
209
    # setting early stopping via global params should be possible
    params = _choose_param_value(
        main_param_name="early_stopping_round",
        params=params,
210
        default_value=None,
211
212
213
    )
    if params["early_stopping_round"] is None:
        params.pop("early_stopping_round")
214
    first_metric_only = params.get("first_metric_only", False)
215

216
    predictor: Optional[_InnerPredictor] = None
217
    if isinstance(init_model, (str, Path)):
218
        predictor = _InnerPredictor.from_model_file(model_file=init_model, pred_parameter=params)
wxchan's avatar
wxchan committed
219
    elif isinstance(init_model, Booster):
220
        predictor = _InnerPredictor.from_booster(booster=init_model, pred_parameter=dict(init_model.params, **params))
221
222
223
224
225

    if predictor is not None:
        init_iteration = predictor.current_iteration()
    else:
        init_iteration = 0
Guolin Ke's avatar
Guolin Ke committed
226

227
228
229
    train_set._update_params(params)._set_predictor(predictor).set_feature_name(feature_name).set_categorical_feature(
        categorical_feature
    )
Guolin Ke's avatar
Guolin Ke committed
230

wxchan's avatar
wxchan committed
231
232
    is_valid_contain_train = False
    train_data_name = "training"
Guolin Ke's avatar
Guolin Ke committed
233
    reduced_valid_sets = []
wxchan's avatar
wxchan committed
234
    name_valid_sets = []
235
    if valid_sets is not None:
Guolin Ke's avatar
Guolin Ke committed
236
237
        if isinstance(valid_sets, Dataset):
            valid_sets = [valid_sets]
238
        if isinstance(valid_names, str):
wxchan's avatar
wxchan committed
239
            valid_names = [valid_names]
Guolin Ke's avatar
Guolin Ke committed
240
        for i, valid_data in enumerate(valid_sets):
241
            # reduce cost for prediction training data
Guolin Ke's avatar
Guolin Ke committed
242
            if valid_data is train_set:
wxchan's avatar
wxchan committed
243
244
245
246
                is_valid_contain_train = True
                if valid_names is not None:
                    train_data_name = valid_names[i]
                continue
Nikita Titov's avatar
Nikita Titov committed
247
            reduced_valid_sets.append(valid_data._update_params(params).set_reference(train_set))
248
            if valid_names is not None and len(valid_names) > i:
wxchan's avatar
wxchan committed
249
250
                name_valid_sets.append(valid_names[i])
            else:
251
                name_valid_sets.append(f"valid_{i}")
252
    # process callbacks
253
    if callbacks is None:
254
        callbacks_set = set()
wxchan's avatar
wxchan committed
255
256
    else:
        for i, cb in enumerate(callbacks):
257
            cb.__dict__.setdefault("order", i - len(callbacks))
258
        callbacks_set = set(callbacks)
wxchan's avatar
wxchan committed
259

260
    if callback._should_enable_early_stopping(params.get("early_stopping_round", 0)):
261
262
        callbacks_set.add(
            callback.early_stopping(
263
                stopping_rounds=params["early_stopping_round"],  # type: ignore[arg-type]
264
                first_metric_only=first_metric_only,
265
                min_delta=params.get("early_stopping_min_delta", 0.0),
266
267
268
                verbose=_choose_param_value(
                    main_param_name="verbosity",
                    params=params,
269
270
271
                    default_value=1,
                ).pop("verbosity")
                > 0,
272
273
            )
        )
274

275
    callbacks_before_iter_set = {cb for cb in callbacks_set if getattr(cb, "before_iteration", False)}
276
    callbacks_after_iter_set = callbacks_set - callbacks_before_iter_set
277
278
    callbacks_before_iter = sorted(callbacks_before_iter_set, key=attrgetter("order"))
    callbacks_after_iter = sorted(callbacks_after_iter_set, key=attrgetter("order"))
wxchan's avatar
wxchan committed
279

280
    # construct booster
281
282
283
284
    try:
        booster = Booster(params=params, train_set=train_set)
        if is_valid_contain_train:
            booster.set_train_data_name(train_data_name)
285
        for valid_set, name_valid_set in zip(reduced_valid_sets, name_valid_sets):
286
287
288
289
290
            booster.add_valid(valid_set, name_valid_set)
    finally:
        train_set._reverse_update_params()
        for valid_set in reduced_valid_sets:
            valid_set._reverse_update_params()
291
    booster.best_iteration = 0
wxchan's avatar
wxchan committed
292

293
    # start training
294
    for i in range(init_iteration, init_iteration + num_boost_round):
wxchan's avatar
wxchan committed
295
        for cb in callbacks_before_iter:
296
297
298
299
300
301
302
303
304
305
            cb(
                callback.CallbackEnv(
                    model=booster,
                    params=params,
                    iteration=i,
                    begin_iteration=init_iteration,
                    end_iteration=init_iteration + num_boost_round,
                    evaluation_result_list=None,
                )
            )
wxchan's avatar
wxchan committed
306
307
308

        booster.update(fobj=fobj)

309
        evaluation_result_list: List[_LGBM_BoosterEvalMethodResultType] = []
wxchan's avatar
wxchan committed
310
        # check evaluation result.
311
        if valid_sets is not None:
wxchan's avatar
wxchan committed
312
313
314
315
316
            if is_valid_contain_train:
                evaluation_result_list.extend(booster.eval_train(feval))
            evaluation_result_list.extend(booster.eval_valid(feval))
        try:
            for cb in callbacks_after_iter:
317
318
319
320
321
322
323
324
325
326
                cb(
                    callback.CallbackEnv(
                        model=booster,
                        params=params,
                        iteration=i,
                        begin_iteration=init_iteration,
                        end_iteration=init_iteration + num_boost_round,
                        evaluation_result_list=evaluation_result_list,
                    )
                )
327
328
        except callback.EarlyStopException as earlyStopException:
            booster.best_iteration = earlyStopException.best_iteration + 1
wxchan's avatar
wxchan committed
329
            evaluation_result_list = earlyStopException.best_score
wxchan's avatar
wxchan committed
330
            break
331
    booster.best_score = defaultdict(OrderedDict)
wxchan's avatar
wxchan committed
332
333
    for dataset_name, eval_name, score, _ in evaluation_result_list:
        booster.best_score[dataset_name][eval_name] = score
334
    if not keep_training_booster:
335
        booster.model_from_string(booster.model_to_string()).free_dataset()
wxchan's avatar
wxchan committed
336
337
338
    return booster


339
class CVBooster:
340
341
    """CVBooster in LightGBM.

342
    Auxiliary data structure to hold and redirect all boosters of ``cv()`` function.
343
    This class has the same methods as Booster class.
344
345
346
347
348
349
    All method calls, except for the following methods, are actually performed for underlying Boosters and
    then all returned results are returned in a list.

    - ``model_from_string()``
    - ``model_to_string()``
    - ``save_model()``
350
351
352
353
354
355
356
357

    Attributes
    ----------
    boosters : list of Booster
        The list of underlying fitted models.
    best_iteration : int
        The best iteration of fitted model.
    """
358

359
360
    def __init__(
        self,
361
        model_file: Optional[Union[str, Path]] = None,
362
    ):
363
364
        """Initialize the CVBooster.

365
366
367
368
        Parameters
        ----------
        model_file : str, pathlib.Path or None, optional (default=None)
            Path to the CVBooster model file.
369
        """
370
        self.boosters: List[Booster] = []
371
        self.best_iteration = -1
372

373
374
375
376
377
378
379
380
381
        if model_file is not None:
            with open(model_file, "r") as file:
                self._from_dict(json.load(file))

    def _from_dict(self, models: Dict[str, Any]) -> None:
        """Load CVBooster from dict."""
        self.best_iteration = models["best_iteration"]
        self.boosters = []
        for model_str in models["boosters"]:
382
            self.boosters.append(Booster(model_str=model_str))
383
384
385
386
387

    def _to_dict(self, num_iteration: Optional[int], start_iteration: int, importance_type: str) -> Dict[str, Any]:
        """Serialize CVBooster to dict."""
        models_str = []
        for booster in self.boosters:
388
389
390
391
392
            models_str.append(
                booster.model_to_string(
                    num_iteration=num_iteration, start_iteration=start_iteration, importance_type=importance_type
                )
            )
393
394
        return {"boosters": models_str, "best_iteration": self.best_iteration}

395
    def __getattr__(self, name: str) -> Callable[[Any, Any], List[Any]]:
396
        """Redirect methods call of CVBooster."""
397

398
        def handler_function(*args: Any, **kwargs: Any) -> List[Any]:
399
            """Call methods with each booster, and concatenate their results."""
400
401
402
403
            ret = []
            for booster in self.boosters:
                ret.append(getattr(booster, name)(*args, **kwargs))
            return ret
404

405
        return handler_function
wxchan's avatar
wxchan committed
406

407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
    def __getstate__(self) -> Dict[str, Any]:
        return vars(self)

    def __setstate__(self, state: Dict[str, Any]) -> None:
        vars(self).update(state)

    def model_from_string(self, model_str: str) -> "CVBooster":
        """Load CVBooster from a string.

        Parameters
        ----------
        model_str : str
            Model will be loaded from this string.

        Returns
        -------
        self : CVBooster
            Loaded CVBooster object.
        """
        self._from_dict(json.loads(model_str))
        return self

    def model_to_string(
        self,
        num_iteration: Optional[int] = None,
        start_iteration: int = 0,
433
        importance_type: str = "split",
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
    ) -> str:
        """Save CVBooster to JSON string.

        Parameters
        ----------
        num_iteration : int or None, optional (default=None)
            Index of the iteration that should be saved.
            If None, if the best iteration exists, it is saved; otherwise, all iterations are saved.
            If <= 0, all iterations are saved.
        start_iteration : int, optional (default=0)
            Start index of the iteration that should be saved.
        importance_type : str, optional (default="split")
            What type of feature importance should be saved.
            If "split", result contains numbers of times the feature is used in a model.
            If "gain", result contains total gains of splits which use the feature.

        Returns
        -------
        str_repr : str
            JSON string representation of CVBooster.
        """
        return json.dumps(self._to_dict(num_iteration, start_iteration, importance_type))

    def save_model(
        self,
        filename: Union[str, Path],
        num_iteration: Optional[int] = None,
        start_iteration: int = 0,
462
        importance_type: str = "split",
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
    ) -> "CVBooster":
        """Save CVBooster to a file as JSON text.

        Parameters
        ----------
        filename : str or pathlib.Path
            Filename to save CVBooster.
        num_iteration : int or None, optional (default=None)
            Index of the iteration that should be saved.
            If None, if the best iteration exists, it is saved; otherwise, all iterations are saved.
            If <= 0, all iterations are saved.
        start_iteration : int, optional (default=0)
            Start index of the iteration that should be saved.
        importance_type : str, optional (default="split")
            What type of feature importance should be saved.
            If "split", result contains numbers of times the feature is used in a model.
            If "gain", result contains total gains of splits which use the feature.

        Returns
        -------
        self : CVBooster
            Returns self.
        """
        with open(filename, "w") as file:
            json.dump(self._to_dict(num_iteration, start_iteration, importance_type), file)

        return self

491

492
493
494
495
496
497
def _make_n_folds(
    full_data: Dataset,
    folds: Optional[Union[Iterable[Tuple[np.ndarray, np.ndarray]], _LGBMBaseCrossValidator]],
    nfold: int,
    params: Dict[str, Any],
    seed: int,
498
499
500
    fpreproc: Optional[_LGBM_PreprocFunction],
    stratified: bool,
    shuffle: bool,
501
    eval_train_metric: bool,
502
) -> CVBooster:
503
    """Make a n-fold list of Booster from random indices."""
wxchan's avatar
wxchan committed
504
505
    full_data = full_data.construct()
    num_data = full_data.num_data()
506
    if folds is not None:
507
508
509
510
511
512
        if not hasattr(folds, "__iter__") and not hasattr(folds, "split"):
            raise AttributeError(
                "folds should be a generator or iterator of (train_idx, test_idx) tuples "
                "or scikit-learn splitter object with split method"
            )
        if hasattr(folds, "split"):
513
514
            group_info = full_data.get_group()
            if group_info is not None:
515
                group_info = np.asarray(group_info, dtype=np.int32)
516
                flatted_group = np.repeat(range(len(group_info)), repeats=group_info)
517
            else:
518
                flatted_group = np.zeros(num_data, dtype=np.int32)
519
            folds = folds.split(X=np.empty(num_data), y=full_data.get_label(), groups=flatted_group)
wxchan's avatar
wxchan committed
520
    else:
521
522
523
524
525
        if any(
            params.get(obj_alias, "")
            in {"lambdarank", "rank_xendcg", "xendcg", "xe_ndcg", "xe_ndcg_mart", "xendcg_mart"}
            for obj_alias in _ConfigAliases.get("objective")
        ):
wxchan's avatar
wxchan committed
526
            if not SKLEARN_INSTALLED:
527
                raise LightGBMError("scikit-learn is required for ranking cv")
528
            # ranking task, split according to groups
529
            group_info = np.asarray(full_data.get_group(), dtype=np.int32)
530
            flatted_group = np.repeat(range(len(group_info)), repeats=group_info)
531
            group_kfold = _LGBMGroupKFold(n_splits=nfold)
532
            folds = group_kfold.split(X=np.empty(num_data), groups=flatted_group)
wxchan's avatar
wxchan committed
533
534
        elif stratified:
            if not SKLEARN_INSTALLED:
535
                raise LightGBMError("scikit-learn is required for stratified cv")
536
            skf = _LGBMStratifiedKFold(n_splits=nfold, shuffle=shuffle, random_state=seed)
537
            folds = skf.split(X=np.empty(num_data), y=full_data.get_label())
extremin's avatar
extremin committed
538
        else:
wxchan's avatar
wxchan committed
539
540
541
542
543
            if shuffle:
                randidx = np.random.RandomState(seed).permutation(num_data)
            else:
                randidx = np.arange(num_data)
            kstep = int(num_data / nfold)
544
            test_id = [randidx[i : i + kstep] for i in range(0, num_data, kstep)]
545
546
            train_id = [np.concatenate([test_id[i] for i in range(nfold) if k != i]) for k in range(nfold)]
            folds = zip(train_id, test_id)
wxchan's avatar
wxchan committed
547

548
    ret = CVBooster()
wxchan's avatar
wxchan committed
549
    for train_idx, test_idx in folds:
550
551
        train_set = full_data.subset(sorted(train_idx))
        valid_set = full_data.subset(sorted(test_idx))
wxchan's avatar
wxchan committed
552
553
        # run preprocessing on the data set if needed
        if fpreproc is not None:
wxchan's avatar
wxchan committed
554
            train_set, valid_set, tparam = fpreproc(train_set, valid_set, params.copy())
wxchan's avatar
wxchan committed
555
        else:
wxchan's avatar
wxchan committed
556
            tparam = params
557
        booster_for_fold = Booster(tparam, train_set)
558
        if eval_train_metric:
559
560
            booster_for_fold.add_valid(train_set, "train")
        booster_for_fold.add_valid(valid_set, "valid")
561
        ret.boosters.append(booster_for_fold)
wxchan's avatar
wxchan committed
562
563
    return ret

wxchan's avatar
wxchan committed
564

565
def _agg_cv_result(
566
    raw_results: List[List[_LGBM_BoosterEvalMethodResultType]],
567
) -> List[_LGBM_BoosterEvalMethodResultWithStandardDeviationType]:
568
    """Aggregate cross-validation results."""
569
    cvmap: Dict[str, List[float]] = OrderedDict()
570
    metric_type: Dict[str, bool] = {}
wxchan's avatar
wxchan committed
571
572
    for one_result in raw_results:
        for one_line in one_result:
573
            key = f"{one_line[0]} {one_line[1]}"
574
            metric_type[key] = one_line[3]
575
            cvmap.setdefault(key, [])
576
            cvmap[key].append(one_line[2])
577
    return [("cv_agg", k, float(np.mean(v)), metric_type[k], float(np.std(v))) for k, v in cvmap.items()]
wxchan's avatar
wxchan committed
578

wxchan's avatar
wxchan committed
579

580
581
582
583
584
585
586
587
588
589
590
def cv(
    params: Dict[str, Any],
    train_set: Dataset,
    num_boost_round: int = 100,
    folds: Optional[Union[Iterable[Tuple[np.ndarray, np.ndarray]], _LGBMBaseCrossValidator]] = None,
    nfold: int = 5,
    stratified: bool = True,
    shuffle: bool = True,
    metrics: Optional[Union[str, List[str]]] = None,
    feval: Optional[Union[_LGBM_CustomMetricFunction, List[_LGBM_CustomMetricFunction]]] = None,
    init_model: Optional[Union[str, Path, Booster]] = None,
591
592
    feature_name: _LGBM_FeatureNameConfiguration = "auto",
    categorical_feature: _LGBM_CategoricalFeatureConfiguration = "auto",
593
594
595
596
    fpreproc: Optional[_LGBM_PreprocFunction] = None,
    seed: int = 0,
    callbacks: Optional[List[Callable]] = None,
    eval_train_metric: bool = False,
597
    return_cvbooster: bool = False,
598
) -> Dict[str, Union[List[float], CVBooster]]:
Andrew Ziem's avatar
Andrew Ziem committed
599
    """Perform the cross-validation with given parameters.
wxchan's avatar
wxchan committed
600
601
602
603

    Parameters
    ----------
    params : dict
604
605
        Parameters for training. Values passed through ``params`` take precedence over those
        supplied via arguments.
Guolin Ke's avatar
Guolin Ke committed
606
    train_set : Dataset
607
        Data to be trained on.
608
    num_boost_round : int, optional (default=100)
wxchan's avatar
wxchan committed
609
        Number of boosting iterations.
610
    folds : generator or iterator of (train_idx, test_idx) tuples, scikit-learn splitter object or None, optional (default=None)
611
        If generator or iterator, it should yield the train and test indices for each fold.
612
        If object, it should be one of the scikit-learn splitter classes
613
        (https://scikit-learn.org/stable/modules/classes.html#splitter-classes)
614
        and have ``split`` method.
615
        This argument has highest priority over other data split arguments.
616
    nfold : int, optional (default=5)
wxchan's avatar
wxchan committed
617
        Number of folds in CV.
618
619
    stratified : bool, optional (default=True)
        Whether to perform stratified sampling.
620
    shuffle : bool, optional (default=True)
621
        Whether to shuffle before splitting data.
622
    metrics : str, list of str, or None, optional (default=None)
623
624
        Evaluation metrics to be monitored while CV.
        If not None, the metric in ``params`` will be overridden.
625
    feval : callable, list of callable, or None, optional (default=None)
626
        Customized evaluation function.
627
        Each evaluation function should accept two parameters: preds, eval_data,
628
        and return (eval_name, eval_result, is_higher_better) or list of such tuples.
629

630
            preds : numpy 1-D array or numpy 2-D array (for multi-class task)
631
                The predicted values.
632
                For multi-class task, preds are numpy 2-D array of shape = [n_samples, n_classes].
633
                If custom objective function is used, predicted values are returned before any transformation,
634
                e.g. they are raw margin instead of probability of positive class for binary task in this case.
635
636
            eval_data : Dataset
                A ``Dataset`` to evaluate.
637
            eval_name : str
Andrew Ziem's avatar
Andrew Ziem committed
638
                The name of evaluation function (without whitespace).
639
640
641
642
643
            eval_result : float
                The eval result.
            is_higher_better : bool
                Is eval result higher better, e.g. AUC is ``is_higher_better``.

644
645
        To ignore the default metric corresponding to the used objective,
        set ``metrics`` to the string ``"None"``.
646
    init_model : str, pathlib.Path, Booster or None, optional (default=None)
647
        Filename of LightGBM model or Booster instance used for continue training.
648
    feature_name : list of str, or 'auto', optional (default="auto")
649
        **Deprecated.** Set ``feature_name`` on ``train_set`` instead.
650
651
        Feature names.
        If 'auto' and data is pandas DataFrame, data columns names are used.
652
    categorical_feature : list of str or int, or 'auto', optional (default="auto")
653
        **Deprecated.** Set ``categorical_feature`` on ``train_set`` instead.
654
655
        Categorical features.
        If list of int, interpreted as indices.
656
        If list of str, interpreted as feature names (need to specify ``feature_name`` as well).
657
        If 'auto' and data is pandas DataFrame, pandas unordered categorical columns are used.
658
        All values in categorical features will be cast to int32 and thus should be less than int32 max value (2147483647).
659
        Large values could be memory consuming. Consider using consecutive integers starting from zero.
660
        All negative values in categorical features will be treated as missing values.
661
        The output cannot be monotonically constrained with respect to a categorical feature.
662
        Floating point numbers in categorical features will be rounded towards 0.
663
664
    fpreproc : callable or None, optional (default=None)
        Preprocessing function that takes (dtrain, dtest, params)
wxchan's avatar
wxchan committed
665
        and returns transformed versions of those.
666
    seed : int, optional (default=0)
wxchan's avatar
wxchan committed
667
        Seed used to generate the folds (passed to numpy.random.seed).
668
    callbacks : list of callable, or None, optional (default=None)
669
        List of callback functions that are applied at each iteration.
670
        See Callbacks in Python API for more information.
671
672
673
    eval_train_metric : bool, optional (default=False)
        Whether to display the train metric in progress.
        The score of the metric is calculated again after each training step, so there is some impact on performance.
674
675
    return_cvbooster : bool, optional (default=False)
        Whether to return Booster models trained on each fold through ``CVBooster``.
wxchan's avatar
wxchan committed
676

677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
    Note
    ----
    A custom objective function can be provided for the ``objective`` parameter.
    It should accept two parameters: preds, train_data and return (grad, hess).

        preds : numpy 1-D array or numpy 2-D array (for multi-class task)
            The predicted values.
            Predicted values are returned before any transformation,
            e.g. they are raw margin instead of probability of positive class for binary task.
        train_data : Dataset
            The training dataset.
        grad : numpy 1-D array or numpy 2-D array (for multi-class task)
            The value of the first order derivative (gradient) of the loss
            with respect to the elements of preds for each sample point.
        hess : numpy 1-D array or numpy 2-D array (for multi-class task)
            The value of the second order derivative (Hessian) of the loss
            with respect to the elements of preds for each sample point.

    For multi-class task, preds are numpy 2-D array of shape = [n_samples, n_classes],
    and grad and hess should be returned in the same format.

wxchan's avatar
wxchan committed
698
699
    Returns
    -------
700
701
    eval_results : dict
        History of evaluation results of each metric.
702
        The dictionary has the following format:
703
704
        {'valid metric1-mean': [values], 'valid metric1-stdv': [values],
        'valid metric2-mean': [values], 'valid metric2-stdv': [values],
705
        ...}.
706
        If ``return_cvbooster=True``, also returns trained boosters wrapped in a ``CVBooster`` object via ``cvbooster`` key.
707
708
709
710
711
        If ``eval_train_metric=True``, also returns the train metric history.
        In this case, the dictionary has the following format:
        {'train metric1-mean': [values], 'valid metric1-mean': [values],
        'train metric2-mean': [values], 'valid metric2-mean': [values],
        ...}.
wxchan's avatar
wxchan committed
712
    """
Guolin Ke's avatar
Guolin Ke committed
713
    if not isinstance(train_set, Dataset):
714
715
716
717
718
        raise TypeError(f"cv() only accepts Dataset object, train_set has type '{type(train_set).__name__}'.")

    if num_boost_round <= 0:
        raise ValueError(f"num_boost_round must be greater than 0. Got {num_boost_round}.")

719
720
721
722
723
724
725
    # raise deprecation warnings if necessary
    # ref: https://github.com/microsoft/LightGBM/issues/6435
    if categorical_feature != "auto":
        _emit_dataset_kwarg_warning("cv", "categorical_feature")
    if feature_name != "auto":
        _emit_dataset_kwarg_warning("cv", "feature_name")

726
    params = copy.deepcopy(params)
727
    params = _choose_param_value(
728
        main_param_name="objective",
729
        params=params,
730
        default_value=None,
731
    )
732
    fobj: Optional[_LGBM_CustomObjectiveFunction] = None
733
734
    if callable(params["objective"]):
        fobj = params["objective"]
735
        params["objective"] = "none"
736
    for alias in _ConfigAliases.get("num_iterations"):
737
        if alias in params:
738
            _log_warning(f"Found '{alias}' in params. Will use it instead of 'num_boost_round' argument")
739
            num_boost_round = params.pop(alias)
740
    params["num_iterations"] = num_boost_round
741
742
743
744
    # setting early stopping via global params should be possible
    params = _choose_param_value(
        main_param_name="early_stopping_round",
        params=params,
745
        default_value=None,
746
747
748
    )
    if params["early_stopping_round"] is None:
        params.pop("early_stopping_round")
749
    first_metric_only = params.get("first_metric_only", False)
750

751
    if isinstance(init_model, (str, Path)):
752
753
        predictor = _InnerPredictor.from_model_file(
            model_file=init_model,
754
            pred_parameter=params,
755
        )
Guolin Ke's avatar
Guolin Ke committed
756
    elif isinstance(init_model, Booster):
757
758
        predictor = _InnerPredictor.from_booster(
            booster=init_model,
759
            pred_parameter=dict(init_model.params, **params),
760
        )
Guolin Ke's avatar
Guolin Ke committed
761
762
763
    else:
        predictor = None

Peter's avatar
Peter committed
764
    if metrics is not None:
765
766
        for metric_alias in _ConfigAliases.get("metric"):
            params.pop(metric_alias, None)
767
        params["metric"] = metrics
wxchan's avatar
wxchan committed
768

769
770
771
    train_set._update_params(params)._set_predictor(predictor).set_feature_name(feature_name).set_categorical_feature(
        categorical_feature
    )
772

773
    results = defaultdict(list)
774
775
776
777
778
779
780
781
782
783
784
    cvfolds = _make_n_folds(
        full_data=train_set,
        folds=folds,
        nfold=nfold,
        params=params,
        seed=seed,
        fpreproc=fpreproc,
        stratified=stratified,
        shuffle=shuffle,
        eval_train_metric=eval_train_metric,
    )
wxchan's avatar
wxchan committed
785
786

    # setup callbacks
787
    if callbacks is None:
788
        callbacks_set = set()
wxchan's avatar
wxchan committed
789
790
    else:
        for i, cb in enumerate(callbacks):
791
            cb.__dict__.setdefault("order", i - len(callbacks))
792
        callbacks_set = set(callbacks)
793

794
    if callback._should_enable_early_stopping(params.get("early_stopping_round", 0)):
795
        callbacks_set.add(
796
            callback.early_stopping(
797
                stopping_rounds=params["early_stopping_round"],  # type: ignore[arg-type]
798
                first_metric_only=first_metric_only,
799
                min_delta=params.get("early_stopping_min_delta", 0.0),
800
801
802
                verbose=_choose_param_value(
                    main_param_name="verbosity",
                    params=params,
803
804
805
                    default_value=1,
                ).pop("verbosity")
                > 0,
806
807
            )
        )
wxchan's avatar
wxchan committed
808

809
    callbacks_before_iter_set = {cb for cb in callbacks_set if getattr(cb, "before_iteration", False)}
810
    callbacks_after_iter_set = callbacks_set - callbacks_before_iter_set
811
812
    callbacks_before_iter = sorted(callbacks_before_iter_set, key=attrgetter("order"))
    callbacks_after_iter = sorted(callbacks_after_iter_set, key=attrgetter("order"))
wxchan's avatar
wxchan committed
813

814
    for i in range(num_boost_round):
wxchan's avatar
wxchan committed
815
        for cb in callbacks_before_iter:
816
817
818
819
820
821
822
823
824
825
            cb(
                callback.CallbackEnv(
                    model=cvfolds,
                    params=params,
                    iteration=i,
                    begin_iteration=0,
                    end_iteration=num_boost_round,
                    evaluation_result_list=None,
                )
            )
826
827
        cvfolds.update(fobj=fobj)  # type: ignore[call-arg]
        res = _agg_cv_result(cvfolds.eval_valid(feval))  # type: ignore[call-arg]
wxchan's avatar
wxchan committed
828
        for _, key, mean, _, std in res:
829
830
            results[f"{key}-mean"].append(mean)
            results[f"{key}-stdv"].append(std)
wxchan's avatar
wxchan committed
831
832
        try:
            for cb in callbacks_after_iter:
833
834
835
836
837
838
839
840
841
842
                cb(
                    callback.CallbackEnv(
                        model=cvfolds,
                        params=params,
                        iteration=i,
                        begin_iteration=0,
                        end_iteration=num_boost_round,
                        evaluation_result_list=res,
                    )
                )
843
844
        except callback.EarlyStopException as earlyStopException:
            cvfolds.best_iteration = earlyStopException.best_iteration + 1
845
846
            for bst in cvfolds.boosters:
                bst.best_iteration = cvfolds.best_iteration
wxchan's avatar
wxchan committed
847
            for k in results:
848
                results[k] = results[k][: cvfolds.best_iteration]
wxchan's avatar
wxchan committed
849
            break
850
851

    if return_cvbooster:
852
        results["cvbooster"] = cvfolds  # type: ignore[assignment]
853

wxchan's avatar
wxchan committed
854
    return dict(results)