engine.py 34.3 KB
Newer Older
wxchan's avatar
wxchan committed
1
# coding: utf-8
2
"""Library with training routines of LightGBM."""
3
import copy
4
import json
5
from collections import OrderedDict, defaultdict
wxchan's avatar
wxchan committed
6
from operator import attrgetter
7
from pathlib import Path
8
from typing import Any, Callable, Dict, Iterable, List, Optional, Tuple, Union
9

wxchan's avatar
wxchan committed
10
import numpy as np
11

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

30
31
32
33
34
35
36
__all__ = [
    'cv',
    'CVBooster',
    'train',
]


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

48
49
50
51
52
_LGBM_PreprocFunction = Callable[
    [Dataset, Dataset, Dict[str, Any]],
    Tuple[Dataset, Dataset, Dict[str, Any]]
]

53
54
55
56
57
58
59
60
61

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,
62
63
    feature_name: _LGBM_FeatureNameConfiguration = 'auto',
    categorical_feature: _LGBM_CategoricalFeatureConfiguration = 'auto',
64
65
66
    keep_training_booster: bool = False,
    callbacks: Optional[List[Callable]] = None
) -> Booster:
67
    """Perform the training with given parameters.
wxchan's avatar
wxchan committed
68
69
70
71

    Parameters
    ----------
    params : dict
72
73
        Parameters for training. Values passed through ``params`` take precedence over those
        supplied via arguments.
Guolin Ke's avatar
Guolin Ke committed
74
    train_set : Dataset
75
76
        Data to be trained on.
    num_boost_round : int, optional (default=100)
wxchan's avatar
wxchan committed
77
        Number of boosting iterations.
78
    valid_sets : list of Dataset, or None, optional (default=None)
79
        List of data to be evaluated on during training.
80
    valid_names : list of str, or None, optional (default=None)
81
        Names of ``valid_sets``.
82
    feval : callable, list of callable, or None, optional (default=None)
wxchan's avatar
wxchan committed
83
        Customized evaluation function.
Akshita Dixit's avatar
Akshita Dixit committed
84
        Each evaluation function should accept two parameters: preds, eval_data,
85
        and return (eval_name, eval_result, is_higher_better) or list of such tuples.
86

87
            preds : numpy 1-D array or numpy 2-D array (for multi-class task)
88
                The predicted values.
89
                For multi-class task, preds are numpy 2-D array of shape = [n_samples, n_classes].
90
                If custom objective function is used, predicted values are returned before any transformation,
91
                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
92
            eval_data : Dataset
93
                A ``Dataset`` to evaluate.
94
            eval_name : str
95
                The name of evaluation function (without whitespaces).
96
97
98
99
100
            eval_result : float
                The eval result.
            is_higher_better : bool
                Is eval result higher better, e.g. AUC is ``is_higher_better``.

101
102
        To ignore the default metric corresponding to the used objective,
        set the ``metric`` parameter to the string ``"None"`` in ``params``.
103
    init_model : str, pathlib.Path, Booster or None, optional (default=None)
104
        Filename of LightGBM model or Booster instance used for continue training.
105
    feature_name : list of str, or 'auto', optional (default="auto")
106
107
        Feature names.
        If 'auto' and data is pandas DataFrame, data columns names are used.
108
    categorical_feature : list of str or int, or 'auto', optional (default="auto")
109
110
        Categorical features.
        If list of int, interpreted as indices.
111
        If list of str, interpreted as feature names (need to specify ``feature_name`` as well).
112
        If 'auto' and data is pandas DataFrame, pandas unordered categorical columns are used.
113
        All values in categorical features will be cast to int32 and thus should be less than int32 max value (2147483647).
114
        Large values could be memory consuming. Consider using consecutive integers starting from zero.
115
        All negative values in categorical features will be treated as missing values.
116
        The output cannot be monotonically constrained with respect to a categorical feature.
117
        Floating point numbers in categorical features will be rounded towards 0.
118
119
120
    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.
121
        This means you won't be able to use ``eval``, ``eval_train`` or ``eval_valid`` methods of the returned Booster.
122
123
        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``.
124
        You can still use _InnerPredictor as ``init_model`` for future continue training.
125
    callbacks : list of callable, or None, optional (default=None)
126
        List of callback functions that are applied at each iteration.
127
        See Callbacks in Python API for more information.
wxchan's avatar
wxchan committed
128

129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
    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
150
151
    Returns
    -------
152
153
    booster : Booster
        The trained Booster model.
wxchan's avatar
wxchan committed
154
    """
