B {hyperdirichlet}R Documentation

Normalizing constant for the hyperdirichlet distribution

Description

Uses numerical techniques for calculating the normalizing constant for the hyperdirichlet distribution

Usage

B(x, ...)
NC(x)
calculate_B(x, disallowed=NULL, give=FALSE, ...)
probability(x, disallowed, ...)
mgf(x, powers, ...)
mean(x, ...)
is.proper(x,irregardless)
validated(x)

Arguments

x Object of class “hyperdirichlet” (or coerced thereto)
powers Vector of length dim(x) whose elements are the powers of the expectation; see details section
irregardless Boolean; see details section
disallowed Function specifying a subset of the simplex over which to integrate; default NULL means to integrate over the whole simplex. The integration proceeds over p with disallowed(p) evaluating to FALSE
give Boolean, with default FALSE meaning to return the value of the integral and TRUE meaning to return the full output of adapt()
... Further arguments passed to adapt()

Details

Author(s)

Robin K. S. Hankin

See Also

hyperdirichlet

Examples


a <- hyperdirichlet(c(4,3,6,5,4,3,2,1))
B(a)                                    # Not recommended
a <- as.hyperdirichlet(a,TRUE)          # Recommended

mgf(a,powers=1:3)    # expectation of p1^1 * p2^2 * p3^3


[Package hyperdirichlet version 1.1-7 Index]