Оценка плотности по одной точке - PullRequest
Новая
       85

Оценка плотности по одной точке

1 голос
/ 20 октября 2019

Я пытаюсь построить плотности с помощью непараметрического метода оценки ядра. Для этой цели я использую асимметричные ядра, такие как Weibull, Gamma, Inverse Gaussian и т. Д. Я успешно построил графики плотности, приведенные в их статьях. Мой запрос касается тех графиков, которые построены в одной точке x , особенно x = 0 , как указано в оценках Бирнбаума-Сондерса и Логнормального ядра для длительностей моделирования на высокой частотеФинансовые данные на странице # 108, а также некоторые другие представили то же самое. Простые графики плотности могут быть построены с использованием следующего кода R. Пожалуйста, помогите мне, как построить на конкретной точке. 

n <- 200
k <- 400
y <- rexp(n, 1)

h <- 0.79 * IQR(y) * length(y) ^ (-1/5)
x <- seq(min(y) + 0.05, max(y), length=k)
Kbs <- matrix(rep(0, k * n), ncol=k)
fhat <- rep(0, k)

########### BS ###########

for (j in 1:k) {
  for (i in 1:n) {
    Kbs[i, j] <- (1 / (2*sqrt(2*h*pi))) *
      ((sqrt(1 / (x[j]*y[i]))) + (sqrt(x[j] / (y[i]^3)))) *
      exp(- (y[i] / (2*h*x[j])) + (1/h) - (x[j]/(2*h*y[i])))
    Kbs[is.nan(Kbs)] <- 0
  }
  fhat[j] <- 1/n * (sum(Kbs[, j]))
}

d1 <- density(y, bw=h)

plot(x, fhat, type="s", ylab="Density Function", lty=1, xlab="Time")
lines(d1, type="p", col="red")
legend("topright", c("Real Density", "Density by Birnbaum-Saunders Kernel"),
       col=c("red", "black"), lty=c(1, 2))

С уважением.

...