immunedeconv估算免疫细胞比例
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估算免疫细胞比例相关推荐
- Python爬取京东回力鞋购买情况看看码数比例
背景 突然想起前阵子看到一篇爬bar并统计女生罩杯情况的文章,心想.我擦,这么好的题材让别人先写了我风暴哭泣.心想,算了,人家本来就比我优秀.然后随便去京东逛了逛.看到回力旗舰店给我推了,那我也来搞搞 ...
- Cibersort免疫浸润的在线分析及R语言代码实现
上期展示了ESITMATE(基于转录组数据)计算免疫得分和肿瘤纯度的一个例子,详见ggplot2实现分半小提琴图绘制基因表达谱和免疫得分.实际上计算肿瘤纯度的方法还有InfiniumPurify(基于 ...
- 政策定价风控审批策略
政策与定价 一.不同应用场景下量化风控政策设定 市场主流产品 无定向用途贷款(信用贷款):2.定向用途贷款(商品贷款) 无定向用途贷款申请流程: 进件 审核 审核结果/授信 贷后管理 政策制定关注点: ...
- 新冠患者样本单细胞测序文献汇总
科研工作者的信仰就是将真相大白于天下 NGS系列文章包括NGS基础.转录组分析 (Nature重磅综述|关于RNA-seq你想知道的全在这).ChIP-seq分析 (ChIP-seq基本分析流程).单 ...
- 2022年考研计算机-数据库原理8-11章
文章链接 https://gitee.com/fakerlove/database-management-system 8. 数据库编程 SQL缺少流程控制能力,难以实现业务应用中的逻辑控制:SQL编 ...
- IF:5+ 7种癌症免疫治疗证明DNA损伤反应通路突变可作为免疫检查点阻断疗效的潜在生物标志物
点击关注,桓峰基因 桓峰基因 生物信息分析,SCI文章撰写及生物信息基础知识学习:R语言学习,perl基础编程,linux系统命令,Python遇见更好的你 105篇原创内容 公众号 桓峰基因的教程不 ...
- 本周最新文献速递20220522
本周最新文献速递20220522 一.精细解读文献 一 文献题目: Single-nucleus chromatin accessibility profiling highlights regula ...
- 免疫浸润计算方法是CIBERSORT和ssgsea 画图
目前主流的免疫浸润计算方法是CIBERSORT和ssgsea,今天介绍CIBERSORT. 1.输入数据要什么 下面这段话摘自CIBERSORT的介绍 Importantly, all express ...
- 撬动PMBOK6运用指南第九章
编者按:本文是撬动PMBOK6运用指南系列第九章,作者栾志强是达人圈资深网友,感谢作者对达人圈的大力支持!关注IT达人圈公众号并输入"PMBOK6",即可下载到PMBOK6管理指南 ...
- 重磅综述:精准精神病学机器学习的现代观点
1. 全局 在美国和世界各地,心理健康问题是一种流行病,给医疗保健系统和社会带来了巨大的负担.迄今为止,仍然缺乏生物标志物和个体化治疗指南.近年来,机器学习(ML)和人工智能(AI)在分析精神病学的神 ...
最新文章
- RPA实施过程中可能会遇到的14个坑
- 软考之CPU的寻址方式
- 从上云到云管,只差一名青云“最佳小助手”
- NuGet Error:Unable to find version ‘0.3.3.0‘ of package ‘glog.overlay-x64_v120_Release_dynamic‘.
- 实现springboot的starter
- Linux开机启动过程(7):内核执行入口点
- angularjs1.x版本,父子组件之间的双向绑定
- 从RDS中获取binlog
- 从文本界面安装RHEL5操作系统详解
- 117 Populating Next Right Pointers in Each Node II
- sql server 标准版升级专业版或者开发版没有sql server profiler
- 自控专业工程设计用标准及规范
- 传智播客reactnative_传智播客黑马前端36期(2018年)
- 【ELIXIR】简单说下elixir的历史
- C语言实现乘法原理,伽罗华乘法原理与C语言实现
- python day2
- 软件实施工作个人看法
- ASP.NET mvcConf Videos Available
- ZYNQ空间地址详细配置
- Cubic spline(三次样条插值)(转载)
热门文章
- 谷歌翻译用不了解决办法
- Java常用设计模式及应用场景介绍
- 干货分享 ▎软考论文怎么写?
- UDS - 诊断协议和协议栈介绍
- 关于阿里云ICON矢量图SVG上传问题
- 账龄分析表excel模板_老会计分享财务报表及EXCEL会计报表分析系统模板!收藏领取!...
- sys caffe教程
- 前端响应式布局与Bootstrap栅格系统快速了解与应用
- 公司用的非标普通自动化用单片机还是plc_合格电气自动化工程师必备十大技能...
- 14_星仔带你学Java之Java编码规范、常用类