155
156
157
158
159
160
161
162
163
164
165
166
167
168
    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__}'."
                )

169
    # create predictor first
170
    params = copy.deepcopy(params)
171
172
173
174
175
    params = _choose_param_value(
        main_param_name='objective',
        params=params,
        default_value=None
    )
176
    fobj: Optional[_LGBM_CustomObjectiveFunction] = None
177
178
179
    if callable(params["objective"]):
        fobj = params["objective"]
        params["objective"] = 'none'
180
    for alias in _ConfigAliases.get("num_iterations"):
181
        if alias in params:
182
            num_boost_round = params.pop(alias)
183
            _log_warning(f"Found `{alias}` in params. Will use it instead of argument")
184
    params["num_iterations"] = num_boost_round
185
186
187
188
189
190
191
192
    # setting early stopping via global params should be possible
    params = _choose_param_value(
        main_param_name="early_stopping_round",
        params=params,
        default_value=None
    )
    if params["early_stopping_round"] is None:
        params.pop("early_stopping_round")
193
    first_metric_only = params.get('first_metric_only', False)
194

195
    predictor: Optional[_InnerPredictor] = None
196
    if isinstance(init_model, (str, Path)):
197
198
199
200
        predictor = _InnerPredictor.from_model_file(
            model_file=init_model,
            pred_parameter=params
        )
wxchan's avatar
wxchan committed
201
    elif isinstance(init_model, Booster):
202
203
204
205
206
207
208
209
210
        predictor = _InnerPredictor.from_booster(
            booster=init_model,
            pred_parameter=dict(init_model.params, **params)
        )

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

212
213
214
215
    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
216

wxchan's avatar
wxchan committed
217
218
    is_valid_contain_train = False
    train_data_name = "training"
Guolin Ke's avatar
Guolin Ke committed
219
    reduced_valid_sets = []
wxchan's avatar
wxchan committed
220
    name_valid_sets = []
221
    if valid_sets is not None:
Guolin Ke's avatar
Guolin Ke committed
222
223
        if isinstance(valid_sets, Dataset):
            valid_sets = [valid_sets]
224
        if isinstance(valid_names, str):
wxchan's avatar
wxchan committed
225
            valid_names = [valid_names]
Guolin Ke's avatar
Guolin Ke committed
226
        for i, valid_data in enumerate(valid_sets):
227
            # reduce cost for prediction training data
Guolin Ke's avatar
Guolin Ke committed
228
            if valid_data is train_set:
wxchan's avatar
wxchan committed
229
230
231
232
                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
233
            reduced_valid_sets.append(valid_data._update_params(params).set_reference(train_set))
234
            if valid_names is not None and len(valid_names) > i:
wxchan's avatar
wxchan committed
235
236
                name_valid_sets.append(valid_names[i])
            else:
237
                name_valid_sets.append(f'valid_{i}')
238
    # process callbacks
239
    if callbacks is None:
240
        callbacks_set = set()
wxchan's avatar
wxchan committed
241
242
243
    else:
        for i, cb in enumerate(callbacks):
            cb.__dict__.setdefault('order', i - len(callbacks))
244
        callbacks_set = set(callbacks)
wxchan's avatar
wxchan committed
245

246
247
248
    if "early_stopping_round" in params:
        callbacks_set.add(
            callback.early_stopping(
249
                stopping_rounds=params["early_stopping_round"],  # type: ignore[arg-type]
250
251
252
253
254
255
256
257
                first_metric_only=first_metric_only,
                verbose=_choose_param_value(
                    main_param_name="verbosity",
                    params=params,
                    default_value=1
                ).pop("verbosity") > 0
            )
        )
258

259
260
261
262
    callbacks_before_iter_set = {cb for cb in callbacks_set if getattr(cb, 'before_iteration', False)}
    callbacks_after_iter_set = callbacks_set - callbacks_before_iter_set
    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
263

264
    # construct booster
265
266
267
268
    try:
        booster = Booster(params=params, train_set=train_set)
        if is_valid_contain_train:
            booster.set_train_data_name(train_data_name)
269
        for valid_set, name_valid_set in zip(reduced_valid_sets, name_valid_sets):
270
271
272
273
274
            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()
275
    booster.best_iteration = 0
wxchan's avatar
wxchan committed
276

277
    # start training
278
    for i in range(init_iteration, init_iteration + num_boost_round):
