linkET包可视化WGCNA模块与代谢物相关性
2022-11-22 本文已影响0人
队长的生物实验室
直接上图
当需要可视化WGCNA结果中某些模块与代谢物的相关性,就需要考虑图形的直观与美观,linkET包比较好的满足了这一要求。
输入文件
df1
各模块的表达量
df2
各代谢物在不同处理下的表达量
devtools::install_github("Hy4m/linkET", force = TRUE)
library(RColorBrewer)
library(linkET)
df1 <- read.csv("cor_GM_exp.csv",header=T,row.names = 1)
df2 <- read.csv("cor_Metab_exp.csv",header=T,row.names = 1)
correlate(df2,df1)%>%
qcorrplot()+
geom_square()+
scale_fill_gradientn(colours=RColorBrewer::brewer.pal(11,"RdBu"))#利用RColorBrewer包自行定义颜色