ry_bmf: Helper Function for Sampling a Bingham-von...

Description Usage Arguments Value Author(s) References Examples

View source: R/ry_bmf.R

Description

C interface to perform a Gibbs update of y with invariant distribution proportional to exp( sum(l*y^2+y*d) with respect to the uniform measure on the sphere.

Usage

1
ry_bmf(y, l, d)

Arguments

y

a normal vector.

l

a vector.

d

a vector.

Value

a normal vector

Author(s)

Peter Hoff

References

Hoff(2009)

Examples

1
2
3
4
5
6
## The function is currently defined as
function (y, l, d) 
{
    .C("ry_bmf", y = as.double(y), l = as.double(l), d = as.double(d), 
        n = as.integer(length(y)))$y
  }

rstiefel documentation built on June 12, 2018, 5:19 p.m.