sklearn.py 64.4 KB
Newer Older
wxchan's avatar
wxchan committed
1
# coding: utf-8
2
"""Scikit-learn wrapper interface for LightGBM."""
3

4
import copy
5
from inspect import signature
6
from pathlib import Path
7
from typing import Any, Callable, Dict, List, Optional, Tuple, Union
8

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

12
from .basic import (
13
    _MULTICLASS_OBJECTIVES,
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
    Booster,
    Dataset,
    LightGBMError,
    _choose_param_value,
    _ConfigAliases,
    _LGBM_BoosterBestScoreType,
    _LGBM_CategoricalFeatureConfiguration,
    _LGBM_EvalFunctionResultType,
    _LGBM_FeatureNameConfiguration,
    _LGBM_GroupType,
    _LGBM_InitScoreType,
    _LGBM_LabelType,
    _LGBM_WeightType,
    _log_warning,
)
29
from .callback import _EvalResultDict, record_evaluation
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
from .compat import (
    SKLEARN_INSTALLED,
    LGBMNotFittedError,
    _LGBMAssertAllFinite,
    _LGBMCheckArray,
    _LGBMCheckClassificationTargets,
    _LGBMCheckSampleWeight,
    _LGBMCheckXY,
    _LGBMClassifierBase,
    _LGBMComputeSampleWeight,
    _LGBMCpuCount,
    _LGBMLabelEncoder,
    _LGBMModelBase,
    _LGBMRegressorBase,
    dt_DataTable,
    pd_DataFrame,
)
wxchan's avatar
wxchan committed
47
from .engine import train
48

49
__all__ = [
50
51
52
53
    "LGBMClassifier",
    "LGBMModel",
    "LGBMRanker",
    "LGBMRegressor",
54
55
]

56
57
58
59
60
_LGBM_ScikitMatrixLike = Union[
    dt_DataTable,
    List[Union[List[float], List[int]]],
    np.ndarray,
    pd_DataFrame,
61
    scipy.sparse.spmatrix,
62
]
63
_LGBM_ScikitCustomObjectiveFunction = Union[
64
    # f(labels, preds)
65
    Callable[
66
        [Optional[np.ndarray], np.ndarray],
67
        Tuple[np.ndarray, np.ndarray],
68
    ],
69
    # f(labels, preds, weights)
70
    Callable[
71
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray]],
72
        Tuple[np.ndarray, np.ndarray],
73
    ],
74
    # f(labels, preds, weights, group)
75
    Callable[
76
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray], Optional[np.ndarray]],
77
        Tuple[np.ndarray, np.ndarray],
78
    ],
79
80
]
_LGBM_ScikitCustomEvalFunction = Union[
81
    # f(labels, preds)
82
    Callable[
83
        [Optional[np.ndarray], np.ndarray],
84
        _LGBM_EvalFunctionResultType,
85
86
    ],
    Callable[
87
        [Optional[np.ndarray], np.ndarray],
88
        List[_LGBM_EvalFunctionResultType],
89
    ],
90
    # f(labels, preds, weights)
91
    Callable[
92
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray]],
93
        _LGBM_EvalFunctionResultType,
94
    ],
95
96
    Callable[
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray]],
97
        List[_LGBM_EvalFunctionResultType],
98
99
100
101
    ],
    # f(labels, preds, weights, group)
    Callable[
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray], Optional[np.ndarray]],
102
        _LGBM_EvalFunctionResultType,
103
104
105
    ],
    Callable[
        [Optional[np.ndarray], np.ndarray, Optional[np.ndarray], Optional[np.ndarray]],
106
107
        List[_LGBM_EvalFunctionResultType],
    ],
108
]
109
110
111
_LGBM_ScikitEvalMetricType = Union[
    str,
    _LGBM_ScikitCustomEvalFunction,
112
    List[Union[str, _LGBM_ScikitCustomEvalFunction]],
113
]
114
_LGBM_ScikitValidSet = Tuple[_LGBM_ScikitMatrixLike, _LGBM_LabelType]
115

wxchan's avatar
wxchan committed
116

117
118
119
120
121
122
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."
    )
123
    assert group is None or isinstance(group, np.ndarray), error_msg
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
    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."
    )
143
    assert weight is None or isinstance(weight, np.ndarray), error_msg
144
145
146
    return weight


147
class _ObjectiveFunctionWrapper:
148
    """Proxy class for objective function."""
149

150
    def __init__(self, func: _LGBM_ScikitCustomObjectiveFunction):
151
        """Construct a proxy class.
152

153
154
        This class transforms objective function to match objective function with signature ``new_func(preds, dataset)``
        as expected by ``lightgbm.engine.train``.
155

156
157
158
        Parameters
        ----------
        func : callable
159
160
161
162
            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)``
163
164
            and returns (grad, hess):

165
                y_true : numpy 1-D array of shape = [n_samples]
166
                    The target values.
167
                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)
168
                    The predicted values.
169
170
                    Predicted values are returned before any transformation,
                    e.g. they are raw margin instead of probability of positive class for binary task.
171
172
                weight : numpy 1-D array of shape = [n_samples]
                    The weight of samples. Weights should be non-negative.
173
                group : numpy 1-D array
174
175
176
                    Group/query data.
                    Only used in the learning-to-rank task.
                    sum(group) = n_samples.
177
178
                    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.
179
                grad : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape [n_samples, n_classes] (for multi-class task)
180
181
                    The value of the first order derivative (gradient) of the loss
                    with respect to the elements of y_pred for each sample point.
182
                hess : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
183
184
                    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
185

Nikita Titov's avatar
Nikita Titov committed
186
187
        .. note::

188
            For multi-class task, y_pred is a numpy 2-D array of shape = [n_samples, n_classes],
189
            and grad and hess should be returned in the same format.
190
191
        """
        self.func = func
wxchan's avatar
wxchan committed
192

