R包estimate评估肿瘤组织中基质及免疫细胞浸润水平

根据表达数据，ESTIMATE为研究人员提供了肿瘤纯度、存在的基质细胞水平和肿瘤组织中免疫细胞浸润水平的分数。https://bioinformatics.mdanderson.org/estimate/index.html

1. estimate安装

install.packages("estimate")rforge <- "http://r-forge.r-project.org"
install.packages("estimate", repos=rforge, dependencies=TRUE)
library(estimate)
ls("package:estimate") in.file <- system.file("extdata", "sample_input.txt", package="estimate")
out.file <- tempfile(pattern="estimate", fileext=".gct")

2. estimate分析

#进行estimate分析
setwd("/Users/test")
filterCommonGenes(in.file, output.f=out.file,id="GeneSymbol")estimateScore(input.ds = out.file,output.ds="OV_estimate_score.gct", platform="affymetrix")
# platform = c("affymetrix", "agilent", "illumina"))plotPurity(scores="OV_estimate_score.gct", platform="affymetrix",output.dir="estimated_purity_plots")
#生成plot文件

3. 累计分布曲线作图

scores=read.table("OV_estimate_score.gct",skip = 2,header = T)
rownames(scores)=scores[,1]
scores=t(scores[,3:ncol(scores)])
View(scores)
scores<- as.data.frame(scores)
scores$SampleID <- rownames(scores)save(scores,file = 'BRCA_estimate_score.rdata')
# m_score <- get(load('BRCA_estimate_score.rdata'))
class(scores)# median(scores[,'ImmuneScore'])## 画单条累计分布曲线
library(ggplot2)
p1<-ggplot(scores,aes(x=ImmuneScore)) +stat_ecdf(color = "green") +labs(y="accumulative propotion")p2<-ggplot(scores,aes(x=StromalScore)) +stat_ecdf(color = "red") +labs(y="accumulative propotion")## 画多条累计分布曲线
# 要先把数据框转化成长数据格式
library(reshape2)
scores2 <- scores[,-4] # 去掉不需要的列
dim(scores2)
long = melt(scores2, id=c("SampleID"),variable.name= 'Class', value.name = 'Value')
# 创建新的变量，含有分组信息
long$group <- rep(c(rep("g1",6),rep("g2",4)),3)
library(dplyr)
long<- mutate(long,group_score = paste(Class,group,sep="_"))## 长格式转成宽格式数据框
#wide <- dcast(long,SampleID~Class,value.var='Value')#作多条累计分布曲线，以不同颜色区分。
ggplot(long,aes(x=Value,color =Class,linetype=group)) +stat_ecdf(size =1) + # 线粗细labs(x="Score",y="Accumulative propotion",title="Accumulative Plotting") +theme(legend.position = c(0.15,0.7), #图例位置 ("none", "left", "right", "bottom", "top", or two-element numeric vector)legend.background =element_rect(fill = "white", colour = "grey50"),#图例背景panel.background=element_rect(fill = "grey95", # grey90colour = "black",size = 1), #画布背景颜色plot.title = element_text(hjust=0.5,size=16,vjust=0.5), #标题位置legend.text=element_text(size=10,colour='black'), #图例文字axis.text=element_text(size=8,colour="black"), #坐标轴文字axis.title.y = element_text(size = rel(1.3), angle = 90),#y坐标轴名称文字axis.title.x = element_text(size = rel(1.3)),#x坐标轴名称文字)

R包estimate评估肿瘤组织中基质及免疫细胞浸润水平相关推荐

Front Immunol 复现 | 4. 使用estimate包评估肿瘤纯度
前几天有同学问了一篇文章里的一个方法的实现,看了一下这篇文章除了qPCR验证基本都是纯生信,今天就试着来复现一下.随缘复现哈,如果阅读数据不好看的话,可能就放弃了,希望大家多多点赞.在看,转发支持. ...
使用EPIC预测肿瘤微环境中免疫细胞构成
欢迎关注"生信修炼手册"! 在传统的RNA_seq测序中,每个样本取样后实际包含了成千上万个细胞,和单细胞测序的single cell相比,这样的样本称之为bulk samples ...
学位论文精读-hBMSCs在肿瘤微环境中分泌IL-6并上调IL-17水平协同促进DLBCL生长的研究
摘要弥漫大B细胞淋巴瘤(DLBCL)是最常见的非霍奇金淋巴瘤亚型,是一种侵袭性淋巴瘤.虽然近年来一线的化疗方案R-CHOP(利妥昔单抗.环磷酰胺.多柔比星.长春新碱和强的松)改善DLBCL患者的预后 ...
生信分析之R语言常用R包一步下载
系列文章目录生信分析第一步:R语言基础应用以及数据前处理文章目录 R包下载使用GEOquery包下载原始数据芯片数据读取 GEOquery 下载并读取数据提取GEO表达矩阵提取GEO注释信 ...
16S预测宏基因组最强R包-Tax4Fun
之前在公众号的文章<根据16S预测微生物群落功能最全攻略>阅读人数近3000人,有需求的用户还是非常多的.其中提到了4个软件,之前已经介绍了其中非常有特点的三种,分别为: - PICRUS ...
R语言McSpatial_R语言中帮助和R包
#####获得帮助##### #关于函数的帮助 help(glm) ?glm #关于操作符的帮助 ?`+` #查看帮助文档中例子 example(glm) #搜索主题 help.search(&quo ...
R语言stringr包str_detect函数检测字符串中模式存在与否实战
R语言stringr包str_detect函数检测字符串中模式存在与否实战目录 R语言stringr包str_detect函数检测字符串中模式存在与否实战 #导入stringr包
R语言ggplot2可视化散点图并添加平滑曲线、WVPlots包BinaryYScatterPlot函数可视化一起中一个变量为二元离散变量的散点图（随着年龄变化是否有健康保险的客户的比例）
R语言ggplot2可视化散点图并添加平滑曲线.WVPlots包BinaryYScatterPlot函数可视化一起中一个变量为二元离散变量的散点图(随着年龄变化是否有健康保险的客户的比例) 目录
R语言观察日志(part13)--从GitHub中下载R包
学习笔记,仅供参考,有错必纠从GitHub中下载R包这里,我们将介绍两种下载并加载R包的方式,即通过在线安装和本地安装. 在线安装我们在使用在线安装的方法时,可能会遇到各种各样的问题,这里,我们 ...

R包estimate评估肿瘤组织中基质及免疫细胞浸润水平

1. estimate安装

2. estimate分析

3. 累计分布曲线作图

R包estimate评估肿瘤组织中基质及免疫细胞浸润水平相关推荐

最新文章

热门文章