calculate common allometric component
calculate common allometric component
CAC(x, size, groups = NULL, log = FALSE)
x |
datamatrix (e.g. with PC-scores) or 3D-array with landmark coordinates |
size |
vector with Centroid sizes |
groups |
grouping variable |
log |
logical: use |
CACscores |
common allometric component scores |
CAC |
common allometric component |
x |
(group-) centered data |
sc |
CAC reprojected into original space by applying |
RSCscores |
residual shape component scores |
RSC |
residual shape components |
gmeans |
groupmeans |
CS |
the centroid sizes (log transformed if |
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.
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 ## Not run: large <- showPC(max(cac$CACscores),cac$CAC,proc$mshape) small <- showPC(min(cac$CACscores),cac$CAC,proc$mshape) deformGrid3d(small,large,ngrid=0) ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.