wxchan's avatar
wxchan committed
279
280
        for cb in callbacks_before_iter:
            cb(callback.CallbackEnv(model=booster,
281
                                    params=params,
wxchan's avatar
wxchan committed
282
                                    iteration=i,
283
284
                                    begin_iteration=init_iteration,
                                    end_iteration=init_iteration + num_boost_round,
wxchan's avatar
wxchan committed
285
286
287
288
                                    evaluation_result_list=None))

        booster.update(fobj=fobj)

289
        evaluation_result_list: List[_LGBM_BoosterEvalMethodResultType] = []
wxchan's avatar
wxchan committed
290
        # check evaluation result.
291
        if valid_sets is not None:
wxchan's avatar
wxchan committed
292
293
294
295
296
297
            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:
                cb(callback.CallbackEnv(model=booster,
298
                                        params=params,
wxchan's avatar
wxchan committed
299
                                        iteration=i,
300
301
                                        begin_iteration=init_iteration,
                                        end_iteration=init_iteration + num_boost_round,
wxchan's avatar
wxchan committed
302
                                        evaluation_result_list=evaluation_result_list))
303
304
        except callback.EarlyStopException as earlyStopException:
            booster.best_iteration = earlyStopException.best_iteration + 1
wxchan's avatar
wxchan committed
305
            evaluation_result_list = earlyStopException.best_score
wxchan's avatar
wxchan committed
306
            break
307
    booster.best_score = defaultdict(OrderedDict)
wxchan's avatar
wxchan committed
308
309
    for dataset_name, eval_name, score, _ in evaluation_result_list:
        booster.best_score[dataset_name][eval_name] = score
310
    if not keep_training_booster:
311
        booster.model_from_string(booster.model_to_string()).free_dataset()
wxchan's avatar
wxchan committed
312
313
314
    return booster


315
class CVBooster:
316
317
    """CVBooster in LightGBM.

318
    Auxiliary data structure to hold and redirect all boosters of ``cv()`` function.
319
    This class has the same methods as Booster class.
320
321
322
323
324
325
    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()``
326
327
328
329
330
331
332
333

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

335
336
337
338
    def __init__(
        self,
        model_file: Optional[Union[str, Path]] = None
    ):
339
340
        """Initialize the CVBooster.

341
342
343
344
        Parameters
        ----------
        model_file : str, pathlib.Path or None, optional (default=None)
            Path to the CVBooster model file.
