dotcvm: Calculate 3D covariance matrix using unscaled congruence...

In timcdlucas/paleomorph: Geometric Morphometric Tools for Paleobiology

Description

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

Usage

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

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

timcdlucas/paleomorph documentation built on May 31, 2019, 1:48 p.m.