肿瘤突变负荷计算

Python037

肿瘤突变负荷计算,第1张

之前看到TMB还以为在骂人呢。后来知道这是一个随着肿瘤免疫火起来的新名词。

定义: 肿瘤基因组每个编码区域的非同义突变的总数。

好吧,那就需要知道什么是非同义突变了。

不导致氨基酸改变的核苷酸变异我们称为同义突变,反之则称为非同义突变

这个应该是可以理解到的,因为密码子的简并性,有的点突变虽然发生,但并没有改变密码子所编码的氨基酸,突变了个寂寞。

TMB的等级划分

maf文件读进R语言后,一个函数即可计算,还出了个图呢。

我用的是卵巢癌的数据。随便拿哪个maf都行。

x第三列就是TMB值了,可以用四分位数看看tmb值的分布情况:

图的意思:

横坐标没写,就是TMB排序后的样本,纵坐标是log后的TMB值,也就是x表格里的第四列。默认取log,是为了画图直观。如果你不想取log,那就设置一下参数:

齐活啦!感谢maftools!

mutation load

英[mju:ˈteɪʃən ləud]美[mjuˈteʃən lod]

[词典]突变负荷

[例句]It is analyzed and compared that the influence on the mutation load of the DC capacitor voltage.

讨论了负载突变对直流电容上电压的影响,用仿真结果对比分析加入负载电流前馈控制前后直流电容电压控制性能的差异。

mutate:变异 突变 改变 数据修改

紧接着创建新的列gain和speed

新创建的列同时也可以使用(但是保留的方法仍然是赋值给某个名称):

由于系统显示限制,最后一列没有展示出来,运行view()函数即可:

如果只想要保留新的变量,那就使用transmute()函数:

此时参与计算的arr_delay、dep_delay、air_time、hours都消失,只有新转换的gain、hours 、gain_per_hour三列。

mutare()函数可以和前面提到的几个函数结合起来使用。其中进行运算的时候,肯定会涉及到R语言的计算语言,以下列出几个常用的:

举例: