从给定X值的密度函数中检索Y值

Question

给出一个简单的密度直方图和如下所示的曲线，如何检索给定X值的Y值。例如mean(dat)处的y值

dat<-c(5,7,4,6,4,3,55,6,7,5,4,3,33,44,5,2,33,22)
hist (dat,freq=F)
lines(density(dat), col="red", lwd=2)

感谢。

Answer 1

您可以approxfun()使用density的结果来获得近似密度的函数

dat<-c(5,7,4,6,4,3,55,6,7,5,4,3,33,44,5,2,33,22)
hist (dat,freq=F)
lines(d<-density(dat), col="red", lwd=2)

#get density function
dd<-approxfun(d$x, d$y)
dd(mean(dat))
# [1] 0.015039

#plot results
abline(v=mean(dat), lty=2)
points(mean(dat), dd(mean(dat)), cex=1.2, pch=20, col="blue")