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

#图的样子:

相关推荐
alsmile4 小时前
什么是组态?不会代码也能搭出工业大屏和数字可视化
低代码·数据可视化
IT研究室5 小时前
大数据毕业设计选题推荐-基于大数据的国家药品采集药品数据可视化分析系统-Spark-Hadoop-Bigdata
大数据·hadoop·信息可视化·spark·毕业设计·数据可视化·bigdata
q5673152310 小时前
告别低效:构建健壮R爬虫的工程思维
开发语言·爬虫·r语言
没有梦想的咸鱼185-1037-166317 小时前
【高分论文密码】大尺度空间模拟预测与数字制图
信息可视化·数据分析·r语言
IT研究室1 天前
大数据毕业设计选题推荐-基于大数据的贵州茅台股票数据分析系统-Spark-Hadoop-Bigdata
大数据·hadoop·spark·毕业设计·源码·数据可视化·bigdata
IT毕设梦工厂1 天前
大数据毕业设计选题推荐-基于大数据的国家基站整点数据分析系统-Hadoop-Spark-数据可视化-BigData
大数据·hadoop·spark·毕业设计·源码·数据可视化
kebeiovo1 天前
项目必备流程图,类图,E-R图实例速通
开发语言·r语言·流程图
@HNUSTer2 天前
Python数据可视化科技图表绘制系列教程(六)
python·数据可视化·科技论文·专业制图·科研图表
界面开发小八哥2 天前
数据可视化图表库LightningChart JS v8.0上线:全新图例系统 + 数据集重构
javascript·信息可视化·数据可视化·lightningchart
zhangfeng11332 天前
BiocManager下载失败 R语言 解决办法
开发语言·r语言