sklearn.py 62.1 KB
Newer Older
wxchan's avatar
wxchan committed
1
# coding: utf-8
2
"""Scikit-learn wrapper interface for LightGBM."""
3
import copy
4
from inspect import signature
5
from pathlib import Path
6
from typing import Any, Callable, Dict, List, Optional, Tuple, Union
7

wxchan's avatar
wxchan committed
8
import numpy as np
9
import scipy.sparse
10

11
from .basic import (Booster, Dataset, LightGBMError, _choose_param_value, _ConfigAliases, _LGBM_BoosterBestScoreType,
12
                    _LGBM_CategoricalFeatureConfiguration, _LGBM_EvalFunctionResultType, _LGBM_FeatureNameConfiguration,
13
                    _LGBM_GroupType, _LGBM_InitScoreType, _LGBM_LabelType, _LGBM_WeightType, _log_warning)
14
from .callback import _EvalResultDict, record_evaluation
15
16
from .compat import (SKLEARN_INSTALLED, LGBMNotFittedError, _LGBMAssertAllFinite, _LGBMCheckArray,
                     _LGBMCheckClassificationTargets, _LGBMCheckSampleWeight, _LGBMCheckXY, _LGBMClassifierBase,
17
18
                     _LGBMComputeSampleWeight, _LGBMCpuCount, _LGBMLabelEncoder, _LGBMModelBase, _LGBMRegressorBase,
                     dt_DataTable, pd_DataFrame)
wxchan's avatar
wxchan committed
19
from .engine import train
20

21
22
23
24
25
26
27
__all__ = [
    'LGBMClassifier',
    'LGBMModel',
    'LGBMRanker',
    'LGBMRegressor',
]

28
29
30
31
32
33
34
_LGBM_ScikitMatrixLike = Union[
    dt_DataTable,
    List[Union[List[float], List[int]]],
    np.ndarray,
    pd_DataFrame,
    scipy.sparse.spmatrix
]
35
_LGBM_ScikitCustomObjectiveFunction = Union[
36
    # f(labels, preds)
37
    Callable[
38
        [Optional[np.ndarray], np.ndarray],
39
        Tuple[np.ndarray, np.ndarray]
40
    ],
41
    # f(labels, preds, weights)
42
    Callable[
43
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray]],
44
        Tuple[np.ndarray, np.ndarray]
45
    ],
46
    # f(labels, preds, weights, group)
47
    Callable[
48
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray], Optional[np.ndarray]],
49
50
        Tuple[np.ndarray, np.ndarray]
    ],
51
52
]
_LGBM_ScikitCustomEvalFunction = Union[
53
    # f(labels, preds)
54
    Callable[
55
56
        [Optional[np.ndarray], np.ndarray],
        _LGBM_EvalFunctionResultType
57
58
    ],
    Callable[
59
60
        [Optional[np.ndarray], np.ndarray],
        List[_LGBM_EvalFunctionResultType]
61
    ],
62
    # f(labels, preds, weights)
63
    Callable[
64
65
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray]],
        _LGBM_EvalFunctionResultType
66
    ],
67
68
69
70
71
72
73
74
75
76
77
78
79
    Callable[
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray]],
        List[_LGBM_EvalFunctionResultType]
    ],
    # f(labels, preds, weights, group)
    Callable[
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray], Optional[np.ndarray]],
        _LGBM_EvalFunctionResultType
    ],
    Callable[
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray], Optional[np.ndarray]],
        List[_LGBM_EvalFunctionResultType]
    ]
80
]
81
82
83
84
85
_LGBM_ScikitEvalMetricType = Union[
    str,
    _LGBM_ScikitCustomEvalFunction,
    List[Union[str, _LGBM_ScikitCustomEvalFunction]]
]
86
_LGBM_ScikitValidSet = Tuple[_LGBM_ScikitMatrixLike, _LGBM_LabelType]
87

wxchan's avatar
wxchan committed
88

89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
def _get_group_from_constructed_dataset(dataset: Dataset) -> Optional[np.ndarray]:
    group = dataset.get_group()
    error_msg = (
        "Estimators in lightgbm.sklearn should only retrieve query groups from a constructed Dataset. "
        "If you're seeing this message, it's a bug in lightgbm. Please report it at https://github.com/microsoft/LightGBM/issues."
    )
    assert (group is None or isinstance(group, np.ndarray)), error_msg
    return group


def _get_label_from_constructed_dataset(dataset: Dataset) -> np.ndarray:
    label = dataset.get_label()
    error_msg = (
        "Estimators in lightgbm.sklearn should only retrieve labels from a constructed Dataset. "
        "If you're seeing this message, it's a bug in lightgbm. Please report it at https://github.com/microsoft/LightGBM/issues."
    )
    assert isinstance(label, np.ndarray), error_msg
    return label


def _get_weight_from_constructed_dataset(dataset: Dataset) -> Optional[np.ndarray]:
    weight = dataset.get_weight()
    error_msg = (
        "Estimators in lightgbm.sklearn should only retrieve weights from a constructed Dataset. "
        "If you're seeing this message, it's a bug in lightgbm. Please report it at https://github.com/microsoft/LightGBM/issues."
    )
    assert (weight is None or isinstance(weight, np.ndarray)), error_msg
    return weight


119
class _ObjectiveFunctionWrapper:
120
    """Proxy class for objective function."""
121

122
    def __init__(self, func: _LGBM_ScikitCustomObjectiveFunction):
123
        """Construct a proxy class.
124

125
126
        This class transforms objective function to match objective function with signature ``new_func(preds, dataset)``
        as expected by ``lightgbm.engine.train``.
127

128
129
130
        Parameters
        ----------
        func : callable
131
132
133
134
            Expects a callable with following signatures:
            ``func(y_true, y_pred)``,
            ``func(y_true, y_pred, weight)``
            or ``func(y_true, y_pred, weight, group)``
135
136
            and returns (grad, hess):

137
                y_true : numpy 1-D array of shape = [n_samples]
138
                    The target values.
139
                y_pred : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
140
                    The predicted values.
141
142
                    Predicted values are returned before any transformation,
                    e.g. they are raw margin instead of probability of positive class for binary task.
143
144
                weight : numpy 1-D array of shape = [n_samples]
                    The weight of samples. Weights should be non-negative.
145
                group : numpy 1-D array
146
147
148
                    Group/query data.
                    Only used in the learning-to-rank task.
                    sum(group) = n_samples.
149
150
                    For example, if you have a 100-document dataset with ``group = [10, 20, 40, 10, 10, 10]``, that means that you have 6 groups,
                    where the first 10 records are in the first group, records 11-30 are in the second group, records 31-70 are in the third group, etc.
151
                grad : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape [n_samples, n_classes] (for multi-class task)
152
153
                    The value of the first order derivative (gradient) of the loss
                    with respect to the elements of y_pred for each sample point.
154
                hess : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
155
156
                    The value of the second order derivative (Hessian) of the loss
                    with respect to the elements of y_pred for each sample point.
wxchan's avatar
wxchan committed
157

Nikita Titov's avatar
Nikita Titov committed
158
159
        .. note::

160
            For multi-class task, y_pred is a numpy 2-D array of shape = [n_samples, n_classes],
161
            and grad and hess should be returned in the same format.
162
163
        """
        self.func = func
wxchan's avatar
wxchan committed
164

165
    def __call__(self, preds: np.ndarray, dataset: Dataset) -> Tuple[np.ndarray, np.ndarray]:
166
167
168
169
        """Call passed function with appropriate arguments.

        Parameters
        ----------
170
        preds : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
171
172
173
174
175
176
            The predicted values.
        dataset : Dataset
            The training dataset.

        Returns
        -------
177
        grad : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
178
179
            The value of the first order derivative (gradient) of the loss
            with respect to the elements of preds for each sample point.
180
        hess : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
181
182
            The value of the second order derivative (Hessian) of the loss
            with respect to the elements of preds for each sample point.
183
        """
184
        labels = _get_label_from_constructed_dataset(dataset)
185
        argc = len(signature(self.func).parameters)
186
        if argc == 2:
187
            grad, hess = self.func(labels, preds)  # type: ignore[call-arg]
