linkET包繪制mantal test 相關性熱圖

今天小雲想學習利用linkET包來繪制mantel test相關性熱圖,代碼如下:

01

安裝需要的R包

install.packages(“ggplot2”)install.packages(“tidyverse”)install.packages(“devtools”)devtools::install_github("Hy4m/linkET", force = TRUE)

02

載入需要的R包

library(devtools)library(ggplot2)library(tidyverse)library(linkET)

03

代碼展示

#讀取數據varechem<-read.csv(file="varechem.csv",row.names=1)

varespec<-read.csv(file="varespec.csv",row.names=1)

b228284c8c6807dac506eab34e7d96da

# mantel testmantel <- mantel_test(varespec, varechem, # 傳入需要檢測的兩個矩陣 mantel_fun = 'mantel', spec_select = list(Spec01 = 1:7, Spec02 = 8:18, Spec03 = 19:37, Spec04 = 38:44)) %>% mutate(r = cut(r, breaks = c(-Inf, 0.2, 0.4, Inf), labels = c("< 0.2", "0.2 - 0.4", ">= 0.4")), p = cut(p, breaks = c(-Inf, 0.01, 0.05, Inf), labels = c("< 0.01", "0.01 - 0.05", ">= 0.05")))#計算環境因子相關性系數:cor2 <- correlate(varechem)

#繪制相關性熱圖p1<- qcorrplot(cor2,grid_col = "grey50",grid_size = 0.2,type = "upper",diag = FALSE) +geom_square() +scale_fill_gradientn(colours = c("#610214", "#d05646", "#f5f4f4", "#569cc7", "#0b3b71"),limits = c(-1, 1))

#添加顯著性標簽:p2 <- p1+geom_mark(size = 4,only_mark = T,sig_level = c(0.05, 0.01, 0.001),sig_thres = 0.05,colour = 'white')

#在相關性熱圖上添加mantel連線:p3<- p2+geom_couple(data = mantel,aes(colour = p, size = r),curvature = nice_curvature())+scale_size_manual(values = c(0.5, 1.2, 2)) + #連線粗細scale_colour_manual(values = c("#d85c01", "#29d300", "#A2A2A288")) + #連線配色#修改圖例:guides(size = guide_legend(title = "Mantel r",override.aes = list(colour = "grey35"),order = 2),colour = guide_legend(title = "Mantel p",override.aes = list(size = 3),order = 1),fill = guide_colorbar(title = "Pearson r", order = 3))

d64555bae25c8f83fa169f08532079a4

小雲感覺效果還不錯,該包與ggplot2語法兼容,在繪圖的同時可以學習相關的ggplot2的繪圖知識,值得推薦,今天小雲的分享就到這裡瞭,有需要的可以借鑒學習,下期再見。

選分析方向可以從最新發表文獻入手,當然最便捷的方法當然還是持續關註小雲的生信思路推薦,想定制創新性思路歡迎直接call小雲,雲生信團隊竭誠為您的科研助力!

发表回复

相关推荐

小米手機怎樣熄屏顯示時間?教你三種玩法,漲知識瞭

小米手機怎樣熄屏顯示時間?現在很多主流的手機都支持“熄屏顯示”功能,不僅可以在熄屏狀態下顯示時間,而且還能熄屏顯示圖片...

· 2秒前

2022年中山市永安中學招生簡章

2022年中山市永安中學招生簡章—、總則招生依據:根據《中山市2022年高中階段學校考試招生工 作方案》(中教體通〔 2022〕84號...

· 8分钟前

常用除草剂注意事项大汇总!

乙草胺 乙草胺内吸性酰胺类除草剂,是选择性芽前除草剂。可被植物幼芽吸收,单子叶植物通过芽鞘吸收,双子叶植物下胚轴吸收 ...

· 10分钟前

股東是否實繳出資如何查詢

如何確認實繳?1.可以去企業所在地工商局網站查詢。2.看章程。3.看看企業銀行的流水。目前一般公司成立時,不需要預先實繳資...

· 16分钟前

筆記本電腦怎麼看參數?小白必看的筆記本參數

如今的筆記本電腦市場,品牌種類繁多,對於一些不懂電腦參數的小白該如何挑選呢?一般萌妹看筆記本電腦時,對於亂七八糟參數...

· 24分钟前