Skip to contents

Mantel test (correlation between two distance matrices (in R).)

mantel.rtest.Rd

Performs a Mantel test between two distance matrices.

Usage

mantel.rtest(m1, m2, nrepet = 99, ...)

Arguments

m1

an object of class dist

m2

an object of class dist

nrepet

the number of permutations

...

further arguments passed to or from other methods

Value

an object of class rtest (randomization tests)

References

Mantel, N. (1967) The detection of disease clustering and a generalized regression approach. Cancer Research, 27, 209--220.

Author

Daniel Chessel
Stéphane Dray stephane.dray@univ-lyon1.fr

Examples

    data(yanomama)
    gen <- quasieuclid(as.dist(yanomama$gen))
    geo <- quasieuclid(as.dist(yanomama$geo))
    plot(r1 <- mantel.rtest(geo,gen), main = "Mantel's test")

    r1
#> Monte-Carlo test
#> Call: mantel.rtest(m1 = geo, m2 = gen)
#> 
#> Observation: 0.5095199 
#> 
#> Based on 99 replicates
#> Simulated p-value: 0.01 
#> Alternative hypothesis: greater 
#> 
#>     Std.Obs Expectation    Variance 
#>  3.36944904  0.01018575  0.02196164