193
194
195
196
197
    def __call__(
        self,
        preds: np.ndarray,
        dataset: Dataset,
    ) -> Tuple[np.ndarray, np.ndarray]:
198
199
200
201
        """Call passed function with appropriate arguments.

        Parameters
        ----------
202
        preds : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
203
204
205
206
207
208
            The predicted values.
        dataset : Dataset
            The training dataset.

        Returns
        -------
209
        grad : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
210
211
            The value of the first order derivative (gradient) of the loss
            with respect to the elements of preds for each sample point.
212
        hess : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
213
214
            The value of the second order derivative (Hessian) of the loss
            with respect to the elements of preds for each sample point.
215
        """
216
        labels = _get_label_from_constructed_dataset(dataset)
217
        argc = len(signature(self.func).parameters)
218
        if argc == 2:
219
            grad, hess = self.func(labels, preds)  # type: ignore[call-arg]
220
221
222
223
224
225
226
227
228
229
230
231
            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
232

wxchan's avatar
wxchan committed
233

234
class _EvalFunctionWrapper:
235
    """Proxy class for evaluation function."""
236

237
    def __init__(self, func: _LGBM_ScikitCustomEvalFunction):
238
        """Construct a proxy class.
239

240
241
        This class transforms evaluation function to match evaluation function with signature ``new_func(preds, dataset)``
        as expected by ``lightgbm.engine.train``.
242

243
244
245
246
247
248
249
250
251
252
        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):

253
                y_true : numpy 1-D array of shape = [n_samples]
254
                    The target values.
255
                y_pred : numpy 1-D array of shape = [n_samples] or numpy 2-D array shape = [n_samples, n_classes] (for multi-class task)
256
                    The predicted values.
257
258
                    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.
259
                weight : numpy 1-D array of shape = [n_samples]
260
                    The weight of samples. Weights should be non-negative.
261
                group : numpy 1-D array
262
263
264
                    Group/query data.
                    Only used in the learning-to-rank task.
                    sum(group) = n_samples.
265
266
                    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.
267
                eval_name : str
Andrew Ziem's avatar
Andrew Ziem committed
268
                    The name of evaluation function (without whitespace).
269
270
271
272
273
274
                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
275

276
277
278
    def __call__(
        self,
        preds: np.ndarray,
279
        dataset: Dataset,
280
    ) -> Union[_LGBM_EvalFunctionResultType, List[_LGBM_EvalFunctionResultType]]:
281
        """Call passed function with appropriate arguments.
282

283
284
        Parameters
        ----------
285
        preds : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
286
287
288
289
290
291
            The predicted values.
        dataset : Dataset
            The training dataset.

        Returns
        -------
292
        eval_name : str
Andrew Ziem's avatar
Andrew Ziem committed
293
            The name of evaluation function (without whitespace).
294
295
296
297
298
        eval_result : float
            The eval result.
        is_higher_better : bool
            Is eval result higher better, e.g. AUC is ``is_higher_better``.
        """
299
        labels = _get_label_from_constructed_dataset(dataset)
300
        argc = len(signature(self.func).parameters)
301
        if argc == 2:
302
            return self.func(labels, preds)  # type: ignore[call-arg]
303
304
305
306
307
308
309
310
311
312

        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}")
313

wxchan's avatar
wxchan committed
314

315
316
317
# documentation templates for LGBMModel methods are shared between the classes in
# this module and those in the ``dask`` module

318
_lgbmmodel_doc_fit = """
319
320
321
322
323
324
325
326
327
    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}
328
        Weights of training data. Weights should be non-negative.
329
    init_score : {init_score_shape}
330
331
332
333
334
335
336
337
338
        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.
339
    eval_names : list of str, or None, optional (default=None)
340
        Names of eval_set.
341
    eval_sample_weight : {eval_sample_weight_shape}
342
        Weights of eval data. Weights should be non-negative.
343
344
    eval_class_weight : list or None, optional (default=None)
        Class weights of eval data.
345
    eval_init_score : {eval_init_score_shape}
346
        Init score of eval data.
347
    eval_group : {eval_group_shape}
348
        Group data of eval data.
349
350
    eval_metric : str, callable, list or None, optional (default=None)
        If str, it should be a built-in evaluation metric to use.
351
352
353
354
        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.
355
    feature_name : list of str, or 'auto', optional (default='auto')
356
357
        Feature names.
        If 'auto' and data is pandas DataFrame, data columns names are used.
358
    categorical_feature : list of str or int, or 'auto', optional (default='auto')
359
360
        Categorical features.
        If list of int, interpreted as indices.
361
        If list of str, interpreted as feature names (need to specify ``feature_name`` as well).
362
        If 'auto' and data is pandas DataFrame, pandas unordered categorical columns are used.
363
        All values in categorical features will be cast to int32 and thus should be less than int32 max value (2147483647).
364
365
366
        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.
367
        Floating point numbers in categorical features will be rounded towards 0.
368
    callbacks : list of callable, or None, optional (default=None)
369
370
        List of callback functions that are applied at each iteration.
        See Callbacks in Python API for more information.
371
    init_model : str, pathlib.Path, Booster, LGBMModel or None, optional (default=None)
372
373
374
375
        Filename of LightGBM model, Booster instance or LGBMModel instance used for continue training.

    Returns
    -------
376
    self : LGBMModel
377
378
379
380
381
382
383
384
385
386
387
388
        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):

389
        y_true : numpy 1-D array of shape = [n_samples]
390
            The target values.
391
        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)
392
            The predicted values.
393
394
            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.
395
        weight : numpy 1-D array of shape = [n_samples]
396
            The weight of samples. Weights should be non-negative.
397
        group : numpy 1-D array
398
399
400
401
402
            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.
403
        eval_name : str
Andrew Ziem's avatar
Andrew Ziem committed
404
            The name of evaluation function (without whitespace).
405
406
407
408
409
410
        eval_result : float
            The eval result.
        is_higher_better : bool
            Is eval result higher better, e.g. AUC is ``is_higher_better``.
