Morpho (version 2.12)

CAC: calculate common allometric component

Description

calculate common allometric component

Usage

CAC(x, size, groups = NULL, log = FALSE)

Value

CACscores

common allometric component scores

CAC

common allometric component

x

(group-) centered data

sc

CAC reprojected into original space by applying CAC %*% x

RSCscores

residual shape component scores

RSC

residual shape components

gmeans

groupmeans

CS

the centroid sizes (log transformed if log = TRUE)

Arguments

x

datamatrix (e.g. with PC-scores) or 3D-array with landmark coordinates

size

vector with Centroid sizes

groups

grouping variable

log

logical: use log(size)

References

Mitteroecker P, Gunz P, Bernhard M, Schaefer K, Bookstein FL. 2004. Comparison of cranial ontogenetic trajectories among great apes and humans. Journal of Human Evolution 46(6):679-97.

Examples

data(boneData)
proc <- procSym(boneLM)
pop.sex <- name2factor(boneLM,which=3:4)
cac <- CAC(proc$rotated,proc$size,pop.sex)
plot(cac$CACscores,cac$size)#plot scores against Centroid size
cor.test(cac$CACscores,cac$size)#check for correlation
#visualize differences between large and small on the sample's consensus
if (FALSE) {
large <- restoreShapes(max(cac$CACscores),cac$CAC,proc$mshape)
small <- restoreShapes(min(cac$CACscores),cac$CAC,proc$mshape)
deformGrid3d(small,large,ngrid=0)
}