Morpho (version 2.1)

CAC: calculate common allometric component

Description

calculate common allometric component

Usage

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

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)

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)

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
large <- showPC(max(cac$CACscores),cac$CAC,proc$mshape)
small <- showPC(min(cac$CACscores),cac$CAC,proc$mshape)
deformGrid3d(small,large,ngrid=0)