Я пытаюсь получить лучшее понимание оценки плотности ядра.
Использование определения из Википедии: https://en.wikipedia.org/wiki/Kernel_density_esvaluation#Definition
Давайте возьмем в качестве прямоугольной функции, которая дает если находится в от до и противном случае, а (размер окна) равен 1.
Я понимаю, что плотность - это свертка двух функций, но я не уверен, что знаю, как определить эти две функции. Одна из них должна (вероятно) быть функцией данных, которая для каждой точки в R сообщает нам, сколько точек данных у нас в этом месте (в основном ). А другой функцией, вероятно, должна быть какая-то модификация функции ядра в сочетании с размером окна. Но я не уверен, как это определить.
Какие-либо предложения?
Ниже приведен пример кода R, который (я подозреваю) копирует настройки, которые я определил выше (со смесью из двух гауссианов и ), на которых я надеюсь увидеть «доказательство» того, что функции, которые должны быть свернуты, являются, как мы подозреваем, ,
# example code:
set.seed(2346639)
x <- c(rnorm(50), rnorm(50,2))
plot(density(x, kernel='rectangular', width=1, n = 10**4))
rug(x)