使用Scikit-Learn计算平均AUC的差异

我有以下代码:

from sklearn import svm, datasetsfrom sklearn.metrics import roc_curve, aucfrom sklearn.model_selection import StratifiedKFold, cross_val_scoreimport numpy as npfrom scipy import interpseed = 7np.random.seed(seed)iris = datasets.load_iris()X = iris.datay = iris.targetX, y = X[y != 2], y[y != 2]n_samples, n_features = X.shape# Add noisy featuresrandom_state = np.random.RandomState(0)X = np.c_[X, random_state.randn(n_samples, 200 * n_features)]cv = StratifiedKFold(n_splits=10)classifier = svm.SVC(kernel='linear', probability=True, random_state=seed)mean_tpr = 0.0mean_fpr = np.linspace(0, 1, 100)i= 0for train, test in cv.split(X, y):    probas_ = classifier.fit(X[train], y[train]).predict_proba(X[test])    # Compute ROC curve and area the curve    fpr, tpr, thresholds = roc_curve(y[test], probas_[:, 1])    mean_tpr += interp(mean_fpr, fpr, tpr)    mean_tpr[0] = 0.0    roc_auc = auc(fpr, tpr)    i += 1mean_tpr /= cv.get_n_splits(X, y)mean_tpr[-1] = 1.0mean_auc_1 = auc(mean_fpr, mean_tpr)print "#--- Method 1 to calculate mean AUC ---"print mean_auc_1print "#--- Method 2 to calculate mean AUC ---"results = cross_val_score(classifier, X, y, cv=cv)mean_auc_2 =  "{:.3f}".format(results.mean())print mean_auc_2

它产生了以下结果:

#--- Method 1 to calculate mean AUC ---0.801818181818#--- Method 2 to calculate mean AUC ---0.700

计算平均AUC的第一种方法是通过循环,如这个Scikit教程所建议的。第二种方法使用Scikit内置的cross_val_score()方法来计算平均AUC。

我的问题是,为什么会有差异?我应该相信哪个平均AUC值?如何修改第二种方法,使其结果与第一种方法相同?

我使用的是以下版本的Scikit-Learn:

In [442]: sklearn.__version__Out[442]: '0.18'

回答:

你的第二个例子中没有AUC计算。你应该添加一个自定义评分函数。请查看cross_val_score的API。

你只是在计算平均准确率。这通常是用于分类器的标准评分函数。你可以在文档中查看svm的标准评分函数

像这样

cross_val_score(classifier, X, y, cv=cv, scoring='roc_auc')

应该可以工作

Related Posts

使用LSTM在Python中预测未来值

这段代码可以预测指定股票的当前日期之前的值,但不能预测…

如何在gensim的word2vec模型中查找双词组的相似性

我有一个word2vec模型,假设我使用的是googl…

dask_xgboost.predict 可以工作但无法显示 – 数据必须是一维的

我试图使用 XGBoost 创建模型。 看起来我成功地…

ML Tuning – Cross Validation in Spark

我在https://spark.apache.org/…

如何在React JS中使用fetch从REST API获取预测

我正在开发一个应用程序,其中Flask REST AP…

如何分析ML.NET中多类分类预测得分数组?

我在ML.NET中创建了一个多类分类项目。该项目可以对…

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注