Alb*_*nto 3 python regression neural-network scikit-learn
我目前有一个包含变量和观测值的数据集。我想预测一个变量(需求),它是一个连续的变量,因此我需要使用回归模型。我尝试过Linear Regression,并使用R2大约的指标对其进行了评估0.85。我想用其他模型评估其性能,其中一个是NNs。我相信神经网络更适合分类之类的其他任务,不过我想尝试一下。
我scikit-learn之所以决定使用它,主要是因为它提供了两种模型(线性回归和多层感知器),但事实是,R2与线性回归相比,该度量标准太过复杂了。因此,我得出的结论是我缺少许多重要的配置。在下面,您可以查看我的代码以及数据的来源。
我的数据有以下几列,只demand(这是我的标签)population,gdp,day和year是数值连续,其余都是绝对的。
['demand','holy','gdp','population', 'day','year', 'f0', 'f1', 'f2', 'f3', 'f4','f5', 'f6', 'f7', 'f8', 'f9', 'f10', 'f11', 'f12', 'f13', 'f14', 'f15', 'f16', 'f17', 'f18', 'f19', 'f20', 'f21', 'f22', 'f23', 'g0', 'g1', 'g2', 'g3', 'g4', 'g5', 'g6', 'g7', 'g8', 'g9', 'g10', 'g11']
Run Code Online (Sandbox Code Playgroud)
这是我实际上所做的,我删除了一些输出。
import pandas as pd
import numpy as np
import math as math
from sklearn.linear_model import LinearRegression
from sklearn.neural_network import MLPRegressor
from sklearn.metrics import r2_score
training_data, validation_data = np.split(data.sample(frac=1), [int(.8*len(data))])
linear_model = LinearRegression().fit(training_data[[c for c in data.columns if c != "demand"]], training_data[["demand"]])
validation_data_predictions = linear_model.predict(validation_data[[c for c in training_data.columns if c != "demand"]])
validation_predictions_pd = pd.DataFrame(data=validation_data_predictions,
index=validation_data.index.values,
columns=["prediction"])
# join both pandas
result_df = validation_data.join(validation_predictions_pd, how="inner")
r2_error = r2_score(y_true=result_df[["demand"]], y_pred=result_df[["prediction"]], multioutput="uniform_average")
print(r2_error) # outputs 0.85
# NN section
clf = MLPRegressor(hidden_layer_sizes=(10,), max_iter=100000)
neural_model = clf.fit(training_data[[c for c in training_data.columns if c != "demand"]], training_data[["demand"]])
validation_data_predictions = neural_model.predict(validation_data[[c for c in training_data.columns if c != "demand"]])
validation_predictions_pd = pd.DataFrame(data=validation_data_predictions,
index=validation_data.index.values,
columns=["prediction"])
result_df = validation_data.join(validation_predictions_pd, how="inner")
r2_error = r2_score(y_true=result_df[["demand"]], y_pred=result_df[["prediction"]], multioutput="uniform_average")
print(r2_error) # outputs 0.23
Run Code Online (Sandbox Code Playgroud)
因此,如您所见,NN的性能非常差。而且我认为它的性能可以提高,有什么暗示吗?
| 归档时间: |
|
| 查看次数: |
11967 次 |
| 最近记录: |