345
        """
346
        self.boosters: List[Booster] = []
347
        self.best_iteration = -1
348

349
350
351
352
353
354
355
356
357
        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"]:
358
            self.boosters.append(Booster(model_str=model_str))
359
360
361
362
363
364
365
366
367

    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:
            models_str.append(booster.model_to_string(num_iteration=num_iteration, start_iteration=start_iteration,
                                                      importance_type=importance_type))
        return {"boosters": models_str, "best_iteration": self.best_iteration}

368
    def __getattr__(self, name: str) -> Callable[[Any, Any], List[Any]]:
369
        """Redirect methods call of CVBooster."""
370
        def handler_function(*args: Any, **kwargs: Any) -> List[Any]:
371
            """Call methods with each booster, and concatenate their results."""
372
373
374
375
            ret = []
            for booster in self.boosters:
                ret.append(getattr(booster, name)(*args, **kwargs))
            return ret
376
        return handler_function
wxchan's avatar
wxchan committed
377

378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
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
433
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
    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,
        importance_type: str = 'split'
    ) -> 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,
        importance_type: str = 'split'
    ) -> "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

462

463
464
465
466
467
468
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,
469
470
471
472
    fpreproc: Optional[_LGBM_PreprocFunction],
    stratified: bool,
    shuffle: bool,
    eval_train_metric: bool
473
) -> CVBooster:
474
    """Make a n-fold list of Booster from random indices."""
wxchan's avatar
wxchan committed
475
476
    full_data = full_data.construct()
    num_data = full_data.num_data()
477
    if folds is not None:
478
479
480
481
482
483
        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'):
            group_info = full_data.get_group()
            if group_info is not None:
484
                group_info = np.array(group_info, dtype=np.int32, copy=False)
485
                flatted_group = np.repeat(range(len(group_info)), repeats=group_info)
486
            else:
487
                flatted_group = np.zeros(num_data, dtype=np.int32)
488
            folds = folds.split(X=np.empty(num_data), y=full_data.get_label(), groups=flatted_group)
wxchan's avatar
wxchan committed
489
    else:
490
491
492
        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
493
            if not SKLEARN_INSTALLED:
494
                raise LightGBMError('scikit-learn is required for ranking cv')
495
            # ranking task, split according to groups
496
            group_info = np.array(full_data.get_group(), dtype=np.int32, copy=False)
497
            flatted_group = np.repeat(range(len(group_info)), repeats=group_info)
498
            group_kfold = _LGBMGroupKFold(n_splits=nfold)
499
            folds = group_kfold.split(X=np.empty(num_data), groups=flatted_group)
wxchan's avatar
wxchan committed
500
501
        elif stratified:
            if not SKLEARN_INSTALLED:
502
                raise LightGBMError('scikit-learn is required for stratified cv')
503
            skf = _LGBMStratifiedKFold(n_splits=nfold, shuffle=shuffle, random_state=seed)
504
            folds = skf.split(X=np.empty(num_data), y=full_data.get_label())
extremin's avatar
extremin committed
505
        else:
wxchan's avatar
wxchan committed
506
507
508
509
510
            if shuffle:
                randidx = np.random.RandomState(seed).permutation(num_data)
            else:
                randidx = np.arange(num_data)
            kstep = int(num_data / nfold)
511
512
513
            test_id = [randidx[i: i + kstep] for i in range(0, num_data, kstep)]
            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
514

515
    ret = CVBooster()
wxchan's avatar
wxchan committed
516
    for train_idx, test_idx in folds:
517
518
        train_set = full_data.subset(sorted(train_idx))
        valid_set = full_data.subset(sorted(test_idx))
wxchan's avatar
wxchan committed
519
520
        # run preprocessing on the data set if needed
        if fpreproc is not None:
wxchan's avatar
wxchan committed
521
            train_set, valid_set, tparam = fpreproc(train_set, valid_set, params.copy())
wxchan's avatar
wxchan committed
522
        else:
wxchan's avatar
wxchan committed
523
            tparam = params
524
        booster_for_fold = Booster(tparam, train_set)
525
        if eval_train_metric:
526
527
528
            booster_for_fold.add_valid(train_set, 'train')
        booster_for_fold.add_valid(valid_set, 'valid')
        ret.boosters.append(booster_for_fold)
wxchan's avatar
wxchan committed
529
530
    return ret

wxchan's avatar
wxchan committed
531

532
def _agg_cv_result(
533
534
    raw_results: List[List[_LGBM_BoosterEvalMethodResultType]]
) -> List[_LGBM_BoosterEvalMethodResultWithStandardDeviationType]:
535
    """Aggregate cross-validation results."""
536
    cvmap: Dict[str, List[float]] = OrderedDict()
537
    metric_type: Dict[str, bool] = {}
wxchan's avatar
wxchan committed
538
539
    for one_result in raw_results:
        for one_line in one_result:
540
            key = f"{one_line[0]} {one_line[1]}"
541
            metric_type[key] = one_line[3]
542
            cvmap.setdefault(key, [])
543
            cvmap[key].append(one_line[2])
544
    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
545

wxchan's avatar
wxchan committed
546

547
548
549
550
551
552
553
554
555
556
557
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,
558
559
    feature_name: _LGBM_FeatureNameConfiguration = 'auto',
    categorical_feature: _LGBM_CategoricalFeatureConfiguration = 'auto',
560
561
562
563
564
    fpreproc: Optional[_LGBM_PreprocFunction] = None,
    seed: int = 0,
    callbacks: Optional[List[Callable]] = None,
    eval_train_metric: bool = False,
    return_cvbooster: bool = False
565
) -> Dict[str, Union[List[float], CVBooster]]:
Andrew Ziem's avatar
Andrew Ziem committed
566
    """Perform the cross-validation with given parameters.
wxchan's avatar
wxchan committed
567
568
569
570

    Parameters
    ----------
    params : dict
571
572
        Parameters for training. Values passed through ``params`` take precedence over those
        supplied via arguments.
Guolin Ke's avatar
Guolin Ke committed
573
    train_set : Dataset
574
        Data to be trained on.
575
    num_boost_round : int, optional (default=100)
wxchan's avatar
wxchan committed
576
        Number of boosting iterations.
577
    folds : generator or iterator of (train_idx, test_idx) tuples, scikit-learn splitter object or None, optional (default=None)
578
        If generator or iterator, it should yield the train and test indices for each fold.
579
        If object, it should be one of the scikit-learn splitter classes
