Как мне добавить медианные значения в этом графике регрессии через децили? - PullRequest
0 голосов
/ 11 июля 2019

У меня линейная регрессия по децилям переменной (сильно искажена), и я нанес линию на децили, чтобы предсказать изменчивость результата. График дает мне линию регрессии через децили, однако я хочу добавить медианные значения каждого дециля на оси X, чтобы лучше объяснить распределение моей переменной.

ggplotRegression <- function (fit) {

require(ggplot2)

ggplot(fit$model, aes_string(x = names(fit$model)[2], y = names(fit$model)[1])) + 
  geom_point() +
  stat_smooth(method = "lm", col = "red") +
  labs(title = paste("Adj R2 = ",signif(summary(fit)$adj.r.squared, 5),
                     "Intercept =",signif(fit$coef[[1]],5 ),
                     " Slope =",signif(fit$coef[[2]], 5),
                     " P =",signif(summary(fit)$coef[2,4], 5)))
}

ggplotRegression(lm(math_score ~ lead_decile, data=mkf1))

...