【Python | 机器学习】Python中进行特征重要性分析的9个常用方法（含源代码）

桃扇骨 2023-10-16 07:12 39阅读 0赞

#### 文章目录 ####

*  一、为什么特征重要性分析很重要?
 *  二、特征重要性分析方法
 *   *  2.1 排列重要性 PermutationImportance
     *  2.2 内置特征重要性(coef\_或feature\_importances\_)
     *  2.3 Leave-one-out
     *  2.4 相关性分析
     *  2.5 递归特征消除 Recursive Feature Elimination
     *  2.6 XGBoost特性重要性
     *  2.7 主成分分析 PCA
     *  2.8 方差分析 ANOVA
     *  2.9 卡方检验
 *  三、为什么不同的方法会检测到不同的特征?
 *  四、选择特征重要性分析方法的一些最佳实践

特征重要性分析用于了解每个特征(变量或输入)对于做出预测的有用性或价值。目标是确定对模型输出影响最大的最重要的特征，它是机器学习中经常使用的一种方法。

## 一、为什么特征重要性分析很重要? ##

特征重要性分析在数据科学和机器学习中扮演着重要的角色，具有以下重要性：

理解数据：特征重要性分析有助于理解数据集中每个特征对目标变量的影响程度。这有助于识别哪些特征在解决问题时可能是关键的，哪些可能是次要的。

特征选择：在建模过程中，特征重要性分析可用于选择最相关的特征，以减少模型的复杂性和提高模型的泛化能力。通过消除不相关或冗余的特征，可以改善模型的性能，减少过拟合的风险。

模型解释性：特征重要性分析有助于解释模型的预测结果。它可以告诉我们哪些特征对于模型的某个预测是最重要的，这对于业务决策和解释模型的工作原理非常有帮助。

异常检测：通过分析特征重要性，可以识别出与预期模式不符的异常数据点。这些异常数据可能是数据质量问题或需要特殊处理的情况，例如欺诈检测或故障检测。

特征工程：特征重要性分析可以为特征工程提供指导。它可以帮助数据科学家和机器学习工程师选择最相关的特征，或者启发他们创建新的特征来提高模型性能。

优化模型：了解特征的重要性可以帮助优化模型的超参数选择，例如决策树中的树的深度或随机森林中的树的数量。这可以提高模型的性能并减少计算成本。

监测模型稳定性：特征重要性的变化可以用于监测模型的稳定性。如果模型的重要特征发生变化，可能意味着数据分布或问题本身发生了变化，需要进一步的调查。

如果有一个包含数十个甚至数百个特征的数据集，每个特征都可能对你的机器学习模型的性能有所贡献。但是并不是所有的特征都是一样的。有些可能是冗余的或不相关的，这会增加建模的复杂性并可能导致过拟合。

特征重要性分析可以识别并关注最具信息量的特征，从而带来以下几个优势:

*  改进的模型性能
 *  减少过度拟合
 *  更快的训练和推理
 *  增强的可解释性

下面我们深入了解在Python中的一些特性重要性分析的方法。

## 二、特征重要性分析方法 ##

### 2.1 排列重要性 PermutationImportance ###

该方法会随机排列每个特征的值，然后监控模型性能下降的程度。如果获得了更大的下降意味着特征更重要：

from sklearn.datasets import load_breast_cancer
     from sklearn.ensemble import RandomForestClassifier
     from sklearn.inspection import permutation_importance
     from sklearn.model_selection import train_test_split
     import matplotlib.pyplot as plt
     
     cancer = load_breast_cancer()
     
     X_train, X_test, y_train, y_test = train_test_split(cancer.data, cancer.target, random_state=1)
     
     rf = RandomForestClassifier(n_estimators=100, random_state=1)
     rf.fit(X_train, y_train)
     
     baseline = rf.score(X_test, y_test)
     result = permutation_importance(rf, X_test, y_test, n_repeats=10, random_state=1, scoring='accuracy')
     
     importances = result.importances_mean
     
     # Visualize permutation importances
     plt.bar(range(len(importances)), importances)
     plt.xlabel('Feature Index')
     plt.ylabel('Permutation Importance')
     plt.show()

![在这里插入图片描述][41946dec8f41492e8dbbe61a87afc7da.png]

### 2.2 内置特征重要性(coef\_或feature\_importances\_) ###

一些模型，如线性回归和随机森林，可以直接输出特征重要性分数。这些显示了每个特征对最终预测的贡献。

