Python社区  »  机器学习算法

TCGA官方数据挖掘文章教你机器学习or深度学习

生信技能树 • 1 年前 • 273 次点击  


学徒和学员已经陆续出师,是时候把生信技能树的舞台交给后辈了!

最近我们又组织了:《机器学习加深度学习资料大放送(附上资料群)》交流群,感觉吧,大家松鼠症发作收集整理了大把资料最后却束之高阁,也不是一个事啊。所以就安排学徒系统性讲解一下机器学习的应用。本次教程参考TCGA 官方的一篇文章Machine Learning Detects Pan-cancer Ras Pathway Activation in The Cancer Genome Atlas文章,文章的源码也可以在官方的 github 上搜到:

  • https://www.cell.com/cell-reports/fulltext/S2211-1247(18)30389-9
  • https://github.com/greenelab/pancancer

 

本笔记来源于福建医科大学研究生

前些日子,老师给了我一份大纲并让我重启一下机器学习的系列,但是我看完之后发现,可能是由于老师在写大纲的时间比较早,并且也长时间没更新,导致内容略微有些过时,于是我打算用我自己的方式重启这个系列。其实并不是我不想写,是因为我看了鲍志炜师兄的机器学习教程,讲实话,真的很不错的。我的目标就是这个系列能够有他一半的阅读量就行。而且我觉得我可能没有资格把我的系列取名后面加上教程,因为我觉得我可能并没有这个资格。。。(也不知道老师为何如此看重我然我重启这个项目。。)

还是深度学习香

  深度学习可能大家都不陌生,不过我感觉大部分的人只停留在不陌生的阶段。虽然大家可能都知道,深度学习的主流应用方向还是图像识别,人脸识别,seq2seq ,或者目前非常火热的自动驾驶等等。这样子的话听起来似乎深度学习和生信一点关系都没有。

  但是事实并非如此。到目前来说其实已经有很多的文章已经把深度学习的方法应用的计算生物学和生物信息学的领域之中了(这篇文章就是关于深度学习在生物信息学领域应用的综述),所以本次系列课程的任务重心旨在把深度学习的各种方法合理的应用在生物信息学的各个领域之中(尤其是各组学领域,但是因为我对基因组的领域不是很熟悉,所以重心可能会放在下游的表达谱这一块)

Why Deep learning?

  为什么选择深度学习而不是机器学习,最重要的原因就是上面提到的鲍志炜师兄已经做过了,我觉得我做肯定不如他,所以就直接跳到深度学习了。

  第二个原因就是,深度学习所需要的大数据刚好在生信领域完美契合,有很多公共的大型数据库供我们使用,最典型的就是 TCGA、GEO这两个我们非常熟悉的数据库

  并且我们这次系列主要还是以实战为主,头篇参考的文章是来自 TCGA  官方的一篇文章Machine Learning Detects Pan-cancer Ras Pathway Activation in The Cancer Genome Atlas文章的源码也可以在官方的 github 上搜到。虽然这篇文章我在我的 github 上发过一次源码教程(在生信菜鸟团也发过一部分,但是由于审核的问题后面就没发了)。但是这一次系列课程是完全不一样的,因为 TCGA 用的是机器学习的方法,我打算使用深度学习的方法重启这篇文章

源码在 Scripts 里面推荐大家先把源码看完,虽然内容是机器学习:https://github.com/greenelab/pancancer

生信相关的知识,在本次系列教程都不会有任何的赘述,如果需要的话可以去菜鸟团补课,或者《机器学习加深度学习资料大放送(附上资料群)》交流群。

python 相关的语法也是不会有任何的讲解的,万能的B站太多视频资料了,所以默认大家都会python哈!

万物皆函数!

  OK,不讲废话了,直接进入系列的第一个篇章。这篇文章的目的很明确——利用基因表达谱+额外的信息从而预测癌样本 RAS 通路激活情况。在多说一句废话:由于 TCGA 上可以利用的资源比较多,也可以尝试用同样的思路预测别想要的信息

  那么,深度学习到底学习的是什么呢?下面直接用文献中的图来解释:

image-20201005201458377

  这张图就基本上解释了机器学习或者深度学习的本质——函数映射。对,就是那个我们在中学的时候就学过的函数 f(x) = y ,这个函数的映射关系可以用到宇宙万物。

  例如,我们可以将 x 记录为:

  1. 早晨起来投的硬币是正面还是反面
  2. 今天天气怎么样
  3. 出门穿的是什么样的衣服
  4. .....

  预测:y 为晚餐吃得是什么?

  或许有的人认为这个是扯淡的预测,但是很遗憾,并不是,因为只要我们能够找到正确的 f(x)—>y 的映射关系,就一定能够根据 x 的信息预测这个人晚餐吃得是啥,因为马克思说过:联系是事务普遍的规律。之所以觉得不可能,是因为光靠我们这个脑子是无法计算出一个正确的映射关系的。因此,对于任何一个机器学习还是深度学习,其本质就是找到一个映射关系将 x 映射到 y。当然,实际情况及其复杂一般都是这样的映射网络 f1(f2(f3(f4(f5(.....fn(X)))))) 就和俄罗斯套娃一样套了好几层的映射关系。因此在很多情况下,即便是机器都很难学习到真正的映射网络。在实际情况中,学习的到的映射模型只是一种和真实模型很像的模型。但即便是这样,这种和很像的模型通常情况下也能够解决很多事情。

学习的过程?

  接下来的问题就在于要怎么让电脑去学习这个映射关系 f()或者说f1(f2(f3(f4(f5(.....fn(X)))))),我们从最简单的模型开始讲吧!

image-20201005201553896

其中输入层就是手头上的数据 X(输入数据) 隐含层就是 f(.) 而输出层就是 y'(称为预测值),那么如何判断隐藏层的f(.)就是我们想要的映射关系呢?当然只需要判断输出的 y'y(真实数据) 是否相同,理想情况是输出的每一个 y' == y,在实际情况中是不可能出现这种情况的,于是我们构造了一个能够判断预测值和真实数据之间的差异的函数,叫做 loss(损失)函数,理论上 loss函数越小,训练出来的模型就越接近真实模型,那么降低loss就是深度学习的核心目的。由此可以得出,深度学习 == 求解最小的loss。而 loss = g(y - f(x)), 其中 g()是 loss 函数映射(这种映射是不需要学习的,因为别人已经根据不同的方案给好了),由于 x 是定值,y 也是定值,g(.)也是别人定义好的。那么我们只需要**“找”出 loss 最小值的时候对应的 f(.)即可。我们将“找”的过程称为迭代,不同的“找法”也不产生不同的结果,我们将找法**称为优化器。至此,深度学习的核心组件全部完成接下来只需要将其组装在一起就行,下面我举一个简单的例子。

  机器的学习过程和人的学习过程其实非常像,回顾一下我们备考四六级的时候的样子我们每天埋头备考,我们买了一沓的星火英语,华研外语等等,里面有模拟题(准备训练集样本)也有答案(准备训练集标签),然后开始做题(f(.)),做完核对,发现做错了几十道题目了甩自己几十耳光(计算loss),不懂的找知识点(优化器),然后回去查看一下到底为什么做错(反向传递),然后继续做题(修正f(.))。然后继续做题继续错继续抽自己继续找知识点。。。。(迭代)。然后上考场(验证集)最后考了666 分。这样就训练出了一个较为完美深度学习框架。

  总结一下到底怎样四六级可以考高分呢(训练一个好的模型)?首先是肯定要有模拟题(训练集),模拟题越多越好,当然这里又要关心模拟题的质量,如果你拿少儿剑桥英语去备考四六级(质量低的数据),那么即便你做再多的题目也没用,此外模拟题的答案也很重要(训练集标签),如果你的答案就是错的,那么在核对的过程你只会加深错误的印象,考试就不可能考高分(模型泛化能力差)。那么题目要怎么做呢?f(.)(隐藏层的设置)我们要针对四六级的题型去选择题型去做(隐藏层类型),每种题型做几道(隐藏层个数),这些都会影响到我们最后考试的成绩。除此之外合理的惩罚制度(loss 函数的选择)也非常重要,错 100 题都不愿意甩自己一巴掌的人最后考试的结果肯定不如错一题甩自己一巴掌的人考的好(当然,如果没做错也乱甩自己到最后自己都不知道自己哪里错了也不行)。最后,如何整理和找到自己真正不懂的知识点也及其的重要(优化器的选择),这一点如果做得好,会让你事半功倍

  综上,机器学习最重要的组件有:

  1. 较为良好的数据和相对正确的标签(这一点在生信领域尤为重要,因为图像识别领域标签一定是对的,但是我们这里的话定义就很模糊,后期会讲)
  2. 隐藏层的设置
  3. loss 函数的选择
  4. 优化器

