Мне интересно, существует ли базовый способ R для получения точной плотности ядра в любой желаемой точке? Например, как я могу получить точную плотность ядра в 3 следующих точках -2 , 0, +2 по оси X на графике, как показано ниже?
set.seed(2937107)
plot( density(rnorm(1e4)) )
Используйте линейную интерполяцию, чтобы найти его.
d <- density(rnorm(10000))
approx(d$x, d$y, xout = c(-2, 0, 2))
Точность интерполяции может быть выше, если вы установите большее значение n в density. По умолчанию n = 512, поэтому интерполяция основана на 512 точках.