188
189
190
191
192
193
194
195
196
197
198
199
            return grad, hess

        weight = _get_weight_from_constructed_dataset(dataset)
        if argc == 3:
            grad, hess = self.func(labels, preds, weight)  # type: ignore[call-arg]
            return grad, hess

        if argc == 4:
            group = _get_group_from_constructed_dataset(dataset)
            return self.func(labels, preds, weight, group)  # type: ignore[call-arg]

        raise TypeError(f"Self-defined objective function should have 2, 3 or 4 arguments, got {argc}")
wxchan's avatar
wxchan committed
200

wxchan's avatar
wxchan committed
201

202
class _EvalFunctionWrapper:
203
    """Proxy class for evaluation function."""
204

205
    def __init__(self, func: _LGBM_ScikitCustomEvalFunction):
206
        """Construct a proxy class.
207

208
209
        This class transforms evaluation function to match evaluation function with signature ``new_func(preds, dataset)``
        as expected by ``lightgbm.engine.train``.
210

211
212
213
214
215
216
217
218
219
220
        Parameters
        ----------
        func : callable
            Expects a callable with following signatures:
            ``func(y_true, y_pred)``,
            ``func(y_true, y_pred, weight)``
            or ``func(y_true, y_pred, weight, group)``
            and returns (eval_name, eval_result, is_higher_better) or
            list of (eval_name, eval_result, is_higher_better):

221
                y_true : numpy 1-D array of shape = [n_samples]
222
                    The target values.
223
                y_pred : numpy 1-D array of shape = [n_samples] or numpy 2-D array shape = [n_samples, n_classes] (for multi-class task)
224
                    The predicted values.
225
226
                    In case of custom ``objective``, predicted values are returned before any transformation,
                    e.g. they are raw margin instead of probability of positive class for binary task in this case.
227
                weight : numpy 1-D array of shape = [n_samples]
228
                    The weight of samples. Weights should be non-negative.
229
                group : numpy 1-D array
230
231
232
                    Group/query data.
                    Only used in the learning-to-rank task.
                    sum(group) = n_samples.
233
234
                    For example, if you have a 100-document dataset with ``group = [10, 20, 40, 10, 10, 10]``, that means that you have 6 groups,
                    where the first 10 records are in the first group, records 11-30 are in the second group, records 31-70 are in the third group, etc.
235
                eval_name : str
Andrew Ziem's avatar
Andrew Ziem committed
236
                    The name of evaluation function (without whitespace).
237
238
239
240
241
242
                eval_result : float
                    The eval result.
                is_higher_better : bool
                    Is eval result higher better, e.g. AUC is ``is_higher_better``.
        """
        self.func = func
243

244
245
246
247
248
    def __call__(
        self,
        preds: np.ndarray,
        dataset: Dataset
    ) -> Union[_LGBM_EvalFunctionResultType, List[_LGBM_EvalFunctionResultType]]:
249
        """Call passed function with appropriate arguments.
250

251
252
        Parameters
        ----------
253
        preds : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
254
255
256
257
258
259
            The predicted values.
        dataset : Dataset
            The training dataset.

        Returns
        -------
260
        eval_name : str
Andrew Ziem's avatar
Andrew Ziem committed
261
            The name of evaluation function (without whitespace).
262
263
264
265
266
        eval_result : float
            The eval result.
        is_higher_better : bool
            Is eval result higher better, e.g. AUC is ``is_higher_better``.
        """
267
        labels = _get_label_from_constructed_dataset(dataset)
268
        argc = len(signature(self.func).parameters)
269
        if argc == 2:
270
            return self.func(labels, preds)  # type: ignore[call-arg]
271
272
273
274
275
276
277
278
279
280

        weight = _get_weight_from_constructed_dataset(dataset)
        if argc == 3:
            return self.func(labels, preds, weight)  # type: ignore[call-arg]

        if argc == 4:
            group = _get_group_from_constructed_dataset(dataset)
            return self.func(labels, preds, weight, group)  # type: ignore[call-arg]

        raise TypeError(f"Self-defined eval function should have 2, 3 or 4 arguments, got {argc}")
281

wxchan's avatar
wxchan committed
282

283
284
285
286
287
288
289
290
291
292
293
294
295
296
# documentation templates for LGBMModel methods are shared between the classes in
# this module and those in the ``dask`` module

_lgbmmodel_doc_fit = (
    """
    Build a gradient boosting model from the training set (X, y).

    Parameters
    ----------
    X : {X_shape}
        Input feature matrix.
    y : {y_shape}
        The target values (class labels in classification, real numbers in regression).
    sample_weight : {sample_weight_shape}
297
        Weights of training data. Weights should be non-negative.
298
    init_score : {init_score_shape}
299
300
301
302
303
304
305
306
307
        Init score of training data.
    group : {group_shape}
        Group/query data.
        Only used in the learning-to-rank task.
        sum(group) = n_samples.
        For example, if you have a 100-document dataset with ``group = [10, 20, 40, 10, 10, 10]``, that means that you have 6 groups,
        where the first 10 records are in the first group, records 11-30 are in the second group, records 31-70 are in the third group, etc.
    eval_set : list or None, optional (default=None)
        A list of (X, y) tuple pairs to use as validation sets.
308
    eval_names : list of str, or None, optional (default=None)
309
        Names of eval_set.
310
    eval_sample_weight : {eval_sample_weight_shape}
311
        Weights of eval data. Weights should be non-negative.
312
313
    eval_class_weight : list or None, optional (default=None)
        Class weights of eval data.
314
    eval_init_score : {eval_init_score_shape}
315
        Init score of eval data.
316
    eval_group : {eval_group_shape}
317
        Group data of eval data.
318
319
    eval_metric : str, callable, list or None, optional (default=None)
        If str, it should be a built-in evaluation metric to use.
320
321
322
323
        If callable, it should be a custom evaluation metric, see note below for more details.
        If list, it can be a list of built-in metrics, a list of custom evaluation metrics, or a mix of both.
        In either case, the ``metric`` from the model parameters will be evaluated and used as well.
        Default: 'l2' for LGBMRegressor, 'logloss' for LGBMClassifier, 'ndcg' for LGBMRanker.
324
    feature_name : list of str, or 'auto', optional (default='auto')
325
326
        Feature names.
        If 'auto' and data is pandas DataFrame, data columns names are used.
327
    categorical_feature : list of str or int, or 'auto', optional (default='auto')
328
329
        Categorical features.
        If list of int, interpreted as indices.
330
        If list of str, interpreted as feature names (need to specify ``feature_name`` as well).
331
        If 'auto' and data is pandas DataFrame, pandas unordered categorical columns are used.
332
        All values in categorical features will be cast to int32 and thus should be less than int32 max value (2147483647).
333
334
335
        Large values could be memory consuming. Consider using consecutive integers starting from zero.
        All negative values in categorical features will be treated as missing values.
        The output cannot be monotonically constrained with respect to a categorical feature.
336
        Floating point numbers in categorical features will be rounded towards 0.
337
    callbacks : list of callable, or None, optional (default=None)
338
339
        List of callback functions that are applied at each iteration.
        See Callbacks in Python API for more information.
340
    init_model : str, pathlib.Path, Booster, LGBMModel or None, optional (default=None)
341
342
343
344
        Filename of LightGBM model, Booster instance or LGBMModel instance used for continue training.

    Returns
    -------
345
    self : LGBMModel
346
347
348
349
350
351
352
353
354
355
356
357
358
        Returns self.
    """
)

_lgbmmodel_doc_custom_eval_note = """
    Note
    ----
    Custom eval function expects a callable with following signatures:
    ``func(y_true, y_pred)``, ``func(y_true, y_pred, weight)`` or
    ``func(y_true, y_pred, weight, group)``
    and returns (eval_name, eval_result, is_higher_better) or
    list of (eval_name, eval_result, is_higher_better):

359
        y_true : numpy 1-D array of shape = [n_samples]
360
            The target values.
361
        y_pred : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
362
            The predicted values.
363
364
            In case of custom ``objective``, predicted values are returned before any transformation,
            e.g. they are raw margin instead of probability of positive class for binary task in this case.
365
        weight : numpy 1-D array of shape = [n_samples]
366
            The weight of samples. Weights should be non-negative.
367
        group : numpy 1-D array
368
369
370
371
372
            Group/query data.
            Only used in the learning-to-rank task.
            sum(group) = n_samples.
            For example, if you have a 100-document dataset with ``group = [10, 20, 40, 10, 10, 10]``, that means that you have 6 groups,
            where the first 10 records are in the first group, records 11-30 are in the second group, records 31-70 are in the third group, etc.
373
        eval_name : str
Andrew Ziem's avatar
Andrew Ziem committed
374
            The name of evaluation function (without whitespace).
375
376
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
        eval_result : float
            The eval result.
        is_higher_better : bool
            Is eval result higher better, e.g. AUC is ``is_higher_better``.
"""

