【R>>IOBR】counts转TPM

分析测序数据时,常常需要将counts数据转换为TPM格式,而这个转变过程就需要涉及每个基因的长度,幸好有专业人士已经帮我们处理好这个东东,我们可以一键进行操作。
首先来认识下这个牛气冲天的R包IOBR(Immuno-Oncology Biological Research):

IOBR is an R package to perform comprehensive analysis of tumor microenvironment and signatures for immuno-oncology.


提供8种计算免疫细胞浸润的方法,收录255种构建的signature,这么强大的功能咱们以后慢慢学,这次先学习下count2tpm功能。

#devtools::install_github("IOBR/IOBR",ref="master")

rm(list = ls())
library(IOBR)
library(UCSCXenaTools)
help("count2tpm")

eset_prad<-XenaGenerate(subset = XenaCohorts =="GDC TCGA Prostate Cancer (PRAD)") %>% 
  XenaFilter(filterDatasets    = "TCGA-PRAD.htseq_counts.tsv") %>% 
  XenaQuery() %>%
  XenaDownload() %>% 
  XenaPrepare()
eset_prad$Ensembl_ID <- substring(eset_prad$Ensembl_ID,1,15)
eset_prad <- column_to_rownames(eset_prad,var = "Ensembl_ID")
eset_prad<-(2^eset_prad)-1
eset_prad <- count2tpm(countMat = eset_prad,idType = "Ensembl",source = "default")


参考链接:
IOBR: 一步完成RNAseq: counts到TPM的转化

最后编辑于
©著作权归作者所有,转载或内容合作请联系作者
【社区内容提示】社区部分内容疑似由AI辅助生成,浏览时请结合常识与多方信息审慎甄别。
平台声明:文章内容(如有图片或视频亦包括在内)由作者上传并发布,文章内容仅代表作者本人观点,简书系信息发布平台,仅提供信息存储服务。

相关阅读更多精彩内容

  • Bioinformatics Abstract The Gene Network Estimation Tool ...
    WH川川阅读 4,160评论 0 0
  • Web Resources Contents 1 Google DataSet Search Engine 2 G...
    生信义博阅读 5,382评论 0 4
  • Abstract MotivationSingle-cell RNA-sequencing (scRNA-seq)...
    丑小鸭_b360阅读 3,051评论 0 0
  • 表情是什么,我认为表情就是表现出来的情绪。表情可以传达很多信息。高兴了当然就笑了,难过就哭了。两者是相互影响密不可...
    Persistenc_6aea阅读 127,936评论 2 7
  • 16宿命:用概率思维提高你的胜算 以前的我是风险厌恶者,不喜欢去冒险,但是人生放弃了冒险,也就放弃了无数的可能。 ...
    yichen大刀阅读 11,318评论 0 4

友情链接更多精彩内容