问题描述
我正在尝试从一个较高的不完全伽马函数生成随机数，并且我意识到SciPy中实现的反不完全伽马函数gammainccinv与gammaincc并不完全对应。以下代码可用于说明差异：
import numpy
from scipy.stats import poisson
from scipy.special import gammainc,gammaincc
from scipy.special import gammaincinv,gammainccinv
import matplotlib.pyplot as plt
lam = 10
xr = numpy.arange(0,20,0.01)
xri = numpy.arange(20)
y = numpy.arange(0,1,0.001)
plt.plot(xri,poisson.pmf(xri,lam),"*",label="Poisson")
plt.plot(xr-0.5,(gammaincc(xr+1e-6,lam)-gammaincc(xr,lam))/1e-6,label="Upper incomplete gamma")
x = gammainccinv(lam,y)
plt.plot(0.5*(x[1:]+x[:-1]),(y[:-1]-y[1:])/(x[1:]-x[:-1]),label="Inverse upper incomplete gamma\n(transformed)")
plt.legend(fontsize=8)
plt.show()
在Python中是否存在逆上不完全伽马函数的替代实现？
解决方法
暂无找到可以解决该程序问题的有效方法，小编努力寻找整理中！
如果你已经找到好的解决方法，欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@）