_lgbmmodel_doc_predict = (
    """
    {description}

    Parameters
    ----------
    X : {X_shape}
        Input features matrix.
    raw_score : bool, optional (default=False)
        Whether to predict raw scores.
    start_iteration : int, optional (default=0)
        Start index of the iteration to predict.
        If <= 0, starts from the first iteration.
    num_iteration : int or None, optional (default=None)
        Total number of iterations used in the prediction.
        If None, if the best iteration exists and start_iteration <= 0, the best iteration is used;
        otherwise, all iterations from ``start_iteration`` are used (no limits).
        If <= 0, all iterations from ``start_iteration`` are used (no limits).
    pred_leaf : bool, optional (default=False)
        Whether to predict leaf index.
    pred_contrib : bool, optional (default=False)
        Whether to predict feature contributions.

        .. note::

            If you want to get more explanations for your model's predictions using SHAP values,
            like SHAP interaction values,
            you can install the shap package (https://github.com/slundberg/shap).
            Note that unlike the shap package, with ``pred_contrib`` we return a matrix with an extra
            column, where the last column is the expected value.

412
413
414
    validate_features : bool, optional (default=False)
        If True, ensure that the features used to predict match the ones used to train.
        Used only if data is pandas DataFrame.
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
    **kwargs
        Other parameters for the prediction.

    Returns
    -------
    {output_name} : {predicted_result_shape}
        The predicted values.
    X_leaves : {X_leaves_shape}
        If ``pred_leaf=True``, the predicted leaf of every tree for each sample.
    X_SHAP_values : {X_SHAP_values_shape}
        If ``pred_contrib=True``, the feature contributions for each sample.
    """
)


430
431
class LGBMModel(_LGBMModelBase):
    """Implementation of the scikit-learn API for LightGBM."""
wxchan's avatar
wxchan committed
432

433
434
435
436
437
438
439
440
    def __init__(
        self,
        boosting_type: str = 'gbdt',
        num_leaves: int = 31,
        max_depth: int = -1,
        learning_rate: float = 0.1,
        n_estimators: int = 100,
        subsample_for_bin: int = 200000,
441
        objective: Optional[Union[str, _LGBM_ScikitCustomObjectiveFunction]] = None,
442
443
444
445
446
447
448
449
450
451
        class_weight: Optional[Union[Dict, str]] = None,
        min_split_gain: float = 0.,
        min_child_weight: float = 1e-3,
        min_child_samples: int = 20,
        subsample: float = 1.,
        subsample_freq: int = 0,
        colsample_bytree: float = 1.,
        reg_alpha: float = 0.,
        reg_lambda: float = 0.,
        random_state: Optional[Union[int, np.random.RandomState]] = None,
452
        n_jobs: Optional[int] = None,
453
454
455
        importance_type: str = 'split',
        **kwargs
    ):
456
        r"""Construct a gradient boosting model.
wxchan's avatar
wxchan committed
457
458
459

        Parameters
        ----------
460
        boosting_type : str, optional (default='gbdt')
461
462
463
464
            'gbdt', traditional Gradient Boosting Decision Tree.
            'dart', Dropouts meet Multiple Additive Regression Trees.
            'rf', Random Forest.
        num_leaves : int, optional (default=31)
wxchan's avatar
wxchan committed
465
            Maximum tree leaves for base learners.
466
        max_depth : int, optional (default=-1)
467
            Maximum tree depth for base learners, <=0 means no limit.
468
        learning_rate : float, optional (default=0.1)
469
            Boosting learning rate.
470
471
472
            You can use ``callbacks`` parameter of ``fit`` method to shrink/adapt learning rate
            in training using ``reset_parameter`` callback.
            Note, that this will ignore the ``learning_rate`` argument in training.
473
        n_estimators : int, optional (default=100)
wxchan's avatar
wxchan committed
474
            Number of boosted trees to fit.
475
        subsample_for_bin : int, optional (default=200000)
wxchan's avatar
wxchan committed
476
            Number of samples for constructing bins.
477
        objective : str, callable or None, optional (default=None)
wxchan's avatar
wxchan committed
478
479
            Specify the learning task and the corresponding learning objective or
            a custom objective function to be used (see note below).
480
            Default: 'regression' for LGBMRegressor, 'binary' or 'multiclass' for LGBMClassifier, 'lambdarank' for LGBMRanker.
481
482
483
484
        class_weight : dict, 'balanced' or None, optional (default=None)
            Weights associated with classes in the form ``{class_label: weight}``.
            Use this parameter only for multi-class classification task;
            for binary classification task you may use ``is_unbalance`` or ``scale_pos_weight`` parameters.
485
486
487
            Note, that the usage of all these parameters will result in poor estimates of the individual class probabilities.
            You may want to consider performing probability calibration
            (https://scikit-learn.org/stable/modules/calibration.html) of your model.
488
489
490
            The 'balanced' mode uses the values of y to automatically adjust weights
            inversely proportional to class frequencies in the input data as ``n_samples / (n_classes * np.bincount(y))``.
            If None, all classes are supposed to have weight one.
491
            Note, that these weights will be multiplied with ``sample_weight`` (passed through the ``fit`` method)
492
            if ``sample_weight`` is specified.
493
        min_split_gain : float, optional (default=0.)
wxchan's avatar
wxchan committed
494
            Minimum loss reduction required to make a further partition on a leaf node of the tree.
495
        min_child_weight : float, optional (default=1e-3)
496
            Minimum sum of instance weight (Hessian) needed in a child (leaf).
497
        min_child_samples : int, optional (default=20)
498
            Minimum number of data needed in a child (leaf).
499
        subsample : float, optional (default=1.)
wxchan's avatar
wxchan committed
500
            Subsample ratio of the training instance.
501
        subsample_freq : int, optional (default=0)
Andrew Ziem's avatar
Andrew Ziem committed
502
            Frequency of subsample, <=0 means no enable.
503
        colsample_bytree : float, optional (default=1.)
wxchan's avatar
wxchan committed
504
            Subsample ratio of columns when constructing each tree.
505
        reg_alpha : float, optional (default=0.)
506
            L1 regularization term on weights.
507
        reg_lambda : float, optional (default=0.)
508
            L2 regularization term on weights.
509
        random_state : int, RandomState object or None, optional (default=None)
wxchan's avatar
wxchan committed
510
            Random number seed.
511
512
513
            If int, this number is used to seed the C++ code.
            If RandomState object (numpy), a random integer is picked based on its state to seed the C++ code.
            If None, default seeds in C++ code are used.
514
515
516
517
518
519
520
521
522
523
524
525
        n_jobs : int or None, optional (default=None)
            Number of parallel threads to use for training (can be changed at prediction time by
            passing it as an extra keyword argument).

            For better performance, it is recommended to set this to the number of physical cores
            in the CPU.

            Negative integers are interpreted as following joblib's formula (n_cpus + 1 + n_jobs), just like
            scikit-learn (so e.g. -1 means using all threads). A value of zero corresponds the default number of
            threads configured for OpenMP in the system. A value of ``None`` (the default) corresponds
            to using the number of physical cores in the system (its correct detection requires
            either the ``joblib`` or the ``psutil`` util libraries to be installed).
526
527
528

            .. versionchanged:: 4.0.0

529
        importance_type : str, optional (default='split')
530
            The type of feature importance to be filled into ``feature_importances_``.
531
532
533
534
            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.
        **kwargs
            Other parameters for the model.
wxchan's avatar
wxchan committed
535
            Check http://lightgbm.readthedocs.io/en/latest/Parameters.html for more parameters.
536

Nikita Titov's avatar
Nikita Titov committed
537
538
539
            .. warning::

                \*\*kwargs is not supported in sklearn, it may cause unexpected issues.
wxchan's avatar
wxchan committed
540
541
542

        Note
        ----
543
544
        A custom objective function can be provided for the ``objective`` parameter.
        In this case, it should have the signature
545
546
547
        ``objective(y_true, y_pred) -> grad, hess``,
        ``objective(y_true, y_pred, weight) -> grad, hess``
        or ``objective(y_true, y_pred, weight, group) -> grad, hess``:
wxchan's avatar
wxchan committed
548

549
            y_true : numpy 1-D array of shape = [n_samples]
550
                The target values.
551
            y_pred : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
552
                The predicted values.
553
554
                Predicted values are returned before any transformation,
                e.g. they are raw margin instead of probability of positive class for binary task.
555
556
            weight : numpy 1-D array of shape = [n_samples]
                The weight of samples. Weights should be non-negative.
557
            group : numpy 1-D array
558
559
560
                Group/query data.
                Only used in the learning-to-rank task.
                sum(group) = n_samples.
561
562
                For example, if you have a 100-document dataset with ``group = [10, 20, 40, 10, 10, 10]``, that means that you have 6 groups,
                where the first 10 records are in the first group, records 11-30 are in the second group, records 31-70 are in the third group, etc.
563
            grad : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
564
565
                The value of the first order derivative (gradient) of the loss
                with respect to the elements of y_pred for each sample point.
566
            hess : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
567
568
                The value of the second order derivative (Hessian) of the loss
                with respect to the elements of y_pred for each sample point.
wxchan's avatar
wxchan committed
569

570
        For multi-class task, y_pred is a numpy 2-D array of shape = [n_samples, n_classes],
571
        and grad and hess should be returned in the same format.
wxchan's avatar
wxchan committed
572
        """
