Python实现贝叶斯优化器(Bayes_opt)优化深度森林(Deep Forest)回归模型(deepforest回归算法)项目实战

说明：这是一个机器学习实战项目（附带数据+代码+文档+视频讲解），如需数据+代码+文档+视频讲解可以直接到文章最后获取。

1.项目背景

深度森林是一种基于集成学习的方法，它受到深度学习成功的启发，但不依赖于梯度下降和反向传播等复杂机制。深度森林通过堆叠多个决策树模型（通常是随机森林）并在层级结构中结合它们的输出来构建最终的模型。这种方法可以在许多情况下提供与深度神经网络相当的性能，同时减少了对计算资源的需求和训练时间。

贝叶斯优化是一种全局优化技术，主要用于黑盒函数的优化问题，这些函数往往计算成本较高且没有解析形式。贝叶斯优化通过构建代理模型（如高斯过程）来近似目标函数，并通过获取代理模型的信息来指导下一步搜索的位置。这种技术特别适用于超参数调优，因为它能够有效地探索参数空间，找到最优解。

在本项目中，我们的目标是应用贝叶斯优化来自动寻找深度森林回归模型的最佳超参数组合。通过这种方式，我们可以提高模型的预测准确性和泛化能力，同时减少手动调整参数所需的时间和精力。

本项目通过Python实现贝叶斯优化器(Bayes_opt)优化深度森林(Deep Forest)回归模型(deepforest回归算法)项目实战。   

2.数据获取

本次建模数据来源于网络(本项目撰写人整理而成)，数据项统计如下：

编号	变量名称	描述
1	x1
2	x2
3	x3
4	x4
5	x5
6	x6
7	x7
8	x8
9	x9
10	x10
11	y	因变量

数据详情如下(部分展示)：

3.数据预处理

3.1 用Pandas工具查看数据

使用Pandas工具的head()方法查看前五行数据：

关键代码：

3.2数据缺失查看

使用Pandas工具的info()方法查看数据信息：

从上图可以看到，总共有11个变量，数据中无缺失值，共2000条数据。

关键代码： 

3.3数据描述性统计

通过Pandas工具的describe()方法来查看数据的平均值、标准差、最小值、分位数、最大值。

关键代码如下：  

4.探索性数据分析

4.1 y变量分布直方图

用Matplotlib工具的hist()方法绘制直方图：

4.2 相关性分析

从上图中可以看到，数值越大相关性越强，正值是正相关、负值是负相关。

5.特征工程

5.1 建立特征数据和标签数据

关键代码如下：

5.2 数据集拆分

通过train_test_split()方法按照80%训练集、20%测试集进行划分，关键代码如下：

6.构建贝叶斯优化算法优化深度森林回归模型 

主要使用贝叶斯优化算法优化深度森林回归算法，用于目标回归。 

6.1 构建调优模型

编号	模型名称	调优参数
1	深度森林回归模型	n_estimators
2		n_trees
3		max_layers

6.2 最优参数展示

最优参数结果展示：

6.3 最优参数构建模型 

编号	模型名称	调优参数
1	深度森林回归模型	n_estimators=int(params_best['n_estimators'])
2		n_trees=int(params_best['n_trees'])
3		max_layers=int(params_best['max_layers']

7.模型评估

7.1评估指标及结果 

评估指标主要包括R方、均方误差、解释性方差、绝对误差等等。

模型名称	指标名称	指标值
测试集
深度森林回归模型	R方	0.9686
	均方误差	1285.5238
	解释方差分	0.9686
	绝对误差	26.7584

从上表可以看出，R方分值为0.9686，说明模型效果比较好。

关键代码如下：  

  

7.2 真实值与预测值对比图

从上图可以看出真实值和预测值波动基本一致，模型效果良好。  

8.结论与展望

综上所述，本文采用了贝叶斯优化算法寻找深度森林回归算法的最优参数值来构建回归模型，最终证明了我们提出的模型效果良好。此模型可用于日常产品的预测。

# 描述性统计分析
print(df.describe())
print('******************************')

# y变量分布直方图
fig = plt.figure(figsize=(8, 5))  # 设置画布大小
plt.rcParams['font.sans-serif'] = 'SimHei'  # 设置中文显示
plt.rcParams['axes.unicode_minus'] = False  # 解决保存图像是负号'-'显示为方块的问题
data_tmp = df['y']  # 过滤出y变量的样本
# 绘制直方图  bins：控制直方图中的区间个数 auto为自动填充个数  color：指定柱子的填充色
plt.hist(data_tmp, bins='auto', color='g')  # 绘图
plt.xlabel('y')  # 设置x轴名称
plt.ylabel('数量')  # 设置y轴名称
plt.title('y变量分布直方图')  # 设置标题名称
plt.show()  # 展示图片

# 获取方式：
 
# 项目实战合集导航：
 
# https://docs.qq.com/sheet/DTVd0Y2NNQUlWcmd6?tab=BB08J2

# https://item.taobao.com/item.htm?ft=t&id=831101293649

param_grid_simple = {'n_estimators': (2.0, 10.0),  # 估计器数量
                     'n_trees': (50.0, 150.0),  # 每个估计器树的数量
                     'max_layers': (10.0, 20.0)  # 深度森林中最大级联层数
                     }

params_best, score_best = param_bayes_opt_DeepForest(5, 3)  # 调用贝叶斯优化器

print('最优参数组合:  ',
      'n_estimators的参数值为：', int(params_best['n_estimators']),
      'n_trees的参数值为：', int(params_best['n_trees']),
      'max_layers的参数值为：', int(params_best['max_layers'])
      )  # 打印最优参数组合
print('最优分数：  ', abs(score_best))  # 打印最优参数评分

'欢迎关注下方公众号，谢谢'