final

  那么,回到我们生信相关的具体的例子来说就是,x 就是患者所有基因的表达值,Y 就是想要预测的特征(某个基因的突变,患者生存风险,某条通路的异常等等),接下来的教程我会直接用 TCGA 的数据进行实战演练,带领大家做完一次完整的深度学习框架

建议

  1. 如果可以的话,希望有兴趣和我一起学习深度学习的小伙伴们可以吧上面 TCGA 官方的文章的内容复现一遍。
  2. 如果不是数学或者物理科班出生(或者说高中以来没有参加过数学竞赛的,考研数一没 135 的)的然后转行到生信的建议在后续的教程中关于一些数学知识不懂的话没有必要死磕,因为并不影响调参
  3. 如果有不信邪的可以试着尝试先看一下 B 站的一个白板系列劝退课程

第二篇笔记

接下来的两篇教程不会涉及到具体的深度学习的代码讲解,不懂没有关系,因为后续会进行讲解,因为最近我自己手头上一下子要参与两个课题所以更新可能会慢一点,如果有需要的朋友可以去菜鸟团补课,或者《机器学习加深度学习资料大放送(附上资料群)》交流群。

  没有意义的开场的废话就不讲了,直接先上手跑一遍大概的流程

(接上一篇的 TCGA 文章)原始数据我们就不用了,直接用他处理好的数据:https://github.com/greenelab/pancancer

image-20201005201639667

Start

环境配置

# 导入模块
import os
import math
import torch  # pythorch 模块
import itertools
import warnings

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

from torch import nn

T = True
F = False
%matplotlib inline
%config InlineBackend.figure_format = 'svg'

warnings.filterwarnings("ignore")
plt.rc('font', family='Times New Roman')
plt.rcParams['figure.figsize'] = (43)
my_colors = ["#1EB2A6""#ffc4a3""#e2979c""#F67575"]

  接下来呢我们就按照文章的代码进行数据预处理,当然,只有数据预处理的这一部分是一样的,后续的建模过程会完全不一样




    
# 数据导入
rnaseq_full_df = pd.read_table('data/pancan_rnaseq_freeze.tsv.gz', index_col=0)
mutation_df = pd.read_table('data/pancan_mutation_freeze.tsv.gz', index_col=0)
sample_freeze = pd.read_table('data/sample_freeze.tsv', index_col=0)
mut_burden = pd.read_table('data/mutation_burden_freeze.tsv')

# 查看数据类型
"""sample_freeze 样本肿瘤信息矩阵"""
sample_freeze.head()
sample_freeze.shape
"""mut_burden 样本肿瘤突变负荷"""
mut_burden.head()
mut_burden.shape
"""rnaseq_full_df rnaseq表达谱"""
rnaseq_full_df.head()
rnaseq_full_df.shape
"""mutation_df mutation(0-1)矩阵"""
mutation_df.head()
mutation_df.shape

这里多说一句,建议大家使用 jupyter,配置什么的可以百度搜一下生信菜鸟团鲍志炜的目录

  接下来就是提取文章中要预测的RAS 通路相关的信息

提取目标基因的突变信息

genes = 'KRAS,HRAS,NRAS' 
genes = genes.split(',')
common_genes = set(mutation_df.columns).intersection(genes) # 目标基因和mutation矩阵基因交集
common_genes = list(common_genes.intersection(rnaseq_full_df.columns)) # 交集后和rnaseq矩阵交集
y = mutation_df[common_genes] # 提取目标基因中和mutation矩阵和rnaseq交集的基因,在样本中的突变情况
"""y"""
y.head()
image-20201005201753798

  不要忘记把表达谱中如果有测到这个基因的话是要删除的,不然的话后续的预测就没有意义了

