RDKit | 基于随机森林(RF)的机器学习模型预测hERG阻断剂活性

从分子相似性评估到使用机器学习技术的定量构效关系分析各种建模方法已应用于不同大小和组成的数据集(阻断剂和非阻滞剂的数量)。本研究中使用从公共生物活性数据开发用于预测hERG阻断剂的稳健分类器。随机森林被用来开发使用不同分子描述符,活性阈值和训练集合成的预测模型。与先前提取数据集的研究报告相比,该模型在外部验证中表现出优异的性能。


#导入依赖库
import pandas as pd
import numpy as np
import warnings; warnings.simplefilter('ignore')

from rdkit import Chem, DataStructs
from rdkit.Chem.Draw import IPythonConsole
from rdkit.Chem import PandasTools
from rdkit.Chem import AllChem, Draw

from sklearn.ensemble import RandomForestClassifier
#from sklearn.model_selection import StratifiedKFold
from imblearn.under_sampling import RandomUnderSampler
from sklearn.metrics import recall_score, roc_auc_score
from sklearn.model_selection import KFold, StratifiedKFold,StratifiedShuffleSplit
from sklearn.model_selection import train_test_split
from matplotlib import cm

import math
import pickle
import os

  • 2
    点赞
  • 13
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
随机森林是一种集成学习方法,它由多个决策树组成。在随机森林中,每个决策树都是由随机选择的样本和特征训练而成的。这种随机性使得随机森林可以减少过拟合,并提高模型的泛化能力。 在sklearn中,可以使用RandomForestClassifier类来构建随机森林分类器。具体的步骤如下: 1. 导入数据 2. 划分数据集为训练集和测试集 3. 实例化RandomForestClassifier类 4. 训练模型 5. 预测测试集 6. 评估模型性能 下面是一个示例代码,用于基于随机森林模型预测hERG阻断剂活性: ```python from sklearn.ensemble import RandomForestClassifier from sklearn.model_selection import train_test_split from sklearn.metrics import accuracy_score import pandas as pd # 导入数据集 data = pd.read_csv('data.csv') # 划分数据集 X = data.drop(['hERG'], axis=1) y = data['hERG'] X_train, X_test, y_train, y_test = train_test_split(X, y, test_size=0.2, random_state=42) # 实例化随机森林分类器 rf = RandomForestClassifier(n_estimators=100, random_state=42) # 训练模型 rf.fit(X_train, y_train) # 预测测试集 y_pred = rf.predict(X_test) # 评估模型性能 print('Accuracy:', accuracy_score(y_test, y_pred)) ``` 在这个例子中,我们首先导入了数据集,然后将其划分为训练集和测试集。接着,我们实例化了一个随机森林分类器,并使用训练集训练了模型。最后,我们使用测试集进行预测,并计算了模型的准确率。 需要注意的是,这只是一个简单的示例代码,对于实际问题,需要根据具体情况进行调参和优化。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

DrugAI

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值