Как можно разбить данные на 2 раздела из пика на графике плотности?

Question

я могу получить пик на графике и нарисовать линию на нем, но теперь мне нужно разделить мои данные на 2 части (до пика) и (после пика)

как это сделать?

это код, который я использовал, чтобы получить пик и нарисовать на нем линию:

library(data.table)
library(ggplot2)

ggplot(DS, aes(DS$loan_amount_ink)) + geom_density()


Ymax <- which.max(density( DS$loan_amount_ink)$y)  

Xmax <- density(DS$loan_amount_ink)$x[Ymax]       


ggplot(DS, aes(DS$loan_amount_ink)) + geom_density() + 
       geom_vline(xintercept = density( DS$loan_amount_ink)$x[Ymax])

Answer 1

library(ggplot2)
library(dplyr)
View(faithful)
ggplot(faithful, aes(waiting)) + geom_density()

which.max(density(faithful$waiting)$y) # Peak
peak <- density(faithful$waiting)$x[326]
peak # 79.96245
ggplot(faithful, aes(waiting)) + 
  geom_density() +
  geom_vline(xintercept = density(faithful$waiting)$x[326])

Разделение набора данных в соответствии с пиковой плотностью.

df1 <- faithful %>% 
  dplyr::filter(waiting <= peak)

df2 <- faithful %>% 
  dplyr::filter(waiting > peak)

Источник - найти пиковую плотность: http://ianmadd.github.io/pages/PeakDensityDistribution.html