"""

411
_lgbmmodel_doc_predict = """
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
    {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.

441
442
443
    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.
444
445
446
447
448
449
450
451
452
453
454
455
456
457
    **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.
    """


458
459
460
461
462
463
464
465
466
467
468
469
470
def _extract_evaluation_meta_data(
    *,
    collection: Optional[Union[Dict[Any, Any], List[Any]]],
    name: str,
    i: int,
) -> Optional[Any]:
    """Try to extract the ith element of one of the ``eval_*`` inputs."""
    if collection is None:
        return None
    elif isinstance(collection, list):
        # It's possible, for example, to pass 3 eval sets through `eval_set`,
        # but only 1 init_score through `eval_init_score`.
        #
471
        # This if-else accounts for that possibility.
472
473
474
475
476
477
478
479
480
481
        if len(collection) > i:
            return collection[i]
        else:
            return None
    elif isinstance(collection, dict):
        return collection.get(i, None)
    else:
        raise TypeError(f"{name} should be dict or list")


482
483
class LGBMModel(_LGBMModelBase):
    """Implementation of the scikit-learn API for LightGBM."""
wxchan's avatar
wxchan committed
484

485
486
    def __init__(
        self,
487
        boosting_type: str = "gbdt",
488
489
490
491
492
        num_leaves: int = 31,
        max_depth: int = -1,
        learning_rate: float = 0.1,
        n_estimators: int = 100,
        subsample_for_bin: int = 200000,
493
        objective: Optional[Union[str, _LGBM_ScikitCustomObjectiveFunction]] = None,
494
        class_weight: Optional[Union[Dict, str]] = None,
495
        min_split_gain: float = 0.0,
496
497
        min_child_weight: float = 1e-3,
        min_child_samples: int = 20,
498
        subsample: float = 1.0,
499
        subsample_freq: int = 0,
500
501
502
        colsample_bytree: float = 1.0,
        reg_alpha: float = 0.0,
        reg_lambda: float = 0.0,
503
        random_state: Optional[Union[int, np.random.RandomState, np.random.Generator]] = None,
504
        n_jobs: Optional[int] = None,
505
        importance_type: str = "split",
506
        **kwargs: Any,
507
    ):
508
        r"""Construct a gradient boosting model.
wxchan's avatar
wxchan committed
509
510
511

        Parameters
        ----------
512
        boosting_type : str, optional (default='gbdt')
513
514
515
516
            '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
517
            Maximum tree leaves for base learners.
518
        max_depth : int, optional (default=-1)
519
            Maximum tree depth for base learners, <=0 means no limit.
520
            If setting this to a positive value, consider also changing ``num_leaves`` to ``<= 2^max_depth``.
521
        learning_rate : float, optional (default=0.1)
522
            Boosting learning rate.
523
524
525
            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.
526
        n_estimators : int, optional (default=100)
wxchan's avatar
wxchan committed
527
            Number of boosted trees to fit.
528
        subsample_for_bin : int, optional (default=200000)
wxchan's avatar
wxchan committed
529
            Number of samples for constructing bins.
530
        objective : str, callable or None, optional (default=None)
wxchan's avatar
wxchan committed
531
532
            Specify the learning task and the corresponding learning objective or
            a custom objective function to be used (see note below).
533
            Default: 'regression' for LGBMRegressor, 'binary' or 'multiclass' for LGBMClassifier, 'lambdarank' for LGBMRanker.
534
535
536
537
        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.
538
539
540
            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.
541
542
543
            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.
544
            Note, that these weights will be multiplied with ``sample_weight`` (passed through the ``fit`` method)
545
            if ``sample_weight`` is specified.
546
        min_split_gain : float, optional (default=0.)
wxchan's avatar
wxchan committed
547
            Minimum loss reduction required to make a further partition on a leaf node of the tree.
548
        min_child_weight : float, optional (default=1e-3)
549
            Minimum sum of instance weight (Hessian) needed in a child (leaf).
550
        min_child_samples : int, optional (default=20)
551
            Minimum number of data needed in a child (leaf).
552
        subsample : float, optional (default=1.)
wxchan's avatar
wxchan committed
553
            Subsample ratio of the training instance.
554
        subsample_freq : int, optional (default=0)
Andrew Ziem's avatar
Andrew Ziem committed
555
            Frequency of subsample, <=0 means no enable.
556
        colsample_bytree : float, optional (default=1.)
wxchan's avatar
wxchan committed
557
            Subsample ratio of columns when constructing each tree.
558
        reg_alpha : float, optional (default=0.)
559
            L1 regularization term on weights.
560
        reg_lambda : float, optional (default=0.)
561
            L2 regularization term on weights.
562
        random_state : int, RandomState object or None, optional (default=None)
wxchan's avatar
wxchan committed
563
            Random number seed.
564
            If int, this number is used to seed the C++ code.
565
            If RandomState or Generator object (numpy), a random integer is picked based on its state to seed the C++ code.
566
            If None, default seeds in C++ code are used.
567
568
569
570
571
572
573
574
575
576
577
578
        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).
579
580
581

            .. versionchanged:: 4.0.0

582
        importance_type : str, optional (default='split')
583
            The type of feature importance to be filled into ``feature_importances_``.
584
585
586
587
            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
588
            Check http://lightgbm.readthedocs.io/en/latest/Parameters.html for more parameters.
589

Nikita Titov's avatar
Nikita Titov committed
590
591
592
            .. warning::

                \*\*kwargs is not supported in sklearn, it may cause unexpected issues.
wxchan's avatar
wxchan committed
593
594
595

        Note
        ----
596
597
        A custom objective function can be provided for the ``objective`` parameter.
        In this case, it should have the signature
598
599
600
        ``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
601

602
            y_true : numpy 1-D array of shape = [n_samples]
603
                The target values.
604
            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)
605
                The predicted values.
