有沒有辦法在源代碼中查看 Scipy.Special.Gamma() function 實現？
[英]Is there a way to see Scipy.Special.Gamma() function implementation in the source?
我嘗試編寫函數來計算特殊函數（例如指數函數、伽瑪函數、erf 函數等），但為了進行求和或乘積運算，我使用了 10k 轉的 while 循環。 這是非常耗時的。 然后，我意識到 scipy.special.gamma() function 比我的 while 循環快得多，我想看看和研究 gamma ...
多重精度 Gamma function in R
[英]Multiple precision Gamma function in R
我需要計算 R 中的總和，其中涉及每個元素的伽瑪 function。 當 gamma 的 arguments 增加時，我得到 NaN 結果，我懷疑問題與 gamma 的評估有關。 我已經閱讀了 Rmpfr 和 gmp 的文檔，但我只發現了整數的階乘。 我也把代碼貼在這里，也許你對錯誤的來源有更好的 ...
將伽馬分布參數擬合到期望和概率
[英]Fit gamma distribution parameters to expectation and probabilty
我的問題如下。 我測量了一堆不同的物理特性，並將方法和測量的不確定性一直傳播到某種效率比。 對於我所有的物理特性，正態分布似乎是一個不錯的選擇，並且對於前幾次計算和相應的傳播，我的不確定性相當低。 我將所有不確定性作為 kp=2 的擴展不確定性轉發，這意味着覆蓋了所有可能值的 95.45%。 但是， ...
R 中的 Gamma 函數
[英]Gamma Function in R
我正在嘗試使用以下規則從頭開始實現伽馬函數： 如果 z 等於 1，我們返回 1； 項目清單 gama(z) (使用函數的遞歸) (z-1)*gamma(z-1); 如果 z 是一個自然數，那么我們返回 (z-1) 的階乘； 如果 z 等於 1/2，我們返回 sqrt(pi)。 ...
scipy 中 t 分布的熵：如何輸入 digamma 和 beta 函數的自由度？
[英]Entropy of t-distribution in scipy: How to input degrees of freedom to digamma and beta functions?
對於隨后的t分布的變量X，衍生的熵的封閉形式的解析解這里，是 看到python有函數scipy.special.digamma和scipy.special.beta ，上面的公式如何在代碼中實現？ 令我困惑的是，根據文檔，剛才提到的函數沒有將自由度參數nu(v)作為輸入。 一個正在運行的例子會 ...
為 Weibull 分銷續訂 Function
[英]Renewal Function for Weibull Distribution
t = 10的 Weibull 分布m(t)的更新 function 如下所示。 我想找到m(t)的值。 我編寫了以下r代碼來計算m(t) output 是m(t) = 2.670408e+93 。 我的迭代程序正確嗎？ 謝謝。 ...
計算負實數值的 Gamma function (C++, Boost)
[英]calculating the Gamma function for negative real values (C++, Boost)
以下無限嚴重需要計算非整數、負數、實數的階乘： （這是一種計算橢圓周長的方法，a 和 b 是半長軸和半短軸，h 定義為： h = (ab)^2/(a+b)^2) 階乘 function 可以通過 Gamma function 擴展到負值，該 Gamma 是為所有非負整數的實數定義的。 在編寫 ...
Bachelier Normal Implicit Vol Python 計算（幫助） Jekel
[英]Bachelier Normal Implied Vol Python Calculation (Help) Jekel
編寫 python 腳本來計算隱含正常體積； 符合 Jekel 文章（行業標准）。 https://jaeckel.000webhostapp.com/ImpliedNormalVolatility.pdf 他們說他們正在使用廣義不完全伽瑪 Function 逆。 來電： F(x)=v/(K ...
如果隨機數超過 0.96，Gamma_inv 返回錯誤
[英]Gamma_inv returns error if random number is over 0.96
如果傳遞給 Gamma_inv 的隨機數大約高於 0.96，則返回錯誤。 這發生在 Excel 和 VBA 中。 我正在運行 Monte Carlo，而這僅在 40,000 多次運行中發生 10 次左右。 這是我使用它的方式。 alpha 和 beta 總是 >0。 如果“gammatrun ...
估計給定分位數的逆伽馬分布的參數
[英]Estimate the parameters of an inverse gamma distribution given quantiles
我試圖估計逆 Gamma 分布的參數，因為它的分位數為 0.025 和 0.975。 目前，我發現rriskDistributions::get.gamma.par可以讓我估計給定 Gamma 分布的分位數的參數。 但是，我無法弄清楚 Gamma 的分位數和逆 Gamma 之間的關系。 我應該如 ...
scipy.special.gammainc 不接受負輸入
[英]scipy.special.gammainc does not accept negative input
我使用 SageMath 象征性地整合了一個表達式。 結果包含一個具有兩個輸入參數的gamma函數： 顯然，這被稱為不完全伽馬函數。 現在我想在 Python 代碼中使用這個表達式。 我已經將不完整的伽馬函數追蹤到scipy.special.gammainc 。 不幸的是，這個函數不允許負輸入參 ...
如何使用 R 類公式在 Python 中運行 GLM 伽馬回歸
[英]How to run a GLM gamma regression in Python with R like formulas
我使用以下代碼使用 statsmodels 在 Python 中運行 GLM 回歸。 我特別想實現一個日志鏈接功能。 我能夠使用 Statsmodels 編寫類似 R 的公式。 以下代碼成功運行，但引發棄用警告。 任何人都可以建議如何擺脫這個警告。 謝謝您的幫助。 代碼： 警告：Deprec ...
對非中心卡方分布求和 Gamma 函數導致“mpf”object 不能解釋為 python 中的 integer
[英]Summing Gamma functions for a non-central chi square distribution leads to 'mpf' object cannot be interpreted as an integer in python
為了對非中心卡方分布進行抽樣以對歐式看漲期權定價，使用 mp.sum 似乎是最佳選擇； 但是，我放在一起的程序在運行時給出了 TypeError: 'mpf' object 不能解釋為 integer 。 我將代碼發送給一個朋友，他能夠運行完全相同的代碼並獲得 output 而沒有錯誤。 我不確定 ...
是否存在計算digamma函數導數的函數？
[英]Is there a function calculating derivative of digamma function?
我想計算R中的digamma函數的導數。 R有一個稱為digamma的現有函數，但我找不到digamma函數的導數的函數。 是否可以為R中的差分digamma函數編寫代碼？ 謝謝。 ...
計算兩個伽瑪分布列表之間的KL散度
[英]Calculate KL Divergence between two gamma distribution list
我有兩個清單。 兩者都包括標准化百分比： actual_population_distribution = [0.2,0.3,0.3,0.2] sample_population_distribution = [0.1,0.4,0.2,0.3] 我希望將這兩個列表擬合 ...
在ggplot中繪制具有伽馬分布的模型
[英]Plotting model with gamma distribution in ggplot
我正在繪制我的物種中雌性和雄性的飛行速度和時間之間的關系。 我的廣義線性混合模型（個體 ID 的隨機截距）表明男性與女性之間存在差異，因此在圖中，我想展示這些差異。 到目前為止，我有以下情節： 但現在我想添加線條來顯示這種關系。 我嘗試了兩種不同的方法： 1) 按物種使用 geom_s ...
為什么對於大型數組，此對數伽馬numba函數比scipy慢，但對於單個值，它快嗎？
[英]Why is this log gamma numba function slower than scipy for large arrays, but faster for single values?
我有一個函數來計算用numba.njit裝飾的對數伽馬函數 。 當我對大型數組使用gammaln_nr ，例如np.linspace(0.001, 100, 10**7) ，我的運行時間比scipy慢大約7倍（請參見下面附錄中的代碼）。 但是，如果我運行任何單個值，我的numba函數總 ...
對數伽瑪函數的快速算法
[英]Fast algorithm for log gamma function
我正在嘗試編寫一種快速算法來計算對數伽馬函數 。 目前，我的實現似乎很幼稚，只是迭代了1000萬次來計算gamma函數的對數（我也使用numba來優化代碼）。 我將代碼計時為scipy.special.gammaln實現，而我的速度實際上慢了100,000倍。 因此，我正在做一些非常 ...
Gamma函數實現不生成正確的值
[英]Gamma function implementation not producing correct values
在Fortran 95中編程的函數用於從數學計算Gamma函數的值不會產生正確的值。 我試圖在Fortran 95中實現一個遞歸函數，使用Lanczos近似計算Gamma函數的值（是的，我知道在2003標准及以后版本中有一個內在函數）。 我非常仔細地遵循標准公式，所以我不確定是什么問題。 ...