问题标签 [gamma-function]

我正在绘制我的物种中雌性和雄性的飞行速度和时间之间的关系。我的广义线性混合模型（个人 ID 的随机截距）表明女性与男性之间存在差异，所以在图中，我想展示这些差异。

到目前为止，我有以下情节：

但现在我想添加线条来显示这种关系。我尝试了两种不同的方法：

1) 按物种使用 geom_smooth 和 facet

2）从我的 GLMM 中获取斜率和截距值并使用 abline

这些似乎都没有按我的意愿工作。所以，我的问题是，我怎样才能以一种对我的模型有意义的方式来展示女性和男性飞行速度的关系？

作为参考，我的模型是：

我的数据是：

我正在推导截断伽马分布的表达式。我推导出了一个包含不完整伽马函数的表达式。当我使用'gammainc（x，alpha）'函数在matlab中实现相同的值时，我得到了复杂的值，因为我的x值是负数。有没有其他方法可以计算负值的不完全伽马函数？

我想查找详细信息，例如随机伽马计算误差及其范围。但目前我还没有找到。如果有人知道它以及由 NetLogo 中的随机伽玛语法引起的一定范围的风险，请告诉我。谢谢你。

我有两个清单。两者都包括标准化百分比：

• 实际人口分布 = [0.2,0.3,0.3,0.2]
• sample_population_distribution = [0.1,0.4,0.2,0.3]

我希望将这两个列表适合伽马分布，然后计算返回的两个列表以获得 KL 值。

我已经能够得到KL。

这是我用来计算 gamma 的函数：

将两个列表拟合到 gamma 分布中：

这是我用来计算 KL 的代码：

上面的代码不会产生任何错误。

但我怀疑我对 gamma 所做的方式是错误的，因为np.mean/var要获得均值和方差。

事实上，这个数字不同于：

如果我用这种方式。

通过使用“ mean, var, skew, kurt = gamma.stats(fit_alpha, loc = fit_loc, scale = fit_beta, moments = 'mvsk')”，我将得到一个大于 1 的 KL 值，因此这两种方法都无法获得正确的 KL。

我想念什么？

我想计算 R 中 digamma 函数的导数。

R 有一个名为 的现有函数digamma，但我找不到 digamma 函数的导数函数。是否可以在 R 中为微分 digamma 函数编写代码？

谢谢。

我正在尝试在 Tensorflow 中计算 gamma 函数。函数中的参数是来自张量和虚部的浮点数之和。通常，我可以对tf.math.lgamma函数取幂，但是，它不接受复数。是否有任何解决方法可以解决此问题？

为了对非中心卡方分布进行抽样以对欧式看涨期权定价，使用 mp.sum 似乎是最佳选择；但是，我放在一起的程序给出了 TypeError: 'mpf' ​​对象在运行时不能被解释为整数。

我将代码发送给了一个朋友，他能够运行完全相同的代码并获得没有错误的输出。我不确定是什么原因造成的。问题应该来自我的 ncx2_samp1 和 ncx2_samp2 中的 mp.gammainc 和 mp.gamma 函数。

关于可能导致问题的任何想法？它应该保存为 mpf 并被解释为字符串；问题是由于内存引起的吗？

谢谢

I am running a GLM regression in Python using statsmodels using the following code. I specifically want to implement a log link function. I am able to write R like formulas using Statsmodels.

The following code successfully runs but throws up a Deprecation Warning. Can anyone suggest how to get rid of this warning. Thanks for the help.

Code:

Warning: DeprecationWarning: Calling Family(..) with a link class as argument is deprecated. Use an instance of a link class instead.

我使用 SageMath 来象征性地整合一个表达式。结果包含一个gamma带有两个输入参数的函数：

显然这被称为不完全伽马函数。现在我想在 Python 代码中使用这个表达式。我已将不完整的 gamma 函数追踪到scipy.special.gammainc. 不幸的是，这个函数不允许负输入参数，我必须使用-1它作为第一个输入参数。我该如何解决这个问题？

我试图估计逆 Gamma 分布的参数，因为它的分位数为 0.025 和 0.975。

目前，我发现rriskDistributions::get.gamma.par它可以让我估计给定 Gamma 分布的分位数的参数。但是，我无法弄清楚 Gamma 的分位数和逆 Gamma 之间的关系。

我应该如何继续，或者是否有可以为我执行此操作的软件包？