将注释/描述性标签成对添加到每个子图

Question

我有几对数据集，我想为每个数据集添加相关系数，以便可以直观地看到相关性。例如：

corr1 <- cor(richstruc[,c(4,10,11)])
pairs(richstruc[,c(4,10,11)], main = "comparison 1")

有没有什么方法可以将注释“p =（来自corr1的相应值）”添加到对中的每个相应子图？

Answer 1

使用参数panel，您可以控制每个面板中显示的内容。您可以自定义添加边缘标签。我增加了边距大小（gap=2）并减小了字体大小（cex=0.8）以使文本适合。由于您不提供数据，我将使用内置的虹膜数据进行说明。

Display = function(x,y) {
    points(x,y, pch=16, col=rainbow(3)[iris$Species])
    CXY = round(cor(x,y), 3)
    mtext(paste("cor =", CXY), 1, cex=0.8)
}
pairs(iris[1:3], panel=Display, gap=2)