Возможно ли иметь перекрывающиеся графики плотности и линию на графике рассеяния?
Справочная информация : у меня есть данные, по которым я делаю регрессионную модель. Я беру уникальные значения одной из переменных в моих данных, и я вижу, как прогнозы меняются между этими уникальными значениями, и я строю эти прогнозы. Мне бы хотелось, чтобы плотность исходных данных, то есть плотность / гистограмма графика рассеяния уникальных значений, лежала на фоне моей предсказанной линии.
Для иллюстрации:
library(dplyr)
iris2 = iris %>% mutate(pw_int=round(Petal.Width))
iris2_summary = iris2 %>% group_by(pw_int) %>% summarize(mean=mean(Sepal.Length))
Итак, я делаю вид, что среднее значение Sepal.Length для каждого округленного значения Petal.Width является моей моделью.
library(ggpubr)
ggscatterhist(
iris, x = "Petal.Width", y = "Sepal.Length",
#color = "Species", # comment out this and last line to remove the split by species
#margin.plot = "histogram", density plots
margin.params = list(color = "red", size = 0.2)
)
Мне бы хотелось иметь только гистограмму / плотность по оси X, но я не знаю, как.
Покушение
lines(iris2_summary$pw_int,iris2_summary$mean,col="green",type="l")
ничего не делает для меня. Кажется, что следующий метод может быть, но он выглядит очень плохо:
par(new=TRUE)
plot( iris2_summary$pw_int, iris2_summary$mean, type="l", col="green" )
Существует также проблема, которую я не знаю, как получить плотности на заднем плане.
Желаемый выход :
