Dav*_*ora 3 python function machine-learning scikit-learn logistic-regression
我构建了一个显示单个模型的一些评估指标的函数,现在我想将此函数应用于我估计的模型池。
旧函数的输入是:
OldFunction(code: str, x, X_train: np.array, X_test: np.array, X:pd.DataFrame)
Run Code Online (Sandbox Code Playgroud)
在哪里:
code 是一个字符串,用于创建数据框的列名
x 是模型名称
X_train 和 X_test 是数据分割器的 np.arrays
X 是整个数据的数据框
为了估计模型池的指标,我尝试通过在函数中添加循环并将模型放入列表中来修改函数。
但这不起作用。
出现问题是因为我无法迭代模型列表,那么我有什么选择呢?你有什么想法吗?
我将新功能留在下面。
import numpy as np
import pandas as pd
from sklearn.metrics import roc_curve, auc
from sklearn.metrics import accuracy_score, recall_score, precision_score
from sklearn.model_selection import cross_val_score
def displaymetrics(code: list, models: list, X_train: np.array, X_test: np.array, X: pd.DataFrame):
for i in models:
y_score = models[i].fit(X_train, y_train).decision_function(X_test)
fpr, tpr, _ = roc_curve(y_test, y_score)
roc_auc = auc(fpr, tpr)
# Traditional Scores
y_pred = pd.DataFrame(model[i].predict(X_train)).reset_index(drop=True)
Recall_Train,Precision_Train, Accuracy_Train = recall_score(y_train, y_pred), precision_score(y_train, y_pred), accuracy_score(y_train, y_pred)
y_pred = pd.DataFrame(model[i].predict(X_test)).reset_index(drop=True)
Recall_Test = recall_score(y_test, y_pred)
Precision_Test = precision_score(y_test, y_pred)
Accuracy_Test = accuracy_score(y_test, y_pred)
#Cross Validation
cv_au = cross_val_score(models[i], X_test, y_test, cv=30, scoring='roc_auc')
cv_f1 = cross_val_score(models[i], X_test, y_test, cv=30, scoring='f1')
cv_pr = cross_val_score(models[i], X_test, y_test, cv=30, scoring='precision')
cv_re = cross_val_score(models[i], X_test, y_test, cv=30, scoring='recall')
cv_ac = cross_val_score(models[i], X_test, y_test, cv=30, scoring='accuracy')
cv_ba = cross_val_score(models[i], X_test, y_test, cv=30, scoring='balanced_accuracy')
cv_au_m, cv_au_std = cv_au.mean() , cv_au.std()
cv_f1_m, cv_f1_std = cv_f1.mean() , cv_f1.std()
cv_pr_m, cv_pr_std = cv_pr.mean() , cv_pr.std()
cv_re_m, cv_re_std= cv_re.mean() , cv_re.std()
cv_ac_m, cv_ac_std = cv_ac.mean() , cv_ac.std()
cv_ba_m, cv_ba_std= cv_ba.mean() , cv_ba.std()
cv_au, cv_f1, cv_pr = (cv_au_m, cv_au_std), (cv_f1_m, cv_f1_std), (cv_pr_m, cv_pr_std)
cv_re, cv_ac, cv_ba = (cv_re_m, cv_re_std), (cv_ac_m, cv_ac_std), (cv_ba_m, cv_ba_std)
tuples = [cv_au, cv_f1, cv_pr, cv_re, cv_ac, cv_ba]
tuplas = [0]*len(tuples)
for i in range(len(tuples)):
tuplas[i] = [round(x,4) for x in tuples[i]]
results = pd.DataFrame()
results['Metrics'] = ['roc_auc', 'Accuracy_Train', 'Precision_Train', 'Recall_Train', 'Accuracy_Test',
'Precision_Test','Recall_Test', 'cv_roc-auc (mean, std)', 'cv_f1score(mean, std)',
'cv_precision (mean, std)', 'cv_recall (mean, std)', 'cv_accuracy (mean, std)',
'cv_bal_accuracy (mean, std)']
results.set_index(['Metrics'], inplace=True)
results['Model_'+code[i]] = [roc_auc, Accuracy_Train, Precision_Train, Recall_Train, Accuracy_Test,
Precision_Test, Recall_Test, tuplas[0], tuplas[1], tuplas[2], tuplas[3],
tuplas[4], tuplas[5]]
return results
Run Code Online (Sandbox Code Playgroud)
输出应该是一个数据框,其中每列代表每个模型,行代表指标。
您应该使用字典而不是列表,如下例所示:
dict_classifiers = {
"Logreg": LogisticRegression(solver='lbfgs'),
"NN": KNeighborsClassifier(),
"LinearSVM": SVC(probability=True, kernel='linear'), #class_weight='balanced'
"GBC": GradientBoostingClassifier(),
"DT": tree.DecisionTreeClassifier(),
"RF": RandomForestClassifier(),
"NB": GaussianNB(),
}
Run Code Online (Sandbox Code Playgroud)
然后使用,例如:
for model, model_instantiation in dict_classifiers.iteritems():
y_score = model_instantiation.fit(X_train, y_train).decision_function(X_test)
...
Run Code Online (Sandbox Code Playgroud)
希望对您有所帮助,让我知道您的进展如何!
| 归档时间: |
|
| 查看次数: |
8537 次 |
| 最近记录: |