immunedeconv:一个R软件包,用于统一访问从大量RNA测序数据估算免疫细胞组分的计算方法。

an R package for unified access to computational methods for estimating immune cell fractions from bulk RNA sequencing data.

1. 安装

conda create -n deconvolution
conda activate deconvolution
conda install -c bioconda -c conda-forge r-immunedeconv## 连接超时
#install.packages("remotes")
#remotes::install_github("icbi-lab/immunedeconv")#devtools::install_github('icbi-lab/immunedeconv', repos='http://cran.rstudio.com')

2. 准备表达谱数据

The input data is a gene ×× sample gene expression matrix. In general values should be

  • TPM-normalized
  • not log-transformed.
gene_expr <- read.table("/data/query_expr_ready.txt",sep='\t',header=T)head(gene_expr)colnames(gene_expr)
rownames(gene_expr)rownames(gene_expr) <- gene_expr$IDgene_expr <- gene_expr[,-1]
exprs_matrix <- as.matrix(gene_expr)
head(exprs_matrix)

3. 估算免疫细胞比例

library(tidyr)
library(immunedeconv)res = deconvolute(exprs_matrix, "quantiseq")res = deconvolute(exprs_matrix, "quantiseq") %>%map_result_to_celltypes(c("T cell CD4+"), "quantiseq")knitr::kable(res, digits=2)res1 <- deconvolute(your_mixture_matrix, "mcp_counter")res2 <- deconvolute(your_mixture_matrix, "xcell")res3 <- deconvolute(your_mixture_matrix, "epic")## method:
#quantiseq
#timer
#cibersort
#cibersort_abs
#mcp_counter
#xcell
#epic

参考:

https://github.com/icbi-lab/immunedeconv

immunedeconv估算免疫细胞比例相关推荐

  1. Python爬取京东回力鞋购买情况看看码数比例

    背景 突然想起前阵子看到一篇爬bar并统计女生罩杯情况的文章,心想.我擦,这么好的题材让别人先写了我风暴哭泣.心想,算了,人家本来就比我优秀.然后随便去京东逛了逛.看到回力旗舰店给我推了,那我也来搞搞 ...

  2. Cibersort免疫浸润的在线分析及R语言代码实现

    上期展示了ESITMATE(基于转录组数据)计算免疫得分和肿瘤纯度的一个例子,详见ggplot2实现分半小提琴图绘制基因表达谱和免疫得分.实际上计算肿瘤纯度的方法还有InfiniumPurify(基于 ...

  3. 政策定价风控审批策略

    政策与定价 一.不同应用场景下量化风控政策设定 市场主流产品 无定向用途贷款(信用贷款):2.定向用途贷款(商品贷款) 无定向用途贷款申请流程: 进件 审核 审核结果/授信 贷后管理 政策制定关注点: ...

  4. 新冠患者样本单细胞测序文献汇总

    科研工作者的信仰就是将真相大白于天下 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单 ...

  5. 2022年考研计算机-数据库原理8-11章

    文章链接 https://gitee.com/fakerlove/database-management-system 8. 数据库编程 SQL缺少流程控制能力,难以实现业务应用中的逻辑控制:SQL编 ...

  6. IF:5+ 7种癌症免疫治疗证明DNA损伤反应通路突变可作为免疫检查点阻断疗效的潜在生物标志物

    点击关注,桓峰基因 桓峰基因 生物信息分析,SCI文章撰写及生物信息基础知识学习:R语言学习,perl基础编程,linux系统命令,Python遇见更好的你 105篇原创内容 公众号 桓峰基因的教程不 ...

  7. 本周最新文献速递20220522

    本周最新文献速递20220522 一.精细解读文献 一 文献题目: Single-nucleus chromatin accessibility profiling highlights regula ...

  8. 免疫浸润计算方法是CIBERSORT和ssgsea 画图

    目前主流的免疫浸润计算方法是CIBERSORT和ssgsea,今天介绍CIBERSORT. 1.输入数据要什么 下面这段话摘自CIBERSORT的介绍 Importantly, all express ...

  9. 撬动PMBOK6运用指南第九章

    编者按:本文是撬动PMBOK6运用指南系列第九章,作者栾志强是达人圈资深网友,感谢作者对达人圈的大力支持!关注IT达人圈公众号并输入"PMBOK6",即可下载到PMBOK6管理指南 ...

  10. 重磅综述:精准精神病学机器学习的现代观点

    1. 全局 在美国和世界各地,心理健康问题是一种流行病,给医疗保健系统和社会带来了巨大的负担.迄今为止,仍然缺乏生物标志物和个体化治疗指南.近年来,机器学习(ML)和人工智能(AI)在分析精神病学的神 ...

最新文章

  1. RPA实施过程中可能会遇到的14个坑
  2. 软考之CPU的寻址方式
  3. 从上云到云管,只差一名青云“最佳小助手”
  4. NuGet Error:Unable to find version ‘0.3.3.0‘ of package ‘glog.overlay-x64_v120_Release_dynamic‘.
  5. 实现springboot的starter
  6. Linux开机启动过程(7):内核执行入口点
  7. angularjs1.x版本,父子组件之间的双向绑定
  8. 从RDS中获取binlog
  9. 从文本界面安装RHEL5操作系统详解
  10. 117 Populating Next Right Pointers in Each Node II
  11. sql server 标准版升级专业版或者开发版没有sql server profiler
  12. 自控专业工程设计用标准及规范
  13. 传智播客reactnative_传智播客黑马前端36期(2018年)
  14. 【ELIXIR】简单说下elixir的历史
  15. C语言实现乘法原理,伽罗华乘法原理与C语言实现
  16. python day2
  17. 软件实施工作个人看法
  18. ASP.NET mvcConf Videos Available
  19. ZYNQ空间地址详细配置
  20. Cubic spline(三次样条插值)(转载)

热门文章

  1. 谷歌翻译用不了解决办法
  2. Java常用设计模式及应用场景介绍
  3. 干货分享 ▎软考论文怎么写?
  4. UDS - 诊断协议和协议栈介绍
  5. 关于阿里云ICON矢量图SVG上传问题
  6. 账龄分析表excel模板_老会计分享财务报表及EXCEL会计报表分析系统模板!收藏领取!...
  7. sys caffe教程
  8. 前端响应式布局与Bootstrap栅格系统快速了解与应用
  9. 公司用的非标普通自动化用单片机还是plc_合格电气自动化工程师必备十大技能...
  10. 14_星仔带你学Java之Java编码规范、常用类