Гладкий график плотности ядра

Я пытаюсь создать график апостериорной вероятности. Код, который я использую:

qplot(N2016, geom="density", xlim=c((n_star_2016 -5), (CIs2016[[3]])),
      xlab="Número estimado del total de asesinatos 2016",
      ylab="Probabilidad de la estimación", adjust=5)

где N2016 — апостериорная вероятность на 100 000 выборок.

Проблема, с которой я сталкиваюсь, заключается в том, что когда я использую график типа «плотность», я получаю некоторые вероятности со значениями, которые не заслуживают доверия, потому что мои наблюдения равны 160. Следовательно, я должен получить вероятность, равную нулю до 160.

График с корректировкой=5

img1

Я играл с параметром настройки, но это не похоже на ответ, так как кривая совсем не гладкая.

График с корректировкой = 0,25

img2

Я также прилагаю график с использованием гистограмм, чтобы показать, что проблема не в данных:

Гистограмма

img3

Кто-нибудь знает, как получить гладкую кривую, которая не начинается до наблюдений? В данном случае 160?

Спасибо!


r graph ggplot2 kernel-density
person ThesisStudent    schedule 15.11.2019    source источник

Ответы (1)


arrow_upward
3
arrow_downward

Используйте параметр trim:

r <- rnorm(100)
ggplot() + geom_density(aes(r), trim = FALSE) + xlim(-4, 4),
ggplot() + geom_density(aes(r), trim = TRUE) + xlim(-4, 4)

введите здесь описание изображения

person Axeman    schedule 15.11.2019
comment
Сработало идеально! Спасибо - person ThesisStudent; 21.11.2019