Постоянная полоса пропускания в оценке плотности ядра (R) - PullRequest
0 голосов
/ 19 августа 2011

У меня сейчас некоторые трудности с R.

Что у меня есть:

  • вектор с минимальными значениями для каждого треугольного ядра
  • вектор с максимальными значениями
  • вектор с "средней" точкой (где должен быть пик каждого треугольника)

Что мне нужно:

  • оценка плотности (треугольные ядра ofc)

1 Ответ

0 голосов
/ 19 августа 2011

Посмотрите, можете ли вы опираться на этот пример изменяющихся спецификаций полосы пропускания для функции базовой плотности с триагулярным ядром:

data("faithful", package = "datasets") 
x <- faithful$waiting 
opar <- par(mfcol= c(2,2))
for (bw in c(3,6,12,18) ) {
  hist(x, xlab = "Waiting interval", ylab = "Eruption Freq", 
  probability = TRUE, main = "Density plot: Triangular kernel", 
  border = "gray", sub= bquote("Bandwidth = "*.(bw) ) ) 
  lines(density(x, width = bw, window = "triangular"), lty=3, col="red")
  }
par(opar)

Four bandwidths to triangular kernel density plot

...