wxchan's avatar
wxchan committed
573
        if not SKLEARN_INSTALLED:
574
575
            raise LightGBMError('scikit-learn is required for lightgbm.sklearn. '
                                'You must install scikit-learn and restart your session to use this module.')
wxchan's avatar
wxchan committed
576

577
        self.boosting_type = boosting_type
578
        self.objective = objective
wxchan's avatar
wxchan committed
579
580
581
582
        self.num_leaves = num_leaves
        self.max_depth = max_depth
        self.learning_rate = learning_rate
        self.n_estimators = n_estimators
wxchan's avatar
wxchan committed
583
        self.subsample_for_bin = subsample_for_bin
wxchan's avatar
wxchan committed
584
585
586
587
588
589
590
591
        self.min_split_gain = min_split_gain
        self.min_child_weight = min_child_weight
        self.min_child_samples = min_child_samples
        self.subsample = subsample
        self.subsample_freq = subsample_freq
        self.colsample_bytree = colsample_bytree
        self.reg_alpha = reg_alpha
        self.reg_lambda = reg_lambda
592
593
        self.random_state = random_state
        self.n_jobs = n_jobs
594
        self.importance_type = importance_type
595
        self._Booster: Optional[Booster] = None
596
597
        self._evals_result: _EvalResultDict = {}
        self._best_score: _LGBM_BoosterBestScoreType = {}
598
        self._best_iteration: int = -1
599
        self._other_params: Dict[str, Any] = {}
600
        self._objective = objective
601
        self.class_weight = class_weight
602
603
        self._class_weight: Optional[Union[Dict, str]] = None
        self._class_map: Optional[Dict[int, int]] = None
604
605
        self._n_features: int = -1
        self._n_features_in: int = -1
606
        self._classes: Optional[np.ndarray] = None
607
        self._n_classes: int = -1
608
        self.set_params(**kwargs)
wxchan's avatar
wxchan committed
609

610
    def _more_tags(self) -> Dict[str, Any]:
611
612
613
614
615
616
617
618
619
620
        return {
            'allow_nan': True,
            'X_types': ['2darray', 'sparse', '1dlabels'],
            '_xfail_checks': {
                'check_no_attributes_set_in_init':
                'scikit-learn incorrectly asserts that private attributes '
                'cannot be set in __init__: '
                '(see https://github.com/microsoft/LightGBM/issues/2628)'
            }
        }
Nikita Titov's avatar
Nikita Titov committed
621

622
623
624
    def __sklearn_is_fitted__(self) -> bool:
        return getattr(self, "fitted_", False)

625
    def get_params(self, deep: bool = True) -> Dict[str, Any]:
626
627
628
629
630
631
632
633
634
635
636
637
638
        """Get parameters for this estimator.

        Parameters
        ----------
        deep : bool, optional (default=True)
            If True, will return the parameters for this estimator and
            contained subobjects that are estimators.

        Returns
        -------
        params : dict
            Parameter names mapped to their values.
        """
639
        params = super().get_params(deep=deep)
640
        params.update(self._other_params)
wxchan's avatar
wxchan committed
641
642
        return params

643
    def set_params(self, **params: Any) -> "LGBMModel":
644
645
646
647
648
649
650
651
652
653
654
655
        """Set the parameters of this estimator.

        Parameters
        ----------
        **params
            Parameter names with their new values.

        Returns
        -------
        self : object
            Returns self.
        """
wxchan's avatar
wxchan committed
656
657
        for key, value in params.items():
            setattr(self, key, value)
658
659
            if hasattr(self, f"_{key}"):
                setattr(self, f"_{key}", value)
660
            self._other_params[key] = value
wxchan's avatar
wxchan committed
661
        return self
wxchan's avatar
wxchan committed
662

663
664
665
666
667
668
669
670
671
672
673
674
675
676
    def _process_params(self, stage: str) -> Dict[str, Any]:
        """Process the parameters of this estimator based on its type, parameter aliases, etc.

        Parameters
        ----------
        stage : str
            Name of the stage (can be ``fit`` or ``predict``) this method is called from.

        Returns
        -------
        processed_params : dict
            Processed parameter names mapped to their values.
        """
        assert stage in {"fit", "predict"}
677
678
679
680
681
        params = self.get_params()

        params.pop('objective', None)
        for alias in _ConfigAliases.get('objective'):
            if alias in params:
682
                obj = params.pop(alias)
683
                _log_warning(f"Found '{alias}' in params. Will use it instead of 'objective' argument")
684
685
686
687
688
689
690
691
692
693
694
695
696
                if stage == "fit":
                    self._objective = obj
        if stage == "fit":
            if self._objective is None:
                if isinstance(self, LGBMRegressor):
                    self._objective = "regression"
                elif isinstance(self, LGBMClassifier):
                    if self._n_classes > 2:
                        self._objective = "multiclass"
                    else:
                        self._objective = "binary"
                elif isinstance(self, LGBMRanker):
                    self._objective = "lambdarank"
697
                else:
698
                    raise ValueError("Unknown LGBMModel type.")
699
        if callable(self._objective):
700
            if stage == "fit":
701
702
703
                params['objective'] = _ObjectiveFunctionWrapper(self._objective)
            else:
                params['objective'] = 'None'
704
        else:
705
            params['objective'] = self._objective
706

707
        params.pop('importance_type', None)
wxchan's avatar
wxchan committed
708
        params.pop('n_estimators', None)
709
        params.pop('class_weight', None)
710

711
712
        if isinstance(params['random_state'], np.random.RandomState):
            params['random_state'] = params['random_state'].randint(np.iinfo(np.int32).max)
713
        if self._n_classes > 2:
714
715
            for alias in _ConfigAliases.get('num_class'):
                params.pop(alias, None)
716
717
            params['num_class'] = self._n_classes
        if hasattr(self, '_eval_at'):
718
            eval_at = self._eval_at
719
            for alias in _ConfigAliases.get('eval_at'):
720
721
722
723
                if alias in params:
                    _log_warning(f"Found '{alias}' in params. Will use it instead of 'eval_at' argument")
                    eval_at = params.pop(alias)
            params['eval_at'] = eval_at
wxchan's avatar
wxchan committed
724

725
        # register default metric for consistency with callable eval_metric case
726
        original_metric = self._objective if isinstance(self._objective, str) else None
