gamma-function

2熱度

1回答

我想繪製一個伽瑪密度函數的參數（形狀）阿爾法= 3457和（速度）貝塔= 84。如果我這樣做使用： curve(dgamma(x,shape=3457,rate=84),from=35,to=50,xlab="posterior theta",ylab="density") 一切正常，我得到的密度集中在〜41。而如果我先建立密度爲： para_density = function(x,s=34

0熱度

3回答

Python中的Gamma函數

我對Python很新，但想在我寫的函數中使用歐拉的伽瑪函數。我不想把它寫成一個整體，並想知道是否有我可以導入的東西，很容易定義gamma函數。謝謝

1熱度

1回答

使用scipy.integrate.quad集成伽瑪函數

當a = 3/2使用scipy.integrate.quad時，如何從0找到Gamma(x) = x^(a-1) * e^-x的積分？到目前爲止，我已經嘗試了一些代碼，但我被告知「quad：第一個參數不可調用」。這是我所使用的代碼： import numpy as np import math as m import matplotlib.pyplot as plt from math

0熱度

1回答

R中的digamma函數的根

我正在研究最大似然估計器，其中一個參數是使用digamma函數估計的。我試圖用uniroot解決方程，但無法做到這一點。這裏是我的代碼： dig = function(alpha){ digamma(2 + alpha) - digamma(alpha) - (1/(2+alpha)) + (2/(2+alpha)) } curve(dig, from = 0, to = 10)

2熱度

2回答

任意精度gamma函數

我在C++中實現了一個任意的精度算術庫，並且在實現gamma函數時我幾乎陷入了困境。通過使用等價gamma(n) = gamma(n - 1) * n和gamma(n) = gamma(n + 1)/n，分別，我可以在對所有實數值x範圍(1; 2]獲得的有理數r。我不知道如何評估gamma(r)。對於Lanczos近似（https://en.wikipedia.org/wiki/Lanczos_

-1熱度

1回答

是否有R A功能，尋找到Z，如果我知道的Y伽馬（Z）= Y

我期待知道是否有在R，它尋找到Z，如果我知道Y和伽瑪（Z任何功能）= y Uniroot可能是有用的，但不知道我該如何使用它。由於

1熱度

1回答

的R - 日誌不完全γ函數的導數

以下功能中的R是可用的： gamma計算伽馬函數 digamma計算衍生物日誌伽馬函數的 pgamma來計算不完全γ函數？來計算對數不完全伽馬函數的導數我很想知道什麼函數可以計算對數不完全伽馬函數的導數。我注意到gsl包有一個函數gamma_inc，但不知道如何計算這個函數的日誌的導數。如果沒有函數存在，是否有一個簡單的方法來近似這個導數R？

0熱度

2回答

不完全γ函數算法

沒有計算完全伽瑪功能非常簡潔的算法： https://people.sc.fsu.edu/~jburkardt/f_src/asa147/asa147.html 我們在C++編碼此。有一件事我不明白這個算法。在一個地方來計算下面的表達式：它被取代：顯然，這是一樣的，但是它爲什麼這樣做呢？計算lgamma的exp是否比計算tgamma函數更有效率（在C++中都可以使用lgamma和tgamma）

0熱度

1回答

Rcpp伽馬積分

我想重寫（R）CPP原始R函數，利用伽馬函數（來自雙輸入）。低於原始來源。當與sourceCpp配合時，出現以下錯誤：「沒有匹配函數調用'gamma（Rcpp :: traits :: storage_type（< 14> :. type）'' gamma函數應該放在sugar 。下面用平均值），所以我想到應該簡單地調 #include <Rcpp.h> #include <math.h> u

2熱度

1回答

R語言源代碼 - C代碼中的NOMORE_FOR_THREADS是什麼

我正在通過指向this C source file的函數gamma()的R源代碼。代碼的片段如下： #ifdef NOMORE_FOR_THREADS static int ngam = 0; static double xmin = 0, xmax = 0., xsml = 0., dxrel = 0.; /* Initialize machine depend