R语言基因表达量转换（TPM、FPKM、RPKM）

基因表达量一般以TPM或FPKM为单位来展示。

TPM，Transcripts Per Kilobase Million
计算公式：
TPMi=(Ni/Li)*1000000/sum(Ni/Li+………+ Nm/Lm)
Ni：mapping到基因i上的read数； Li：基因i的外显子长度的总和

计算方法也同RPKM/FPKM类似，
先对每个基因的read数用基因的长度进行校正（标准化）
RPK <- count/(efflength/1000)

校正后的这个基因read数(Ni/Li)与校正后的这个样本的所有read数（sum(Ni/Li+………+ Nm/Lm)）求商
RPK/sum(RPK)

最后再乘以10^6
RPK/sum(RPK)*1e6

TPM可以看作是RPKM/FPKM值的百分比。

counts2TPM <- function(count=count, efflength=efflen){RPK <- count/(efflength/1000)   #每千碱基reads (Reads Per Kilobase) 长度标准化PMSC_rpk <- sum(RPK)/1e6        #RPK的每百万缩放因子 (“per million” scaling factor ) 深度标准化RPK/PMSC_rpk
}  tpm <- as.data.frame(apply(counts,2,counts2TPM))

R语言基因表达量转换（TPM、FPKM、RPKM）相关推荐

r语言入门——颜色转换
r语言提供了许多种设置颜色的函数,其中rgb函数颜色强度要求取值范围在[0,1],而通常rgb是用3个取值在[0,255]的数表示,转换方法很简单,只需要自定义一个新的函数即可通过输入普通的rgb值来 ...
R语言countrycode包转换国家名字和代码
本文首发于公众号:医学和生信笔记,完美观看体验请至公众号查看本文. 不同的数据源使用不同的编码方案来表示国家(例如CoW或ISO).这带来了两个主要问题:(1)其中一些编码方案不够直观,(2)合并这些 ...
西门子 SCL 语言模拟量转换小例程
软硬件平台软件是博图 TIA v15,硬件是西门子300系列 PLC 简介使用scl语言编写一个FC子程序,对模拟量进行处理,对模拟数据进行转换. 详细介绍添加新块,添加一个FC块,语言选择sc ...
【生信】R语言进行id转换的方法（附可直接使用代码）
本文我都默认已经下载好了表达矩阵exp了哦代码都是直接给出来了,需要修改的地方我进行了标记一般只要修改一下都能直接用了方法一:下载平台数据以得到对应信息然后进入官网https://www.nc ...
R语言基础知识入门学习（一）
目录系列文章目录一.软件下载二.基本知识 1. 对象 2. 向量 3. 向量化 4. 因子总结系列文章目录 R语言基础知识入门学习(一) 一.软件下载我们可以通过这个网址对R语言软件进行下 ...
r语言boxcox异方差_基于R语言进行Box-Cox变换
原标题:基于R语言进行Box-Cox变换作者简介作者:吴健中国科学院大学 R语言.统计学爱好者,尤其擅长R语言和Arcgis在生态领域的应用分享个人公众号:统计与编程语言 Q: 为什么要进行Bo ...
R语言隐马尔可夫模型HMM识别股市变化分析报告
了解不同的市场状况如何影响您的策略表现可能会对您的收益产生巨大的影响.最近我们被客户要求撰写关于HMM的研究报告,包括一些图形和统计输出. 某些策略在波动剧烈的市场中表现良好,而其他策略则需要强劲而平 ...
R语言中实现马尔可夫链蒙特卡罗MCMC模型
什么是MCMC,什么时候使用它? MCMC只是一个从分布抽样的算法. 这只是众多算法之一.这个术语代表"马尔可夫链蒙特卡洛",因为它是一种使用"马尔可夫链"(我 ...
R语言与数据分析-01-数据类型与数据输入-01-基础
一. R语言的数据类型原文链接 : r语言与数据分析实战1 1.1 知识描述搞清楚我们面对的数据类型是什么,并能实现数据类型之间的转换 1.2 相关知识需要认真学习以下函数: mode(): 用 ...

R语言基因表达量转换（TPM、FPKM、RPKM）

R语言基因表达量转换（TPM、FPKM、RPKM）相关推荐

最新文章

热门文章