r - 在 R 中绘制伽马分布

Question

我正在尝试使用 R 绘制伽马分布直方图

所以我有

gam(10, 0.5)

我之前计算的平均值为

10* 0.5 = 5

所以我应该用 100 个观测值绘制直方图scale = 10 and shape = 0.5

所以我试过了

x <- round(rgamma(100,shape = 0.5,rate = 10),1)
hist(x)

我明白了

这是错误的，因为平均值应该是 5，但我的情节没有产生 5

我哪里错了？

Answer 1

使用scale而不是rate. 所以模拟你的数据如下：

rgamma(100,shape = 0.5,scale = 10)

阅读文档以获取更多信息。

Answer 2

这是 Alex 建议的输出：

 x <- round(rgamma(100,shape = 0.5,rate = 10),1)
 hist(x)