polykde (version 1.1.7)

d_unif_polysph: Density of the uniform distribution on the polysphere

Description

Computes the density of the uniform distribution on the polysphere.

Usage

d_unif_polysph(x, d, log = FALSE)

Value

A vector of size nx with the evaluated density.

Arguments

x

a matrix of size c(nx, sum(d) + r) with the evaluation points.

d

vector of size r with dimensions.

log

compute the logarithm of the density? Defaults to FALSE.

Examples

# Simple check of integration on S^1 x S^2
d <- c(1, 2)
x <- r_unif_polysph(n = 1e4, d = d)
mean(1 / d_unif_polysph(x = x, d = d)) / prod(rotasym::w_p(p = d + 1))