Directional (version 4.0)

Simulation from a Matrix Fisher distribution on SO(3): Simulation from a Matrix Fisher distribution on SO(3)

Description

It simulates random samples (rotation matrices) from a Matrix Fisher distribution with any given parameter matrix, F (3x3).

Usage

rmatrixfisher(n, F)

Arguments

n

the sample size.

F

An arbitrary 3x3 matrix.

Value

An array with simulated rotation matrices.

Details

Firstly corresponding Bingham parameter A is determined for a given Matrix Fisher parameter F using John Kent (2013) algorithm and then Bingham samples for parameter A are generated using rbingham code. Finally convert Bingham samples to Matrix Fisher samples according to the Kent (2013) transformation.

References

Kent J.T., Ganeiber A.M. and Mardia K.V. (2013). A new method to simulate the Bingham and related distributions in directional data analysis with applications. http://arxiv.org/pdf/1310.8110v1.pdf

Examples

# NOT RUN {
F <- matrix( c(85, 11, 41, 78, 39, 60, 43, 64, 48), ncol = 3) / 10   ### An arbitrary F matrix
a <- rmatrixfisher(10, F)
# }