R corrplot：Plot 相关系数以及显着性星？

Question

使用 R corrplot，我还没有找到将方框中的相关系数与其重要性一起绘制的解决方案，即 0.84*** 这是仅绘制重要性星的代码。 如何在那里添加相关系数？

M<-cor(mtcars)
res1 <- cor.mtest(mtcars, conf.level = .95)
corrplot(cor(mtcars),
     method="square",
     type="lower",
     p.mat = res1$p,
     insig = "label_sig",
     sig.level = c(.001, .01, .05),
     pch.cex = 0.8,
     pch.col = "red",
     tl.col="black",
     tl.cex=1,
     outline=TRUE)

如果我按照第一个答案的建议添加 addCoef.col = "black"，则文本会覆盖重要星，因此它们不再被真正看到：

Answer 1

您只需将选项addCoef.col = "black"添加到corrplot

Answer 2

为了解决星星的颜色冲突，我采用了这种方法来制作更大的白色星星。 我将它与以下颜色和布局一起使用。 我无法避免重叠，但它足够清晰，我可以使用：

cex.before <- par("cex")
par(cex = 0.7)
col <- colorRampPalette(c("#BB4444", "#EE9988", "#FFFFFF", "#77AADD", "#4477AA"))
corrplot::corrplot(cor(mtcars), 
    method="color", 
    col=col(200),  
    type="lower", 
    # Combine with significance
    p.mat = res1$p, 
    insig = "label_sig",
    sig.level = c(.001, .01, .05), 
    pch.cex = 3, # Increase size of stars
    pch.col = "white", # Colour of stars
    # hide correlation coefficient on the principal diagonal
    diag=FALSE,
    addCoef.col = "black", # Add coefficient of correlation
    tl.col="black", tl.srt=45, #Text label color and rotation
    tl.cex = 1/par("cex"), cl.cex = 1/par("cex") #Reduce text size of coefficients               
)
par(cex = cex.before)

Answer 3

重要星的 position 由place_points function 中的corrplot function 定义。

问题：

如果两者都应该显示相关系数和显着性水平，它们就会重叠（我把星星选为黄色，因为我有一些色觉问题......）。

library(corrplot)
#> corrplot 0.90 loaded

M<-cor(mtcars)
res1 <- cor.mtest(mtcars, conf.level = .95)

corrplot(cor(mtcars),
         method="square",
         type="lower",
         p.mat = res1$p,
         insig = "label_sig",
         sig.level = c(.001, .01, .05),
         pch.cex = 0.8,
         pch.col = "yellow",
         tl.col="black",
         tl.cex=1,
         addCoef.col = "black",
         tl.pos="n",
         outline=TRUE)

^{由代表 package (v2.0.1) 于 2021 年 10 月 13 日创建}

快速且临时（每次新加载 corrplot 包时都必须重新执行此步骤）解决方案：

更改place_points function 中的corrplot function。 为此，请运行：

trace(corrplot, edit=TRUE)

然后在第 443 行替换

place_points = function(sig.locs, point) {
  text(pos.pNew[, 1][sig.locs], pos.pNew[, 2][sig.locs], 
       labels = point, col = pch.col, cex = pch.cex, 
       lwd = 2)

和：

# adjust text(X,Y ...) according to your needs, here +0.25 is added to the Y-position    
place_points = function(sig.locs, point) {
      text(pos.pNew[, 1][sig.locs], (pos.pNew[, 2][sig.locs])+0.25, 
           labels = point, col = pch.col, cex = pch.cex, 
           lwd = 2)

然后点击“保存”按钮。

结果：

library(corrplot)
#> corrplot 0.90 loaded

#change the corrplot function as described above 
trace(corrplot, edit=TRUE)
#> Tracing function "corrplot" in package "corrplot"
#> [1] "corrplot"

M<-cor(mtcars)
res1 <- cor.mtest(mtcars, conf.level = .95)

corrplot(cor(mtcars),
         method="square",
         type="lower",
         p.mat = res1$p,
         insig = "label_sig",
         sig.level = c(.001, .01, .05),
         pch.cex = 0.8,
         pch.col = "yellow",
         tl.col="black",
         tl.cex=1,
         addCoef.col = "black",
         tl.pos="n",
         outline=TRUE)

^{由代表 package (v2.0.1) 于 2021 年 10 月 13 日创建}