在R中规范化2F1超几何

时间:2016-01-24 11:32:52

标签: r wolfram-mathematica

在Mathematica wolfram中,有一个计算超几何2F1()函数的函数。例如在Mathematica中,两个计算如下:

in:   Hypergeometric2F1Regularized[15, 8, 10, -0.651542]
out:  7.86105*10^-9
in:   Hypergeometric2F1Regularized[31, 13, 18, -2.651542]
out:  5.04779*10^-26

在R中我尝试了以下等效,但结果与上述不相似:

library("BAS")
> hypergeometric2F1(15, 8 ,10, -0.651542,log=F )
[1] 2.441767e-05
> hypergeometric2F1(31, 13, 18, -2.651542,log=F)
[1] Inf
Warning message:
In hypergeometric2F1(31, 13, 18, -2.651542, log = F) :
  integral in 2F1 diverges

我想知道R中Mathematica的Hypergeometric2F1Regularized []函数的确切等价。

感谢

1 个答案:

答案 0 :(得分:5)

hypergeometric2F1(15, 8, 10, -0.651542, log = FALSE) / factorial(10 - 1)
[1] 7.861054e-09

因为正则化版本为2F1(a,b,c;z) / Gamma(c)