geomorph (version 4.0.8)

compare.ZVrel: Comparisons of Effect Sizes from Overall Integration Analyses

Description

Function performs an analysis to compare the effect sizes of two or more Vrel effects

Usage

compare.ZVrel(..., two.tailed = TRUE)

Value

An object of class compare.rel.eig, returns a list of the following

sample.Re.obs

A vector of observed Vrel for each sample.

sample.Z.obs

A vector of effect sizes for each sample.

sample.Z.var

A vector of variances for each effect size.

pairwise.z

A matrix of pairwise, two-sample z scores between all pairs of effect sizes.

pairwise.p

A matrix of corresponding P-values.

Arguments

...

saved analyses of class rel.eig

two.tailed

A logical value to indicate whether a two-tailed test (typical and default) should be performed.

Author

Dean Adams

Details

The function statistically compares the effect sizes of two or more Vrel analyses. Typically, this function might be used to compare the strength of integration in one dataset as compared with another (see Conaway and Adams 2022).

The analysis performs two-sample z-tests based on effect sizes (Z-scores) of Vrel. The method follows that of Conaway Adams (2022) used to compare the strength of integration across datasets.

To use this function, simply perform integration.Vrel on as many samples or as desired. Any number of objects of class rel.eig can be input. One may perform the comparison as either a one-tailed or a two-tailed (default) test.

References

Conaway, M.A., and D.C. Adams. 2022. An effect size for comparing the strength of morphological integration across studies. Evolution. 76: 2244-2259.

Examples

if (FALSE) {
 data("plethodon")
 Y.gpa <- gpagen(plethodon$land)
 
 coords.gp <- coords.subset(Y.gpa$coords, plethodon$species)
 Vrel.gp <- Map(function(x) integration.Vrel(x), coords.gp) 
 
 out <- compare.ZVrel(Vrel.gp$Jord, Vrel.gp$Teyah)
 
 summary(out)
 }