scRNA-data中的R值

愿武艺晴小朋友一定得每天都开心


当我们测序拿得到各个样本中基因的表达值,就可以用基因表达值来表征样本间的相关性

代码如下:

#样本间相似性:R值 相关性 捕获到的基因在两个样本间表达趋势一致性

exp_RNA <- AverageExpression(fasting_memory,

group.by = "Sample",layer = "data") #CPM值来自data图层

exp_RNA <- as.data.frame(exp_RNA)

colnames(exp_RNA) <- c("fed","health","memory_10d","memory_35d","memory_66d")

library(ArchR)

library(viridis)

head(exp_RNA)

df<-exp_RNA[,c(1,5)] #依次计算各个组

head(df)

#为了提高数据质量和准确性,使用两组间表达值都非0的基因用于R值的计算

df<-subset.data.frame(df,df$fed!=0)

df<-subset.data.frame(df,df$memory_66d!=0)

cor(df[,2],df[,1])

library(ggrepel)

df$gene <- rownames(df)

dfslope \<- dfmemory_66d/df$fed #斜率代表在66d组中跟fed组间的表达差别很大

head(df)

label <- subset.data.frame(df,df$slope>1000)

head(label)

ggPoint(x = dffed,y = dfmemory_66d,size=1,

title = "r=0.41",

colorDensity = TRUE,

continuousSet = "solarExtra",

ylabel = "memory_66d:log2(CPM+1)",

xlabel = "fed:log2(CPM+1)",

xlim = c(0,170),

ylim = c(0,170))+ mytheme+

geom_hline(yintercept = 40, lty = "dashed")+

geom_vline(xintercept = 40, lty = "dashed")

#图的样子:

相关推荐
爱分享的飘哥20 小时前
第十九篇 自动化报表生成:Python一键生成可视化Excel图表与专业PDF报告,老板看了都点赞!
自动化·办公自动化·数据可视化·excel自动化·python报表·pdf报告
医工交叉实验工坊20 小时前
R 语言绘制 10 种精美火山图:转录组差异基因可视化
python·信息可视化·r语言
银之夏雪丶20 小时前
AntV G6 基础元素详解(React版)
前端·javascript·数据可视化
智算菩萨20 小时前
传统机器学习在信用卡交易预测中的卓越表现:从R²=-0.0075到1.0000的华丽转身
人工智能·机器学习·r语言
大数据魔法师1 天前
Matplotlib(一)- 数据可视化与Matplotlib
matplotlib·数据可视化
2201_753054892 天前
应用回归分析,R语言,多元线性回归总结(下)
回归·r语言·线性回归
李昊哲小课3 天前
销售数据可视化分析项目
python·信息可视化·数据分析·matplotlib·数据可视化·seaborn
LabEx5 天前
科研数据可视化核心技术:基于 AI 与 R 语言的热图、火山图及网络图绘制实践指南
人工智能·信息可视化·r语言·r语言绘图·乐备实·labex·科研数据绘图
Jet45055 天前
第100+43步 ChatGPT学习:R语言实现特征选择曲线图
学习·chatgpt·r语言
Chef_Chen5 天前
从0开始学习R语言--Day40--Kruskal-Wallis检验
开发语言·学习·r语言