Р: Остановка опорной линии на перехват в плотности участка

Question

У меня есть следующий график плотности:

Мои входные данные выглядят следующим образом:

PA <- c(3028.0,3083.0,2958.0,2889.0,2758.0,2815.0,2877.0,3145.0,3072.0,3056.0,2856.0,2934.0,3112.0,2913.0,2837.0)

Я использовал следующие команды для построения графика:

plot(density(PA), main="Avidin PA")
abline(v=2913, col="red")

Я добавил базовую линию в красном, как я делаю упор опорной линии, когда он перехватывает участок

Answer 1

Нам нужно будет немного поработать, чтобы выяснить, какова плотность приблизительно при данном значении x. Как только мы получим это, мы можем использовать команду линий для построения линии интереса. Это должно дать начало:

PA <- c(3028.0,3083.0,2958.0,2889.0,2758.0,2815.0,2877.0,3145.0,3072.0,3056.0,2856.0,2934.0,3112.0,2913.0,2837.0)

den.PA <- density(PA)
lineat = 2913

# Find which x value that the density was computed
# for is the closest to the x value we want a line at
lineheight <- den.PA$y[which.min(abs(den.PA$x - lineat))]
plot(den.PA, main="Avidin PA",)
lines(c(lineat, lineat), c(0, lineheight), col = "red")