, Robert A LaBudde wrote:
At 05:10 PM 7/5/2010, Gregory Gentlemen wrote:
Dear R-users,
Is there an R function to compute the multinomial beta function? That
is, the normalizing constant that arises in a Dirichlet distribution.
For example, with three parameters the beta function is
Beta(n1,n2,n2
Dear R-users,
Is there an R function to compute the multinomial beta function? That is, the
normalizing constant that arises in a Dirichlet distribution. For example, with
three parameters the beta function is Beta(n1,n2,n2) =
Gamma(n1)*Gamma(n2)*Gamma(n3)/Gamma(n1+n2+n3)
Thanks in advance
How about this?
mbeta - function(...) {
exp(sum(lgamma(c(...)))-lgamma(sum(c(...
}
gamma(5)*gamma(6)*gamma(7)/gamma(18)
[1] 5.829838e-09
mbeta(5,6,7)
[1] 5.829838e-09
On Mon, 2010-07-05 at 17:10 -0400, Gregory Gentlemen wrote:
Dear R-users,
Is there an R function to compute
At 05:10 PM 7/5/2010, Gregory Gentlemen wrote:
Dear R-users,
Is there an R function to compute the multinomial beta function?
That is, the normalizing constant that arises in a Dirichlet
distribution. For example, with three parameters the beta function
is Beta(n1,n2,n2) = Gamma(n1)*Gamma(n2
4 matches
Mail list logo