corrplot：R 面板圖中的 corrplot 函數

Question

我想創建一個包含由 corrplot() 創建的不同相關圖的圖形。 為此，一旦我創建了圖，我將它們分配給一個變量並將它們與 ggpubr:ggarrange() 對齊，但是創建的單個變量的內容為 NULL，並且當齒輪范圍（正確）創建圖形時，內容為空.

library(corrplot)
library(psych)
library(ggpubr)

data(iris)

res_pearson.c_setosa<-iris%>%
  filter(Species=="setosa")%>%
  select(Sepal.Length:Petal.Width)%>%
  corr.test(., y = NULL, use = "complete",method="pearson",adjust="bonferroni", alpha=.05,ci=TRUE,minlength=5)

corr.a<-corrplot(res_pearson.c_setosa$r[,1:3],
         type="lower", 
         order="original", 
         p.mat = res_pearson.c_setosa$p[,1:3], 
         sig.level = 0.05, 
         insig = "blank", 
         col=col4(10), 
         tl.pos = "ld",
         tl.cex = .8, 
         tl.srt=45, 
         tl.col = "black",
         cl.cex = .8)+
  my.theme


res_pearson.c_virginica<-iris%>%
  filter(Species=="virginica")%>%
  select(Sepal.Length:Petal.Width)%>%
  corr.test(., y = NULL, use = "complete",method="pearson",adjust="bonferroni", alpha=.05,ci=TRUE,minlength=5)

corr.b<-corrplot(res_pearson.c_virginica$r[,1:3],
         type="lower", 
         order="original", 
         p.mat = res_pearson.c_virginica$p[,1:3], 
         sig.level = 0.05, 
         insig = "blank", 
         col=col4(10), 
         tl.pos = "ld",
         tl.cex = .8, 
         tl.srt=45, 
         tl.col = "black",
         cl.cex = .8)+
  my.theme

ggarrange(corr.a, corr.b,
          common.legend = TRUE,
          legend = "bottom",
          ncol = 9, nrow = 1)

在網上搜索了幾次之后，似乎這是第一次有人嘗試創建具有多個相關圖的圖形。

Answer 1

您可以使用此代碼par(mfrow=c(1,2))並排繪制圖形，如下所示：

library(corrplot)
library(psych)
library(ggpubr)
library(dplyr)

data(iris)

res_pearson.c_setosa<-iris%>%
  filter(Species=="setosa")%>%
  select(Sepal.Length:Petal.Width)%>%
  corr.test(., y = NULL, use = "complete",method="pearson",adjust="bonferroni", alpha=.05,ci=TRUE,minlength=5)

par(mfrow=c(1,2))

corr.a<-corrplot(res_pearson.c_setosa$r[,1:3],
                 type="lower", 
                 order="original", 
                 p.mat = res_pearson.c_setosa$p[,1:3], 
                 sig.level = 0.05, 
                 insig = "blank", 
                 #col=col4(10), 
                 tl.pos = "ld",
                 tl.cex = .8, 
                 tl.srt=45, 
                 tl.col = "black",
                 cl.cex = .8)


res_pearson.c_virginica<-iris%>%
  filter(Species=="virginica")%>%
  select(Sepal.Length:Petal.Width)%>%
  corr.test(., y = NULL, use = "complete",method="pearson",adjust="bonferroni", alpha=.05,ci=TRUE,minlength=5)

corr.b<-corrplot(res_pearson.c_virginica$r[,1:3],
                 type="lower", 
                 order="original", 
                 p.mat = res_pearson.c_virginica$p[,1:3], 
                 sig.level = 0.05, 
                 insig = "blank", 
                 #col=col4(10), 
                 tl.pos = "ld",
                 tl.cex = .8, 
                 tl.srt=45, 
                 tl.col = "black",
                 cl.cex = .8)

^{由reprex 包(v2.0.1) 於 2022-07-13 創建}