r - R中的伽马分布(qgamma)
问题描述
这可能是一个非常愚蠢的问题,但有谁知道为什么我没有得到第二位等于平均值(100)?
#beta=4, alpha=5, mean=20
qgamma(0.5, 5, 1/4)
# 18.68364
#beta=2500, alpha=0.04, mean=100
qgamma(0.5,0.04,1/2500)
# 0.00004320412
解决方案
这是因为您使用的是分位数函数,并且qgamma(0.5, shape, scale)
对应于中位数 - 而不是您期望的平均值。
请参见下面的示例;
x <- rgamma(50000, shape = 0.04, scale = 2500)
mean(x)
# [1] 98.82911
median(x)
# [1] 3.700012e-05
推荐阅读
- c# - c#标签文本数据绑定不刷新文本
- responsive-design - 是否可以在 Semantic-UI-React 中创建一个在移动设备上看起来很棒的响应式表格?
- javascript - Rails 选择表单:在提交之前显示和使用选择表单值?
- oracle - 在不同的循环中填充相同的 sys_refcursor
- javascript - 从 http 网页请求本地资源时响应为 401 错误
- netsuite - 如何在 Netsuite 中的高级 pdf 检查模板中获取采购订单项目
- linux - 如何从用户 1 切换到用户 2 并在 linux 中运行一些命令?
- wordpress - WooCommerce 中的产品彼此过于接近
- excel - 停止用目标文件名更改源文件名
- c++ - 为什么 QFileDialog::selectFile 不起作用?