606
607
                Predicted values are returned before any transformation,
                e.g. they are raw margin instead of probability of positive class for binary task.
608
609
            weight : numpy 1-D array of shape = [n_samples]
                The weight of samples. Weights should be non-negative.
610
            group : numpy 1-D array
611
612
613
                Group/query data.
                Only used in the learning-to-rank task.
                sum(group) = n_samples.
614
615
                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.
616
            grad : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
617
618
                The value of the first order derivative (gradient) of the loss
                with respect to the elements of y_pred for each sample point.
619
            hess : numpy 1-D array of shape = [n_samples] or numpy 2-D array of shape = [n_samples, n_classes] (for multi-class task)
620
621
                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
622

623
        For multi-class task, y_pred is a numpy 2-D array of shape = [n_samples, n_classes],
624
        and grad and hess should be returned in the same format.
wxchan's avatar
wxchan committed
625
        """
wxchan's avatar
wxchan committed
626
        if not SKLEARN_INSTALLED:
627
628
629
630
            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
631

632
        self.boosting_type = boosting_type
633
        self.objective = objective
wxchan's avatar
wxchan committed
634
635
636
637
        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
638
        self.subsample_for_bin = subsample_for_bin
wxchan's avatar
wxchan committed
639
640
641
642
643
644
645
646
        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
647
648
        self.random_state = random_state
        self.n_jobs = n_jobs
649
        self.importance_type = importance_type
650
        self._Booster: Optional[Booster] = None
651
652
        self._evals_result: _EvalResultDict = {}
        self._best_score: _LGBM_BoosterBestScoreType = {}
653
        self._best_iteration: int = -1
654
        self._other_params: Dict[str, Any] = {}
655
        self._objective = objective
656
        self.class_weight = class_weight
657
658
        self._class_weight: Optional[Union[Dict, str]] = None
        self._class_map: Optional[Dict[int, int]] = None
659
660
        self._n_features: int = -1
        self._n_features_in: int = -1
661
        self._classes: Optional[np.ndarray] = None
662
        self._n_classes: int = -1
663
        self.set_params(**kwargs)
wxchan's avatar
wxchan committed
664

665
    def _more_tags(self) -> Dict[str, Any]:
666
        return {
667
668
669
670
671
672
673
            "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)"
            },
674
        }
Nikita Titov's avatar
Nikita Titov committed
675

676
677
678
    def __sklearn_is_fitted__(self) -> bool:
        return getattr(self, "fitted_", False)

679
    def get_params(self, deep: bool = True) -> Dict[str, Any]:
680
681
682
683
684
685
686
687
688
689
690
691
692
        """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.
        """
693
        params = super().get_params(deep=deep)
694
        params.update(self._other_params)
wxchan's avatar
wxchan committed
695
696
        return params

697
    def set_params(self, **params: Any) -> "LGBMModel":
698
699
700
701
702
703
704
705
706
707
708
709
        """Set the parameters of this estimator.

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

        Returns
        -------
        self : object
            Returns self.
        """
wxchan's avatar
wxchan committed
710
711
        for key, value in params.items():
            setattr(self, key, value)
712
713
            if hasattr(self, f"_{key}"):
                setattr(self, f"_{key}", value)
714
            self._other_params[key] = value
wxchan's avatar
wxchan committed
715
        return self
wxchan's avatar
wxchan committed
716

717
718
719
720
721
722
723
724
725
726
727
728
729
730
    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"}
731
732
        params = self.get_params()

733
734
        params.pop("objective", None)
        for alias in _ConfigAliases.get("objective"):
735
            if alias in params:
736
                obj = params.pop(alias)
737
                _log_warning(f"Found '{alias}' in params. Will use it instead of 'objective' argument")
738
739
740
741
742
743
744
745
746
747
748
749
750
                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"
751
                else:
752
                    raise ValueError("Unknown LGBMModel type.")
753
        if callable(self._objective):
754
            if stage == "fit":
755
                params["objective"] = _ObjectiveFunctionWrapper(self._objective)
756
            else:
757
                params["objective"] = "None"
758
        else:
759
            params["objective"] = self._objective
760

761
762
763
        params.pop("importance_type", None)
        params.pop("n_estimators", None)
        params.pop("class_weight", None)
764

765
766
        if isinstance(params["random_state"], np.random.RandomState):
            params["random_state"] = params["random_state"].randint(np.iinfo(np.int32).max)
767
        elif isinstance(params["random_state"], np.random.Generator):
768
            params["random_state"] = int(params["random_state"].integers(np.iinfo(np.int32).max))
769
        if self._n_classes > 2:
770
            for alias in _ConfigAliases.get("num_class"):
771
                params.pop(alias, None)
772
773
            params["num_class"] = self._n_classes
        if hasattr(self, "_eval_at"):
774
            eval_at = self._eval_at
775
            for alias in _ConfigAliases.get("eval_at"):
776
777
778
                if alias in params:
                    _log_warning(f"Found '{alias}' in params. Will use it instead of 'eval_at' argument")
                    eval_at = params.pop(alias)
779
            params["eval_at"] = eval_at
wxchan's avatar
wxchan committed
780

781
        # register default metric for consistency with callable eval_metric case
782
        original_metric = self._objective if isinstance(self._objective, str) else None
783
784
785
786
787
788
789
790
791
792
        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
793
        params = _choose_param_value("metric", params, original_metric)
794

795
796
797
798
799
800
        # 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"])

801
802
        return params

803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
    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

823
824
    def fit(
        self,
825
826
        X: _LGBM_ScikitMatrixLike,
        y: _LGBM_LabelType,
827
828
        sample_weight: Optional[_LGBM_WeightType] = None,
        init_score: Optional[_LGBM_InitScoreType] = None,
829
        group: Optional[_LGBM_GroupType] = None,
830
        eval_set: Optional[List[_LGBM_ScikitValidSet]] = None,
831
        eval_names: Optional[List[str]] = None,
832
833
834
835
        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,
836
        eval_metric: Optional[_LGBM_ScikitEvalMetricType] = None,
837
838
        feature_name: _LGBM_FeatureNameConfiguration = "auto",
        categorical_feature: _LGBM_CategoricalFeatureConfiguration = "auto",
839
        callbacks: Optional[List[Callable]] = None,
840
        init_model: Optional[Union[str, Path, Booster, "LGBMModel"]] = None,
841
    ) -> "LGBMModel":
842
843
844
845
846
        """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