580
        (https://scikit-learn.org/stable/modules/classes.html#splitter-classes)
581
        and have ``split`` method.
582
        This argument has highest priority over other data split arguments.
583
    nfold : int, optional (default=5)
wxchan's avatar
wxchan committed
584
        Number of folds in CV.
585
586
    stratified : bool, optional (default=True)
        Whether to perform stratified sampling.
587
    shuffle : bool, optional (default=True)
588
        Whether to shuffle before splitting data.
589
    metrics : str, list of str, or None, optional (default=None)
590
591
        Evaluation metrics to be monitored while CV.
        If not None, the metric in ``params`` will be overridden.
592
    feval : callable, list of callable, or None, optional (default=None)
593
        Customized evaluation function.
594
        Each evaluation function should accept two parameters: preds, eval_data,
595
        and return (eval_name, eval_result, is_higher_better) or list of such tuples.
596

597
            preds : numpy 1-D array or numpy 2-D array (for multi-class task)
598
                The predicted values.
599
                For multi-class task, preds are numpy 2-D array of shape = [n_samples, n_classes].
600
                If custom objective function is used, predicted values are returned before any transformation,
601
                e.g. they are raw margin instead of probability of positive class for binary task in this case.
602
603
            eval_data : Dataset
                A ``Dataset`` to evaluate.
604
            eval_name : str
Andrew Ziem's avatar
Andrew Ziem committed
605
                The name of evaluation function (without whitespace).
606
607
608
609
610
            eval_result : float
                The eval result.
            is_higher_better : bool
                Is eval result higher better, e.g. AUC is ``is_higher_better``.

611
612
        To ignore the default metric corresponding to the used objective,
        set ``metrics`` to the string ``"None"``.
613
    init_model : str, pathlib.Path, Booster or None, optional (default=None)
614
        Filename of LightGBM model or Booster instance used for continue training.
615
    feature_name : list of str, or 'auto', optional (default="auto")
616
617
        Feature names.
        If 'auto' and data is pandas DataFrame, data columns names are used.
618
    categorical_feature : list of str or int, or 'auto', optional (default="auto")
619
620
        Categorical features.
        If list of int, interpreted as indices.
621
        If list of str, interpreted as feature names (need to specify ``feature_name`` as well).
622
        If 'auto' and data is pandas DataFrame, pandas unordered categorical columns are used.
623
        All values in categorical features will be cast to int32 and thus should be less than int32 max value (2147483647).
624
        Large values could be memory consuming. Consider using consecutive integers starting from zero.
625
        All negative values in categorical features will be treated as missing values.
626
        The output cannot be monotonically constrained with respect to a categorical feature.
627
        Floating point numbers in categorical features will be rounded towards 0.
628
629
    fpreproc : callable or None, optional (default=None)
        Preprocessing function that takes (dtrain, dtest, params)
wxchan's avatar
wxchan committed
630
        and returns transformed versions of those.
631
    seed : int, optional (default=0)
wxchan's avatar
wxchan committed
632
        Seed used to generate the folds (passed to numpy.random.seed).
633
    callbacks : list of callable, or None, optional (default=None)
634
        List of callback functions that are applied at each iteration.
635
        See Callbacks in Python API for more information.
636
637
638
    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.
639
640
    return_cvbooster : bool, optional (default=False)
        Whether to return Booster models trained on each fold through ``CVBooster``.
wxchan's avatar
wxchan committed
641

642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
    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
663
664
    Returns
    -------
665
666
    eval_results : dict
        History of evaluation results of each metric.
667
        The dictionary has the following format:
668
669
        {'valid metric1-mean': [values], 'valid metric1-stdv': [values],
        'valid metric2-mean': [values], 'valid metric2-stdv': [values],
670
        ...}.
671
        If ``return_cvbooster=True``, also returns trained boosters wrapped in a ``CVBooster`` object via ``cvbooster`` key.
672
673
674
675
676
        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
677
    """
Guolin Ke's avatar
Guolin Ke committed
678
    if not isinstance(train_set, Dataset):
679
680
681
682
683
        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}.")

684
    params = copy.deepcopy(params)
685
686
687
688
689
    params = _choose_param_value(
        main_param_name='objective',
        params=params,
        default_value=None
    )
690
    fobj: Optional[_LGBM_CustomObjectiveFunction] = None
691
692
693
    if callable(params["objective"]):
        fobj = params["objective"]
        params["objective"] = 'none'
694
    for alias in _ConfigAliases.get("num_iterations"):
695
        if alias in params:
696
            _log_warning(f"Found '{alias}' in params. Will use it instead of 'num_boost_round' argument")
697
            num_boost_round = params.pop(alias)
698
    params["num_iterations"] = num_boost_round
699
700
701
702
703
704
705
706
    # setting early stopping via global params should be possible
    params = _choose_param_value(
        main_param_name="early_stopping_round",
        params=params,
        default_value=None
    )
    if params["early_stopping_round"] is None:
        params.pop("early_stopping_round")
707
    first_metric_only = params.get('first_metric_only', False)
708

709
    if isinstance(init_model, (str, Path)):
710
711
712
713
        predictor = _InnerPredictor.from_model_file(
            model_file=init_model,
            pred_parameter=params
        )
Guolin Ke's avatar
Guolin Ke committed
714
    elif isinstance(init_model, Booster):
715
716
717
718
        predictor = _InnerPredictor.from_booster(
            booster=init_model,
            pred_parameter=dict(init_model.params, **params)
        )
Guolin Ke's avatar
Guolin Ke committed
719
720
721
    else:
        predictor = None

Peter's avatar
Peter committed
722
    if metrics is not None:
723
724
        for metric_alias in _ConfigAliases.get("metric"):
            params.pop(metric_alias, None)
Peter's avatar
Peter committed
725
        params['metric'] = metrics
wxchan's avatar
wxchan committed
726

727
728
729
730
731
    train_set._update_params(params) \
             ._set_predictor(predictor) \
             .set_feature_name(feature_name) \
             .set_categorical_feature(categorical_feature)

732
    results = defaultdict(list)
733
    cvfolds = _make_n_folds(full_data=train_set, folds=folds, nfold=nfold,
734
                            params=params, seed=seed, fpreproc=fpreproc,
735
736
                            stratified=stratified, shuffle=shuffle,
                            eval_train_metric=eval_train_metric)
wxchan's avatar
wxchan committed
737
738

    # setup callbacks
739
    if callbacks is None:
740
        callbacks_set = set()
wxchan's avatar
wxchan committed
741
742
743
    else:
        for i, cb in enumerate(callbacks):
            cb.__dict__.setdefault('order', i - len(callbacks))
744
        callbacks_set = set(callbacks)
745
746

    if "early_stopping_round" in params:
747
        callbacks_set.add(
748
            callback.early_stopping(
749
                stopping_rounds=params["early_stopping_round"],  # type: ignore[arg-type]
750
751
752
753
754
755
756
757
                first_metric_only=first_metric_only,
                verbose=_choose_param_value(
                    main_param_name="verbosity",
                    params=params,
                    default_value=1
                ).pop("verbosity") > 0
            )
        )
wxchan's avatar
wxchan committed
758

759
760
761
762
    callbacks_before_iter_set = {cb for cb in callbacks_set if getattr(cb, 'before_iteration', False)}
    callbacks_after_iter_set = callbacks_set - callbacks_before_iter_set
    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
763

764
    for i in range(num_boost_round):
wxchan's avatar
wxchan committed
765
        for cb in callbacks_before_iter:
766
767
            cb(callback.CallbackEnv(model=cvfolds,
                                    params=params,
wxchan's avatar
wxchan committed
768
769
770
771
                                    iteration=i,
                                    begin_iteration=0,
                                    end_iteration=num_boost_round,
                                    evaluation_result_list=None))
772
773
        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
774
        for _, key, mean, _, std in res:
775
776
            results[f'{key}-mean'].append(mean)
            results[f'{key}-stdv'].append(std)
wxchan's avatar
wxchan committed
777
778
        try:
            for cb in callbacks_after_iter:
779
780
                cb(callback.CallbackEnv(model=cvfolds,
                                        params=params,
wxchan's avatar
wxchan committed
781
782
783
784
                                        iteration=i,
                                        begin_iteration=0,
                                        end_iteration=num_boost_round,
                                        evaluation_result_list=res))
785
786
        except callback.EarlyStopException as earlyStopException:
            cvfolds.best_iteration = earlyStopException.best_iteration + 1
787
788
            for bst in cvfolds.boosters:
                bst.best_iteration = cvfolds.best_iteration
wxchan's avatar
wxchan committed
789
            for k in results:
790
                results[k] = results[k][:cvfolds.best_iteration]
wxchan's avatar
wxchan committed
791
            break
792
793

    if return_cvbooster:
794
        results['cvbooster'] = cvfolds  # type: ignore[assignment]
795

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