добавить звезду значения p в матрицу корреляции R

Question

Я использую коды с этого сайта для построения матрицы рассеяния и корреляции: http://www.sthda.com/english/wiki/scatter-plot-matrices-r-base-graphs

panel.cor <- function(x, y){
    usr <- par("usr"); on.exit(par(usr))
    par(usr = c(0, 1, 0, 1))
    r <- round(cor(x, y), digits=2)
    txt <- paste0("R = ", r)
    cex.cor <- 0.8/strwidth(txt)
    text(0.5, 0.5, txt, cex = cex.cor * r)
}
# Customize upper panel
upper.panel<-function(x, y){
  points(x,y, pch = 19, col = my_cols[iris$Species])
}
# Create the plots
pairs(iris[,1:4], 
      lower.panel = panel.cor,
      upper.panel = upper.panel)

Как рассчитать значение p и добавить звезды рядом с коэффициентом корреляции? Я считаю, что пакет «PerformanceAnalytics» может это сделать, но я хочу использовать функцию smoothScatter на диаграмме рассеяния на верхней панели. Я не знаю, как использовать функцию smoothScatter в PerformanceAnalytics. В настоящее время моя функция верхней панели выглядит следующим образом. Копируется с другого сайта

pairs(df, lower.panel = panel.cor,
        upper.panel = function(...) smoothScatter(..., nrpoints = 0, add = TRUE), gap = 0.2)

Answer 1

 r2 <- cor.test(x, y)$p.value
 r3 <- symnum(r2, cutpoints = c(0, 0.001, 0.01, 0.05, 1), 
                   symbols = c("***","**","*",""))

может использоваться для получения значения p.

, поэтому вся функция panel.cor выглядит следующим образом:

panel.cor <- function(x, y){
  usr <- par("usr"); on.exit(par(usr))
  par(usr = c(0, 1, 0, 1))
  r <- round(cor(x, y), digits=2)
  r2 <- cor.test(x, y)$p.value
  r3 <- symnum(r2, cutpoints = c(0, 0.001, 0.01, 0.05, 1), 
                   symbols = c("***","**","*",""))
  txt <- paste(r, r3, sep = " ")
  cex.cor <- 0.8/strwidth(txt)
  text(0.5, 0.5, txt, cex = cex.cor*r)
}