Gibbs Sampling for the Vector-variate Bingham Distribution
吉布斯的矢量变量宾厄姆分布的抽样
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Simulate a random normal vector from the Bingham distribution using Gibbs sampling.
随机模拟法向量的Bingham分布Gibbs抽样。
用法----------Usage----------
rbing.vector.gibbs(A, x)
参数----------Arguments----------
参数:A
a symmetric matrix.
一个对称矩阵。
参数:x
the current value of the random normal vector.
随机法线向量的当前值。
值----------Value----------
a new value of the vector x obtained by Gibbs sampling.
一个新的矢量值xGibbs抽样获得。
注意----------Note----------
This provides one Gibbs scan. The function should be used iteratively.
这提供了1吉布斯扫描。应反复使用的功能。
(作者)----------Author(s)----------
Peter Hoff
参考文献----------References----------
实例----------Examples----------
## The function is currently defined as[#功能目前被定义为]
function (A, x)
{
evdA <- eigen(A)
E <- evdA$vec
l <- evdA$val
y <- t(E) %*% x
x <- E %*% ry_bing(y, l)
x/sqrt(sum(x^2))
}