R/mantel.randtest.R

In ade4: Analysis of Ecological Data : Exploratory and Euclidean Methods in Environmental Sciences

"mantel.randtest" <- function(m1, m2, nrepet=999) {
    if (!inherits(m1, "dist")) 
        stop("Object of class 'dist' expected")
    if (!inherits(m2, "dist")) 
        stop("Object of class 'dist' expected")
    n <- attr(m1, "Size")
    if (n != attr(m2, "Size")) 
        stop("Non convenient dimension")
    m1 <- as.matrix(m1)
    m2 <- as.matrix(m2)
    col <- ncol(m1)
    isim<-testmantel(nrepet, col, as.matrix(m1), as.matrix(m2))
    obs<-isim[1]
    return(as.randtest(isim[-1],obs,call=match.call()))
}