from sklearn.datasets import load_breast_cancer
     from sklearn.ensemble import RandomForestClassifier
     
     X, y = load_breast_cancer(return_X_y=True)
     
     rf = RandomForestClassifier(n_estimators=100, random_state=1)
     rf.fit(X, y)
     
     importances = rf.feature_importances_
     
     # Plot importances
     plt.bar(range(X.shape[1]), importances)
     plt.xlabel('Feature Index')
     plt.ylabel('Feature Importance')
     plt.show()

![在这里插入图片描述][d8878eb848a94116ad915c509cf7968b.png]

### 2.3 Leave-one-out ###

迭代地每次删除一个特征并评估准确性。

from sklearn.datasets import load_breast_cancer
     from sklearn.model_selection import train_test_split
     from sklearn.ensemble import RandomForestClassifier
     from sklearn.metrics import accuracy_score
     import matplotlib.pyplot as plt
     import numpy as np
     
     # Load sample data
     X, y = load_breast_cancer(return_X_y=True)
     
     # Split data into train and test sets
     X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.3, random_state=1)
     
     # Train a random forest model
     rf = RandomForestClassifier(n_estimators=100, random_state=1)
     rf.fit(X_train, y_train)
     
     # Get baseline accuracy on test data
     base_acc = accuracy_score(y_test, rf.predict(X_test))
     
     # Initialize empty list to store importances
     importances = []
     
     # Iterate over all columns and remove one at a time
     for i in range(X_train.shape[1]):
        X_temp = np.delete(X_train, i, axis=1)
        rf.fit(X_temp, y_train)
        acc = accuracy_score(y_test, rf.predict(np.delete(X_test, i, axis=1)))
        importances.append(base_acc - acc)
         
     # Plot importance scores    
     plt.bar(range(len(importances)), importances)
     plt.show()

![在这里插入图片描述][1e30939092e4436c95d91562dd2e1162.png]

### 2.4 相关性分析 ###

计算各特征与目标变量之间的相关性。相关性越高的特征越重要。

import pandas as pd
     from sklearn.datasets import load_breast_cancer
     
     X, y = load_breast_cancer(return_X_y=True)
     df = pd.DataFrame(X, columns=range(30))
     df['y'] = y
     
     correlations = df.corrwith(df.y).abs()
     correlations.sort_values(ascending=False, inplace=True)
     
     correlations.plot.bar()

![在这里插入图片描述][6a0fef5e08fc4f3ab4b68837cbc8394d.png]

### 2.5 递归特征消除 Recursive Feature Elimination ###

递归地删除特征并查看它如何影响模型性能。删除时会导致更大下降的特征更重要。

from sklearn.ensemble import RandomForestClassifier
     from sklearn.feature_selection import RFE
     import pandas as pd
     from sklearn.datasets import load_breast_cancer
     import matplotlib.pyplot as plt
     
     X, y = load_breast_cancer(return_X_y=True)
     df = pd.DataFrame(X, columns=range(30))
     df['y'] = y
     
     rf = RandomForestClassifier()
     
     rfe = RFE(rf, n_features_to_select=10)
     rfe.fit(X, y)
     
     print(rfe.ranking_)

输出结果为：

输出为[6 4 11 12 7 11 18 21 8 16 10 3 15 14 19 17 20 13 11 11 12 9 11 5 11]

### 2.6 XGBoost特性重要性 ###

计算一个特性用于跨所有树拆分数据的次数。更多的分裂意味着更重要。

import xgboost as xgb
     import pandas as pd
     from sklearn.datasets import load_breast_cancer
     import matplotlib.pyplot as plt
     
     X, y = load_breast_cancer(return_X_y=True)
     df = pd.DataFrame(X, columns=range(30))
     df['y'] = y
     
     model = xgb.XGBClassifier()
     model.fit(X, y)
     
     importances = model.feature_importances_
     importances = pd.Series(importances, index=range(X.shape[1]))
     importances.plot.bar()

![在这里插入图片描述][aaa15a02e7ff433cbbb088e30e8ed86f.png]

### 2.7 主成分分析 PCA ###

对特征进行主成分分析，并查看每个主成分的解释方差比。在前几个组件上具有较高负载的特性更为重要。

