公众号“生信小课堂”
TCGA数据分析课程:生物信息学教学
相关性分析是我们生信分析中必不可少的技能,单基因的批量相关性分析,可以用于做单基因的GO,KEGG富集分析和GSEA分析,也有多基因之间的相关性分析,或者多个基因集之间的相关性分析。今天我们利用TCGA肝癌数据,挑选一些基因做相关性分析,并用不同的R包展示。
load(file = "mRNAdata.Rda")#加载数据
library(tidyverse)
library(corrplot)
library(circlize)TCGA数据下载与ID转换mp.weixin.qq.com
mRNAdata
mRNAdata" />
选择列,我们随便选一些基因
new
new
计算相关性系数
cor_new
cor_new
corrplot包画图
method:指定形状,可以是circle圆形(默认),square方形,ellipse,椭圆形,number数值,shade阴影,color颜色,pie饼图。type:指定显示范围,full完全(默认),lower下三角,upper上三角
corrplot(cor_new, method = "circle")
corrplot(cor_new, method = "square")
corrplot(cor_new, method = "pie")
corrplot(cor_new, method = "color")
corrplot(cor_new, method = "number")
更改颜色
corrplot(cor_new, method = "color", col = colorRampPalette(c("