Описание тега gamma-function

Описание тега Вопросы с тегом
Все, что связано с математической гамма-функцией, также известной как обобщенная факториальная функция.
3 ответа

Быстрый алгоритм для лог-гамма-функции

Я пытаюсь написать быстрый алгоритм для вычисления лог-гамма-функции. В настоящее время моя реализация кажется наивной и просто повторяется 10 миллионов раз для вычисления журнала гамма-функции (я также использую numba для оптимизации кода). import …
24 фев '19 в 10:33
1 ответ

Реализация гамма-функции не дает правильных значений

Функция, запрограммированная в Fortran 95 для вычисления значений гамма-функции из математики, не дает правильных значений. Я пытаюсь реализовать рекурсивную функцию в Fortran 95, которая вычисляет значения гамма-функции, используя приближение Ланцо…
23 дек '18 в 17:44
2 ответа

Как я могу оптимизировать эту функцию S-кривой?

Я работаю над гамма-функцией, которая генерирует "S-кривую". Мне нужно запустить его в среде реального времени, поэтому мне нужно максимально ускорить его. Код выглядит следующим образом: float Gamma = 2.0f; //Input Variable float GammaMult = pow(0.…
18 янв '16 в 11:52
1 ответ

Числовые задачи точности в R?

У меня проблема со следующей функцией в R: test <- function(alpha, beta, n){ result <- exp(lgamma(alpha) + lgamma(n + beta) - lgamma(alpha + beta + n) - (lgamma(alpha) + lgamma(beta) - lgamma(alpha + beta))) return(result) } Теперь, если вы вс…
06 авг '15 в 10:59
2 ответа

Неполный алгоритм гамма-функции

Существует очень краткий алгоритм для вычисления нижней неполной гамма-функции: https://people.sc.fsu.edu/~jburkardt/f_src/asa147/asa147.html Мы закодировали это в C++. В этом алгоритме есть одна вещь, которую я не понимаю. В одном месте вычислить с…
14 сен '17 в 17:17
1 ответ

Верхняя Неполная Гамма-функция порядка 0 в scipy

Я пытаюсь реализовать верхнюю неполную гамма-функцию нулевого порядка в Python. Обычно мы используем gammaincc функция, но в соответствии с документами, он определен только для положительного а. Есть ли способ реализовать это в python для случая =0?…
05 ноя '18 в 03:15
1 ответ

В R неполная гамма-функция со сложным вводом?

Неполные гамма-функции могут быть вычислены в R с помощью pgamma, или с помощью gamma_inc_Q из библиотеки (gsl), или с помощью gammainc из библиотеки (expint). Однако все эти функции требуют только реального ввода. Мне нужна реализация неполной гамм…
06 дек '17 в 00:31
0 ответов

Деление двух неполных гамма-функций

Я работаю над уравнением, которое необходимо разделить две неполные гамма-функции. Они просто отличаются по первому члену, поскольку они являются IG(a,b) и IG(2a,b). Есть ли замкнутая форма или простое приближение для деления?
15 ноя '15 в 16:08
1 ответ

Гамма-распределение scipy не соответствует формуле в Википедии

Согласно Википедии, ожидаемое значение гамма-распределения k * theta, Однако, когда я подключаю это к CDF для гамма-распределения, я не получаю 0.5, как ожидалось. k = 1.5 theta = 2.1 expected_value = k * theta scipy.stats.gamma.cdf(expected_value, …
25 окт '17 в 01:02
3 ответа

Как эффективно кешировать и применить scipy gammaln к массиву numpy?

Я ищу способ кешировать результаты функции gammaln, применить ее к массиву и затем суммировать результаты. В настоящее время я могу кэшировать результаты gammln с помощью простого словарного подхода, но, похоже, я не могу эффективно итерировать масс…
04 мар '17 в 20:53
1 ответ

rcpp: как применить гамма-функцию к скаляру?

Я думаю, функция gamma only works for a vector as the input. Is there a way to apply it to a scalar, say,гамма (3)`? На самом деле, я бы получил правильный вывод, если я включу gamma(3) как часть моего кода, но есть предупреждение.... #include <R…
16 апр '18 в 18:16
1 ответ

График гамма-плотности: dgamma работает, а собственная функция возвращает ошибку

Я пытаюсь построить функцию плотности гамма для параметров (форма) альфа =3457 и (скорость) бета =84. Если я сделаю это с помощью: curve(dgamma(x,shape=3457,rate=84),from=35,to=50,xlab="posterior theta",ylab="density") все работает просто отлично, и…
18 фев '16 в 18:03
3 ответа

Гамма-функция в Python

Я довольно плохо знаком с Python, но хотел бы использовать гамма-функцию Эйлера в функции, которую я пишу. Я бы предпочел не записывать это как интеграл, и мне было интересно, есть ли что-то, что я могу импортировать, чтобы легко определить гамма-фу…
04 ноя '16 в 17:04
1 ответ

Накопительная функция распределения с гамма-функцией в Python

Я имею дело с функцией яркости Шехтера, которая выглядит следующим образом: phi(L)dL = norm. Factor * (L/Lstar)^(a) * exp (L/Lstar) d(L/Lstar) Скажем, L/Lstar это л. Аналитическое решение его кумулятивной функции распределения дается гамма-функцией:…
20 ноя '13 в 19:32
1 ответ

Rcpp гамма-интеграл

Я пытаюсь переписать в (R) CPP оригинальную функцию R, которая использует гамма-функцию (от двойного ввода). Ниже первоисточник. При компиляции с sourceCpp возникает следующая ошибка: "нет подходящей функции для вызова" gamma(Rcpp::traits::storage_t…
27 апр '17 в 12:05
0 ответов

Matlab: ошибка при использовании gammaln... при построении психометрических функций

Я использую набор инструментов psignifit для построения психометрических функций из набора данных. Мой код выглядит по существу следующим образом: load data.mat options = struct; options.sigmoidname = 'logistic'; result = psignifit(data,options); % …
09 янв '18 в 19:38
1 ответ

Декодировать гамма-кодированное изображение

Я работаю над проектом, в котором пользователь может указать гамма-функцию на изображении, например: y=(pow(10,(1023*x-681)/444)-.0408)/(1-.0408)) Затем изображение кодируется с помощью этой гамма-кривой, и тогда мой код должен иметь возможность вер…
25 июл '15 в 18:48
1 ответ

Экспоненциальный интеграл и гамма-функция

Чтобы решить эту проблему, без итерации: рекурсия: сумма серий из n членов расчет для заданного n: 1 + 2 * 3 + 3 * 4 * 5 + 4 * 5 * 6 * 7 +... + n * (n + 1)...(2n-1) с этим математическим ответом: https://math.stackexchange.com/questions/1590673/form…
27 дек '15 в 16:44
1 ответ

Исходный код языка R - что такое `NOMORE_FOR_THREADS` в коде C

Я шел через исходный код R для функции gamma() который указывает на этот исходный файл C. Фрагмент кода выглядит следующим образом: #ifdef NOMORE_FOR_THREADS static int ngam = 0; static double xmin = 0, xmax = 0., xsml = 0., dxrel = 0.; /* Initializ…
03 ноя '15 в 00:14
1 ответ

Функции тригаммы и дигаммы для переменной bigfloat в python

Мне нужно вычислить scipy.special Тригамма и дигамма функции переменной с типом bigfloat в Python, но я получаю следующее сообщение об ошибке: TypeError: ufunc 'psi' not supported for the input types, and the inputs could not be safely coerced to an…
02 авг '16 в 02:38