ry_bmf {rstiefel} | R Documentation |
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.
ry_bmf(y, l, d)
y |
a normal vector. |
l |
a vector. |
d |
a vector. |
a normal vector
Peter Hoff
Hoff(2009)
## 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
}