>>>>> David Winsemius <

[hidden email]>

>>>>> on Thu, 7 May 2015 12:10:25 -0700 writes:

> On May 6, 2015, at 7:00 PM, Benjamin Tyner wrote:

>> Hi

>>

>> I'm wondering if anyone is aware of an R package implementing (i.e.,

>> providing a pdf, cdf, and/or quantile function) for the continuous

>> binomial distribution? Specifically the one characterized here:

>>

>>

http://www2.math.uni-paderborn.de/fileadmin/Mathematik/AG-Indlekofer/Workshop/Satellite_meeting/ilenko.pdf >>

>> Figured I would check here first, before attempting to code it up myself.

> I found that reading the ArXiv version of that material was easier to understand:

>

http://arxiv.org/abs/1303.5990 > zipfR package has an implementation of the incomplete beta function that might make some of the coding of the pdf and cdf more simple.

Well.

To: David Winsemius <

[hidden email]>

Subject: Re: [R] package implementing continuous binomial?

In-Reply-To: <

[hidden email]>

References: <

[hidden email]>

<

[hidden email]>

X-Mailer: VM 8.2.0b under 24.3.1 (x86_64-redhat-linux-gnu)

Reply-To: Martin Maechler <

[hidden email]>

CC: maechler

--text follows this line--

Dear David,

>>>>> David Winsemius <

[hidden email]>

>>>>> on Thu, 7 May 2015 12:10:25 -0700 writes:

> On May 6, 2015, at 7:00 PM, Benjamin Tyner wrote:

>> Hi

>>

>> I'm wondering if anyone is aware of an R package implementing (i.e.,

>> providing a pdf, cdf, and/or quantile function) for the continuous

>> binomial distribution? Specifically the one characterized here:

>>

>>

http://www2.math.uni-paderborn.de/fileadmin/Mathematik/AG-Indlekofer/Workshop/Satellite_meeting/ilenko.pdf >>

>> Figured I would check here first, before attempting to code it up myself.

> I found that reading the ArXiv version of that material was easier to understand:

>

http://arxiv.org/abs/1303.5990 > zipfR package has an implementation of the incomplete beta function that might make some of the coding of the pdf and cdf more simple.

> Searching done with Graves' very useful utility package:

> library('sos')

> findFn("incomplete beta function")

Hmm... R's pbeta() function is a pretty good implementation

of the incomplete beta function ... as is tries to say on its

help page.

If you look closely, these functions (for incomplete gamma,

incomplete beta, and their inverses) are simple wrappers to

pgamma() and qgamma() --- sometimes "regularizing" and sometimes

not -- where regularization is simply a multiplication/division

with gamma() or beta().

I wonder why you did not find R's help page about the beta distribution

{ -> functions dgamma, pgamma, qgamma, rgamma }

which does mention the "incomplete beta function" prominently.

I don't think Benjamin should use the zipfR package just for

these functions [and even the zipfR package help page on these

can be read as saying so .. ]

In the end I wonder if the "continuous Binomial" is not just a

version of the good old Beta distribution... as indeed the

Binomial and the Beta are related in the same way

that the Gamma and the Poisson are.

Martin Maechler

ETH Zurich

> (I did't think that doing a search on "continuous Binomial" was likely to be helpful, but I tried it anyway and did not find any functions named "continuous binomial" in their help page titles.)

> --

> David Winsemius

> Alameda, CA, USA

______________________________________________

[hidden email] mailing list -- To UNSUBSCRIBE and more, see

https://stat.ethz.ch/mailman/listinfo/r-helpPLEASE do read the posting guide

http://www.R-project.org/posting-guide.htmland provide commented, minimal, self-contained, reproducible code.