R/eigenQsphere.R
eigenQsphere.RdThe eigenvalues of the precision matrix Q with dimension (mm, nn)
and polar angle limits phimin, phimax are calculated.
eigenQsphere(phimin, phimax, mm, nn)
| phimin | Polar angle minimum. |
|---|---|
| phimax | Polar angle maximum. |
| mm | Number of rows of precision matrix Q. |
| nn | Number of columns of precision matrix Q. |
A list containing 2 elements:
eigval Row vector containing the eigenvalues of Q.
eigvec Matrix containing the eigenvectors of Q as columns.
The corresponding function for data on a grid is eigenLaplace.
eig_out <- eigenQsphere(phimin = 180/10, phimax = 180 - 180/10, mm = 10, nn = 20)