如何在相关矩阵图中显示相关系数和显著性（用星号表示）？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

R代码：

1. panel.conf<-function (x, y, corr = NULL, col.regions, digits = 2, cex.cor,
   
2.     ...)
   
3. {
   
4.     auto <- missing(cex.cor)
   
5.     usr <- par("usr")
   
6.     on.exit(par(usr))
   
7.     par(usr = c(0, 1, 0, 1))
   
8.     if (!is.null(corr)) {
   
9.         est <- corr
   
10. 
    
11.         est <- formatC(est*100, digits = digits, format = "d")
    
12. 
    
13.         if (auto)
    
14.             cex.cor <- 0.7/strwidth(est)
    
15.         text(0.5, 0.6, est, cex = cex.cor)
    
16.     }
    
17.     else {
    
18.         results <- cor.test(x, y, alternative = "two.sided")
    
19.         est <- results$estimate
    
20.         est <- formatC(est*100, digits = digits, format = "d")
    
21.         Signif <- symnum(results$p.value, corr = FALSE, na = FALSE, cutpoints = c(0, 0.001, 0.01, 0.05, 0.1, 1), symbols = c("***", "**", "*", ".", " "))
    
22.         if (auto)
    
23.             cex.cor <- 0.7/strwidth(est)
    
24.         text(0.5, 0.6, est, cex = cex.cor)
    
25. 
    
26.         text(0.8, 0.8, Signif, cex = cex.cor, col = 2)
    
27. 
    
28.     }
    
29. }
    
30. corrgram(mtcars,order=TRUE,lower.panel=panel.shade,upper.panel=panel.conf,text.panel=panel.txt)

复制代码

以上代码运行仅得到图1，如何达到图2的效果呢？（图2源于http://www.talyarkoni.org/blog/2 ... s-are-correlicious/）
程序编写的具体问题就是，如何将panel.shade和panel.conf合并成一个函数，或者是否还有新的函数？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：相关矩阵 相关系数 Panel pane con 如何

who can help me?

感觉很复杂，顶一个

运行后
     Error in lower.panel(...) : object 'resultsestimate' not found

yijiaqin 发表于 2014-8-26 14:13
运行后
     Error in lower.panel(...) : object 'resultsestimate' not found

你再试试...

你好，如果我只想在右上角显示相关系数怎么弄呢？这个太复杂了，我运行完，右上角不显示相关系数

你用ggplot2试试 geom_tile()和geom_text()这两条函数就可以了

楼主，为什么我运行完之后得到的图里面上半部分数字和星号都巨小，要怎么调整参数才能有你上传的这种效果。

感觉很复杂，学习下

您好，请问您做出图2的效果了吗？