Rでは、複雑な入力を伴う不完全ガンマ関数ですか？

Question

不完全なガンマ関数は、pgammaの場合はR、ライブラリ（gsl）の場合はgamma_inc_Q、ライブラリ（expint）の場合はgammaincで計算できます。しかし、これらの関数はすべて実際の入力のみを取ります。

複雑な入力を行う不完全なガンマ関数の実装が必要です。具体的には、私は第1引数の整数と第2引数の複素数（積分の限界）を持っています。

この関数は複雑な入力（Wikipediaを参照）に対してよく定義されており、Mathematicaで計算しています。 Rには組み込まれていないようだが、どのライブラリでもそれを見ることはできない。

アルゴリズムを調べてCで実装し、Rインタフェースを記述するよりも、これらの計算を行う方が簡単な方法を提案できますか？

（私はそれを自分自身を実装する必要がない場合は、ここで私が見つけた、複雑な入力のための唯一のアルゴリズムです：Kostlan & Gokhman 1987）は

Answer 1

ここでは、あなたが下不完全ガンマ関数をしたいと仮定すると、実装したものです。私はWolframといくつかの値を比較し、それらが一致しています。

library(CharFun) 

incgamma <- function(s,z){ 
    z^s * exp(-z) * hypergeom1F1(z, 1, s+1)/s 
} 

多分sの評価が失敗します。