847
848
849
850
851
852
853
        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)]
854
855
856
857
858

        # 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)]

859
        # concatenate metric from params (or default if not provided in params) and eval_metric
860
861
862
        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"]
        params["metric"] = [metric for metric in params["metric"] if metric is not None]
wxchan's avatar
wxchan committed
863

864
        if not isinstance(X, (pd_DataFrame, dt_DataTable)):
865
            _X, _y = _LGBMCheckXY(X, y, accept_sparse=True, force_all_finite=False, ensure_min_samples=2)
866
867
            if sample_weight is not None:
                sample_weight = _LGBMCheckSampleWeight(sample_weight, _X)
868
869
        else:
            _X, _y = X, y
870

871
872
873
874
        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)
875
876
877
878
            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)
879

880
        self._n_features = _X.shape[1]
881
882
        # copy for consistency
        self._n_features_in = self._n_features
883

884
885
886
887
888
889
890
        train_set = Dataset(
            data=_X,
            label=_y,
            weight=sample_weight,
            group=group,
            init_score=init_score,
            categorical_feature=categorical_feature,
891
            feature_name=feature_name,
892
893
            params=params,
        )
Guolin Ke's avatar
Guolin Ke committed
894

895
        valid_sets: List[Dataset] = []
Guolin Ke's avatar
Guolin Ke committed
896
897
898
899
        if eval_set is not None:
            if isinstance(eval_set, tuple):
                eval_set = [eval_set]
            for i, valid_data in enumerate(eval_set):
900
                # reduce cost for prediction training data
Guolin Ke's avatar
Guolin Ke committed
901
902
903
                if valid_data[0] is X and valid_data[1] is y:
                    valid_set = train_set
                else:
904
905
906
907
908
909
910
911
912
913
                    valid_weight = _extract_evaluation_meta_data(
                        collection=eval_sample_weight,
                        name="eval_sample_weight",
                        i=i,
                    )
                    valid_class_weight = _extract_evaluation_meta_data(
                        collection=eval_class_weight,
                        name="eval_class_weight",
                        i=i,
                    )
914
915
916
917
                    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])
918
919
920
921
                        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)
922
923
924
925
926
927
928
929
930
931
                    valid_init_score = _extract_evaluation_meta_data(
                        collection=eval_init_score,
                        name="eval_init_score",
                        i=i,
                    )
                    valid_group = _extract_evaluation_meta_data(
                        collection=eval_group,
                        name="eval_group",
                        i=i,
                    )
932
933
934
935
936
937
938
939
940
                    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,
                    )
941

Guolin Ke's avatar
Guolin Ke committed
942
943
                valid_sets.append(valid_set)

944
945
946
        if isinstance(init_model, LGBMModel):
            init_model = init_model.booster_

947
948
949
        if callbacks is None:
            callbacks = []
        else:
950
            callbacks = copy.copy(callbacks)  # don't use deepcopy here to allow non-serializable objects
951

952
        evals_result: _EvalResultDict = {}
953
954
955
956
957
958
959
960
        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,
961
            feval=eval_metrics_callable,  # type: ignore[arg-type]
962
            init_model=init_model,
963
            callbacks=callbacks,
964
        )
wxchan's avatar
wxchan committed
965

966
        self._evals_result = evals_result
967
        self._best_iteration = self._Booster.best_iteration
968
        self._best_score = self._Booster.best_score
wxchan's avatar
wxchan committed
969

970
971
        self.fitted_ = True

wxchan's avatar
wxchan committed
972
        # free dataset
973
        self._Booster.free_dataset()
wxchan's avatar
wxchan committed
974
        del train_set, valid_sets
wxchan's avatar
wxchan committed
975
976
        return self

977
978
979
980
981
982
983
984
985
986
987
988
989
990
    fit.__doc__ = (
        _lgbmmodel_doc_fit.format(
            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]",
            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)",
            group_shape="numpy array, pandas Series, list of int or float, or None, optional (default=None)",
            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)",
        )
        + "\n\n"
        + _lgbmmodel_doc_custom_eval_note
    )
991

992
993
    def predict(
        self,
994
        X: _LGBM_ScikitMatrixLike,
995
996
997
998
999
1000
        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,
1001
        **kwargs: Any,
1002
    ):
1003
        """Docstring is set after definition, using a template."""
1004
        if not self.__sklearn_is_fitted__():
1005
            raise LGBMNotFittedError("Estimator not fitted, call fit before exploiting the model.")
1006
        if not isinstance(X, (pd_DataFrame, dt_DataTable)):
1007
            X = _LGBMCheckArray(X, accept_sparse=True, force_all_finite=False)
1008
1009
        n_features = X.shape[1]
        if self._n_features != n_features:
1010
1011
1012
1013
1014
            raise ValueError(
                "Number of features of the model must "
                f"match the input. Model n_features_ is {self._n_features} and "
                f"input n_features is {n_features}"
            )
1015
        # retrieve original params that possibly can be used in both training and prediction
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
        # 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",
1026
            *kwargs.keys(),
1027
1028
1029
        ):
            predict_params.pop(alias, None)
        predict_params.update(kwargs)
1030
1031
1032

        # 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
1033
1034
        predict_params = _choose_param_value("num_threads", predict_params, self.n_jobs)
        predict_params["num_threads"] = self._process_n_jobs(predict_params["num_threads"])
1035

1036
        return self._Booster.predict(  # type: ignore[union-attr]
1037
1038
1039
1040
1041
1042
1043
1044
            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,
1045
        )
