ggcor包做相关性分析

ggcor包是很难安装的,直接使用install.package()和使用Bioc是无法安装的,在R的镜像中也是不存在的,很多人都是将上传到GitHub中。如果你运气好,那么一次就可以安装,但如果你没有一点运气,那就需要花费很长的时间。小果自己安转就花费了很多的时间,如果你按照小果的方法安装,就会省去你很多不必要的麻烦哦。

ggcor包介绍:ggcorr: correlation matrixes with ggplot2 (briatte.github.io)

下载安装

首先下载ggcor包到本地,

本文相关代码及数据都已传到百度云网盘:

链接:https://pan.baidu.com/s/1K9XGIaIiRfTdj2dJgLaMyw?pwd=ea0l

提取码:ea0l

###解压ggcor_master.zip,切换到文件夹目录,安装

devtools::install_local(“ggcor_master”, force = TRUE)

下面是ggcor安装的具体实列及代码:

#使用ggcor绘制图形

library(vegan)#加载vegan包

library(ggcor)

library(ggplot2)

data(“varechem”)#加载vegan内置数据-理化数据

data(“varespec”)#加载vegan内置数据-物种数据

#数据准备

head(varechem)###查看数据

head(varespec)

write.csv(varechem,file=”varechem.csv”) ##保存数据

write.csv(varespec,file=”varespec.csv”)

mantel <- mantel_test(varespec, varechem,

spec.select = list(Spec01 = 1:7,#依次定义四种物种作为Mantel的分析对象

Spec02 = 8:18,

Spec03 = 19:37,

Spec04 = 38:44)) %>%

mutate(rd = cut(r, breaks = c(-Inf, 0.2, 0.4, Inf),

labels = c(“< 0.2”, “0.2 – 0.4”, “>= 0.4”)),#定义Mantel的R值范围标签,便于出图

pd = cut(p.value, breaks = c(-Inf, 0.01, 0.05, Inf),

labels = c(“< 0.01”, “0.01 – 0.05”, “>= 0.05”)))#定义Mantel检验的p值范围标签,便于出图

pdf(“corrplot.pdf”,10,10)

quickcor(varechem, type = “upper”) +#绘制理化数据热图

geom_square() +#定义成方块状

anno_link(aes(colour = pd, size = rd), data = mantel) +#定义连线

scale_size_manual(values = c(0.5, 1, 2))+

guides(size = guide_legend(title = “Mantel’s r”,#定义图例

order = 2),

colour = guide_legend(title = “Mantel’s p”,

order = 3),

fill = guide_colorbar(title = “Pearson’s r”, order = 4))

dev.off()

结果图片:

​欢迎关注”生信果”、“云生信学生物信息学”、“尔云间meta分析”公众号,以及知乎、BiliBili等主要发表或收录生物信息学的教程,以及基于R的分析和可视化等内容,一起见证小白和大佬的成长。