from sklearn.decomposition import PCA
     import pandas as pd
     from sklearn.datasets import load_breast_cancer
     import matplotlib.pyplot as plt
     
     X, y = load_breast_cancer(return_X_y=True)
     df = pd.DataFrame(X, columns=range(30))
     df['y'] = y
     
     pca = PCA()
     pca.fit(X)
     
     plt.bar(range(pca.n_components_), pca.explained_variance_ratio_)
     plt.xlabel('PCA components')
     plt.ylabel('Explained Variance')

![在这里插入图片描述][d7be946b62ca45c08108cea2169cf62c.png]

### 2.8 方差分析 ANOVA ###

使用f\_classif()获得每个特征的方差分析f值。f值越高，表明特征与目标的相关性越强。

from sklearn.feature_selection import f_classif
     import pandas as pd
     from sklearn.datasets import load_breast_cancer
     import matplotlib.pyplot as plt
     
     X, y = load_breast_cancer(return_X_y=True)
     df = pd.DataFrame(X, columns=range(30))
     df['y'] = y
     
     fval = f_classif(X, y)
     fval = pd.Series(fval[0], index=range(X.shape[1]))
     fval.plot.bar()

![在这里插入图片描述][52b9bcff2ca14b21ae215e035b79013a.png]

### 2.9 卡方检验 ###

使用chi2()获得每个特征的卡方统计信息。得分越高的特征越有可能独立于目标。

from sklearn.feature_selection import chi2
     import pandas as pd
     from sklearn.datasets import load_breast_cancer
     import matplotlib.pyplot as plt
     
     X, y = load_breast_cancer(return_X_y=True)
     df = pd.DataFrame(X, columns=range(30))
     df['y'] = y
     
     chi_scores = chi2(X, y)
     chi_scores = pd.Series(chi_scores[0], index=range(X.shape[1]))
     chi_scores.plot.bar()

![在这里插入图片描述][9551d951a4bc4569976d158edb721a04.png]

## 三、为什么不同的方法会检测到不同的特征? ##

不同的特征重要性方法有时可以识别出不同的特征是最重要的，这是因为：

1、他们用不同的方式衡量重要性:

有的使用不同特特征进行预测，监控精度下降

像XGBOOST或者回国模型使用内置重要性来进行特征的重要性排列

而PCA着眼于方差解释

2、不同模型有不同模型的方法：

线性模型倾向于线性关系、树模型倾向于接近根的特征

3、交互作用:

有的方法可以获取特征之间的相互左右，而有一些则不行，这就会导致结果的差异

3、不稳定:

使用不同的数据子集，重要性值可能在同一方法的不同运行中有所不同，这是因为数据差异决定的

4、Hyperparameters:

通过调整超参数，如PCA组件或树深度，也会影响结果

所以不同的假设、偏差、数据处理和方法的可变性意味着它们并不总是在最重要的特征上保持一致。

## 四、选择特征重要性分析方法的一些最佳实践 ##

*  尝试多种方法以获得更健壮的视图
 *  聚合结果的集成方法
 *  更多地关注相对顺序，而不是绝对值
 *  差异并不一定意味着有问题，检查差异的原因会对数据和模型有更深入的了解

[41946dec8f41492e8dbbe61a87afc7da.png]: https://img-blog.csdnimg.cn/41946dec8f41492e8dbbe61a87afc7da.png
[d8878eb848a94116ad915c509cf7968b.png]: https://img-blog.csdnimg.cn/d8878eb848a94116ad915c509cf7968b.png
[1e30939092e4436c95d91562dd2e1162.png]: https://img-blog.csdnimg.cn/1e30939092e4436c95d91562dd2e1162.png
[6a0fef5e08fc4f3ab4b68837cbc8394d.png]: https://img-blog.csdnimg.cn/6a0fef5e08fc4f3ab4b68837cbc8394d.png
[aaa15a02e7ff433cbbb088e30e8ed86f.png]: https://img-blog.csdnimg.cn/aaa15a02e7ff433cbbb088e30e8ed86f.png
[d7be946b62ca45c08108cea2169cf62c.png]: https://img-blog.csdnimg.cn/d7be946b62ca45c08108cea2169cf62c.png
[52b9bcff2ca14b21ae215e035b79013a.png]: https://img-blog.csdnimg.cn/52b9bcff2ca14b21ae215e035b79013a.png
[9551d951a4bc4569976d158edb721a04.png]: https://img-blog.csdnimg.cn/9551d951a4bc4569976d158edb721a04.png