727
728
729
730
731
732
733
734
735
736
        if original_metric is None:
            # try to deduce from class instance
            if isinstance(self, LGBMRegressor):
                original_metric = "l2"
            elif isinstance(self, LGBMClassifier):
                original_metric = "multi_logloss" if self._n_classes > 2 else "binary_logloss"
            elif isinstance(self, LGBMRanker):
                original_metric = "ndcg"

        # overwrite default metric by explicitly set metric
737
        params = _choose_param_value("metric", params, original_metric)
738

739
740
741
742
743
744
        # use joblib conventions for negative n_jobs, just like scikit-learn
        # at predict time, this is handled later due to the order of parameter updates
        if stage == "fit":
            params = _choose_param_value("num_threads", params, self.n_jobs)
            params["num_threads"] = self._process_n_jobs(params["num_threads"])

745
746
        return params

747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
    def _process_n_jobs(self, n_jobs: Optional[int]) -> int:
        """Convert special values of n_jobs to their actual values according to the formulas that apply.

        Parameters
        ----------
        n_jobs : int or None
            The original value of n_jobs, potentially having special values such as 'None' or
            negative integers.

        Returns
        -------
        n_jobs : int
            The value of n_jobs with special values converted to actual number of threads.
        """
        if n_jobs is None:
            n_jobs = _LGBMCpuCount(only_physical_cores=True)
        elif n_jobs < 0:
            n_jobs = max(_LGBMCpuCount(only_physical_cores=False) + 1 + n_jobs, 1)
        return n_jobs

767
768
    def fit(
        self,
769
770
        X: _LGBM_ScikitMatrixLike,
        y: _LGBM_LabelType,
771
772
        sample_weight: Optional[_LGBM_WeightType] = None,
        init_score: Optional[_LGBM_InitScoreType] = None,
773
        group: Optional[_LGBM_GroupType] = None,
774
        eval_set: Optional[List[_LGBM_ScikitValidSet]] = None,
775
        eval_names: Optional[List[str]] = None,
776
777
778
779
        eval_sample_weight: Optional[List[_LGBM_WeightType]] = None,
        eval_class_weight: Optional[List[float]] = None,
        eval_init_score: Optional[List[_LGBM_InitScoreType]] = None,
        eval_group: Optional[List[_LGBM_GroupType]] = None,
780
        eval_metric: Optional[_LGBM_ScikitEvalMetricType] = None,
781
782
        feature_name: _LGBM_FeatureNameConfiguration = 'auto',
        categorical_feature: _LGBM_CategoricalFeatureConfiguration = 'auto',
783
        callbacks: Optional[List[Callable]] = None,
784
        init_model: Optional[Union[str, Path, Booster, "LGBMModel"]] = None
785
    ) -> "LGBMModel":
786
787
788
789
790
        """Docstring is set after definition, using a template."""
        params = self._process_params(stage="fit")

        # Do not modify original args in fit function
        # Refer to https://github.com/microsoft/LightGBM/pull/2619
791
792
793
794
795
796
797
        eval_metric_list: List[Union[str, _LGBM_ScikitCustomEvalFunction]]
        if eval_metric is None:
            eval_metric_list = []
        elif isinstance(eval_metric, list):
            eval_metric_list = copy.deepcopy(eval_metric)
        else:
            eval_metric_list = [copy.deepcopy(eval_metric)]
798
799
800
801
802

        # Separate built-in from callable evaluation metrics
        eval_metrics_callable = [_EvalFunctionWrapper(f) for f in eval_metric_list if callable(f)]
        eval_metrics_builtin = [m for m in eval_metric_list if isinstance(m, str)]

803
        # concatenate metric from params (or default if not provided in params) and eval_metric
804
805
        params['metric'] = [params['metric']] if isinstance(params['metric'], (str, type(None))) else params['metric']
        params['metric'] = [e for e in eval_metrics_builtin if e not in params['metric']] + params['metric']
806
        params['metric'] = [metric for metric in params['metric'] if metric is not None]
wxchan's avatar
wxchan committed
807

808
        if not isinstance(X, (pd_DataFrame, dt_DataTable)):
809
            _X, _y = _LGBMCheckXY(X, y, accept_sparse=True, force_all_finite=False, ensure_min_samples=2)
810
811
            if sample_weight is not None:
                sample_weight = _LGBMCheckSampleWeight(sample_weight, _X)
812
813
        else:
            _X, _y = X, y
814

815
816
817
818
        if self._class_weight is None:
            self._class_weight = self.class_weight
        if self._class_weight is not None:
            class_sample_weight = _LGBMComputeSampleWeight(self._class_weight, y)
819
820
821
822
            if sample_weight is None or len(sample_weight) == 0:
                sample_weight = class_sample_weight
            else:
                sample_weight = np.multiply(sample_weight, class_sample_weight)
823

824
        self._n_features = _X.shape[1]
825
826
        # copy for consistency
        self._n_features_in = self._n_features
827

828
829
830
        train_set = Dataset(data=_X, label=_y, weight=sample_weight, group=group,
                            init_score=init_score, categorical_feature=categorical_feature,
                            params=params)
Guolin Ke's avatar
Guolin Ke committed
831

832
        valid_sets: List[Dataset] = []
Guolin Ke's avatar
Guolin Ke committed
833
        if eval_set is not None:
834

835
            def _get_meta_data(collection, name, i):
836
837
838
839
840
841
842
                if collection is None:
                    return None
                elif isinstance(collection, list):
                    return collection[i] if len(collection) > i else None
                elif isinstance(collection, dict):
                    return collection.get(i, None)
                else:
843
                    raise TypeError(f"{name} should be dict or list")
844

Guolin Ke's avatar
Guolin Ke committed
845
846
847
            if isinstance(eval_set, tuple):
                eval_set = [eval_set]
            for i, valid_data in enumerate(eval_set):
848
                # reduce cost for prediction training data
Guolin Ke's avatar
Guolin Ke committed
849
850
851
                if valid_data[0] is X and valid_data[1] is y:
                    valid_set = train_set
                else:
852
853
854
855
856
857
                    valid_weight = _get_meta_data(eval_sample_weight, 'eval_sample_weight', i)
                    valid_class_weight = _get_meta_data(eval_class_weight, 'eval_class_weight', i)
                    if valid_class_weight is not None:
                        if isinstance(valid_class_weight, dict) and self._class_map is not None:
                            valid_class_weight = {self._class_map[k]: v for k, v in valid_class_weight.items()}
                        valid_class_sample_weight = _LGBMComputeSampleWeight(valid_class_weight, valid_data[1])
858
859
860
861
                        if valid_weight is None or len(valid_weight) == 0:
                            valid_weight = valid_class_sample_weight
                        else:
                            valid_weight = np.multiply(valid_weight, valid_class_sample_weight)
862
863
                    valid_init_score = _get_meta_data(eval_init_score, 'eval_init_score', i)
                    valid_group = _get_meta_data(eval_group, 'eval_group', i)
864
865
866
867
                    valid_set = Dataset(data=valid_data[0], label=valid_data[1], weight=valid_weight,
                                        group=valid_group, init_score=valid_init_score,
                                        categorical_feature='auto', params=params)

Guolin Ke's avatar
Guolin Ke committed
868
869
                valid_sets.append(valid_set)

870
871
872
        if isinstance(init_model, LGBMModel):
            init_model = init_model.booster_

873
874
875
        if callbacks is None:
            callbacks = []
        else:
876
            callbacks = copy.copy(callbacks)  # don't use deepcopy here to allow non-serializable objects
877

878
        evals_result: _EvalResultDict = {}
879
880
881
882
883
884
885
886
        callbacks.append(record_evaluation(evals_result))

        self._Booster = train(
            params=params,
            train_set=train_set,
            num_boost_round=self.n_estimators,
            valid_sets=valid_sets,
            valid_names=eval_names,
887
            feval=eval_metrics_callable,  # type: ignore[arg-type]
888
889
890
891
            init_model=init_model,
            feature_name=feature_name,
            callbacks=callbacks
        )
wxchan's avatar
wxchan committed
892

893
        self._evals_result = evals_result
894
        self._best_iteration = self._Booster.best_iteration
895
        self._best_score = self._Booster.best_score
wxchan's avatar
wxchan committed
896

897
898
        self.fitted_ = True

