关于Gamma密度模拟

Question

我正在模拟一些伽马随机数

plot(density(rgamma(10000,8.1,rate=0.00510)),lwd=2,las=1,cex.axis=0.75,
main=expression(paste("Gamma Distribution with",' scale ',alpha," and  rate 
",beta)))

plot(density(rgamma(10000,2.1,rate=0.00110)),lwd=2,las=1,cex.axis=0.75,
 main=expression(paste("Gamma Distribution with",' scale ',alpha," and  rate 
",beta)))


plot(density(rgamma(10000,2.1,rate=110)),lwd=2,las=1,cex.axis=0.75,
 main=expression(paste("Gamma Distribution with",' scale ',alpha," and  rate 
",beta)))

我需要用一些尾部模拟这种伽玛，平均值大约为1200.我一直在选择随机数，以便在考虑伽马分布的期望和方差的定义时得到这些值，但在第一种情况下，我得到负数，我不想要那个。在第二种情况下相同但在两个图中y轴的概率都很低，我想增加这个概率，但我不知道如何选择足够的参数来获得它。

另一方面，第三个图中给出的参数给出了一个奇怪的密度，因为y轴的概率大于1我可以得到大于1的值。我不明白这一点。

Answer 1

在第一种情况下你没有得到负数。

dgamma

密度超过1时，没有任何问题。应该是曲线下面积恰好为1的区域。

为了绘制分布图，您也可以使用rgamma代替x <- seq(0,.12, length = 10000) y <- dgamma(x,2.1,rate=110) plot(x,y, type = "l") ，如下所示：

ggplot2

以下是使用ggplot(data.frame(x,y), aes(x,y)) + geom_line() 的情节。

{{1}}