Мне нужно как можно точнее найти пик оценки плотности ядра (модальное значение непрерывной случайной величины). Я могу найти приблизительное значение:

x<-rlnorm(100)
d<-density(x)
plot(d)
i<-which.max(d$y)
d$y[i]
d$x[i]

Но при расчетеd$y точная функция известна. Как я могу найти точное значение режима?