r - R计算的伽马分布密度错误？

Question

我正在尝试以数值方式计算边际可能性（边缘化一个正参数）。我使用 Gamma 分布作为该参数的先验。在这里，我查看了两个特定参数设置的 Gamma 分布行为：

s = 28.4; r = 17000
plot(x, dgamma(x, shape=s, rate = r), type = 'l', ylab = 'density')
abline(v = s/r, col = 'red')

我得到以下结果：

然后我尝试了以下方法来获得更紧密的 Gamma 分布：

lines(x, dgamma(x, shape=s*1000, rate = r*1000), col = 'blue')

结果：

我很困惑。随着分布越来越紧密，高度应该越来越高，否则面积不会整合到 1。我错过了什么吗？还是有任何数字问题？谢谢！

Answer 1

您的x变量需要有更多样本才能捕获第二个密度函数中的窄峰：

x = seq(0, .01, .000001)
s = 28.4; r = 17000
plot(x, dgamma(x, shape=s, rate = r), type = 'l', ylab = 'density')
abline(v = s/r, col = 'red')
lines(x, dgamma(x, shape=s*1000, rate = r*1000), col = 'blue')