wxchan's avatar
wxchan committed
899
        # free dataset
900
        self._Booster.free_dataset()
wxchan's avatar
wxchan committed
901
        del train_set, valid_sets
wxchan's avatar
wxchan committed
902
903
        return self

904
    fit.__doc__ = _lgbmmodel_doc_fit.format(
905
906
        X_shape="numpy array, pandas DataFrame, H2O DataTable's Frame , scipy.sparse, list of lists of int or float of shape = [n_samples, n_features]",
        y_shape="numpy array, pandas DataFrame, pandas Series, list of int or float of shape = [n_samples]",
907
908
        sample_weight_shape="numpy array, pandas Series, list of int or float of shape = [n_samples] or None, optional (default=None)",
        init_score_shape="numpy array, pandas DataFrame, pandas Series, list of int or float of shape = [n_samples] or shape = [n_samples * n_classes] (for multi-class task) or shape = [n_samples, n_classes] (for multi-class task) or None, optional (default=None)",
909
        group_shape="numpy array, pandas Series, list of int or float, or None, optional (default=None)",
910
911
912
        eval_sample_weight_shape="list of array (same types as ``sample_weight`` supports), or None, optional (default=None)",
        eval_init_score_shape="list of array (same types as ``init_score`` supports), or None, optional (default=None)",
        eval_group_shape="list of array (same types as ``group`` supports), or None, optional (default=None)"
913
914
    ) + "\n\n" + _lgbmmodel_doc_custom_eval_note

915
916
    def predict(
        self,
917
        X: _LGBM_ScikitMatrixLike,
918
919
920
921
922
923
924
925
        raw_score: bool = False,
        start_iteration: int = 0,
        num_iteration: Optional[int] = None,
        pred_leaf: bool = False,
        pred_contrib: bool = False,
        validate_features: bool = False,
        **kwargs: Any
    ):
926
        """Docstring is set after definition, using a template."""
927
        if not self.__sklearn_is_fitted__():
928
            raise LGBMNotFittedError("Estimator not fitted, call fit before exploiting the model.")
929
        if not isinstance(X, (pd_DataFrame, dt_DataTable)):
930
            X = _LGBMCheckArray(X, accept_sparse=True, force_all_finite=False)
931
932
933
        n_features = X.shape[1]
        if self._n_features != n_features:
            raise ValueError("Number of features of the model must "
934
935
                             f"match the input. Model n_features_ is {self._n_features} and "
                             f"input n_features is {n_features}")
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
        # retrive original params that possibly can be used in both training and prediction
        # and then overwrite them (considering aliases) with params that were passed directly in prediction
        predict_params = self._process_params(stage="predict")
        for alias in _ConfigAliases.get_by_alias(
            "data",
            "X",
            "raw_score",
            "start_iteration",
            "num_iteration",
            "pred_leaf",
            "pred_contrib",
            *kwargs.keys()
        ):
            predict_params.pop(alias, None)
        predict_params.update(kwargs)
951
952
953

        # number of threads can have values with special meaning which is only applied
        # in the scikit-learn interface, these should not reach the c++ side as-is
954
955
        predict_params = _choose_param_value("num_threads", predict_params, self.n_jobs)
        predict_params["num_threads"] = self._process_n_jobs(predict_params["num_threads"])
956

957
958
959
960
961
        return self._Booster.predict(  # type: ignore[union-attr]
            X, raw_score=raw_score, start_iteration=start_iteration, num_iteration=num_iteration,
            pred_leaf=pred_leaf, pred_contrib=pred_contrib, validate_features=validate_features,
            **predict_params
        )
wxchan's avatar
wxchan committed
962

963
964
    predict.__doc__ = _lgbmmodel_doc_predict.format(
        description="Return the predicted value for each sample.",
965
        X_shape="numpy array, pandas DataFrame, H2O DataTable's Frame , scipy.sparse, list of lists of int or float of shape = [n_samples, n_features]",
966
967
968
969
970
971
        output_name="predicted_result",
        predicted_result_shape="array-like of shape = [n_samples] or shape = [n_samples, n_classes]",
        X_leaves_shape="array-like of shape = [n_samples, n_trees] or shape = [n_samples, n_trees * n_classes]",
        X_SHAP_values_shape="array-like of shape = [n_samples, n_features + 1] or shape = [n_samples, (n_features + 1) * n_classes] or list with n_classes length of such objects"
    )

972
    @property
973
    def n_features_(self) -> int:
974
        """:obj:`int`: The number of features of fitted model."""
975
        if not self.__sklearn_is_fitted__():
976
977
978
            raise LGBMNotFittedError('No n_features found. Need to call fit beforehand.')
        return self._n_features

979
    @property
980
    def n_features_in_(self) -> int:
981
        """:obj:`int`: The number of features of fitted model."""
982
        if not self.__sklearn_is_fitted__():
983
984
985
            raise LGBMNotFittedError('No n_features_in found. Need to call fit beforehand.')
        return self._n_features_in

986
    @property
987
    def best_score_(self) -> _LGBM_BoosterBestScoreType:
988
        """:obj:`dict`: The best score of fitted model."""
989
        if not self.__sklearn_is_fitted__():
990
991
992
993
            raise LGBMNotFittedError('No best_score found. Need to call fit beforehand.')
        return self._best_score

    @property
994
    def best_iteration_(self) -> int:
995
        """:obj:`int`: The best iteration of fitted model if ``early_stopping()`` callback has been specified."""
996
        if not self.__sklearn_is_fitted__():
997
            raise LGBMNotFittedError('No best_iteration found. Need to call fit with early_stopping callback beforehand.')
998
999
1000
        return self._best_iteration

    @property
1001
    def objective_(self) -> Union[str, _LGBM_ScikitCustomObjectiveFunction]:
1002
        """:obj:`str` or :obj:`callable`: The concrete objective used while fitting this model."""
1003
        if not self.__sklearn_is_fitted__():
1004
            raise LGBMNotFittedError('No objective found. Need to call fit beforehand.')
1005
        return self._objective  # type: ignore[return-value]
1006

1007
1008
1009
1010
1011
1012
    @property
    def n_estimators_(self) -> int:
        """:obj:`int`: True number of boosting iterations performed.

        This might be less than parameter ``n_estimators`` if early stopping was enabled or
        if boosting stopped early due to limits on complexity like ``min_gain_to_split``.
1013
1014
        
        .. versionadded:: 4.0.0
1015
1016
1017
        """
        if not self.__sklearn_is_fitted__():
            raise LGBMNotFittedError('No n_estimators found. Need to call fit beforehand.')
1018
        return self._Booster.current_iteration()  # type: ignore
1019
1020
1021
1022
1023
1024
1025

    @property
    def n_iter_(self) -> int:
        """:obj:`int`: True number of boosting iterations performed.

        This might be less than parameter ``n_estimators`` if early stopping was enabled or
        if boosting stopped early due to limits on complexity like ``min_gain_to_split``.
1026
1027
        
        .. versionadded:: 4.0.0
1028
1029
1030
        """
        if not self.__sklearn_is_fitted__():
            raise LGBMNotFittedError('No n_iter found. Need to call fit beforehand.')
1031
        return self._Booster.current_iteration()  # type: ignore
1032

1033
    @property
1034
    def booster_(self) -> Booster:
1035
        """Booster: The underlying Booster of this model."""
1036
        if not self.__sklearn_is_fitted__():
1037
            raise LGBMNotFittedError('No booster found. Need to call fit beforehand.')
1038
        return self._Booster  # type: ignore[return-value]
wxchan's avatar
wxchan committed
1039

1040
    @property
1041
    def evals_result_(self) -> _EvalResultDict:
1042
        """:obj:`dict`: The evaluation results if validation sets have been specified."""
1043
        if not self.__sklearn_is_fitted__():
1044
1045
            raise LGBMNotFittedError('No results found. Need to call fit with eval_set beforehand.')
        return self._evals_result
1046
1047

    @property
1048
    def feature_importances_(self) -> np.ndarray:
1049
        """:obj:`array` of shape = [n_features]: The feature importances (the higher, the more important).
1050

Nikita Titov's avatar
Nikita Titov committed
1051
1052
1053
1054
        .. note::

            ``importance_type`` attribute is passed to the function
            to configure the type of importance values to be extracted.
1055
        """
