Power of Mantel's Test

Power calculation of Mantel's permutation test.

mantelPower(obj, effect.size = seq(0, 1, length.out = 50), alpha = 0.05)

Arguments

obj	an object of class "mantelTest". See mantelTest.
effect.size	numeric; the effect size specifying the alternative hypothesis.
alpha	numeric; the significance level at which to compute the power level.

Value

A data frame containing the effect size and its respective power level.

References

Silva, A.R.; Dias, C.T.S.; Cecon, P.R.; Rego, E.R. (2015). An alternative procedure for performing a power analysis of Mantel's test. Journal of Applied Statistics, 42(9): 1984-1992.

Author

Anderson Rodrigo da Silva <anderson.agro@hotmail.com>

See also

mantelTest

Examples

# Mantel test
data(garlicdist)
garlic <- tocher(garlicdist)
coph <- cophenetic(garlic)
mt1 <- mantelTest(garlicdist, coph, xlim = c(-1, 1))

# Power calculation, H1: rho = 0.3
mantelPower(mt1, effect.size = 0.3)
#>   effect.size     power
#> 1         0.3 0.7597598

# Power calculation, multiple H1s and different alphas
p01 <- mantelPower(mt1, alpha = 0.01)
p05 <- mantelPower(mt1, alpha = 0.05)
p10 <- mantelPower(mt1, alpha = 0.10)
plot(p01, type = "l", col = 4)
lines(p05, lty = 2, col = 4)
lines(p10, lty = 3, col = 4)
legend("bottomright", c("0.10", "0.05", "0.01"), 
  title = expression(alpha), col = 4, lty = 3:1, cex = 0.8)

# End (Not run)