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")

#图的样子:

相关推荐
源猿人2 天前
企业级文件浏览系统的Vue实现:架构设计与最佳实践
前端·javascript·数据可视化
CodeCraft Studio2 天前
【案例分享】TeeChart 助力 Softdrill 提升油气钻井数据可视化能力
信息可视化·数据可视化·teechart·油气钻井·石油勘探数据·测井数据
lssjzmn3 天前
基于Spring Boot与Micrometer的系统参数监控指南
java·spring boot·数据可视化
czhc11400756633 天前
LINUX913 shell:set ip [lindex $argv 0],\r,send_user,spawn ssh root@ip “cat “
tcp/ip·r语言·ssh
zhangfeng11333 天前
win7 R 4.4.0和RStudio1.25的版本兼容性以及系统区域设置有关 导致Plots绘图面板被禁用,但是单独页面显示
开发语言·人工智能·r语言·生物信息
zhangfeng11333 天前
在 R 语言里,`$` 只有一个作用 按名字提取“列表型”对象里的单个元素 对象 $ 名字
开发语言·windows·r语言
高-老师3 天前
R语言生物群落(生态)数据统计分析与绘图实践技术应用
开发语言·r语言·生物群落
IT研究室4 天前
大数据毕业设计选题推荐-基于大数据的健康与生活方式数据可视化分析系统-Spark-Hadoop-Bigdata
大数据·hadoop·spark·毕业设计·源码·数据可视化·bigdata
@HNUSTer4 天前
Python数据可视化科技图表绘制系列教程(七)
python·数据可视化·科技论文·专业制图·科研图表
WangYan20224 天前
R语言:数据读取与重构、试验设计(RCB/BIB/正交/析因)、ggplot2高级绘图与统计检验(t检验/方差分析/PCA/聚类)
r语言·ggplot2·dplyr