project: Orthogonal projection onto diffusion signal space

In astamm/bayesr: An implementation of Bayes Space for Gaussians in R

Description

Orthogonal projection onto diffusion signal space

Usage

project(x, reference_tensor, distance = "bayes")

Arguments

x	A symmetric 3x3 matrix resulting from a series of operations in Bayes space.
reference_tensor	The SDP tensor defining the reference measure of the Bayes space.
distance	A string indicating whether the projection should be done in Bayes geometry (default) or in Euclidean geometry.

Value

An SDP tensor.

Examples

R <- diag(c(1.7, 0.3, 0.1), 3)
Rn <- matrix(mvtnorm::rmvnorm(1, c(R)), 3, 3)
det_matrix(Rn)
Rfro <- project(Rn, R, distance = "frobenius")
Rbay <- project(Rn, R, distance = "bayes")
eigen(Rn, TRUE, TRUE)$values
eigen(Rfro, TRUE, TRUE)$values
eigen(Rbay, TRUE, TRUE)$values

astamm/bayesr documentation built on May 10, 2019, 2:05 p.m.