wxchan's avatar
wxchan committed
1046

1047
1048
    predict.__doc__ = _lgbmmodel_doc_predict.format(
        description="Return the predicted value for each sample.",
1049
        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]",
1050
1051
1052
        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]",
1053
        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",
1054
1055
    )

1056
    @property
1057
    def n_features_(self) -> int:
1058
        """:obj:`int`: The number of features of fitted model."""
1059
        if not self.__sklearn_is_fitted__():
1060
            raise LGBMNotFittedError("No n_features found. Need to call fit beforehand.")
1061
1062
        return self._n_features

1063
    @property
1064
    def n_features_in_(self) -> int:
1065
        """:obj:`int`: The number of features of fitted model."""
1066
        if not self.__sklearn_is_fitted__():
1067
            raise LGBMNotFittedError("No n_features_in found. Need to call fit beforehand.")
1068
1069
        return self._n_features_in

1070
    @property
1071
    def best_score_(self) -> _LGBM_BoosterBestScoreType:
1072
        """:obj:`dict`: The best score of fitted model."""
1073
        if not self.__sklearn_is_fitted__():
1074
            raise LGBMNotFittedError("No best_score found. Need to call fit beforehand.")
1075
1076
1077
        return self._best_score

    @property
1078
    def best_iteration_(self) -> int:
1079
        """:obj:`int`: The best iteration of fitted model if ``early_stopping()`` callback has been specified."""
1080
        if not self.__sklearn_is_fitted__():
1081
1082
1083
            raise LGBMNotFittedError(
                "No best_iteration found. Need to call fit with early_stopping callback beforehand."
            )
1084
1085
1086
        return self._best_iteration

    @property
1087
    def objective_(self) -> Union[str, _LGBM_ScikitCustomObjectiveFunction]:
1088
        """:obj:`str` or :obj:`callable`: The concrete objective used while fitting this model."""
1089
        if not self.__sklearn_is_fitted__():
1090
            raise LGBMNotFittedError("No objective found. Need to call fit beforehand.")
1091
        return self._objective  # type: ignore[return-value]
1092

1093
1094
1095
1096
1097
1098
    @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``.
1099

1100
        .. versionadded:: 4.0.0
1101
1102
        """
        if not self.__sklearn_is_fitted__():
1103
            raise LGBMNotFittedError("No n_estimators found. Need to call fit beforehand.")
1104
        return self._Booster.current_iteration()  # type: ignore
1105
1106
1107
1108
1109
1110
1111

    @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``.
1112

1113
        .. versionadded:: 4.0.0
1114
1115
        """
        if not self.__sklearn_is_fitted__():
1116
            raise LGBMNotFittedError("No n_iter found. Need to call fit beforehand.")
1117
        return self._Booster.current_iteration()  # type: ignore
1118

1119
    @property
1120
    def booster_(self) -> Booster:
1121
        """Booster: The underlying Booster of this model."""
1122
        if not self.__sklearn_is_fitted__():
1123
            raise LGBMNotFittedError("No booster found. Need to call fit beforehand.")
1124
        return self._Booster  # type: ignore[return-value]
wxchan's avatar
wxchan committed
1125

1126
    @property
1127
    def evals_result_(self) -> _EvalResultDict:
1128
        """:obj:`dict`: The evaluation results if validation sets have been specified."""
1129
        if not self.__sklearn_is_fitted__():
1130
            raise LGBMNotFittedError("No results found. Need to call fit with eval_set beforehand.")
1131
        return self._evals_result
1132
1133

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

Nikita Titov's avatar
Nikita Titov committed
1137
1138
1139
1140
        .. note::

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

1146
    @property
1147
    def feature_name_(self) -> List[str]:
1148
1149
1150
1151
1152
1153
        """:obj:`list` of shape = [n_features]: The names of features.

        .. note::

            If input does not contain feature names, they will be added during fitting in the format ``Column_0``, ``Column_1``, ..., ``Column_N``.
        """
1154
        if not self.__sklearn_is_fitted__():
1155
            raise LGBMNotFittedError("No feature_name found. Need to call fit beforehand.")
1156
        return self._Booster.feature_name()  # type: ignore[union-attr]
1157

1158
1159
1160
1161
1162
1163
1164
    @property
    def feature_names_in_(self) -> np.ndarray:
        """:obj:`array` of shape = [n_features]: scikit-learn compatible version of ``.feature_name_``."""
        if not self.__sklearn_is_fitted__():
            raise LGBMNotFittedError("No feature_names_in_ found. Need to call fit beforehand.")
        return np.array(self.feature_name_)

wxchan's avatar
wxchan committed
1165

1166
class LGBMRegressor(_LGBMRegressorBase, LGBMModel):
1167
    """LightGBM regressor."""
wxchan's avatar
wxchan committed
1168

1169
    def fit(  # type: ignore[override]
1170
        self,
1171
1172
        X: _LGBM_ScikitMatrixLike,
        y: _LGBM_LabelType,
1173
1174
1175
        sample_weight: Optional[_LGBM_WeightType] = None,
        init_score: Optional[_LGBM_InitScoreType] = None,
        eval_set: Optional[List[_LGBM_ScikitValidSet]] = None,
1176
        eval_names: Optional[List[str]] = None,
1177
1178
        eval_sample_weight: Optional[List[_LGBM_WeightType]] = None,
        eval_init_score: Optional[List[_LGBM_InitScoreType]] = None,
1179
        eval_metric: Optional[_LGBM_ScikitEvalMetricType] = None,
1180
1181
        feature_name: _LGBM_FeatureNameConfiguration = "auto",
        categorical_feature: _LGBM_CategoricalFeatureConfiguration = "auto",
1182
        callbacks: Optional[List[Callable]] = None,
1183
        init_model: Optional[Union[str, Path, Booster, LGBMModel]] = None,
1184
    ) -> "LGBMRegressor":
