R中归一化的低不完全伽玛

normalized lower incomplete gamma in R

我试图在 R 中重现 this 论文的结果,但我被困在两个部分。

  1. 绘制归一化下不完全伽马函数

  2. 计算损失函数

如果 R 中存在计算 normalized lower incomplete gamma function 的函数,则此问题将得到解决。所以我的问题是,如何计算 normalized lower incomplete gamma function

pgamma是归一化下不完全伽马函数,见?pgamma:

pgamma is closely related to the incomplete gamma function. As defined by Abramowitz and Stegun 6.5.1 (and by ‘Numerical Recipes’) this is

P(a,x) = 1/Gamma(a) integral_0^x t^(a-1) exp(-t) dt