Hi everyone, I was comparing a colleague’s MATLAB code to my version (we are coding in parallel, and I am using Julia. I am also trying to convince him to move to Julia. Making progress… :slight_smile: ), and I noticed that the beta function returns something different for large x or y. Because of …

To build on other’s answers that express beta() as \Gamma(p)\Gamma(q)/\Gamma(p+q), it seems you can trivially reproduce Matlab’s behavior in Julia using loggamma() : julia> p = beta1 3.6e13 julia> q = 1/alpha1 + 1 2.7746751786287764 julia> exp(loggamma(p) + loggamma(q) - loggamma(p+q)) 4.473779306…

Differences in MATLAB's versus Julia's beta(x, y) for large x or y?

General Usage

mbauman Split this topic April 28, 2021, 2:15pm 8

4 posts were split to a new topic: gamma(::BigFloat) overflows on Windows

Topic		Replies	Views
gamma(::BigFloat) overflows on Windows General Usage specialfunctions	7	1022	April 29, 2021
Result of exponentials are (sometimes) different from Matlab Numerics	5	386	February 15, 2024
Wrong results in calculation Performance	21	3229	March 12, 2019
Beta Distribution - full precision may not have been achieved Statistics	7	1893	June 13, 2018
Unexpected result about integer General Usage	13	1359	August 17, 2018

Differences in MATLAB's versus Julia's beta(x, y) for large x or y?

Related topics