1185
        """Docstring is inherited from the LGBMModel."""
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=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,
1199
            init_model=init_model,
1200
        )
Guolin Ke's avatar
Guolin Ke committed
1201
1202
        return self

1203
    _base_doc = LGBMModel.fit.__doc__.replace("self : LGBMModel", "self : LGBMRegressor")  # type: ignore
1204
1205
1206
1207
1208
1209
    _base_doc = (
        _base_doc[: _base_doc.find("group :")]  # type: ignore
        + _base_doc[_base_doc.find("eval_set :") :]
    )  # type: ignore
    _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
1210

1211

1212
class LGBMClassifier(_LGBMClassifierBase, LGBMModel):
1213
    """LightGBM classifier."""
wxchan's avatar
wxchan committed
1214

1215
    def fit(  # type: ignore[override]
1216
        self,
1217
1218
        X: _LGBM_ScikitMatrixLike,
        y: _LGBM_LabelType,
1219
1220
1221
        sample_weight: Optional[_LGBM_WeightType] = None,
        init_score: Optional[_LGBM_InitScoreType] = None,
        eval_set: Optional[List[_LGBM_ScikitValidSet]] = None,
1222
        eval_names: Optional[List[str]] = None,
1223
1224
1225
        eval_sample_weight: Optional[List[_LGBM_WeightType]] = None,
        eval_class_weight: Optional[List[float]] = None,
        eval_init_score: Optional[List[_LGBM_InitScoreType]] = None,
1226
        eval_metric: Optional[_LGBM_ScikitEvalMetricType] = None,
1227
1228
        feature_name: _LGBM_FeatureNameConfiguration = "auto",
        categorical_feature: _LGBM_CategoricalFeatureConfiguration = "auto",
1229
        callbacks: Optional[List[Callable]] = None,
1230
        init_model: Optional[Union[str, Path, Booster, LGBMModel]] = None,
1231
    ) -> "LGBMClassifier":
1232
        """Docstring is inherited from the LGBMModel."""
1233
        _LGBMAssertAllFinite(y)
1234
1235
        _LGBMCheckClassificationTargets(y)
        self._le = _LGBMLabelEncoder().fit(y)
1236
        _y = self._le.transform(y)
1237
        self._class_map = dict(zip(self._le.classes_, self._le.transform(self._le.classes_)))
1238
1239
        if isinstance(self.class_weight, dict):
            self._class_weight = {self._class_map[k]: v for k, v in self.class_weight.items()}
1240

1241
        self._classes = self._le.classes_
1242
        self._n_classes = len(self._classes)  # type: ignore[arg-type]
1243
1244
        if self.objective is None:
            self._objective = None
1245

1246
1247
        # adjust eval metrics to match whether binary or multiclass
        # classification is being performed
1248
        if not callable(eval_metric):
1249
1250
1251
1252
1253
1254
            if isinstance(eval_metric, list):
                eval_metric_list = eval_metric
            elif isinstance(eval_metric, str):
                eval_metric_list = [eval_metric]
            else:
                eval_metric_list = []
1255
            if self.__is_multiclass:
1256
                for index, metric in enumerate(eval_metric_list):
1257
                    if metric in {"logloss", "binary_logloss"}:
1258
                        eval_metric_list[index] = "multi_logloss"
1259
                    elif metric in {"error", "binary_error"}:
1260
                        eval_metric_list[index] = "multi_error"
1261
            else:
1262
                for index, metric in enumerate(eval_metric_list):
1263
1264
1265
1266
                    if metric in {"logloss", "multi_logloss"}:
                        eval_metric_list[index] = "binary_logloss"
                    elif metric in {"error", "multi_error"}:
                        eval_metric_list[index] = "binary_error"
1267
            eval_metric = eval_metric_list
wxchan's avatar
wxchan committed
1268

1269
        # do not modify args, as it causes errors in model selection tools
1270
        valid_sets: Optional[List[_LGBM_ScikitValidSet]] = None
wxchan's avatar
wxchan committed
1271
        if eval_set is not None:
1272
1273
            if isinstance(eval_set, tuple):
                eval_set = [eval_set]
1274
1275
            valid_sets = []
            for valid_x, valid_y in eval_set:
1276
                if valid_x is X and valid_y is y:
1277
                    valid_sets.append((valid_x, _y))
1278
                else:
1279
                    valid_sets.append((valid_x, self._le.transform(valid_y)))
1280

1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
1293
1294
        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,
1295
            init_model=init_model,
1296
        )
wxchan's avatar
wxchan committed
1297
1298
        return self

1299
    _base_doc = LGBMModel.fit.__doc__.replace("self : LGBMModel", "self : LGBMClassifier")  # type: ignore
1300
1301
1302
1303
1304
    _base_doc = (
        _base_doc[: _base_doc.find("group :")]  # type: ignore
        + _base_doc[_base_doc.find("eval_set :") :]
    )  # type: ignore
    fit.__doc__ = _base_doc[: _base_doc.find("eval_group :")] + _base_doc[_base_doc.find("eval_metric :") :]
1305

1306
1307
    def predict(
        self,
1308
        X: _LGBM_ScikitMatrixLike,
1309
1310
1311
1312
1313
1314
        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,
1315
        **kwargs: Any,
1316
    ):
1317
        """Docstring is inherited from the LGBMModel."""
1318
1319
1320
1321
1322
1323
1324
1325
        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,
1326
            **kwargs,
1327
        )
1328
        if callable(self._objective) or raw_score or pred_leaf or pred_contrib:
1329
1330
1331
1332
            return result
        else:
            class_index = np.argmax(result, axis=1)
            return self._le.inverse_transform(class_index)
wxchan's avatar
wxchan committed
1333

1334
1335
    predict.__doc__ = LGBMModel.predict.__doc__

1336
1337
    def predict_proba(
        self,
1338
        X: _LGBM_ScikitMatrixLike,
1339
1340
1341
1342
1343
1344
        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,
1345
        **kwargs: Any,
1346
    ):
