上回咱们介绍了多基因之间的相关性分析及绘图
R绘图:相关性分析与作图
今天我们来讲一下,单基因之间得相关性分析及作图,这个在实验性得文章中用到得更广泛
首先还是加载数据,数据来源于TCGA数据下载与ID转换
rm(list=ls())library(ggstatsplot)#加载包load(file = "mRNAdata.Rda")#加载数据
将数据转置,并转换成数据框,取log
dat<-as.data.frame(t(mRNAdata))dat<-log2(dat+1)
我们用LDHA和PKM来示范绘图,方法选择pearson,显示平均值
ggscatterstats(dat, y =LDHA, x =PKM,