R语言中已知Gamma分布累积概率0.95求解对应y值的方法

你当前的需求是计算普通Gamma分布的下侧分位数，不需要用到invgamma相关工具，R基础环境内置的分位数函数就可以直接实现。

• 核心调用函数为qgamma()，是R原生的Gamma分布分位数计算函数，不需要额外安装任何包
• 调用时注意参数匹配：你给出的形状参数r=3对应函数的shape参数，尺度参数ρ=6对应函数的scale参数，不要误用函数默认的速率（rate）参数

代码示例 #

# 直接计算满足P(X < y) = 0.95的y值
y <- qgamma(p = 0.95, shape = 3, scale = 6)
# 输出结果查看
print(y)

结果验证 #

你可以通过Gamma分布的累积分布函数验证计算结果是否正确：

# 验证累积概率是否为0.95
pgamma(q = y, shape = 3, scale = 6)

运行后输出值为0.95即可确认计算无误。

内容的提问来源于stack exchange，提问作者Stacy