如何求出density核密度图的概率

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

plot（density（））做出和密度图，但是如何求图中某一个区间的概率，也就是如何对一个区间进行积分啊？

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏3 回帖

关键词：Density 核密度图 sit 核密度 plot 如何

给你个例子你就明白了：



n <- 100
y <- rnorm(n, 5, 1)

width <- 0.01
scal <- seq(-2, 12, by = 0.01)

ll <- length(scal)

den_y <- density(y, n = ll, from = -2, to = 12)  # here defines the range and width

plot(scal, den_y$y, type = "l")

# use numerical integration
# for example, you need pr(a < x < b)
a <- 2
b <- 4

begin <- which.max(den_y$x[den_y$x < a])
end <- which.max(den_y$x[den_y$x < b])


prob <- sum(den_y$y[begin:end]*width)

# confirm

true_prob <- pnorm(4, 5, 1) - pnorm(2, 5, 1)  # increase n to 100000, true_prob = prob. In small sample, they are different

谢谢楼上回答，对我很有启发，谢谢。

kingswordsc 发表于 2015-7-16 20:16
谢谢楼上回答，对我很有启发，谢谢。

记得加分啊！

不太明白是什么意思，可以解释一下为什么这么求么？