1347
        """Docstring is set after definition, using a template."""
1348
1349
1350
1351
1352
1353
1354
1355
        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,
1356
            **kwargs,
1357
        )
1358
        if callable(self._objective) and not (raw_score or pred_leaf or pred_contrib):
1359
1360
1361
1362
1363
            _log_warning(
                "Cannot compute class probabilities or labels "
                "due to the usage of customized objective function.\n"
                "Returning raw scores instead."
            )
1364
            return result
1365
        elif self.__is_multiclass or raw_score or pred_leaf or pred_contrib:  # type: ignore [operator]
1366
            return result
wxchan's avatar
wxchan committed
1367
        else:
1368
            return np.vstack((1.0 - result, result)).transpose()
1369

1370
1371
    predict_proba.__doc__ = _lgbmmodel_doc_predict.format(
        description="Return the predicted probability for each class for each sample.",
1372
        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]",
1373
        output_name="predicted_probability",
1374
        predicted_result_shape="array-like of shape = [n_samples] or shape = [n_samples, n_classes]",
1375
        X_leaves_shape="array-like of shape = [n_samples, n_trees] or shape = [n_samples, n_trees * n_classes]",
1376
        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",
1377
1378
    )

1379
    @property
1380
    def classes_(self) -> np.ndarray:
1381
        """:obj:`array` of shape = [n_classes]: The class label array."""
1382
        if not self.__sklearn_is_fitted__():
1383
            raise LGBMNotFittedError("No classes found. Need to call fit beforehand.")
1384
        return self._classes  # type: ignore[return-value]
1385
1386

    @property
1387
    def n_classes_(self) -> int:
1388
        """:obj:`int`: The number of classes."""
1389
        if not self.__sklearn_is_fitted__():
1390
            raise LGBMNotFittedError("No classes found. Need to call fit beforehand.")
1391
        return self._n_classes
wxchan's avatar
wxchan committed
1392

1393
1394
1395
1396
1397
    @property
    def __is_multiclass(self) -> bool:
        """:obj:`bool`:  Indicator of whether the classifier is used for multiclass."""
        return self._n_classes > 2 or (isinstance(self._objective, str) and self._objective in _MULTICLASS_OBJECTIVES)

wxchan's avatar
wxchan committed
1398

wxchan's avatar
wxchan committed
1399
class LGBMRanker(LGBMModel):
1400
1401
1402
1403
1404
1405
1406
1407
    """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
1408

1409
    def fit(  # type: ignore[override]
1410
        self,
1411
1412
        X: _LGBM_ScikitMatrixLike,
        y: _LGBM_LabelType,
1413
1414
        sample_weight: Optional[_LGBM_WeightType] = None,
        init_score: Optional[_LGBM_InitScoreType] = None,
1415
        group: Optional[_LGBM_GroupType] = None,
1416
        eval_set: Optional[List[_LGBM_ScikitValidSet]] = None,
1417
        eval_names: Optional[List[str]] = None,
1418
1419
1420
        eval_sample_weight: Optional[List[_LGBM_WeightType]] = None,
        eval_init_score: Optional[List[_LGBM_InitScoreType]] = None,
        eval_group: Optional[List[_LGBM_GroupType]] = None,
1421
        eval_metric: Optional[_LGBM_ScikitEvalMetricType] = None,
1422
        eval_at: Union[List[int], Tuple[int, ...]] = (1, 2, 3, 4, 5),
1423
1424
        feature_name: _LGBM_FeatureNameConfiguration = "auto",
        categorical_feature: _LGBM_CategoricalFeatureConfiguration = "auto",
1425
        callbacks: Optional[List[Callable]] = None,
1426
        init_model: Optional[Union[str, Path, Booster, LGBMModel]] = None,
1427
    ) -> "LGBMRanker":
1428
        """Docstring is inherited from the LGBMModel."""
1429
        # check group data
Guolin Ke's avatar
Guolin Ke committed
1430
        if group is None:
1431
            raise ValueError("Should set group for ranking task")
wxchan's avatar
wxchan committed
1432
1433

        if eval_set is not None:
Guolin Ke's avatar
Guolin Ke committed
1434
            if eval_group is None:
1435
                raise ValueError("Eval_group cannot be None when eval_set is not None")
Guolin Ke's avatar
Guolin Ke committed
1436
            elif len(eval_group) != len(eval_set):
1437
                raise ValueError("Length of eval_group should be equal to eval_set")
1438
1439
1440
1441
1442
1443
1444
1445
1446
1447
            elif (
                isinstance(eval_group, dict)
                and any(i not in eval_group or eval_group[i] is None for i in range(len(eval_group)))
                or isinstance(eval_group, list)
                and any(group is None for group in eval_group)
            ):
                raise ValueError(
                    "Should set group for all eval datasets for ranking task; "
                    "if you use dict, the index should start from 0"
                )
1448

1449
        self._eval_at = eval_at
1450
1451
1452
1453
1454
1455
1456
1457
1458
1459
1460
1461
1462
1463
1464
        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,
1465
            init_model=init_model,
1466
        )
wxchan's avatar
wxchan committed
1467
        return self
1468

1469
    _base_doc = LGBMModel.fit.__doc__.replace("self : LGBMModel", "self : LGBMRanker")  # type: ignore
1470
1471
1472
1473
    fit.__doc__ = (
        _base_doc[: _base_doc.find("eval_class_weight :")]  # type: ignore
        + _base_doc[_base_doc.find("eval_init_score :") :]
    )  # type: ignore
1474
    _base_doc = fit.__doc__
1475
    _before_feature_name, _feature_name, _after_feature_name = _base_doc.partition("feature_name :")
1476
    fit.__doc__ = f"""{_before_feature_name}eval_at : list or tuple of int, optional (default=(1, 2, 3, 4, 5))
1477
        The evaluation positions of the specified metric.
1478
    {_feature_name}{_after_feature_name}"""