Sampling on the sphere from f(x) ∝ x^TAx
A method for drawing random samples of unit vectors x in R^p with density proportional to x^TAx where A is a symmetric, positive definite matrix. Includes an R function which implements the method.
READ FULL TEXT