一种饮用水源抗生素抗性相对风险快速识别方法

1.本发明涉及抗生素抗性识别技术领域，尤其涉及一种饮用水源抗生素抗性相对风险快速识别方法。


背景技术：

2.抗生素的大量使用使得自然水环境中抗生素抗性问题日益严峻。由携带抗性基因的抗性细菌及引起的感染治疗难度更大、治疗时间更长、治疗成本更高，甚至也会造成更高的死亡率。然而，由于缺乏抗性细菌致病及抗性基因污染的剂量效应模型及标准数据，以及获取相关数据存在较大的不确定性，对自然水环境中抗生素抗性风险的识别和定量评估也成了当今一大挑战。有学者从定性和定性-定量结合的层面提出了评估抗生素抗性风险的概念框架，但是目前还没有量化抗生素抗性风险的模型。
3.在此背景下，goh等提出了一个评估饮用水源水环境中抗生素抗性相对风险的框架，以样本间的抗生素抗性相对风险替代难以量化的绝对风险，在该框架中，goh等基于dalys数据库和args排序数据库，根据每种arg和arb的排名，计算了每种arg和arb的相对危害性权重，并结合多个自然流域的不同arg和arb的丰度量化了抗生素抗性的相对风险，该框架基于arb和args的真实检测数据计算抗生素抗性的相对风险。然而在实际应用中，arb和args的检测和提取是一个繁琐复杂且耗时的过程，使用该方法并不能快速识别出饮用水源水体中抗生素抗性的相对风险，基于众多易测量的潜在影响因素来执行抗生素抗性相对风险推断是一种值得尝试的途径，所以亟需一种有效的方法来快速识别饮用水源抗生素抗性相对风险。


技术实现要素：

