Front Immunol 复现 | 4. 使用estimate包评估肿瘤纯度
前几天有同学问了一篇文章里的一个方法的实现,看了一下这篇文章除了
qPCR
验证基本都是纯生信,今天就试着来复现一下。随缘复现哈,如果阅读数据不好看的话,可能就放弃了,希望大家多多点赞
、在看
,转发
支持。
文章标题:Investigation of a Hypoxia-Immune-Related Microenvironment Gene Signature and Prediction Model for Idiopathic Pulmonary Fibrosis.
doi: 10.3389/fimmu.2021.629854
注:因为只是浅复现一下,个人水平也有限,过程中可能曲解作者的原意,难免存在纰漏或瑕疵,希望大家批评指正。
流程
本篇推文所在分析流程位置如图。
示例数据和代码领取
点赞
、在看
本文,分享至朋友圈集赞20个
并保留30分钟
,截图发至微信mzbj0002
领取。
木舟笔记2022年度VIP可免费领取。
木舟笔记2022年度VIP企划
权益:
2022年度木舟笔记所有推文示例数据及代码(在VIP群里实时更新)。
木舟笔记科研交流群。
半价购买
跟着Cell学作图系列合集
(免费教程+代码领取)|跟着Cell学作图系列合集。
收费:
99¥/人。可添加微信:mzbj0002
转账,或直接在文末打赏。
方法简析
ESTIMATE为研究人员提供了基于表达数据的肿瘤纯度、基质细胞存在水平和肿瘤组织中免疫细胞浸润水平的评分。
浸润基质细胞和免疫细胞是肿瘤组织中正常细胞的主要组成部分,不仅在分子研究中干扰肿瘤信号,而且在肿瘤生物学中具有重要作用。使用表达数据估计恶性肿瘤中的基质细胞和免疫细胞(Estimation of STromal and Immune cells in MAlignant Tumours using Expression data,ESTIMATE)是一种使用基因表达特征来推断肿瘤样本中基质细胞和免疫细胞的比例的方法。也就是说,基质细胞和免疫细胞含量多,那肿瘤纯度低,反之肿瘤纯度就高了。
开始分析
数据格式
# estimate下载
getOption('timeout')
options(timeout=1000)
library(utils)
rforge <- "http://r-forge.r-project.org"
install.packages("estimate", repos=rforge, dependencies=TRUE)
library(estimate)
?estimate
library(help="estimate")
# 评分
filterCommonGenes(input.f="combat_exp.txt", output.f="exp_estimate_score.gct", id = "GeneSymbol")
estimateScore('exp_estimate_score.gct','exp_purify.gct' )
result <- read.table('exp_purify.gct',header = T,skip = 2)
head(result)
> head(result)NAME Description GSM1820739 GSM1820740 GSM1820741 GSM1820742
1 StromalScore StromalScore -475.990263 -345.9680632 -390.8492694 -223.2208710
2 ImmuneScore ImmuneScore 2543.250642 2473.3975436 2732.0745355 2760.4281562
3 ESTIMATEScore ESTIMATEScore 2067.260379 2127.4294804 2341.2252661 2537.2072852
4 TumorPurity TumorPurity 0.614979 0.6079905 0.5827791 0.5591635
系列合集
Front Immunol 复现 | 1. GEO数据下载及sva批次校正(PCA可视化)
Front Immunol 复现 | 2. 一个基于缺氧基因集的数据降维聚类分组方法(umap,MSigDB)
Front Immunol 复现 | 3. 基因表达差异分析(limma)
往期内容
(免费教程+代码领取)|跟着Cell学作图系列合集
Q&A | 如何在论文中画出漂亮的插图?
R绘图 | 气泡散点图+拟合曲线
跟着 Cell 学作图 | 桑葚图(ggalluvial)
R绘图 | 对比条形图+连线
R绘图 | 一幅小提琴图的美化之旅
Front Immunol 复现 | 4. 使用estimate包评估肿瘤纯度相关推荐
- ESTIMATE包计算肿瘤纯度
介绍 肿瘤组织中的正常细胞不仅在分子研究中影响肿瘤信号,而且在癌症生物学中也起着重要作用. 估计包使用基因表达数据预测肿瘤组织中基质细胞和免疫细胞的存在. 示例 首先使用通过Affymetrix U13 ...
- Front Immunol 复现 | 1. GEO数据下载及sva批次校正(PCA可视化)
FII 前几天有同学问了一篇文章里的一个方法的实现,看了一下这篇文章除了qPCR验证基本都是纯生信,今天就试着来复现一下.随缘复现哈,如果阅读数据不好看的话,可能就放弃了,希望大家多多点赞.在看,转发 ...
- ESTIMATE 包 error/报错 无法匹配基因
使用 R包 ESTIMATE 在进行肿瘤纯度分析时 出现: [1] "Merged dataset includes 0 genes (10412 mismatched)." 只要 ...
- 【生物信息】ESTIMATE 分析免疫评分和肿瘤纯度
ESTIMATE 分析免疫评分和肿瘤纯度 背景: 从病人身上获取的肿瘤组织不仅仅包括包括肿瘤细胞,还包括与肿瘤微环境相关的很多其他细胞,比如正常的基质细胞,免疫细胞.血管内皮和血管内的血细胞. 其中比 ...
- R包estimate评估肿瘤组织中基质及免疫细胞浸润水平
根据表达数据,ESTIMATE为研究人员提供了肿瘤纯度.存在的基质细胞水平和肿瘤组织中免疫细胞浸润水平的分数.https://bioinformatics.mdanderson.org/estimat ...
- python 复现java中jasypt包的 PBEWITHSHA256AND256BITAES-CBC-BC加解密算法(可java python双向加解密)
运行环境:python 3.6 依赖包Crypto的安装: pip3 install pycryptodome # Make coding more python3-ish from __future ...
- 数据挖掘肿瘤预测_Nature Medicine封面文章:利用单核细胞数量预测及评估肿瘤免疫治疗效果...
免疫检查点阻断已经彻底改变了癌症治疗,临床实验数据表明PD-1.PD-L1抗体的免疫治疗可以有效应对转移性黑色素瘤和多种其他类型的癌症.尽管无进展生存期显着增加,临床治疗结果也呈现很大的个体差异,只有 ...
- 国庆特惠 !| CNS图表复现|生信分析|R绘图 资源分享讨论群!
cover ❝ Q:群里有哪些资源? A:2022.12.31前木舟笔记公众号更新的所有资源.(具体目录详见下方) Q:2022年都快结束了,现在加群不是亏了? A:无论什么时候加群,拿到的资源都是一 ...
- R实战 | NGS数据时间序列分析(maSigPro)
masigpro 跟着Cell学作图 | 6.时间序列分析(Mfuzz包) 一个答疑教程. maSigPro 流程 示例数据 #BiocManager::install('maSigPro') lib ...
最新文章
- 用T-SQL得到数据库的可视化结构
- 关于mtl_transactions_interface表序列
- 每日一天:对象数组去重
- js打开本地文件夹_vue + ArcGIS 地图应用系列一:arcgis api本地部署(开发环境)
- No.7 nginx、nfs、集群
- 前端学习(571):margin负值下的两栏自适应
- LeetCode 1169. 查询无效交易
- bzoj:1692 [Usaco2007 Dec]队列变换1640 [Usaco2007 Nov]Best Cow Line 队列变换
- 窗口操作-关闭,最小化
- 如果计算机原理程序设计,计算机组成原理程序设计.doc
- e5cc温控仪通讯参数设定_应用 | 如何实现S7300与S7200smart通讯?
- 【NLP开发】Python实现聊天机器人(微软小冰)
- 与其埋头啃文献不如关注这些公众号
- OpenWrt路由器通过Aria2开启迅雷离线下载[图文教程]
- C语言学习--------修改控制台输出文字颜色
- 数码照片像素与可冲印照片尺寸对照表
- 【全栈编程系列】SpringBoot整合Shiro(含KickoutSessionControlFilter并发在线人数控制以及不生效问题、配置启动异常No SecurityManager...)
- 利用WireShark分析H264码流
- GBase 8c技术特点
- nginx配置域名访问/禁止ip访问
热门文章
- mysql 给表添加唯一约束、联合唯一约束,指定唯一约束的名字
- Jmeter 4.0自带代理接口录制脚本
- rf 遍历列表_RF的变量list在For循环的用法,试错中学习
- 敏捷遇上UML(深圳 2014.3.15)- 活动报道
- mysql unix_timestamp now() dfdfd f_数据库函数lpad的搜索结果-阿里云开发者社区
- 法规遵从战略:晋升内部IT审计师的策略
- superIO通过PS2接口扩展键盘
- 松下6轴程序模板 plc采用FP-XHC60T ,标准可带6轴程序
- DBA平均月薪17000,入职3年感叹这份工作实在是太难了!
- 流氓软件卸载了又偷偷冒出来,dllhost.exe暗藏安装玄机