如何在ggplot转plotly图表中正确显示stat_cor计算的R和P值

问题核心在于stat_cor生成的是ggplot专属的表达式文本（包含italic()这类语法标签），而plotly无法解析ggplot的表达式语法，导致标签代码直接显示在图表上。以下是两种保留自动计算逻辑的解决方法：

方法一：直接修改stat_cor输出为HTML格式文本 #

通过设置stat_cor的output.type为"text"，并手动构造带HTML斜体标签的文本，让plotly能正确渲染格式：

library(ggplot2)
library(plotly)
library(ggpubr)

p1 <- ggplot(iris) +
  geom_point(aes(Sepal.Length, Petal.Length)) +
  stat_cor(
    mapping = aes(Sepal.Length, Petal.Length),
    output.type = "text",
    label.sep = "\n",
    # 用HTML <em>标签替代ggplot的italic()语法
    aes(label = paste0(
      "<em>R</em> = ", round(r, 2), "\n",
      "<em>P</em> = ", format.pval(p.value, digits = 2)
    ))
  )

# 转换时让plotly解析HTML文本
p2 <- ggplotly(p1) %>%
  style(textposition = "top right")
p2

方法二：提取stat_cor计算值后添加HTML注释 #

先通过ggplot_build提取stat_cor自动计算的相关系数和p值，再用annotate添加HTML格式的文本：

library(ggplot2)
library(plotly)
library(ggpubr)

# 生成原始图，用于提取计算值
temp_plot <- ggplot(iris) +
  geom_point(aes(Sepal.Length, Petal.Length)) +
  stat_cor(mapping = aes(Sepal.Length, Petal.Length))

# 从ggplot对象中提取stat_cor计算的r和p值
cor_results <- ggplot_build(temp_plot)$data[[2]]
r_value <- round(cor_results$r, 2)
p_value <- cor_results$p.value

# 重新绘制图表，添加HTML格式的注释文本
p1_fixed <- ggplot(iris) +
  geom_point(aes(Sepal.Length, Petal.Length)) +
  annotate(
    "text",
    x = max(iris$Sepal.Length) - 1,
    y = max(iris$Petal.Length) - 0.5,
    label = paste0("<em>R</em> = ", r_value, "<br><em>P</em> = ", format.pval(p_value, digits = 2)),
    parse = FALSE # 关闭ggplot的表达式解析，保留HTML标签
  )

# 转换为plotly并确保注释正确渲染
p2 <- ggplotly(p1_fixed) %>%
  layout(annotations = list(
    text = paste0("<em>R</em> = ", r_value, "<br><em>P</em> = ", format.pval(p_value, digits = 2)),
    x = max(iris$Sepal.Length) - 1,
    y = max(iris$Petal.Length) - 0.5,
    showarrow = FALSE,
    xref = "x", yref = "y"
  ))
p2

注意事项 #

• 确保安装并加载最新版本的ggpubr包，避免参数兼容性问题
• 两种方法都保留了stat_cor的自动计算逻辑，无需手动输入相关系数和p值
• HTML标签<em>用于实现斜体效果，plotly能直接识别并渲染

内容的提问来源于stack exchange，提问作者Arda Askin