4.本发明的目的是为了解决现有技术中存在的缺点，而提出的一种饮用水源抗生素抗性相对风险快速识别方法。其优点在于能够根据易测量的环境因素和土地使用类型等众多影响因素快速且准确地识别饮用水源地抗生素抗性的相对风险，从而大大节约了时间成本。
5.为了实现上述目的，本发明采用了如下技术方案：
6.一种饮用水源抗生素抗性相对风险快速识别方法，包括以下步骤：
7.步骤一：选取相关环境因素和土地使用类型因素为初始特征变量，其中环境因素包括：温度(temp)、ph值(ph)、氧化还原电位(orp)、电导率(ec)、电阻率(res)、总溶解固体(tds)、盐度(sal)、压强(press)、溶解氧(do)、浊度(turb)、24小时累计降雨量(24-hr rainfall)；土地使用类型因素包括：居民区比例(residential)、农业区比例(agriculture)、绿化区比例(green)、城市区比例(urban)，对各特征变量进行相关性分析，计算变量之间的spearman相关系数及显著性水平(p值)，剔除相关关系显著的变量(p《0.05)，通过计算发现residential、urban、green、agriculture四个变量之间相关关系显著，且green与amr相对风险值之间相关关系显著，因此剔除residential、urban和
agriculture三个特征变量，同理，ec、res、tds、sal四个变量之间相关关系显著，又res与amr相对风险值之间的相关关系更显著，因此剔除ec、tds、sal三个特征变量，考虑到现场数据获取的不确定性以及相对风险与多个影响变量之间的非线性关系，利用弹性网络对特征变量进行进一步筛选；
8.步骤二：采用goh等提出的amr相对风险评估框架计算出的amr相对风险得分为0～1范围内的连续值，将问题转换为一个二分类问题，相对风险值在中位数以上的样本的目标值编码为1，表示相对来说存在风险，相对风险值在中位数以下的样本的目标值编码为0，表示相对来说风险较小或没有风险；
9.步骤三：将样本数据按照7:3的比例划分训练集和测试集，将训练集数据输入随机森林分类器进行学习训练，用测试集数据测试训练得到的模型的性能；
10.步骤四：结合实际样本数据及遗传算法对随机森林算法的超参数进行自适应校准，在此方法中主要利用遗传算法的全局搜索能力同时对多个超参数进行协同优化，以找到最优超参数组合，从而形成最佳的风险识别机器学习模型。
11.优选地，所述弹性网络是一种使用l1和l2范数作为先验正则项训练的线性回归模型，它是岭回归和套索回归的结合，既具有岭回归的稳定性，又具有套索回归特征选择的功能。
12.优选地，所述弹性网络的目标函数如下：通过目标函数筛选出权重值不为0的特征变量，最终的界定出的模型特征变量包括green、orp、res、do、turb、24-hr rainfall。
13.优选地，所述随机森林分类器是一种以cart决策树为基学习器构建的bagging集成器，由多棵相对独立的决策树组合而成，模型最终结果由多棵决策树的结果投票决定，随机森林算法的随机性主要体现在对总体样本的bootstrap抽样和总体特征的随机选择。
14.设训练集为d，样本数为n，特征数为m，h
t
表示第t棵决策树，其中t＝1，2，
…
，t，则随机森林算法的步骤如下：
15.(1)通过bootstrap方法从数据集d中有放回地随机抽取n个样本，共进行t次抽样，得到训练子集d
t
，t＝1，2，
…
，t；
16.(2)利用训练子集d
t
构建t棵决策树；
17.(3)对决策树的每一个节点，采用无放回随机抽样的方式从m个特征中抽取m个特征(m＜m)，遍历这m个特征，基于“基尼指数最小”原则来确定最优划分属性。对于训练子集d
t
，设样本类别数为n，pk为数据集d
t
中第k类样本所占的比例，其中k＝1，2，
…
，n，则数据集d
t
的基尼指数定义为：
18.如果特征a的某一取值a将数据集d
t
分割成和两部分，则在特征a的条件下，数据集d
t
的基尼指数为：
[0019][0020]
gini(d
t
，a)表示经a＝a划分后数据集d
t
的不确定性，在基于基尼指数选择划分属性时，遍历所有特征和每个特征对应的值，以gini(d
t
，a)最小的特征和特征值为最优划分
特征和最优分割点对结点进行分裂，直至停止分裂；
[0021]
(4)记录每棵决策树的分类结果，采用投票法，以得票最高的类别作为样本的预测结果。
[0022]
优选地，所述超参数是机器学习算法在开始学习之前需要设置的参数，不同的超参数组合也会造成模型的性能不同，通过对超参数进行优化，选择最优的超参数组合可以提高学习的性能和效果，随机森林算法中树的数量和树的深度等都是超参数，以树的数量为例，nestimators太小容易造成欠拟合的问题，太大运算量太大，不能显著提升模型性能，需要对超参数进行调优以选择一个适中的数值。
[0023]
优选地，所述超参数优化主要对随机森林算法的n_estimators、max_depth、min_samples_split、max_leaf_nodes四个超参数进行调优，利用遗传算法进行超参数优化的具体步骤如下：
[0024]
(1)种群初始化：给定每个超参数值的范围，让每个超参数在取值范围内随机取值，并通过排列组合生成多个超参数组合，每个超参数组合为一条染色体。
[0025]
(2)计算适应度值。对于二分类问题，样本类别有正负两类，正类是饮用水源存在抗生素抗性风险，负类是饮用水源不存在抗生素抗性风险，可能出现的预测情况有以下四种：
[0026]
①
真阳性(tp)：实际上有风险，模型预测的结果也是有风险；
[0027]
②
假阳性(fp)：实际上没有风险，但是模型预测的结果是有风险；
[0028]
③
真阴性(tn)：实际上没有风险，模型预测的结果也是没有风险；
[0029]
④
假阴性(fn)：实际上有风险，但是模型预测的结果是没有风险。
[0030]
在二分类问题中，常用的评价指标有准确率(accuracy)、精确率(precision)、召回率(recall)、f1-socre、auc值等，其中auc值即roc曲线下的面积，roc曲线是一条以真阳性率(tpr)为纵轴，以假阳性率(fpr)为横轴的曲线。tpr和fpr的计算公式如下：
[0031][0032]
auc值指随机给定一个正样本和负样本，分类器输出该正样本为正的概率值比上分类器输出该负样本为正的那个概率要大的可能性，可用于评估模型的整体能力，auc值越大，则模型分类正确率越高，以10折交叉验证的auc值的均值为适应度函数，即
[0033]
(3)选择。采用轮盘赌选择法对种群中的个体进行选择，个体i被选中的概率为
[0034][0035]
(4)交叉。采用单点交叉法，随机生成一个0～1之间的小数r，如果r小于交叉概率pc，则在个体编码串中随机生成一个交叉点，将两个个体在交叉点前后的基因互换；否则，不进行交叉。
[0036]
(5)变异。采用基本位变异法，随机生成一个0～1之间的小数r，如果r小于变异概率，则在个体编码串中随机生成一个变异点，通过对变异点的值加减1来进行变异，否则，不
进行变异，直至达到算法收敛或者至最大允许的迭代次数。
[0037]
本发明的有益效果为：
[0038]
1、本发明能够根据易测量的环境因素和土地使用类型等众多影响因素快速且准确地识别饮用水源地抗生素抗性的相对风险，从而大大节约了时间成本。
[0039]
2、本发明结合真实数据样本及遗传算法对随机森林算法的超参数进行全局优化，能够快速有效地搜索出最优超参数组合，使模型能够自适应地针对不同数据集选择最优的超参数组合。
[0040]
3、本发明对真实饮用水源数据样本地进行案例分析研究，结果显示模型在测试数据集及检验数据集上的准确率、f1-score和auc值均超过了90％，且与未进行调参的随机森林模型相比，准确率和auc值分别提升了25％和18.2％。
附图说明
[0041]
图1为本发明提出的一种饮用水源抗生素抗性相对风险快速识别方法的工作流程结构示意图。
具体实施方式
[0042]
下面结合具体实施方式对本专利的技术方案作进一步详细地说明。
[0043]
参照图1，一种饮用水源抗生素抗性相对风险快速识别方法，包括以下步骤：
[0044]
步骤一：选取相关环境因素和土地使用类型因素为初始特征变量，其中环境因素包括：温度(temp)、ph值(ph)、氧化还原电位(orp)、电导率(ec)、电阻率(res)、总溶解固体(tds)、盐度(sal)、压强(press)、溶解氧(do)、浊度(turb)、24小时累计降雨量(24-hr rainfall)；土地使用类型因素包括：居民区比例(residential)、农业区比例(agriculture)、绿化区比例(green)、城市区比例(urban)，对各特征变量进行相关性分析，计算变量之间的spearman相关系数及显著性水平(p值)，剔除相关关系显著的变量(p《0.05)，通过计算发现residential、urban、green、agriculture四个变量之间相关关系显著，且green与amr相对风险值之间相关关系显著，因此剔除residential、urban和agriculture三个特征变量，同理，ec、res、tds、sal四个变量之间相关关系显著，又res与amr相对风险值之间的相关关系更显著，因此剔除ec、tds、sal三个特征变量，考虑到现场数据获取的不确定性以及相对风险与多个影响变量之间的非线性关系，利用弹性网络对特征变量进行进一步筛选；
[0045]
步骤二：采用goh等提出的amr相对风险评估框架计算出的amr相对风险得分为0～1范围内的连续值，将问题转换为一个二分类问题，相对风险值在中位数以上的样本的目标值编码为1，表示相对来说存在风险，相对风险值在中位数以下的样本的目标值编码为0，表示相对来说风险较小或没有风险；
[0046]
步骤三：将样本数据按照7:3的比例划分训练集和测试集，将训练集数据输入随机森林分类器进行学习训练，用测试集数据测试训练得到的模型的性能；
[0047]
步骤四：结合实际样本数据及遗传算法对随机森林算法的超参数进行自适应校准，在此方法中主要利用遗传算法的全局搜索能力同时对多个超参数进行协同优化，以找到最优超参数组合，从而形成最佳的风险识别机器学习模型。
[0048]
本实施例中，弹性网络是一种使用l1和l2范数作为先验正则项训练的线性回归模型，它是岭回归和套索回归的结合，既具有岭回归的稳定性，又具有套索回归特征选择的功能。
[0049]
本实施例中，弹性网络的目标函数如下：通过目标函数筛选出权重值不为0的特征变量，最终的界定出的模型特征变量包括green、orp、res、do、turb、24-hr rainfall。
[0050]
本实施例中，随机森林分类器是一种以cart决策树为基学习器构建的bagging集成器，由多棵相对独立的决策树组合而成，模型最终结果由多棵决策树的结果投票决定，随机森林算法的随机性主要体现在对总体样本的bootstrap抽样和总体特征的随机选择。由于其具有训练速度快、泛化能力强等特点，随机森林算法一直被认为是强大的机器学习算法，并被广泛用于风险识别领域。
[0051]
设训练集为d，样本数为n，特征数为m，h
t
表示第t棵决策树，其中t＝1，2，
…
，t，则随机森林算法的步骤如下：
[0052]
(1)通过bootstrap方法从数据集d中有放回地随机抽取n个样本，共进行t次抽样，得到训练子集d
t
，t＝1，2，
…
，t；
[0053]
(2)利用训练子集d
t
构建t棵决策树；
[0054]
(3)对决策树的每一个节点，采用无放回随机抽样的方式从m个特征中抽取m个特征(m＜m)，遍历这m个特征，基于“基尼指数最小”原则来确定最优划分属性。对于训练子集d
t
，设样本类别数为n，pk为数据集d
t
中第k类样本所占的比例，其中k＝1，2，
…
，n，则数据集d
t
的基尼指数定义为：
[0055]
如果特征a的某一取值a将数据集d
t
分割成和两部分，则在特征a的条件下，数据集d
t
的基尼指数为：
[0056][0057]
gini(d
t
，a)表示经a＝a划分后数据集d
t
的不确定性，在基于基尼指数选择划分属性时，遍历所有特征和每个特征对应的值，以gini(d
t
，a)最小的特征和特征值为最优划分特征和最优分割点对结点进行分裂，直至停止分裂；
[0058]
(4)记录每棵决策树的分类结果，采用投票法，以得票最高的类别作为样本的预测结果。
[0059]
本实施例中，超参数是机器学习算法在开始学习之前需要设置的参数，不同的超参数组合也会造成模型的性能不同，通过对超参数进行优化，选择最优的超参数组合可以提高学习的性能和效果，随机森林算法中树的数量和树的深度等都是超参数，以树的数量为例，nestimators太小容易造成欠拟合的问题，太大运算量太大，不能显著提升模型性能，需要对超参数进行调优以选择一个适中的数值，利用遗传算法对随机森林算法的nestimators，maxdepth，min_samples_split，max_leaf_nodes四个超参数调优，使模型能够自适应地针对不同数据集选择最优的超参数组合。
[0060]
本实施例中，超参数优化主要对随机森林算法的n_estimators、max_depth、min_
samples_split、max_leaf_nodes四个超参数进行调优，利用遗传算法进行超参数优化的具体步骤如下：
[0061]
(1)种群初始化：给定每个超参数值的范围，让每个超参数在取值范围内随机取值，并通过排列组合生成多个超参数组合，每个超参数组合为一条染色体。
[0062]
(2)计算适应度值。对于二分类问题，样本类别有正负两类，正类是饮用水源存在抗生素抗性风险，负类是饮用水源不存在抗生素抗性风险，可能出现的预测情况有以下四种：
[0063]
①
真阳性(tp)：实际上有风险，模型预测的结果也是有风险；
[0064]
②
假阳性(fp)：实际上没有风险，但是模型预测的结果是有风险；
[0065]
③
真阴性(tn)：实际上没有风险，模型预测的结果也是没有风险；
[0066]
④
假阴性(fn)：实际上有风险，但是模型预测的结果是没有风险。
[0067]
在二分类问题中，常用的评价指标有准确率(accuracy)、精确率(precision)、召回率(recall)、f1-socre、auc值等，其中auc值即roc曲线下的面积，roc曲线是一条以真阳性率(tpr)为纵轴，以假阳性率(fpr)为横轴的曲线。tpr和fpr的计算公式如下：
[0068][0069]
auc值指随机给定一个正样本和负样本，分类器输出该正样本为正的概率值比上分类器输出该负样本为正的那个概率要大的可能性，可用于评估模型的整体能力，auc值越大，则模型分类正确率越高，以10折交叉验证的auc值的均值为适应度函数，即
[0070]
(3)选择。采用轮盘赌选择法对种群中的个体进行选择，个体i被选中的概率为
[0071][0072]
(4)交叉。采用单点交叉法，随机生成一个0～1之间的小数r，如果r小于交叉概率pc，则在个体编码串中随机生成一个交叉点，将两个个体在交叉点前后的基因互换；否则，不进行交叉。
[0073]
(5)变异。采用基本位变异法，随机生成一个0～1之间的小数r，如果r小于变异概率，则在个体编码串中随机生成一个变异点，通过对变异点的值加减1来进行变异，否则，不进行变异，直至达到算法收敛或者至最大允许的迭代次数。
[0074]
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，根据本发明的技术方案及其发明构思加以等同替换或改变，都应涵盖在本发明的保护范围之内。