Есть ли утечка данных в моем коде (кривая RO C дает оценку 1,00 AU C)? - PullRequest
Новая
       58

Есть ли утечка данных в моем коде (кривая RO C дает оценку 1,00 AU C)?

0 голосов
/ 03 марта 2020

Люди говорили мне, что произошла утечка данных при разделении моих данных на набор для тестирования и обучения. Моя цель - построить кривую рабочих характеристик c приемника с использованием логарифмической регрессии как исходных данных (после стандартизации), так и модели анализа главных компонентов.

Однако кривая рабочих характеристик приемника c просто дает мне 1,00 область под характеристикой работы приемника c показатель кривой для каждого графика.

Редактировать: Как отметил Дэвид, мне удалось получить показатель RO C AU C <1 для модель спс. Это половина проблемы решена! Модель log reg в обычном наборе данных по-прежнему дает мне идеальный балл 1,00. </p>

Вот код для назначения моих функций и меток .:

# Import wine dataset.
df_1 = pd.read_csv ('./wine.data', delimiter =',', header = None, nrows=200)

# Feature Selection (dropping strongly correlated features):
df_1 = df_1.drop(df_1.columns[7], axis='columns')
df_1

# Separating out the features, columns 1 to 13.
x = df_1.iloc[:, 1:13].values 

# Separating out the target, column 0, with classes 1, 2, and 3.
y = df_1.iloc[:, 0].values 

# Standardizing the features (centering and scaling).
x = StandardScaler().fit_transform(x)

# Dataset is split into training set (70%) and testing set (30%).
X_train,X_test,y_train,y_test=train_test_split(x,y,test_size=0.30,random_state=123)

#PCA is performed for 2 components
pca_2 = PCA(n_components=2)
X_train_pca_2 = pca_2.fit_transform(X_train)
X_test_pca_2 = pca_2.fit_transform(X_test)

А вот остальные мой код, на случай, если утечка данных на самом деле не проблема: 

# Create logistic regression function for one vs all classifier.
logreg = LogisticRegression(random_state=0, multi_class='ovr')

# Train the model using the original pre-processed dataset.
model = logreg.fit(X_train, y_train)

# Predict the class of each wine in the testing data.
y_pred = model.predict(X_test)
y_pred

array([3, 2, 3, 2, 2, 3, 1, 3, 3, 2, 3, 3, 3, 1, 1, 3, 2, 2, 1, 2, 3, 3,
       3, 3, 2, 3, 3, 2, 1, 1, 1, 1, 2, 2, 3, 2, 3, 1, 2, 2, 3, 3, 1, 1,
       2, 1, 1, 2, 1, 2, 2, 3, 3, 2], dtype=int64)

# Train a second model using the PCA.
model_PCA = logreg.fit(X_train_pca_2, y_train)

# Predict the class of each wine in the testing data.
y_pred_PCA = model_PCA.predict(X_test_pca_2)
y_pred_PCA

array([3, 2, 2, 2, 2, 3, 1, 2, 3, 2, 3, 3, 3, 1, 1, 3, 2, 2, 1, 2, 2, 3,
       3, 2, 2, 3, 2, 1, 1, 1, 1, 1, 2, 1, 3, 2, 2, 1, 2, 2, 3, 3, 1, 1,
       2, 1, 1, 2, 1, 1, 2, 2, 3, 2], dtype=int64)

# Binarize the labels
y_train_binary = label_binarize(y_train, classes=[1, 2, 3])
y_test_binary = label_binarize(y_test, classes=[1, 2, 3])
n_classes = y_train_binary.shape[1]

# Score for One vs Rest Classifier
y_score = model.fit(X_train, y_train).decision_function(X_test)

# Score for PCA model
y_score_PCA = model_PCA.fit(X_train_pca_2, y_train).decision_function(X_test_pca_2)

fpr1 = dict()
tpr1 = dict()
roc_auc1 = dict()
for i in range(n_classes):
    fpr1[i], tpr1[i], _ = roc_curve(y_test_binary[:, i], y_score[:, i])
    roc_auc1[i] = auc(fpr1[i], tpr1[i])

fpr2 = dict()
tpr2 = dict()
roc_auc2 = dict()
for j in range(n_classes):
    fpr2[j], tpr2[j], _ = roc_curve(y_test_binary[:, j], y_score_PCA[:, j])
    roc_auc2[j] = auc(fpr2[j], tpr2[j])

# Plotting ROC Curve of all 3 classes using logistic regression on the orginal data set and the PCA model
plt.figure()
lw = 1
plt.plot(fpr1[0], tpr1[0], color='purple',
         lw=lw, label='LogReg Class 1 (area = %0.2f)' % roc_auc1[0])
plt.plot(fpr1[1], tpr1[1], color='blue',
         lw=lw, label='LogReg Class 2 (area = %0.2f)' % roc_auc1[1])
plt.plot(fpr1[2], tpr1[2], color='aqua',
         lw=lw, label='LogReg Class 3 (area = %0.2f)' % roc_auc1[2])
plt.plot(fpr2[0], tpr2[0], color='orange',
         lw=lw, label='PCA Class 1 (area = %0.2f)' % roc_auc2[0])
plt.plot(fpr2[1], tpr2[1], color='green',
         lw=lw, label='PCA Class 2 (area = %0.2f)' % roc_auc2[1])
plt.plot(fpr2[2], tpr2[2], color='red',
         lw=lw, label='PCA Class 3 (area = %0.2f)' % roc_auc2[2])
plt.plot([0, 1], [0, 1], color='navy', lw=lw, linestyle='--')
plt.xlim([0.0, 1.0])
plt.ylim([0.0, 1.05])
plt.xlabel('False Positive Rate')
plt.ylabel('True Positive Rate')
plt.title('ROC Curve')
plt.legend(loc="lower right")
plt.show()

1 Ответ

0 голосов
/ 03 марта 2020

Возможно, вы действительно получаете идеальные результаты. Я не имел дело с набором данных, поэтому я не знаю, насколько это возможно с помощью функции регрессии logisti c. Однако в вашем коде есть один существенный недостаток: вы перезаписываете свою модель pca моделью, обученной во всем пространстве признаков. см .: 

y_score_PCA = model_PCA.fit(X_train, y_train).decision_function(X_test)

Вы используете весь вектор X_train, а не то, что я ожидал бы быть правильным, 

y_score_PCA = model_PCA.fit(X_train_pca_2, y_train).decision_function(X_test_pca_2)

или, в качестве альтернативы, не обновляйте модель 

y_score_PCA = model_PCA.decision_function(X_test_pca_2)

Это делает ваши оценки одинаковыми для обеих моделей. Затем я ожидаю, что ваш AU C для модели PCA будет <1,0. Если это не так, вероятно, возникнет еще больше проблем, поскольку у вас, похоже, разные прогнозы для каждой модели. </p>

...