joint_angle_cluster: Function that estimates joint rank by method of profile...
In justicesuker/DMMD_Package: Double Matched Matrix Decomposition

joint_angle_cluster

R Documentation

Function that estimates joint rank by method of profile likelihood.

Function that estimates joint rank by method of profile likelihood.

joint_angle_cluster(
  angle_vec,
  angle_threshold = 90 * pi/180,
  variance = c("equal", "unequal")
)

`angle_vec`	A vector of principal angles
`angle_threshold`	Optional threshold (radians) that is used to truncate principal angles. Default is 90 degrees, which means all the principal angles are taken into account.
`variance`	Either "equal" or "unequal", i.e., whether the assumption is equal variance or unequal variance. If no input is given, equal variance is assumed.

A list with the following elements:

`joint_rank`	Estimated joint rank
`profileloglikvec`	Profile log likelihood calculated at each index. Function returns NA with a warning message, if less or equal to 2 principal angles are smaller than the threshold.