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	1024	April 29, 2021
Beta Distribution - full precision may not have been achieved Statistics	7	1900	June 13, 2018
Wrong results in calculation Performance	21	3241	March 12, 2019
Unexpected result about integer General Usage	13	1371	August 17, 2018
Why is this Julia code considerably slower than Matlab New to Julia performance	64	8461	March 5, 2017

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

Related topics