Regression plot with R-sqr value

regression_plot.R

# Adding regression line to a plot
set.seed(1234)


# Creating a random dataset
x <- rnorm(n = 1000, mean = 10, sd = 3)
y <- 2 * x + 10 + rnorm(n = 1000, mean = 0, sd = 4)

d <- data.frame(y = y,
                x = x)

head(d)


# fitting regression model
fit <- lm(y ~ x, data = d)

summary(fit)

# The R-sqr valu
summary(fit)$adj.r.squared

library(plotly)


# Plotting the data
plot_ly(data = d,
        x = ~ x,
        y = ~ y,
        type = "scatter",
        mode = "markers",
        name = "Data") %>%
  add_lines(x = ~ x,
            y = fitted(fit),
            name = "Regression Line",
            line = list(color = "red", width = 3, dash = "dash")) %>%
  # Add the R^2 as text annotation
  add_annotations(text = paste("R^2 = ", round(summary(fit)$adj.r.squared, 2), sep = ""),
                  xref = "paper",
                  yref = "paper",
                  x = 0.1,
                  y = 0.9,
                  showarrow = FALSE) %>%
  layout(title = "My Title...",
         yaxis = list(title = "Y axis"),
         xaxis = list(title = "X axis"),
         margin = list(t = 60, b = 60))