R plotly（）：将回归线添加到相关散点图

Question

我想将回归线添加到我的相关散点图中。不幸的是，这并不真正适用于 plot_ly()。我已经尝试过本论坛其他帖子的一些解决方案，但没有用。

我的数据框如下所示（只是其中的一部分）：

我的绘图代码和实际绘图输出如下所示：

CorrelationPlot <- plot_ly(data = df.dataCorrelation, x = ~df.dataCorrelation$prod1, 
                           y = ~df.dataCorrelation$prod2, type = 'scatter', mode = 'markers',
                           marker = list(size = 7, color = "#FF9999", line = list(color = "#CC0000", width = 2))) %>%
                    layout(title = "<b> Correlation Scatter Plot", xaxis = list(title = product1), 
                           yaxis = list(title = product2), showlegend = FALSE)

我想要的是这样的：

这是我用 ggscatter() 函数生成的：

library(ggpubr)
  ggscatter(df.dataCorrelation, x = "prod1", y = "prod2", color = "#CC0000", shape = 21, size = 2,
            add = "reg.line", add.params = list(color = "#CC0000", size = 2), conf.int = TRUE, 
            cor.coef = TRUE, cor.method = "pearson", xlab = product1, ylab = product2)

如何使用 plot_ly() 获得回归线？？

代码编辑：

CorrelationPlot <- plot_ly(data = df.dataCorrelation, x = ~df.dataCorrelation$prod1, 
                           y = ~df.dataCorrelation$prod2, type = 'scatter', mode = 'markers',
                           marker = list(size = 7, color = "#FF9999",
                             line = list(color = "#CC0000", width = 2))) %>%
                   add_trace(x = ~df.dataCorrelation$fitted_values, mode = "lines", type = 'scatter',
                             line = list(color = "black")) %>%
                   layout(title = "<b> Correlation Scatter Plot", xaxis = list(title = product1), 
                           yaxis = list(title = product2), showlegend = FALSE)

给予：

如何在这里找到回归线的线？？

Answer 1

我不认为有像 ggscatter 这样的现成函数，很可能你必须手动完成，比如首先拟合线性模型并将值添加到 data.frame。

我做了一个 data.frame 就像你的数据：

set.seed(111)
df.dataCorrelation = data.frame(prod1=runif(50,20,60))
df.dataCorrelation$prod2 = df.dataCorrelation$prod1 + rnorm(50,10,5)

fit = lm(prod2 ~ prod1,data=df.dataCorrelation)
fitdata = data.frame(prod1=20:60)
prediction = predict(fit,fitdata,se.fit=TRUE)
fitdata$fitted = prediction$fit

线的上下界只是预测的1.96*标准误差：

fitdata$ymin = fitdata$fitted - 1.96*prediction$se.fit
fitdata$ymax = fitdata$fitted + 1.96*prediction$se.fit

我们计算相关性：

COR = cor.test(df.dataCorrelation$prod1,df.dataCorrelation$prod2)[c("estimate","p.value")]
COR_text = paste(c("R=","p="),signif(as.numeric(COR,3),3),collapse=" ")

并将其放入 plotly:

library(plotly)

df.dataCorrelation %>%
plot_ly(x = ~prod1) %>%
add_markers(x=~prod1, y = ~prod2) %>%
add_trace(data=fitdata,x= ~prod1, y = ~fitted, 
mode = "lines",type="scatter",line=list(color="#8d93ab")) %>%
add_ribbons(data=fitdata, ymin = ~ ymin, ymax = ~ ymax,
line=list(color="#F1F3F8E6"),fillcolor ="#F1F3F880" ) %>%
layout(
    showlegend = F,
    annotations = list(x = 50, y = 50,
    text = COR_text,showarrow =FALSE)
)

Answer 2

另一种选择是使用 ggplotly 作为

library(plotly)
ggplotly(
ggplot(iris, aes(x = Sepal.Length, y = Petal.Length))+
  geom_point(color = "#CC0000", shape = 21, size = 2) +
  geom_smooth(method = 'lm') +
  annotate("text", label=paste0("R = ", round(with(iris, cor.test(Sepal.Length, Petal.Length))$estimate, 2),
                                ", p = ", with(iris, cor.test(Sepal.Length, Petal.Length))$p.value), 
x = min(iris$Sepal.Length) + 1, y = max(iris$Petal.Length) + 1, color="steelblue", size=5)+
  theme_classic()
)

R plotly（）：将回归线添加到相关散点图

R plotly(): Adding regression line to a correlation scatter plot

r

plotly

regression

scatter