print("移除的基因是 : %s"  %common_genes )
rnaseq_full_df.drop(common_genes, axis=1, inplace=True)
# 移除的基因是 : ['HRAS', 'NRAS', 'KRAS']

然后,因为 TCGA 的数据比较大,所以不仅有基因的突变信息,还有拷贝数变异信息 所以要将拷贝数变异信息和突变谱信息进行整合。无论是拷贝数还是突变谱,只要有一个发生突变了就当做突变

整个的流程是这样的:

  1. 提取突变信息(上面已经提取过)
  2. 提取基因的拷贝数变异信息
  3. 如果是抑癌基因,则考虑是否有拷贝数删失,如果是原癌基因则考虑是否有拷贝数增加
  4. 整合拷贝数变异信息和突变信息,只要有一个突变则认为样本 RAS 通路异常

这里我不知道我是否有表述清楚(如果看懂了这句话直接无视),因为文章研究的是 ras 通路异常,因为有的样本只测了突变谱,有的样本又只测了拷贝数,有的样本又两个都测,所以为了让样本尽可能的多,作者就进行了这个操作(如果还是没讲清楚还是直接看文章把)

loss_df = pd.read_table('data/copy_number_loss_status.tsv.gz', index_col=0# 缺失
gain_df = pd.read_table('data/copy_number_gain_status.tsv.gz', index_col=0)# 拷贝
cancer_genes_df = pd.read_table('data/vogelstein_cancergenes.tsv'# 癌基因分类表
'''癌基因分类表'''
cancer_genes_df[['Gene Symbol','Ocogene score* ','Tumor Suppressor Gene score* ']].iloc[0:5,:]
genes_sub = cancer_genes_df[cancer_genes_df['Gene Symbol'].isin(common_genes)]# 获取基因信息

tumor_suppressor = genes_sub[genes_sub['Classification*'] == 'TSG'# 提取抑癌基因
oncogene = genes_sub[genes_sub['Classification*'] == 'Oncogene'# 提取原癌基因
copy_loss_sub = loss_df[tumor_suppressor['Gene Symbol']] # 如果是抑癌基因则获取拷贝数缺失信息
copy_gain_sub = gain_df[oncogene['Gene Symbol']]  # 如果是原癌基因则获取拷贝数增加信息

genes_sub.head()

copy_loss_sub.columns = [col + '_loss' for col in copy_loss_sub.columns] # 修改列名
copy_gain_sub.columns = [col + '_gain' for col in copy_gain_sub.columns]

# Add columns to y matrix
y = y.join(copy_loss_sub) # 合并拷贝数变异信息
y = y.join(copy_gain_sub)

# 数据清洗
y = y.fillna(0
y = y.astype(int)
 # 三个目标基因中只要有一种基因突变则人为患者为突变型 
y['total_status'] = y.max(axis=1)
y = y.reset_index().merge(sample_freeze,
                          how='left').set_index('SAMPLE_BARCODE'# 合并样本疾病类型
'''y'''
y.head()


count_df = y.groupby('DISEASE').sum() # 计算不同癌的型基因突变频率
prop_df = count_df.divide(y['DISEASE'].value_counts(sort=False).sort_index(),
                          axis=0)



count_table = count_df.merge(prop_df, left_index=True, right_index=True,
                             suffixes=('_count''_proportion'))
count_table.to_csv("count_table_file.csv")
prop_df.head()

  移除肿瘤突变负荷过高的样本

ras_diseases='BLCA,CESC,COAD,ESCA,HNSC,LUAD,LUSC,OV,PAAD,PCPG,READ,SKCM,STAD,TGCT,THCA,UCEC'
diseases = ras_diseases
diseases = diseases.split(',')  
y_df = y[y.DISEASE.isin(diseases)].total_status # 提取指定疾病的基因患者突变信息
common_samples = list(set(y_df.index) & set(rnaseq_full_df.index)) # 提取指定疾病的基因表达谱
y_df = y_df.loc[common_samples]
rnaseq_df = rnaseq_full_df.loc[y_df.index, :]
burden_filter = mut_burden['log10_mut'] 5 * mut_burden['log10_mut'].std()
mut_burden = mut_burden[burden_filter]

y_matrix = mut_burden.merge(pd.DataFrame(y_df), right_index=True,
                            left_on='SAMPLE_BARCODE').set_index('SAMPLE_BARCODE')

  添加协变量信息

这里简单解释一下什么叫协变量信息,协变量信息是一种独立的信息(医统里面好像有)一般是不受实验操作影响但是又会对结果产生影响的信息。举个最简单的例子就是性别,如果在分析的时候需要用到性别信息的时候,这时候是不能简单把性别为男设为 1 女设为 0(因为你无论怎么设置都会有人 bb),这时候就出现协变量矩阵,直接把性别这一列拆成两列,一列“男”,一列“女”,如果一个样本性别是男,则在“男”这一列设为 1 ,“女”这一列设为 0(也叫 one-hot)。

y_sub = y.loc[y_matrix.index]['DISEASE']
covar_dummy = pd.get_dummies(sample_freeze['DISEASE']).astype(int)
covar_dummy.index = sample_freeze['SAMPLE_BARCODE']
covar = covar_dummy.merge(y_matrix, right_index=True, left_index=True)
covar = covar.drop('total_status', axis=1)


'''covar 疾病协变量矩阵'''
covar.head()
"""y_sub 样本肿瘤信息"""
y_sub.head()
"""y_df 目标基因患者突变状态"""
y_df.head()

就像这样我们把疾病类型变成了协变量矩阵(也叫 one-hot)

最后提炼一下疾病的分组信息然后输出,在输出之前记得归一化一下(以后会讲为什么这么做)

y_df = y_df.loc[y_sub.index] # 基因突变矩阵
strat = y_sub.str.cat(y_df.astype(str))  # 疾病类型
x_df = rnaseq_df.loc[y_df.index, :] # 基因表达矩阵
'''疾病类型+0/1'''
'''例如: CESC0'''
'''表示患者属于cesc 且目标基因未突变'''
strat.head()
from sklearn.preprocessing import StandardScaler 
fitted_scaler = StandardScaler().fit(x_df)
x_df_update = pd.DataFrame(fitted_scaler.transform(x_df),
                           columns=x_df.columns)
x_df_update.index = x_df.index
x_df = x_df_update.merge(covar, left_index=True, right_index=True)
x_df.head()
x_df.shape
x_df.to_csv("raw_exp.csv")
y_df.to_csv("raw_clinic.csv")
y_sub.to_csv("strat.csv")

Final

下一章就和这篇文章完全没关系了,所以后续的内容得自己看哈~~~

第三篇笔记突然累了

三篇笔记合起来没几个字的教程,我却写了四天,加上课题进度又赶已经连续四天晚上 3 点才睡了。读书真的太累了,突然又不想读书想考公务员了。。

正文

  关于深度学习的内容,我会以 pytorch 为基本开始展开这个系列当然官方也给出了教科书,不过是纯英语的(该练英语了)

image-20201005201917479

  包括在网上也有很多入门教程,所以我这个系列的定位是实战,虽然网上也有很多实战的内容,但是我这个系列是主要在表达谱这个版块的应用不要总是天天 limma,deseq2,edgr,时代变了赶紧跟着我的系列进入新时代。

强调!!!!

这次的内容的代码,看不懂没关系,因为后续都会展开讲解。所以只要跟着看就行。目的是让大家知道四个组间分别在哪一个步骤设置

回顾一下深度学习最重要的四个组件

  1. 较为良好的数据和相对正确的标签
  2. 隐藏层的设置
  3. loss 函数的选择
  4. 优化器

Start

环境配置

import os
import math
import torch  # pythorch 模块
import itertools
import warnings

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt
import torch.utils.data as Data
from torch import nn
from sklearn.model_selection import train_test_split
from sklearn.metrics import roc_curve,precision_recall_curve ,roc_auc_score,auc,average_precision_score 

T = True
F = False
%matplotlib inline
%config InlineBackend.figure_format = 'svg'

warnings.filterwarnings("ignore")
plt.rc('font', family='Times New Roman')
plt.rcParams['figure.figsize'] = (43)
my_colors = ["#1EB2A6""#ffc4a3""#e2979c""#F67575"]
def roc_plot(label=None,scores=None,title=None,file_name=None, save=False):
    fpr, tpr, threshold = roc_curve(label, scores)  ###计算真正率和假正率
    roc_auc = roc_auc_score(label, scores)  ###计算auc的值
    plt.figure(figsize=figsize)
    plt.plot(fpr, tpr, color=my_colors[3], lw=lw,
             label='AUC: %0.4f' % roc_auc)  ###假正率为横坐标,真正率为纵坐标做曲线
    plt.plot([01], [01], color=my_colors[0], lw=lw, linestyle='--')
    plt.xlim([-0.011.01])
    plt.ylim([-0.011.01])
    plt.xlabel('1 - Sepcificity', size=label_size)
    plt.ylabel('Sensitivity', size=label_size)
    plt.yticks(fontsize=ticks_size)
    plt.xticks(fontsize=ticks_size)
    plt.title(title, size=title_size)
    #plt.grid(True)
    plt.legend(bbox_to_anchor=legend_sit,
               fontsize=legend_size,
               borderaxespad=0.)
    if save:
        plt.savefig(file_name)
    plt.show()
    plt.clf()

  把上一次的数据导入

第一个组件:数据和标签

Expr = pd.read_csv('raw' + '_exp.csv', index_col=0# 表达谱
Clinic = pd.read_csv( 'raw' + '_clinic.csv', index_col=0# label
strat = pd.read_csv('strat.csv', index_col=0# 分层信息
strat = strat['DISEASE'].str.cat(Clinic.astype(str))

Expr.head()
Clinic.head()
strat.head()

  拆分样本

分为训练集和验证集,一般比例为 9:1(我个人喜欢),验证集是不可以在训练模型中的任何一个步骤出现,只能在验证的时候出现 这就好比我们考研的时候,总是会把最近一年的真题留在靠前的最后一天做,除了近一年以外的真题都叫训练集,只有最后一套才能叫做验证集,因为只有这样才能检验我们一年备考的学习成果,上考场才能更加的有信息。

Samples = Expr.index  # 样本名
Genes = Clinic.columns  # 基因名
Expr = np.array(Expr) # 转换为 numpy.array
Clinic = np.array(Clinic["total_status"]) # 这里注意下,因为我们导入的形式是data.frame的虽然只有一列,但是如果直接转为 numpy 会出现问题,你们可以尝试一下
x_train, x_test, y_train, y_test = train_test_split(
    Expr,
    Clinic,
    test_size=0.1,  # 验证集比例
    random_state=0,  # 随机种子
    stratify=strat  #分层抽样
# 这里用到了 sklearn 的库

del Expr, Clinic ## 删除变量释放一下内存,不然会卡死

# 将数据转为 torch
x_train = torch.from_numpy(x_train)
x_test = torch.from_numpy(x_test)
y_train = torch.from_numpy(y_train)
y_test = torch.from_numpy(y_test)


x_train = x_train.float()
x_test = x_test.float()
# 将 label 转为 tensor 类型
y_train = torch.LongTensor(y_train) 
y_test = torch.LongTensor(y_test)
batch_size = 512
# 将训练集分割成小 batch,当然也可以不分割,但是我劝你最好还是做,当然你也可以不听劝
train_iter = Data.TensorDataset(x_train, y_train)
train_iter = Data.DataLoader(
    dataset=train_iter,      # 数据,封装进Data.TensorDataset()类的数据
    batch_size=batch_size,      # 每块的大小
    shuffle=False,               # 要不要打乱数据 (这里不用打乱顺序)
    num_workers=8,              # 多进程(multiprocess)来读数据
)

# 验证集不分割
test_iter = Data.TensorDataset(x_test, y_test)
test_iter = Data.DataLoader(
    dataset=test_iter,      # 数据,封装进Data.TensorDataset()类的数据
    batch_size=x_test.shape[0],      # 每块的大小
    shuffle=False,               # 要不要打乱数据 
    num_workers=8,              # 多进程(multiprocess)来读数据
)

  构建神经网络

# 组件 2 隐藏的设置
# 这里我们定义一个最最最最最简单的网络
class MyNet(nn.Module):
    def __init__(self):
        super(MyNet, self).__init__()
        self.mlp = nn.Sequential(nn.Linear(x_train.shape[1], 512), nn.ReLU(),
                                 nn.Linear(512, 256), nn.ReLU(),
                                 nn.Linear(256, 2), nn.Sigmoid())
    def forward(self, x):
        return self.mlp(x)

    

Net = MyNet()
def Trainning(Net=None,
              train_iter=None,
              test=None,
              loss=None,
              num_epochs=None,
              optimizer=None):
    for epoch in range(num_epochs):
        train_loss, train_acc,n = 0.0,0.0,0
        for X,y in train_iter:
            y_pred = Net(X)
            l = loss(y_pred,y).sum()
            
            optimizer.zero_grad()
            l.backward()
            optimizer.step()
            
            train_loss += l.item()
            train_acc += (y_pred.argmax(dim=1) == y).sum().item()
            n += y.shape[0]
        with torch.no_grad():
            for test_X ,test_y in test_iter:
                test_acc = (Net(test_X).argmax(dim=1) == test_y).sum().item()/len(y_test)
        print('epoch %d, loss %.4f, train acc %.3f, test acc %.3f'
              % (epoch + 1, train_loss / n, train_acc / n, test_acc))

  训练模型

# 组件 3loss
# 选择分类交叉熵作为损失函数
loss = nn.CrossEntropyLoss()
num_epochs = 5 # 先训练 5 轮看看
# 组件 4 优化器
optimizer = torch.optim.Adamax(params=Net.parameters()) # 优化器,这个就没必要懂了(不仅现在不用懂以后也不用懂),因为解释起来真的很麻烦
Trainning(Net,
              train_iter,
              test_iter,
              loss,
              num_epochs,
              optimizer)
 你看随随便便训练一下就高达 85%的精度,是不是发现 cell 也不过如此

接下来在看看 auc,文章的 auc 是 84,我们用深度学习随随便泡一下也 78% 了

  AUC




    
with torch.no_grad():
    prob= Net(x_test).numpy()
    pridict = Net(x_test).argmax(dim=1)
    scores = [i[1] for i in prob]
label = y_test


# 定义函数
lw = 3.5
ticks_size = 12
label_size = 16
title_size = 18
legend_size = 14
legend_sit = (0.85, 0.3)
figsize = (5, 5)
roc_plot(label,
         scores,
         title="validation"
         #file_name='AUC.pdf',
         #save=True)
        )

Final

  不知道看完大家有什么感受,可能是云里雾里一脸懵逼(虽然我觉得这可能是我教程写的有点太差了,这一点我会尽可能的改进,所以不要喷,毕竟写系列教程我也是第一次)。代码这种东西确实比较枯燥,我也想在别写代码的时候尝试添加一下自己的话,这样有点互动的感觉,而不是整片文章都是代码,因为我最近事情比较多,所以学习的过程可能比较就,但是我也有抽时间把老师写的 R 语言入门教程看了一下,模仿一下老师学习如何写一个比较优秀的入门教程

  希望之后的教程能写好一点吧。。。

肯定是有后续的,基础知识大家可以去菜鸟团补课,或者《机器学习加深度学习资料大放送(附上资料群)》交流群跟作者交流哈。

文末友情推荐

要想真正入门生物信息学建议务必购买全套书籍,一点一滴攻克计算机基础知识,书单在:什么,生信入门全套书籍仅需160 。如果大家没有时间自行慢慢摸索着学习,可以考虑我们生信技能树官方举办的学习班:

如果你课题涉及到转录组,欢迎添加一对一客服:详见:你还在花三五万做一个单细胞转录组吗?

号外:生信技能树知识整理实习生招募,长期招募,也可以简单参与软件测评笔记撰写,开启你的分享人生!另外,:绝大部分生信技能树粉丝都没有机会加我微信,已经多次满了5000好友,所以我开通了一个微信好友,前100名添加我,仅需150元即可,3折优惠期机会不容错过哈。我的微信小号二维码在:0元,10小时教学视频直播《跟着百度李彦宏学习肿瘤基因组测序数据分析》

Python社区是高质量的Python/Django开发社区
本文地址:http://www.python88.com/topic/74096
 
273 次点击  
分享到微博