compare.ZVrel {geomorph} | R Documentation |
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)
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. |
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.
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. |
Author(s)
Dean Adams
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
## Not run:
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)
## End(Not run)