dotcvm: Calculate 3D covariance matrix using unscaled congruence coefficient. Skips any missing values in computation of covariance matrix

Description

Calculate 3D covariance matrix using unscaled congruence coefficient. Skips any missing values in computation of covariance matrix

Usage

1
dotcvm(A)

Arguments

A

An N x 3 x M array where N is the number of landmarks, 3 is the number of dimensions, and M is the number of specimens.

Details

This function does not guarantee that the returned matrix is positive definite. If the covariance matrix is not positive definite a warning is given and the matrix can be bent to create the closest positive definite matrix with as.matrix(Matrix::nearPD(mat)$mat).

Value

N x N covariance matrix

Examples

1
2
A <- array(rnorm(4 * 2 * 3), dim = c(2, 3, 4)) 
A.cvm <- dotcvm(A)

Questions? Problems? Suggestions? or email at ian@mutexlabs.com.

All documentation is copyright its authors; we didn't write any of that.