View source: R/compare.ZVrel.R
compare.ZVrel | R Documentation |
Function performs an analysis to compare the effect sizes of two or more Vrel effects
compare.ZVrel(..., two.tailed = TRUE)
... |
saved analyses of class rel.eig |
two.tailed |
A logical value to indicate whether a two-tailed test (typical and default) should be performed. |
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.
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. |
Dean Adams
Conaway, M.A., and D.C. Adams. 2022. An effect size for comparing the strength of morphological integration across studies. Evolution. 76: 2244-2259.
## 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.