powered by
mantel.rtest(m1, m2, nrepet = 99)
if (require(mva, quietly=TRUE)) { 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 }
Run the code above in your browser using DataLab