1056
        if not self.__sklearn_is_fitted__():
1057
            raise LGBMNotFittedError('No feature_importances found. Need to call fit beforehand.')
1058
        return self._Booster.feature_importance(importance_type=self.importance_type)  # type: ignore[union-attr]
wxchan's avatar
wxchan committed
1059

1060
    @property
1061
1062
    def feature_name_(self) -> List[str]:
        """:obj:`list` of shape = [n_features]: The names of features."""
1063
        if not self.__sklearn_is_fitted__():
1064
            raise LGBMNotFittedError('No feature_name found. Need to call fit beforehand.')
1065
        return self._Booster.feature_name()  # type: ignore[union-attr]
1066

wxchan's avatar
wxchan committed
1067

1068
class LGBMRegressor(_LGBMRegressorBase, LGBMModel):
1069
    """LightGBM regressor."""
wxchan's avatar
wxchan committed
1070

1071
    def fit(  # type: ignore[override]
1072
        self,
1073
1074
        X: _LGBM_ScikitMatrixLike,
        y: _LGBM_LabelType,
1075
1076
1077
        sample_weight: Optional[_LGBM_WeightType] = None,
        init_score: Optional[_LGBM_InitScoreType] = None,
        eval_set: Optional[List[_LGBM_ScikitValidSet]] = None,
1078
        eval_names: Optional[List[str]] = None,
1079
1080
        eval_sample_weight: Optional[List[_LGBM_WeightType]] = None,
        eval_init_score: Optional[List[_LGBM_InitScoreType]] = None,
1081
        eval_metric: Optional[_LGBM_ScikitEvalMetricType] = None,
1082
1083
        feature_name: _LGBM_FeatureNameConfiguration = 'auto',
        categorical_feature: _LGBM_CategoricalFeatureConfiguration = 'auto',
1084
        callbacks: Optional[List[Callable]] = None,
1085
        init_model: Optional[Union[str, Path, Booster, LGBMModel]] = None
1086
    ) -> "LGBMRegressor":
1087
        """Docstring is inherited from the LGBMModel."""
1088
1089
1090
1091
1092
1093
1094
1095
1096
1097
1098
1099
1100
1101
1102
        super().fit(
            X,
            y,
            sample_weight=sample_weight,
            init_score=init_score,
            eval_set=eval_set,
            eval_names=eval_names,
            eval_sample_weight=eval_sample_weight,
            eval_init_score=eval_init_score,
            eval_metric=eval_metric,
            feature_name=feature_name,
            categorical_feature=categorical_feature,
            callbacks=callbacks,
            init_model=init_model
        )
Guolin Ke's avatar
Guolin Ke committed
1103
1104
        return self

1105
    _base_doc = LGBMModel.fit.__doc__.replace("self : LGBMModel", "self : LGBMRegressor")  # type: ignore
1106
1107
    _base_doc = (_base_doc[:_base_doc.find('group :')]  # type: ignore
                 + _base_doc[_base_doc.find('eval_set :'):])  # type: ignore
1108
1109
1110
1111
    _base_doc = (_base_doc[:_base_doc.find('eval_class_weight :')]
                 + _base_doc[_base_doc.find('eval_init_score :'):])
    fit.__doc__ = (_base_doc[:_base_doc.find('eval_group :')]
                   + _base_doc[_base_doc.find('eval_metric :'):])
wxchan's avatar
wxchan committed
1112

1113

1114
class LGBMClassifier(_LGBMClassifierBase, LGBMModel):
1115
    """LightGBM classifier."""
wxchan's avatar
wxchan committed
1116

1117
    def fit(  # type: ignore[override]
1118
        self,
1119
1120
        X: _LGBM_ScikitMatrixLike,
        y: _LGBM_LabelType,
1121
1122
1123
        sample_weight: Optional[_LGBM_WeightType] = None,
        init_score: Optional[_LGBM_InitScoreType] = None,
        eval_set: Optional[List[_LGBM_ScikitValidSet]] = None,
1124
        eval_names: Optional[List[str]] = None,
1125
1126
1127
        eval_sample_weight: Optional[List[_LGBM_WeightType]] = None,
        eval_class_weight: Optional[List[float]] = None,
        eval_init_score: Optional[List[_LGBM_InitScoreType]] = None,
1128
        eval_metric: Optional[_LGBM_ScikitEvalMetricType] = None,
1129
1130
        feature_name: _LGBM_FeatureNameConfiguration = 'auto',
        categorical_feature: _LGBM_CategoricalFeatureConfiguration = 'auto',
1131
        callbacks: Optional[List[Callable]] = None,
1132
        init_model: Optional[Union[str, Path, Booster, LGBMModel]] = None
1133
    ) -> "LGBMClassifier":
1134
        """Docstring is inherited from the LGBMModel."""
1135
        _LGBMAssertAllFinite(y)
1136
1137
        _LGBMCheckClassificationTargets(y)
        self._le = _LGBMLabelEncoder().fit(y)
1138
        _y = self._le.transform(y)
1139
        self._class_map = dict(zip(self._le.classes_, self._le.transform(self._le.classes_)))
1140
1141
        if isinstance(self.class_weight, dict):
            self._class_weight = {self._class_map[k]: v for k, v in self.class_weight.items()}
1142

1143
        self._classes = self._le.classes_
1144
        self._n_classes = len(self._classes)  # type: ignore[arg-type]
1145
1146
        if self.objective is None:
            self._objective = None
1147

1148
1149
        # adjust eval metrics to match whether binary or multiclass
        # classification is being performed
1150
        if not callable(eval_metric):
1151
1152
1153
1154
1155
1156
            if isinstance(eval_metric, list):
                eval_metric_list = eval_metric
            elif isinstance(eval_metric, str):
                eval_metric_list = [eval_metric]
            else:
                eval_metric_list = []
1157
            if self._n_classes > 2:
1158
                for index, metric in enumerate(eval_metric_list):
1159
                    if metric in {'logloss', 'binary_logloss'}:
1160
                        eval_metric_list[index] = "multi_logloss"
1161
                    elif metric in {'error', 'binary_error'}:
1162
                        eval_metric_list[index] = "multi_error"
1163
            else:
1164
                for index, metric in enumerate(eval_metric_list):
1165
                    if metric in {'logloss', 'multi_logloss'}:
1166
                        eval_metric_list[index] = 'binary_logloss'
1167
                    elif metric in {'error', 'multi_error'}:
1168
1169
                        eval_metric_list[index] = 'binary_error'
            eval_metric = eval_metric_list
wxchan's avatar
wxchan committed
1170

1171
        # do not modify args, as it causes errors in model selection tools
1172
        valid_sets: Optional[List[_LGBM_ScikitValidSet]] = None
wxchan's avatar
wxchan committed
1173
        if eval_set is not None:
1174
1175
            if isinstance(eval_set, tuple):
                eval_set = [eval_set]
1176
1177
            valid_sets = []
            for valid_x, valid_y in eval_set:
1178
                if valid_x is X and valid_y is y:
1179
                    valid_sets.append((valid_x, _y))
1180
                else:
1181
                    valid_sets.append((valid_x, self._le.transform(valid_y)))
1182

1183
1184
1185
1186
1187
1188
1189
1190
1191
1192
1193
1194
1195
1196
1197
1198
        super().fit(
            X,
            _y,
            sample_weight=sample_weight,
            init_score=init_score,
            eval_set=valid_sets,
            eval_names=eval_names,
            eval_sample_weight=eval_sample_weight,
            eval_class_weight=eval_class_weight,
            eval_init_score=eval_init_score,
            eval_metric=eval_metric,
            feature_name=feature_name,
            categorical_feature=categorical_feature,
            callbacks=callbacks,
            init_model=init_model
        )
wxchan's avatar
wxchan committed
1199
1200
        return self

1201
    _base_doc = LGBMModel.fit.__doc__.replace("self : LGBMModel", "self : LGBMClassifier")  # type: ignore
1202
1203
    _base_doc = (_base_doc[:_base_doc.find('group :')]  # type: ignore
                 + _base_doc[_base_doc.find('eval_set :'):])  # type: ignore
