Apollo Fitted Line:
\(\text{Y-intercept} = 4.4944\)
\(\text{Slope} = -1.6430\)
\(\text{p-value} = 0.0144\)
\(\text{Residual standard error} = 1.386 \text{ on 103 degrees of freedom}\)
\(\text{Multiple }R^2 = 0.0568\)
\(\text{Adjusted } R^2 = 0.0476\)

Amor Fitted Line:
\(\text{Y-intercept} = -13.809\)
\(\text{Slope} = 15.989\)
\(\text{p-value} = 0.0629\)
\(\text{Residual standard error} = 4.9 \text{ on 59 degrees of freedom}\)
\(\text{Multiple }R^2 = 0.0574\)
\(\text{Adjusted } R^2 = 0.0414\)

Aten Fitted Line:
\(\text{Y-intercept} = 1.7814\)
\(\text{Slope} = -0.8847\)
\(\text{p-value} = 0.00085\)
\(\text{Residual standard error} = 0.127 \text{ on 13 degrees of freedom}\)
\(\text{Multiple }R^2 = 0.5882\)
\(\text{Adjusted } R^2 = 0.5566\)

# Plots three scatter plots based on the orbit class of the object
fig1 <- plot_ly(x = x1,
        y = y1,
        type = "scatter", mode = "markers",
        name = unique(df$Orbit.Class)[1],
        width = 800, height = 430) %>%
  add_lines(x = x1, y = fitted(mod1), name = "Apollo Fitted") %>%
  add_markers(x = x2, y = y2, name = unique(df$Orbit.Class)[2]) %>%
  add_lines(x = x2, y = fitted(mod2), name = "Amor Fitted") %>%
  add_markers(x = x3, y = y3, name = unique(df$Orbit.Class)[3]) %>%
  add_lines(x = x3, y = fitted(mod3), name = "Aten Fitted")