# trying to prove out negbinomdist

Nicholas Ferguson nicholasferguson at wingarch.com
Thu Jun 4 01:44:49 UTC 2020

```Hi  Eike,

Thank you for taking time to give me answer.  Much appreciated.
Regards
Nick

-----Original Message-----
From: LibreOffice [mailto:libreoffice-bounces at lists.freedesktop.org] On Behalf Of Eike Rathke
Sent: Wednesday, June 03, 2020 3:15 PM
To: Nicholas Ferguson
Cc: LibreOffice at lists.freedesktop.org
Subject: Re: trying to prove out negbinomdist

Hi Nicholas,

On Wednesday, 2020-06-03 08:49:01 -0400, Nicholas Ferguson wrote:

> This might be an unusual question for this site.  But I cannot prove out the negative binomial distribution, by so called hard coding the function per LibreOffice docs….a jpg of that function is attached.  Does anyone have a clue. Though I will add…Microsoft Excel has the same problem.   I did search for its corresponding code in C++/H in libreoffice source…but that didn't help either…
>
> =negbinomdist(1;14.4;.92)
>
> =(FACT(1+14.4-1)/(FACT(1)*FACT(14.4-1)))*(power(.92;14.4)*power(1-.92;
> 1))

You are overlooking that X and R are defined to be integer values, so
=negbinomdist(1;14.4;.92)
effectively is calculated as
=negbinomdist(1;14;.92)
and if you substitute 14.4 with 14 in your manual approach to form
=(FACT(1+14-1)/(FACT(1)*FACT(14-1)))*(POWER(0.92;14)*POWER(1-0.92;1))
you'll discover that both results are almost equal, apart from some precision error of 5.55111512312578E-17