1204
1205
    fit.__doc__ = (_base_doc[:_base_doc.find('eval_group :')]
                   + _base_doc[_base_doc.find('eval_metric :'):])
1206

1207
1208
    def predict(
        self,
1209
        X: _LGBM_ScikitMatrixLike,
1210
1211
1212
1213
1214
1215
1216
1217
        raw_score: bool = False,
        start_iteration: int = 0,
        num_iteration: Optional[int] = None,
        pred_leaf: bool = False,
        pred_contrib: bool = False,
        validate_features: bool = False,
        **kwargs: Any
    ):
1218
        """Docstring is inherited from the LGBMModel."""
1219
1220
1221
1222
1223
1224
1225
1226
1227
1228
        result = self.predict_proba(
            X=X,
            raw_score=raw_score,
            start_iteration=start_iteration,
            num_iteration=num_iteration,
            pred_leaf=pred_leaf,
            pred_contrib=pred_contrib,
            validate_features=validate_features,
            **kwargs
        )
1229
        if callable(self._objective) or raw_score or pred_leaf or pred_contrib:
1230
1231
1232
1233
            return result
        else:
            class_index = np.argmax(result, axis=1)
            return self._le.inverse_transform(class_index)
wxchan's avatar
wxchan committed
1234

1235
1236
    predict.__doc__ = LGBMModel.predict.__doc__

1237
1238
    def predict_proba(
        self,
1239
        X: _LGBM_ScikitMatrixLike,
1240
1241
1242
1243
1244
1245
1246
1247
        raw_score: bool = False,
        start_iteration: int = 0,
        num_iteration: Optional[int] = None,
        pred_leaf: bool = False,
        pred_contrib: bool = False,
        validate_features: bool = False,
        **kwargs: Any
    ):
1248
        """Docstring is set after definition, using a template."""
1249
1250
1251
1252
1253
1254
1255
1256
1257
1258
        result = super().predict(
            X=X,
            raw_score=raw_score,
            start_iteration=start_iteration,
            num_iteration=num_iteration,
            pred_leaf=pred_leaf,
            pred_contrib=pred_contrib,
            validate_features=validate_features,
            **kwargs
        )
1259
        if callable(self._objective) and not (raw_score or pred_leaf or pred_contrib):
1260
1261
1262
            _log_warning("Cannot compute class probabilities or labels "
                         "due to the usage of customized objective function.\n"
                         "Returning raw scores instead.")
1263
            return result
1264
        elif self._n_classes > 2 or raw_score or pred_leaf or pred_contrib:  # type: ignore [operator]
1265
            return result
wxchan's avatar
wxchan committed
1266
        else:
1267
            return np.vstack((1. - result, result)).transpose()
1268

1269
1270
    predict_proba.__doc__ = _lgbmmodel_doc_predict.format(
        description="Return the predicted probability for each class for each sample.",
1271
        X_shape="numpy array, pandas DataFrame, H2O DataTable's Frame , scipy.sparse, list of lists of int or float of shape = [n_samples, n_features]",
1272
        output_name="predicted_probability",
1273
        predicted_result_shape="array-like of shape = [n_samples] or shape = [n_samples, n_classes]",
1274
1275
1276
1277
        X_leaves_shape="array-like of shape = [n_samples, n_trees] or shape = [n_samples, n_trees * n_classes]",
        X_SHAP_values_shape="array-like of shape = [n_samples, n_features + 1] or shape = [n_samples, (n_features + 1) * n_classes] or list with n_classes length of such objects"
    )

1278
    @property
1279
    def classes_(self) -> np.ndarray:
1280
        """:obj:`array` of shape = [n_classes]: The class label array."""
1281
        if not self.__sklearn_is_fitted__():
1282
            raise LGBMNotFittedError('No classes found. Need to call fit beforehand.')
1283
        return self._classes  # type: ignore[return-value]
1284
1285

    @property
1286
    def n_classes_(self) -> int:
1287
        """:obj:`int`: The number of classes."""
1288
        if not self.__sklearn_is_fitted__():
1289
1290
            raise LGBMNotFittedError('No classes found. Need to call fit beforehand.')
        return self._n_classes
wxchan's avatar
wxchan committed
1291

wxchan's avatar
wxchan committed
1292

wxchan's avatar
wxchan committed
1293
class LGBMRanker(LGBMModel):
1294
1295
1296
1297
1298
1299
1300
1301
    """LightGBM ranker.

    .. warning::

        scikit-learn doesn't support ranking applications yet,
        therefore this class is not really compatible with the sklearn ecosystem.
        Please use this class mainly for training and applying ranking models in common sklearnish way.
    """
wxchan's avatar
wxchan committed
1302

1303
    def fit(  # type: ignore[override]
1304
        self,
1305
1306
        X: _LGBM_ScikitMatrixLike,
        y: _LGBM_LabelType,
1307
1308
        sample_weight: Optional[_LGBM_WeightType] = None,
        init_score: Optional[_LGBM_InitScoreType] = None,
1309
        group: Optional[_LGBM_GroupType] = None,
1310
        eval_set: Optional[List[_LGBM_ScikitValidSet]] = None,
1311
        eval_names: Optional[List[str]] = None,
1312
1313
1314
        eval_sample_weight: Optional[List[_LGBM_WeightType]] = None,
        eval_init_score: Optional[List[_LGBM_InitScoreType]] = None,
        eval_group: Optional[List[_LGBM_GroupType]] = None,
1315
        eval_metric: Optional[_LGBM_ScikitEvalMetricType] = None,
1316
        eval_at: Union[List[int], Tuple[int, ...]] = (1, 2, 3, 4, 5),
1317
1318
        feature_name: _LGBM_FeatureNameConfiguration = 'auto',
        categorical_feature: _LGBM_CategoricalFeatureConfiguration = 'auto',
1319
        callbacks: Optional[List[Callable]] = None,
1320
        init_model: Optional[Union[str, Path, Booster, LGBMModel]] = None
1321
    ) -> "LGBMRanker":
1322
        """Docstring is inherited from the LGBMModel."""
1323
        # check group data
Guolin Ke's avatar
Guolin Ke committed
1324
        if group is None:
1325
            raise ValueError("Should set group for ranking task")
wxchan's avatar
wxchan committed
1326
1327

        if eval_set is not None:
Guolin Ke's avatar
Guolin Ke committed
1328
            if eval_group is None:
1329
                raise ValueError("Eval_group cannot be None when eval_set is not None")
Guolin Ke's avatar
Guolin Ke committed
1330
            elif len(eval_group) != len(eval_set):
1331
                raise ValueError("Length of eval_group should be equal to eval_set")
1332
            elif (isinstance(eval_group, dict)
1333
                  and any(i not in eval_group or eval_group[i] is None for i in range(len(eval_group)))
1334
1335
                  or isinstance(eval_group, list)
                  and any(group is None for group in eval_group)):
1336
1337
                raise ValueError("Should set group for all eval datasets for ranking task; "
                                 "if you use dict, the index should start from 0")
1338

1339
        self._eval_at = eval_at
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
        super().fit(
            X,
            y,
            sample_weight=sample_weight,
            init_score=init_score,
            group=group,
            eval_set=eval_set,
            eval_names=eval_names,
            eval_sample_weight=eval_sample_weight,
            eval_init_score=eval_init_score,
            eval_group=eval_group,
            eval_metric=eval_metric,
            feature_name=feature_name,
            categorical_feature=categorical_feature,
            callbacks=callbacks,
            init_model=init_model
        )
wxchan's avatar
wxchan committed
1357
        return self
1358

1359
    _base_doc = LGBMModel.fit.__doc__.replace("self : LGBMModel", "self : LGBMRanker")  # type: ignore
1360
1361
    fit.__doc__ = (_base_doc[:_base_doc.find('eval_class_weight :')]  # type: ignore
                   + _base_doc[_base_doc.find('eval_init_score :'):])  # type: ignore
1362
    _base_doc = fit.__doc__
1363
    _before_feature_name, _feature_name, _after_feature_name = _base_doc.partition('feature_name :')
1364
    fit.__doc__ = f"""{_before_feature_name}eval_at : list or tuple of int, optional (default=(1, 2, 3, 4, 5))
1365
        The evaluation positions of the specified metric.
1366
    {_feature_name}{_after_feature_name}"""